50 CAL 



IF( 



PRI 



940 FOR 



IIHU 



c 



c 



60 

 930 



IS 

 REA 

 FOR 

 IF( 

 IF( 

 CON 

 CON 

 NDE 

 NDE 

 VX = 

 IF( 

 NDE 

 NDE 

 VY = 

 IF( 

 VEL 

 ANG 

 IF( 

 PRI 

 FOR 

 GO 

 RET 

 END 

 SUB 

 1 

 DESCRIP 

 LIBRARY 



950 



70 



L SP 

 lERR 

 NT 9 

 MAT ( 

 »9X» 

 READ 

 TO B 

 D 93 

 MAT( 

 X(2) 

 lERR 



2 = 2. 

 RIV( 

 R I V ( 

 CONl 

 lERR 

 RIV( 

 RIV( 

 C0N2 

 lERR 

 = (VX 

 = 57. 

 ANG. 

 NT 9 

 MAT ( 

 TO 6 

 URN 



LCW (N 



OR.GE 



40 



1H1,2 



1HV,6 



THE 

 E COM 

 0,X(2 

 F6.2. 

 .EQ.9 

 OR.GE 

 .2/SI 

 4/SIN 

 1)=0 

 2)=I 

 <>SPlO 

 OR.GE 

 1)=1 

 2) =0 

 »SPLD 

 OR.GF 

 *VX+V 

 296>A 

 LT.O. 

 50»X( 

 1H0,F 

 



DIM»NDATA»XMIN,XMAX»N0DES»XTRAP»C0EF*NCF,NWRK, I ERROR) 

 .32) GO TO 60 



X,3HLAT»8X»4HL0NG,7X»3HDAY,6X»3HDIR,7X«3HSPD,7X» 



X»5HNDATA) 



LATITUDE AND LONGITUDE AT WHICH THE SURFACE CURRENT 



PUTED. 



) »X(1) 



1X,F7.2) 



9.) GO TO 70 



.32) GO TO 60 



N(X(2)/57.296) 



(2.*X(2)/57.296) 



E (NDIM,X»NDERIV.COEF,XMIN»XMAX»NODES* I ERROR) 

 .32) GO TO 70 



£(NDIM,X»NDERIV»COEF»XMIN,XMAX.NODES»IERROR) 



.32) GO TO 70 



Y«vY)**0.5 



TAN2(VX.VY) 



) ANG=ANG+360. 



2) .X ( 1) »CTIMEf ANG,VEL»VX«VY,NDATA 



6.2.5X.F7.2»5X,F4.0*5X,F4.0,3(5X,F5.2).5X,T3) 



ROUTINE SPLCW(NDIM,NDATA*XMIN,XMAX* 



NODES,XTRAP,COEF,NCF,NWRK.IERROR) 

 TION OF SPLCW IS AVAILABLE IN NCAR SOFTWARE SUPPORT 

 , NCAR TECHNICAL NOTE IA-105, MARCH 197S. 



