Changeset 168 for trunk/tools


Ignore:
Timestamp:
11/13/09 17:04:09 (15 years ago)
Author:
pinsard
Message:

syntax errors fixing

Location:
trunk/tools
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • trunk/tools/density_binning/binning_neutral_and_co/bn2.pro

    r165 r168  
    33; Compute the local Brunt-Vaisala frequency  
    44;        
    5 ;  Method : 
    6 ;  ------- 
     5; Method : 
     6; ------- 
    77;       neos = 0  : UNESCO sea water properties 
    88;          The brunt-vaisala frequency is computed using the polynomial 
    … …  
    2323;      and is never used at this level. 
    2424; 
     25; References : 
     26; ----------- 
     27; McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
     28; 
    2529; @param TN {in}  
    2630; potential temperature 
    … …  
    3135; @returns 
    3236; brunt-vaisala frequency 
    33 ; 
    34 ;   References : 
    35 ;   ----------- 
    36 ;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
    3737; 
    3838; @history 

  • trunk/tools/density_binning/binning_neutral_and_co/eos.pro

    r165 r168  
    22; 
    33; computes rho  (in situ volumic mass)  
    4 ; 
    5 ; @categories 
    6 ; calculation 
    7 ; 
    8 ; @examples  
    9 ;       tableau=eos(t,s,sigma) 
    10 ; 
    11 ; @param T {in}{required}  
    12 ; temperature 
    13 ; 
    14 ; @param S  {in}{required} 
    15 ; salinity 
    16 ; 
    17 ; @param SIGMA {in}{required} 
    18 ; reference level (real number belonging to [0,5]) 
    19 ; for instance, sigma=3. means that rho will be referenced to 3000m   
    20 ; 
    21 ; @uses 
    22 ;       common.pro       
    234; 
    245; REFERENCE: 
    … …  
    4223; 
    4324;       Jackett, D.R., and T.J. McDougall. J. Atmos. Ocean. Tech., 1994 
     25; 
     26; @categories 
     27; calculation 
     28; 
     29; @examples  
     30;       tableau=eos(t,s,sigma) 
     31; 
     32; @param T {in}{required}  
     33; temperature 
     34; 
     35; @param S  {in}{required} 
     36; salinity 
     37; 
     38; @param SIGMA {in}{required} 
     39; reference level (real number belonging to [0,5]) 
     40; for instance, sigma=3. means that rho will be referenced to 3000m   
     41; 
     42; @uses 
     43; common.pro     
    4444; 
    4545; @history 

  • trunk/tools/density_binning/binning_neutral_and_co/eos_neutral.pro

    r164 r168  
    22; 
    33; computes rho  (in situ volumic mass)  
     4; 
     5; REFERENCE: 
     6;       Compute the neutral volumic mass (Kg/m3) from known potential  
     7;       temperature and salinity fields using McDougall and Jackett 2005 
     8;       equation of state. 
     9;              potential temperature         t        deg celsius 
     10;              salinity                      s        psu 
     11;              nutral density                rho      kg/m**3 
     12; 
     13;         Check value: rho(35,20) = 1024.59416751197 kg/m**3  
     14;          t = 20 deg celcius, s=35 psu 
     15; 
     16;       McDougall and Jackett, J. Mar Res., 2005 
    417; 
    518; @categories 
    … …  
    1528; @param S {in} 
    1629; salinity 
    17 ; 
    18 ; REFERENCE: 
    19 ;       Compute the neutral volumic mass (Kg/m3) from known potential  
    20 ;       temperature and salinity fields using McDougall and Jackett 2005 
    21 ;       equation of state. 
    22 ;              potential temperature         t        deg celsius 
    23 ;              salinity                      s        psu 
    24 ;              nutral density                rho      kg/m**3 
    25 ; 
    26 ;         Check value: rho(35,20) = 1024.59416751197 kg/m**3  
    27 ;          t = 20 deg celcius, s=35 psu 
    28 ; 
    29 ;       McDougall and Jackett, J. Mar Res., 2005 
    3030; 
    3131; @history 

  • trunk/tools/density_binning/binning_neutral_and_co/fig_msfs_new.pro

    r166 r168  
    11;+ 
     2; 
     3; @todo 
     4; missing 0_initorca2_NDV1 module 
    25; 
    36; @version 

  • trunk/tools/density_binning/binning_neutral_and_co/fsalbt.pro

    r165 r168  
    44; used for the computation of the the local Brunt-Vaisala frequency  
    55;        
     6;   References : 
     7;   ----------- 
     8;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
     9; 
    610; @param PFT{in} 
    711; potential temperature in degrees celcius 
    … …  
    1519; @returns 
    1620; ratio of thermal over saline expension coefficients 
    17 ; 
    18 ;   References : 
    19 ;   ----------- 
    20 ;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
    2121; 
    2222; @history 

  • trunk/tools/density_binning/binning_neutral_and_co/fsbeta.pro

    r165 r168  
    44; used for the computation of the the local Brunt-Vaisala frequency  
    55;        
     6;   References : 
     7;   ----------- 
     8;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
     9; 
    610; @param PFT {in} {required} 
    711; potential temperature in degrees celcius 
    … …  
    1620; saline expension coefficients 
    1721;  
    18 ;   References : 
    19 ;   ----------- 
    20 ;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 
    21 ; 
    2222; @history 
    2323; Original  : 94-07 (G. Madec, M. Imbard) 

  • trunk/tools/density_binning/binning_neutral_and_co/npc.pro

    r165 r168  
    11;+ 
    2 ;     Non penetrative convective adjustment scheme. solve the static 
    3 ;     instability of the water column. 
    4 ;        
    5 ;  Method : 
    6 ;  ------- 
    7 ;       The algorithm used converges in a maximium of jpk iterations. 
    8 ;      instabilities are treated when the vertical density gradient 
    9 ;      is less than 1.e-5. 
     2; 
     3; Non penetrative convective adjustment scheme. solve the static 
     4; instability of the water column. 
     5; 
     6; Method 
     7; ------ 
     8; The algorithm used converges in a maximium of jpk iterations. 
     9; instabilities are treated when the vertical density gradient 
     10; is less than 1.e-5. 
     11; 
     12; References : 
     13; ----------- 
     14; Madec et al., 1991, JPO, 21, 9, 1349-1371. 
    1015; 
    1116; @param S3D {in} 
    … …  
    1419; @returns 
    1520; adjusted potential density field 
    16 ;               
    17 ;   References : 
    18 ;   ----------- 
    19 ;       Madec et al., 1991, JPO, 21, 9, 1349-1371. 
    20 ; 
     21 
    2122; @history 
    2223;      original : 90-09 (G. Madec) 

  • trunk/tools/density_binning/density_bin_IDL_gm/bin_sigma.pro

    r165 r168  
    55; Method : 
    66; -------- 
    7 ;   a- read density 
    8 ;   b- define density grid 
    9 ;   c- project field       
     7; a- read density 
     8; b- define density grid 
     9; c- project field 
    1010; 
    1111; @param CMD {in} 
    … …  
    3535   vargrid = 'T' 
    3636   sig = make_eos(file_nam, '', time_1 = time1_r, time_2 = time2_r, BOXZOOM = boxzoom, ALL_DATA = all_data) 
    37    IF sig_bowl EQ 1 THEN BEGIN  
     37   IF sig_bowl EQ 1 THEN BEGIN 
    3838      sobwlmax = make_sobwlmax(file_nam, '', time_1 = time1_r, time_2 = time2_r) 
    3939   ENDIF ELSE BEGIN 
    4040      sobwlmax = {data:0} 
    41    ENDELSE  
     41   ENDELSE 
    4242 
    4343   mask = tmask[firstxt:lastxt, firstyt:lastyt, firstzt:lastzt] 
    … …  
    5050; c- project field 
    5151 
    52    CASE pfildi.name OF  
     52   CASE pfildi.name OF 
    5353      '@@vodeptht': data_s = sig.data 
    5454      '@@vosigvol': data_s = sig.data 
    5555      '@@vosigthi': data_s = sig.data 
    5656      ELSE: data_s = pfildi.data 
    57    ENDCASE  
     57   ENDCASE 
    5858   final_dim = pfildi.dim 
    5959   IF strpos(cmd.plt, "_") NE -1 OR strpos(cmd.plt, "#") NE -1 THEN BEGIN 
    … …  
    6161      IF cmd.var EQ '@@vosigvol' THEN BEGIN 
    6262         print, '    (using bining3)...' 
    63          bining3, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, vol_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E1T = e1t, E2T = e2t, E3T = e3t, TMASK = mask  
     63         bining3, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, vol_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E1T = e1t, E2T = e2t, E3T = e3t, TMASK = mask 
    6464      ENDIF ELSE BEGIN 
    6565         print, '    (using bining2)...' 
    6666         bining2, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E3T = e3t, E3W = e3w, TMASK = mask 
    67       ENDELSE  
     67      ENDELSE 
    6868   ENDIF ELSE BEGIN 
    6969      print, '    Performing density bining (assuming global domain)...' 
    7070      IF cmd.var EQ '@@vosigvol' THEN BEGIN 
    7171         print, '    (using bining3)...' 
    72          bining3, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, vol_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E1T = e1t, E2T = e2t, E3T = e3t, TMASK = tmask  
     72         bining3, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, vol_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E1T = e1t, E2T = e2t, E3T = e3t, TMASK = tmask 
    7373      ENDIF ELSE BEGIN 
    7474         print, '    (using bining2)...' 
    7575         bining2, sig.data-1000., data_s, sobwlmax.data, sig_bowl, depth_bin, thick_bin, f_bin, bowl_bin, SIGMA = z_sig, DEPTH_T = gdept, DEPTH_W = gdepw, E3T = e3t, E3W = e3w, TMASK = tmask 
    76       ENDELSE  
    77    ENDELSE  
     76      ENDELSE 
     77   ENDELSE 
    7878 
    7979   ; output : depth_bin = depth of isopycnal surfaces (3D array) 
    … …  
    8484;   print, '    Density domain / intervals : ', sig_min, sig_max, sig_del, n_sig 
    8585 
    86    IF sig_surf EQ 0. THEN BEGIN  
     86   IF sig_surf EQ 0. THEN BEGIN 
    8787      ; vertical section (pltz) 
    8888 
    89       CASE pfildi.name OF  
     89      CASE pfildi.name OF 
    9090         '@@vodeptht': nfild = depth_bin(*, *, 0:n_sig-1) 
    9191         '@@vosigthi': nfild = thick_bin(*, *, 0:n_sig-1) 
    … …  
    9797;            nfild = (thick_bin*e12t3)(*, *, 0:n_sig-1) 
    9898 
    99          END  
     99         END 
    100100         ELSE: nfild = f_bin(*, *, 0:n_sig-1) 
    101       ENDCASE  
     101      ENDCASE 
    102102 
    103    ENDIF ELSE BEGIN  
     103   ENDIF ELSE BEGIN 
    104104      ; horizontal projection on sig_surf (plt) example : votemper@s25 
    105105      ; find index of sig_surf in z_sig 
    … …  
    107107      index = index[n_elements(index)-1] 
    108108 
    109       CASE pfildi.name OF  
     109      CASE pfildi.name OF 
    110110         '@@vodeptht': nfild = depth_bin(*, *, index) 
    111111         '@@vosigthi': nfild = thick_bin(*, *, index) 
    112          '@@vosigvol': BEGIN  
     112         '@@vosigvol': BEGIN 
    113113            nfild = vol_bin 
    114114            final_dim = 1 
    … …  
    116116;            e12t3 = reform((reform(e12t(premierx:dernierx, premiery:derniery), nx*ny))#replicate(1,n_sig +2), nx, ny, n_sig+2) 
    117117;            nfild = ((thick_bin*e12t3)(*, *, index)) 
    118          END  
     118         END 
    119119         ELSE: nfild = f_bin(*, *, index) 
    120       ENDCASE  
     120      ENDCASE 
    121121 
    122    ENDELSE  
     122   ENDELSE 
    123123 
    124    IF final_dim eq 3 THEN BEGIN  
     124   IF final_dim eq 3 THEN BEGIN 
    125125      indx = where(finite(nfild, /NAN)) 
    126126      nfild(indx) = -1.e20 
    … …  
    129129      nfild(indx) = !VALUES.F_NAN 
    130130   ENDIF ELSE BEGIN 
    131   
    132    ENDELSE  
     131 
     132   ENDELSE 
    133133 
    134134   pfild2 = {name: pfildi.name, data: nfild, legend: pfildi.legend, units: pfildi.units, origin: pfildi.origin, dim: final_dim, bowl: bowl_bin} 
    … …  
    136136 
    137137 
    138 END  
     138END 

  • trunk/tools/density_binning/density_bin_IDL_gm/bining2.pro

    r165 r168  
    310310      print, '    Computing spiciness...' 
    311311      FOR i = 0, (size(x1_bin))[3]-1 DO BEGIN  
    312          x1_bin(*,*,i) = x1_bin(*,*,i) - mean(x1_bin(*,*,i),/NaN) 
     312         x1_bin(*,*,i) = x1_bin(*,*,i) - mean(x1_bin(*,*,i),/NAN) 
    313313      ENDFOR 
    314314   ENDIF 
    … …  
    317317   IF 0 THEN BEGIN 
    318318      print, '    Removing domain mean...' 
    319       x1_bin(*,*,*) = x1_bin(*,*,*) - mean(x1_bin(*,*,*),/NaN) 
     319      x1_bin(*,*,*) = x1_bin(*,*,*) - mean(x1_bin(*,*,*),/NAN) 
    320320   ENDIF 
    321321 

  • trunk/tools/density_binning/density_bin_IDL_gm/eos.pro

    r165 r168  
    22; 
    33; computes rho  (in situ volumic mass)  
    4 ; 
    5 ; @categories 
    6 ;       calculation 
    7 ; 
    8 ; @examples  
    9 ;       tableau=eos(t,s,sigma) 
    10 ; 
    11 ; @param T {in}{required} 
    12 ; temperature 
    13 ; 
    14 ; @param S  {in}{required} 
    15 ; salinity 
    16 ; 
    17 ; @param SIGMA {in}{required} 
    18 ; reference level (real number belonging to [0,5]) 
    19 ; for instance, sigma=3. means that rho will be referenced to 3000m   
    20 ; 
    21 ; @uses 
    22 ;       common.pro       
    234; 
    245; REFERENCE: 
    … …  
    4223; 
    4324;       Jackett, D.R., and T.J. McDougall. J. Atmos. Ocean. Tech., 1994 
     25; 
     26; @categories 
     27;       calculation 
     28; 
     29; @examples  
     30;       tableau=eos(t,s,sigma) 
     31; 
     32; @param T {in}{required} 
     33; temperature 
     34; 
     35; @param S  {in}{required} 
     36; salinity 
     37; 
     38; @param SIGMA {in}{required} 
     39; reference level (real number belonging to [0,5]) 
     40; for instance, sigma=3. means that rho will be referenced to 3000m   
     41; 
     42; @uses 
     43;       common.pro       
    4444; 
    4545; @history 

  • trunk/tools/density_binning/density_bin_IDL_gm/npc.pro

    r165 r168  
    44; instability of the water column. 
    55;        
    6 ;  Method : 
    7 ;  ------- 
    8 ;       The algorithm used converges in a maximium of jpk iterations. 
    9 ;      instabilities are treated when the vertical density gradient 
    10 ;      is less than 1.e-5. 
     6; Method : 
     7; ------- 
     8; The algorithm used converges in a maximium of jpk iterations. 
     9; instabilities are treated when the vertical density gradient 
     10; is less than 1.e-5. 
     11; 
     12; References : 
     13; ----------- 
     14; Madec et al., 1991, JPO, 21, 9, 1349-1371. 
    1115; 
    1216; @param S3D {in} 
    … …  
    1620; adjusted potential density field 
    1721;               
    18 ;   References : 
    19 ;   ----------- 
    20 ;       Madec et al., 1991, JPO, 21, 9, 1349-1371. 
    21 ; 
    2222; @history 
    2323;      original : 90-09 (G. Madec) 

  • trunk/tools/lma/slec.pro

    r165 r168  
    1414   if strlowcase(name) EQ 'vn' then name = 'AM_win1r' 
    1515;------------------------------------------------------------  
    16    return, read_data(name,debut,fin, PARENT = parent, BOITE=boite, _EXTRA=extra) 
     16   return, read_data(name,debut,fin, PARENT=parent, BOITE=boite, _EXTRA=extra) 
    1717; 
    1818end 

  • trunk/tools/transfo/grads_read.pro

    r165 r168  
    1212FUNCTION grads_read, file_name, var_name $ 
    1313         , NCDF_DB=ncdf_db $ 
    14          , TIME_1=time_1$ 
    15          , TIME_2=time_2$ 
    16          , ALL_DATA=all_data$ 
    17          , NO_KEY_SHIFT=no_key_shift$ 
     14         , TIME_1=time_1 $ 
     15         , TIME_2=time_2 $ 
     16         , ALL_DATA=all_data $ 
     17         , NO_KEY_SHIFT=no_key_shift $ 
    1818         , LEVEL=level 
    1919; 
Note: See TracChangeset for help on using the changeset viewer.