Changeset 48 for trunk/angle.pro

Timestamp:

03/16/14 20:38:39 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules

File:

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

trunk/angle.pro

@@ -16 +16 @@
 ;;      common
 ;;            /comcoh/ glamf,   : longitudes and latitudes at F-points
-;;                     gphif   
+;;                     gphif
 ;;
 ;;   Output :
 ;;   -------
 ;;      common
-;;            /comcoh/ gsinu,gcosu : sinus and cosinus of the angle 
-;;                     gsinv,gcosv   between north-south direction 
-;;                                   and the j-direction of the mesh
+;;            /comcoh/ gsinu,gcosu : sinus and cosinus of the angle
+;;                     gsinv,gcosv   between north-south direction
+;;                                   and the j-direction of the mesh
 ;;
 ;;   Modifications:
@@ -37 +37 @@
 ;;;---------------------------------------------------------------------
 ;
-;
 ; I. Compute the cosinus and sinus
 ; ================================
-;     
+;
 ;     ... north pole direction & modulous (at u-point)
       zxnpu = 0. - fsxspp( glamu, gphiu )
       zynpu = 0. - fsyspp( glamu, gphiu )
       znnpu = zxnpu*zxnpu + zynpu*zynpu
-;     
+;
 ;     ... north pole direction & modulous (at v-point)
       zxnpv = 0. - fsxspp( glamv, gphiv )
       zynpv = 0. - fsyspp( glamv, gphiv )
       znnpv = zxnpv*zxnpv + zynpv*zynpv
-;     
+;
 ;     ... j-direction: f-point segment direction (u-point)
       zxffu=  fsxspp( glamf, gphif ) - fsxspp( shift(glamf, 0, 1), shift(gphif, 0, 1) )
@@ -56 +55 @@
       zmnpfu= sqrt ( znnpu * ( zxffu*zxffu + zyffu*zyffu )  )
 
-     
-;     
+
+;
 ;     ... i-direction: f-point segment direction (v-point)
       zxffv=  fsxspp( glamf, gphif ) - fsxspp( shift(glamf, 1, 0), shift(gphif, 1, 0) )
       zyffv=  fsyspp( glamf, gphif ) - fsyspp( shift(glamf, 1, 0), shift(gphif, 1, 0) )
       zmnpfv= sqrt ( znnpv * ( zxffv*zxffv + zyffv*zyffv )  )
-     
-;     
+
+;
 ;     ... cosinus and sinus using scalar and vectorial products
       gsinu = ( zxnpu*zyffu - zynpu*zxffu ) / zmnpfu
@@ -69 +68 @@
 
 
-;     
+;
 ;     ... cosinus and sinus using scalar and vectorial products
 ;     (caution, rotation of 90 degres)
@@ -75 +74 @@
       gcosv =-( zxnpv*zyffv - zynpv*zxffv ) / zmnpfv
 
-     
-;     
-;     
+
+;
+;
 ;     II. Geographic mesh
 ;     ===================
-;     
-;     
+;
+;
       ind = where(abs(glamf-shift(glamf, 0, 1)) LT 1.e-8)
       gsinu(ind) = 0.d
       gcosu(ind) = 1.d
-;     
+;
       ind = where(abs(gphif-shift(gphif, 1, 0)) LT 1.e-8)
       gsinv(ind) = 0.d
@@ -91 +90 @@
 ;
 return
-END 
+END