external REG  "./regcoef_graffino.so"

; INPUT
      x = (/ 1190.,1455.,1550.,1730.,1745.,1770. \
          ,  1900.,1920.,1960.,2295.,2335.,2490. \
          ,  2720.,2710.,2530.,2900.,2760.,3010. /)

      y = (/ 1115.,1425.,1515.,1795.,1715.,1710. \
          ,  1830.,1920.,1970.,2300.,2280.,2520. \
          ,  2630.,2740.,2390.,2800.,2630.,2970. /)

      xmsg =  -999.0
      ymsg =  -999.0
      dmsg =  -999.0        

; create (preallocate)  return argumnts
; MUST be same type (float<->real)); (double<->double)
; MUST be same shape (scalar or array) 

      rcoef= xmsg    ; defined in fortran sub as 'real'
      tval = xmsg

      xave = dmsg    ; defined in fortran sub as double 
      yave = dmsg 
      rstd = dmsg 
      yint = dmsg 
      
      nptxy= 0       ; defined in fortran sub as integer
      ier  = -9

      N    = dimsizes(y)

; Shared Object

      REG::rcgraf (x,y,N,xmsg,ymsg       \    ; input
                  ,rcoef,tval,nptxy      \    ; output
                  ,xave,yave,rstd,yint,ier)   ; output

; Output
      print ("rcoef= "+ rcoef)
      print ("tval = "+ tval ) 
      print ("yint =" + yint )