Changeset 125 for trunk/SRC/Interpolation/angle.pro

Timestamp:

07/06/06 16:10:25 (18 years ago)

Author:

pinsard

Message:

improvements of Interpolation/*.pro header

File:

• trunk/SRC/Interpolation/angle.pro (modified) (6 diffs)

trunk/SRC/Interpolation/angle.pro

@@ -1 +1 @@
 ;---------
 ;+
-; @file_comments north stereographic polar projection
+; @file_comments 
+; north stereographic polar projection
 ;
 ; @param plam {in}{required}
 ;
-; @param pphi {in}{required} 
+; @param pphi {in}{required}
 ;
 ; @keyword DOUBLE {default=0} use double precision (default is float)
 ;
 ; @returns
-;       gsinu,gcosu : sinus and cosinus of the angle 
-;       gsinv,gcosv   between north-south direction 
+;       gsinu,gcosu : sinus and cosinus of the angle
+;       gsinv,gcosv   between north-south direction
 ;       gsint,gcost   and the j-direction of the mesh
 ;
-; @restrictions to compute the lateral boundary conditions, we assume
-; that: 
+; @restrictions
+; to compute the lateral boundary conditions, we assume that:
 ;     (1) the first line is similar to the second line
-;       =>    gcosu[*, 0] = gcosu[*, 1] 
-;       =>    gsinu[*, 0] = gsinu[*, 1] 
+;       =>    gcosu[*, 0] = gcosu[*, 1]
+;       =>    gsinu[*, 0] = gsinu[*, 1]
 ;     (2) the grid follows OPA x periodicity rule, first column is
 ;     equal to the next to last column
@@ -26 +27 @@
 ;
 ; @history
-;   --------------
 ;       Original :  96-07 (O. Marti)
 ;                   98-06 (G. Madec)
-;       Feb 2005: IDL adaptation S. Masson 
+;       Feb 2005: IDL adaptation S. Masson
 ;
 ; @version $Id$
 ;
 ;-
-;---------
 ;
 FUNCTION fsnspp, plam, pphi, DOUBLE = double
@@ -47 +46 @@
     a = 2. * tan( !pi/4. - !pi/180.*float(pphi)/2. )
     x = cos( !pi/180.*float(plam) ) * a
-    y = sin( !pi/180.*float(plam) ) * a    
+    y = sin( !pi/180.*float(plam) ) * a
   ENDELSE
   RETURN, {x:x, y:y}
@@ -150 +149 @@
   znpv01 = -1; free memory
   IF keyword_set(double) THEN zmnpvt = 1.e-14 > zmnpvt $
-  ELSE zmnpvt = 1.e-6 > zmnpvt 
+  ELSE zmnpvt = 1.e-6 > zmnpvt
 ;   ... j-direction: f-point segment direction (u-point)
   zxffu = znpf00.x - znpf01.x
@@ -157 +156 @@
   znpf01 = -1; free memory
   IF keyword_set(double) THEN zmnpfu = 1.e-14 > zmnpfu $
-  ELSE zmnpfu = 1.e-6 > zmnpfu 
+  ELSE zmnpfu = 1.e-6 > zmnpfu
 ;   ... i-direction: f-point segment direction (v-point)
   zxffv = znpf00.x - znpf10.x
-  zyffv = znpf00.y - znpf10.y 
+  zyffv = znpf00.y - znpf10.y
   znpf00 = -1 &  znpf10 = -1; free memory
   zmnpfv = sqrt ( temporary(znnpv) * ( zxffv*zxffv + zyffv*zyffv )  )
   IF keyword_set(double) THEN zmnpfv = 1.e-14 > zmnpfv $
-  ELSE zmnpfv = 1.e-6 > zmnpfv 
+  ELSE zmnpfv = 1.e-6 > zmnpfv
 ;   ... cosinus and sinus using scalar and vectorial products
   gsint = ( znpt.x*zyvvt - znpt.y*zxvvt ) / zmnpvt
@@ -193 +192 @@
 ; ================================
 ;
-  gcost[*, 0] = gcost[*, 1] 
-  gsint[*, 0] = gsint[*, 1] 
-  gcosu[*, 0] = gcosu[*, 1] 
-  gsinu[*, 0] = gsinu[*, 1] 
+  gcost[*, 0] = gcost[*, 1]
+  gsint[*, 0] = gsint[*, 1]
+  gcosu[*, 0] = gcosu[*, 1]
+  gsinu[*, 0] = gsinu[*, 1]
   gcosv[0, *] = gcosv[jpj-2, *]
   gsinv[0, *] = gsinv[jpj-2, *]