Changeset 161 for trunk/tools/density_binning/binning_neutral_and_co/bin_velocity_New.pro

Timestamp:

11/12/09 17:15:09 (15 years ago)

Author:

pinsard

Message:

conversion from IDL help header syntax to IDLDOC 2. header syntax

File:

• trunk/tools/density_binning/binning_neutral_and_co/bin_velocity_New.pro (modified) (10 diffs)

trunk/tools/density_binning/binning_neutral_and_co/bin_velocity_New.pro

@@ -13 +13 @@
       xmask = umask()
    ENDIF 
-
    IF keyword_set(vv) THEN BEGIN
       density = boundperio(vmask()*(dens+shift(dens, 0, -1, 0))/2, /vv, /orca2, /pscal)
       xmask = vmask()
    ENDIF 
-
       IF keyword_set(tt) THEN BEGIN
       density = dens
@@ -31 +29 @@
       c1_z = fltarr(N_z)        ; profil du contenu vertical de x
       s_z = fltarr(N_z)         ; profil de la densite
-
 ; vertical coordinate: in z-coord only
       z_zt = gdept              ; profondeur au point T  (k=0 -> 5m)
@@ -73 +70 @@
 ;
          i_ocean = where(xmask(i, j, *) EQ 1)
-
          z_s = z_s*0.
          c1_s = c1_s*0.
          x1_s = x1_s*0.
-
          IF (i_ocean[0] NE -1) THEN BEGIN ; on n'entre que si il y a des points ocean
  
             s_z(*) = density(i, j, *)
             c1_z(*) = x1_content(i, j, *)
-
 ; critere supplementaire a imposer sur le profil pour eviter les cas
 ; pathologiques en attendant d'ecrire une vraie commande d'extraction
 ; de progils stritement croissant. Il s'agit donc d'un test adhoc.
             IF s_z(0) LT s_z(i_ocean(n_elements(i_ocean)-1)) THEN BEGIN 
-
 ;------------------------------------------------------------------------
 ; controle si le profil est bien strictement croissant (pour l'instant
@@ -97 +90 @@
 ;            croissant =  ind(0) EQ -1
 ;------------------------------------------------------------------------
-
             i_bottom = i_ocean(n_elements(i_ocean)-1)
-
             z_s(N_s) = z_zw(i_bottom)
             c1_s(N_s) = x1_content(i, j, jpk-1)
@@ -106 +97 @@
             mini = min(s_z(i_ocean))
             maxi = max(s_z(i_ocean))            
-
             i_min = where(s_z(i_ocean) EQ mini)
             i_max = where(s_z(i_ocean) EQ maxi)
@@ -115 +105 @@
              i_min = i_min[0]
              i_max = i_max(n_elements(i_min)-1)
-
 ;            IF i_max GE jpk-1 THEN print, i, j, i_max
 ;            IF i_min LE 1 THEN print, i, j, i_min
@@ -127 +116 @@
                c1_s(ind) = 0                    &$
             ENDIF 
-
 ; Si la valeur du niveau (s_s) est plus elevee que la densite du fond, 
 ; l isopycne est mise au fond (z_s=z_zw(i_bottom))
@@ -154 +142 @@
                
                    i_profil = i_ocean(i_min:i_max)
-
                    z_s(ind) = interpol(z_zt(i_profil), s_z(i_profil), s_s(ind))
                
@@ -183 +170 @@
          ENDIF 
          ENDIF
-
          x1_bin(i, j, *) = x1_s
          
@@ -190 +176 @@
    
    return, x1_bin
-
 END