New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3999 – NEMO

Changeset 3999


Ignore:
Timestamp:
2013-08-06T13:39:29+02:00 (11 years ago)
Author:
flavoni
Message:

update idl scripts for ORCA2 configurations, see ticket #1115

Location:
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
6 added
21 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.pro

    r2751 r3999  
    2626  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    2727;  
    28   domdef, 20, 380, 60, 90 
     28  domdef, 20, 380, 50, 90 
    2929; 
    3030  varunit = S1.unit 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_GlobMeanTS.pro

    r2751 r3999  
    2525  title = 'Salinity Global mean!C'+std_file1_T+' - Levitus (Black)' 
    2626  if std_file1_T NE std_file2_T THEN title = title+'!C'+std_file2_T+' - Levitus (Red)' 
    27   plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = .2, boxzoom = [4., 5300.], /KEEPBOTTOM $ 
     27  plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = .4, boxzoom = [4., 5300.], /KEEPBOTTOM $ 
    2828         , small = [1, 2, 2], XGRIDSTYLE = 2, YGRIDSTYLE = 2, TITLE = title, /NOERASE, _extra = ex 
    2929  if std_file1_T NE std_file2_T then begin 
    30     plt1d, S2.arr - SLev.arr, typein = 'z',  ticklen = 1, MIN = -2., MAX = 2., boxzoom = [4., 5300.], /KEEPBOTTOM $ 
     30    plt1d, S2.arr - SLev.arr, typein = 'z',  ticklen = 1, MIN = -2., MAX = 4., boxzoom = [4., 5300.], /KEEPBOTTOM $ 
    3131           ,  /ov1d, COLOR = 250, TITLE = title, _extra = ex  
    3232  endif 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r3927 r3999  
    4040         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
    4141         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
    42  endif else begin 
    43     title = title+ ' - Observations' 
    44     plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $ 
    45          , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $ 
    46          , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
    47   endelse   
     42;;SF  finire, aggiungere lettura obs 
     43;;SF  endif else begin 
     44;;SF    title = title+ ' - Observations' 
     45;;SF    plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $ 
     46;;SF         , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $ 
     47;;SF         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
     48;;SF  endelse   
     49endif 
    4850 
    4951  domdef 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro

    r3927 r3999  
    2222    domdef, 20, 380, 50, 90 
    2323    vmin = 0. 
    24     vmax = 10. 
     24    vmax = 6. 
    2525    vint = 0.5 
    2626    fmt = '(i2)' 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro

    r3927 r3999  
    5555  if std_file1_I NE std_file2_I then begin             
    5656    title = title + std_file2_I 
    57     plt, Isnow1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
     57    plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    5858         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
    5959         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r3927 r3999  
    199199  ENDELSE 
    200200; 
    201   Isnow_1 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file1_I )  
    202   help, Isnow_1.arr, jpi,jpj 
    203   caldat, time, mm 
    204   april = where(mm EQ 4, cnt) 
    205   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    206   sept = where(mm EQ 9, cnt) 
    207   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
    208   ;undefine, Isnow_1 
    209 ; 
    210   IF std_file2_I NE std_file1_I THEN BEGIN 
    211     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    212     caldat, time, mm 
    213     april = where(mm EQ 4, cnt) 
    214     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    215     sept = where(mm EQ 9, cnt) 
    216     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    217   ;  undefine, Isnow_2 
    218   ENDIF ELSE BEGIN  
    219     Isnow_april_2 = {arr:-1} 
    220     Isnow_sept_2 = {arr:-1} 
    221   ENDELSE 
     201;; Sf commentato perché non funziona snow 
     202;; Sf   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
     203;; Sf   help, Isnow_1.arr, jpi,jpj 
     204;; Sf   caldat, time, mm 
     205;; Sf   april = where(mm EQ 4, cnt) 
     206;; Sf   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     207;; Sf   sept = where(mm EQ 9, cnt) 
     208;; Sf   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
     209;; Sf   ;undefine, Isnow_1 
     210;; Sf ; 
     211;; Sf   IF std_file2_I NE std_file1_I THEN BEGIN 
     212;; Sf     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
     213;; Sf     caldat, time, mm 
     214;; Sf     april = where(mm EQ 4, cnt) 
     215;; Sf     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     216;; Sf     sept = where(mm EQ 9, cnt) 
     217;; Sf     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     218;; Sf   ;  undefine, Isnow_2 
     219;; Sf   ENDIF ELSE BEGIN  
     220;; Sf     Isnow_april_2 = {arr:-1} 
     221;; Sf     Isnow_sept_2 = {arr:-1} 
     222;; Sf   ENDELSE 
    222223; 
    223224  jpt = 1 
     
    325326  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex  
    326327   
    327 ;;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    328 ;;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
     328;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
     329;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    329330 
    330331  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
     
    334335  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    335336   
    336 ;;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    337 ;;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
     337;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
     338;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    338339 
    339340  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
     
    354355  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    355356  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    356  
    357   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
    358   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
    359  
    360   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    361   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    362   
    363   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
    364   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
    365   
    366   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    367   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    368   
    369   cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
    370   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
    371    
    372   cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
    373   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
     357;; SF non funziona, dice che deve essere una struttura e non e!!!! 
     358 cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
     359 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
     360;; SF  
     361;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
     362;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     363;; SF  
     364;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
     365;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
     366;; SF  
     367;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
     368;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     369;  
     370;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
     371;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
     372   
     373;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
     374;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
    374375 
    375376  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_bsf.pro

    r2751 r3999  
    1414  domdef, 0, 6000 
    1515 
    16   bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0]) 
     16  bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0], /transport ) 
    1717  title = 'Barotropic Stream Function!C'+std_file1_U 
    1818  plt, bb1, min = -200., max = 200., int = 10., /portrait, FORMAT = '(I4)', STYLE = 'so0so' $ 
     
    2020 
    2121  if std_file1_U NE std_file2_U then BEGIN 
    22     bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0]) 
     22    bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0], /transport ) 
    2323    title = 'Barotropic Stream Function!C'+std_file1_U+' - '+std_file2_U 
    2424    plt, bb1.arr - bb2.arr, min = -20., max = 20., int = 2.,  FORMAT = '(I2)', STYLE = 'so0so' $ 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mld.pro

    r2751 r3999  
    3232  ENDIF ELSE BEGIN  
    3333    title = titleorg+std_file1_T 
    34     plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 
    35          , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 
     34    plt, MLD1,label=4, cb_label=[0, 15, 25, 50, 75, 100, 125, 150, 250, 400, 600, 850],/NOCONTOUR,FORMAT = '(I3)' $ 
     35            , small = [1, 2, 1],COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 
     36;ORI SF    plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 
     37;ORI SF         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 
    3638  ENDELSE  
    3739;  
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro

    r2760 r3999  
    3333  domdef, 0, 6000 
    3434 
    35   mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma) 
     35  mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma, /transport ) 
    3636  title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V 
    3737  pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $ 
     
    3939 
    4040  if std_file1_V NE std_file2_V then begin 
    41     mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma) 
     41    mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma, /transport) 
    4242    title = title+' - '+std_file2_V 
    4343    pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $ 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r3927 r3999  
    1616#===================== User PATHS ===================== 
    1717# 
    18 idl_command=/Applications/itt/idl64/bin/idl 
     18#idl_command=/Applications/itt/idl64/bin/idl 
     19idl_command=/Applications/itt/idl71/bin/idl 
     20#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl 
     21#idl_command=idl71 
    1922# 
    20 PS_DIR=$( pwd )/ps_plot 
    21 PDF_DIR=$( pwd )/pdf_plot 
     23PS_DIR=$( pwd )/500yfwb0_ps_plot 
     24PDF_DIR=$( pwd )/500yfwb0_pdf_plot 
    2225HTML_DIR=$( pwd )/html_plot 
    2326SAXO_DIR=/Users/sflod/SAXO_DIR 
    2427# 
    25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2629DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    2730DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     
    3437#===================== DATA ===================== 
    3538# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     39#VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
     40#VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     41VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     42VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3843VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3944VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
     
    5257#===================== EXP1 ===================== 
    5358# 
    54 FILE1_T=LIM2_00010101_00011231_1m_grid_T.nc   # exp1 grid T input file 
    55 FILE1_U=LIM2_00010101_00011231_1m_grid_U.nc   # exp1 grid U input file 
    56 FILE1_V=LIM2_00010101_00011231_1m_grid_V.nc   # exp1 grid V input file 
    57 FILE1_I=LIM2_00010101_00011231_1m_icemod.nc   # exp1 ice    input file 
    58 VAR1_T=votemper 
    59 VAR1_S=vosaline 
    60 VAR1_QNET=sohefldo 
    61 VAR1_ERP=sowaflcd 
    62 VAR1_EMP=sowaflup 
     59#FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
     60#FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
     61#FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
     62#FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
     63FILE1_T=500yfwb0_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file 
     64FILE1_U=500yfwb0_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file 
     65FILE1_V=500yfwb0_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file 
     66FILE1_I=500yfwb0_20010101_21001231_1M_icemod.nc   # exp1 ice    input file 
     67VAR1_T=votemper          
     68VAR1_S=vosaline          
     69VAR1_QNET=sohflxdo 
     70VAR1_ERP=sowafldp       
     71VAR1_EMP=sowaflup        
    6372VAR1_MLD=somxl010 
    64 VAR1_U=uocetr_eff 
    65 VAR1_V=vocetr_eff 
    66 VAR1_Ithick=iicethic 
     73VAR1_U=vozoeftr 
     74#VAR1_U=uoce 
     75#VAR1_U=vozocrtx 
     76VAR1_V=vomeeftr 
     77#VAR1_V=voce 
     78#VAR1_V=vomecrty 
     79VAR1_Ithick=iicethic     
    6780VAR1_Ifrac=iiceconc 
    68 VAR1_Isnow=isnowthi 
     81VAR1_SNOW=isnowthi 
    6982# 
    7083#===================== EXP2 ===================== 
    7184# 
    72 FILE2_T=LIM2_00010101_00011231_1m_grid_T.nc   # exp1 grid T input file 
    73 FILE2_U=LIM2_00010101_00011231_1m_grid_U.nc   # exp1 grid U input file 
    74 FILE2_V=LIM2_00010101_00011231_1m_grid_V.nc   # exp1 grid V input file 
    75 FILE2_I=LIM2_00010101_00011231_1m_icemod.nc 
     85#FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
     86#FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
     87#FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
     88#FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
     89FILE2_T=500yfwb0_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file 
     90FILE2_U=500yfwb0_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file 
     91FILE2_V=500yfwb0_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file 
     92FILE2_I=500yfwb0_20010101_21001231_1M_icemod.nc   # exp1 ice    input file 
    7693VAR2_T=votemper 
    7794VAR2_S=vosaline 
    78 VAR2_QNET=sohefldo 
    79 VAR2_ERP=sowaflcd 
     95VAR2_QNET=sohflxdo 
     96VAR2_ERP=sowafldp 
    8097VAR2_EMP=sowaflup 
    8198VAR2_MLD=somxl010 
    82 VAR2_U=uocetr_eff 
    83 VAR2_V=vocetr_eff 
     99VAR2_U=vozoeftr 
     100#VAR2_U=uoce 
     101#VAR2_U=vozocrtx 
     102VAR2_V=vomeeftr 
     103#VAR2_V=voce 
     104#VAR2_V=vomecrty 
    84105VAR2_Ithick=iicethic 
    85106VAR2_Ifrac=iiceconc 
    86 VAR2_Isnow=isnowthi 
     107VAR2_SNOW=isnowthi 
    87108# 
    88109######################### Export Variables ############################### 
     
    100121export FILE_FLUX      VAR_FLUX 
    101122export FILE_MLD       VAR_MLD 
    102 export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH, VAR_ICE_area_NH, VAR_ICE_area_SH 
     123export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 
    103124export FILE_SNOW_ARC  VAR_SNOW_NH 
    104125export FILE_SNOW_ANT  VAR_SNOW_SH 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro

    r2751 r3999  
    2727  ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce )   &   ts_EMP.unit = 'Sv' 
    2828  title = prefix+' '+d1_d2+'!C'+blabla 
    29   pltt, ts_EMP, 't', -.001, .001, date1, date2, /REMPLI $ 
     29  pltt, ts_EMP, 't', -2., 2., date1, date2, /REMPLI $ 
    3030        , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 
    3131   
     
    3939 
    4040    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    41     pltt, ts_EMP.arr - ts_EMP2.arr, 't', -.001, .001, date1, date2, /REMPLI $ 
     41    pltt, ts_EMP.arr - ts_EMP2.arr, 't', -2., 2., date1, date2, /REMPLI $ 
    4242          , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 
    4343 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r3927 r3999  
    8888  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
    8989  ; 
    90   title = 'Northern Hemisphere !C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 
     90  title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 
    9191  jpt=12 
    9292  time=julday(1,15,1900)+30*lindgen(12) 
     
    127127    ICE_N2 = {arr:ICE_N2 * 1.e-12, unit : '10^12 m^2'} 
    128128    ;ICE Extent minus 15% in NORTH Hemisphere 
    129     ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     129    ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
    130130    msk = ICE_N2_15.arr gt 0.15  ; remove 0.15% for observations 
    131131    ICE_N2_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 
     
    137137    ;ICE Extent in SOUTH Hemisphere 
    138138    domdef, 0, jpi-1, -90, -30, /xindex  
    139     ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     139    ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
    140140    ICE_S2 = grossemoyenne(ICE_S2.arr, 'xy', /integration, mask2d = masknp) 
    141141    if jpt mod 12 ne 0 then stop 
     
    145145    ICE_S2 = {arr:ICE_S2 * 1.e-12, unit : '10^12 m^2'} 
    146146    ;ICE Extent minus 15% in SOUTH Hemisphere 
    147     ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     147    ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
    148148    msk = ICE_S2_15.arr gt 0.15  ; remove 0.15% for observations 
    149149    ICE_S2_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 
     
    154154    ICE_S2_15 = {arr:ICE_S2_15 * 1.e-12, unit : '10^12 m^2'} 
    155155  ; 
    156     time = tsave   &   IF n_elements(time) NE jpt THEN stop 
     156 ;   time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    157157 
    158     if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     158    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    159159 
    160     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area  (Red SOLID simulation)  and Extend minus 15% (Blue SOLID simulation)' 
    161     pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 
    162           , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    163     pltt, ICE_N_15.arr -  ICE_N2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    164       , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    165      
    166     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.) (Red SOLID simulation)  and Extend minus 15% (Blue SOLID simulation)' 
    167     pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 
    168           , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    169     pltt, ICE_S_15.arr -  ICE_S2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    170       , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    171      
    172160 
    173     htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     161  title = 'Northern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+'!C'+' Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 
     162  jpt=12 
     163  time=julday(1,15,1900)+30*lindgen(12) 
     164  pltt, ICE_N.arr - ICE_N2.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 
     165        , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     166  pltt, ICE_N_15.arr -  ICE_N2_15.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     167        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     168;SF  pltt, vICE_area_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     169;SF        , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     170;SF  pltt, vICE_ext_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     171;SF        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     172; 
     173;  title ='Southern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+' - '+'!C'+'Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 
     174  title ='Southern Hemisphere'+'!C' 
     175  pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE $ 
     176         , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     177  pltt, ICE_S_15.arr -  ICE_S2_15.arr , 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     178        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     179;SF  pltt,  vICE_area_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     180;SF         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     181;SF  pltt,  vICE_ext_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     182;SF        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     183; 
     184 
     185    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
    174186    if KEYWORD_SET(postscript) then closeps 
    175187 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

    r3929 r3999  
    2626  v2_Ithick = getenv('VAR2_Ithick')  &   prefix2 = getenv('V2It_PREF')    &   suffix2 = getenv('V2It_SUFF') 
    2727 
    28 ; get ice VOLUME climatology info 
    29 ;; 
    30 ;  std_file_ice =  isafile(getenv('FILE_ICE'), title = 'ICE Extent Climatology', iodir = std_iodir_climato) 
    31 ; 
    32 ;  time_ice = ncdf_lec( std_file_ice, VAR='time' ) 
    33 ;  time_ice = (time_ice - FLOOR(time_ice) ) * 12 
    34 ;  time_ice = round(time_ice) ; round to nearest integer  
    35 ;  t1 = where(time_ice eq 0) 
    36 ;  t1 = t1[0] ;  jannuary 
    37 ;  t2 = where(time_ice eq 11, count) 
    38 ;  t2 = t2[count-1] ; last day of december 
    39 ;  nyear = (t2-t1+1)/12 
    40 ;  vICE_Ithick_NH = read_arr2d(std_file_ice, getenv('VAR1_Ithick'), t1, t2 ) 
    41 ;  vICE_Ithick_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_SH'), t1, t2 ) 
    42 ; 
    43 ;  vICE_area_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_NH'), t1, t2 ) 
    44 ;  vICE_area_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_SH'), t1, t2 ) 
    45  
    4628  cdti3 = string(cnt, format = '(i3.3)') 
    4729  print, cdti3 + ') ' + blabla 
    48   filename = cdti3 + '_ts_ICE_'+prefix 
     30  filename = cdti3 + '_ts_ICE_Vol'+prefix 
    4931  if prefix NE prefix2 then filename = filename + '_'+prefix2 
    5032  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     33; 
    5134  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
     35; 
    5236  iodir = std_iodir_data 
    5337  ; ICE Area(=Surface) in NORTH Hemisphere 
     
    5842  ICE_vol_N = (ICE_N.arr < 1.e10 ) * ( ICE_thick.arr < 1.e10)  ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl 
    5943  ICE_vol_N = grossemoyenne(ICE_vol_N, 'xy', /integration, mask2d = masknp) 
     44  ; 
    6045  if jpt mod 12 ne 0 then stop 
    6146  nyr=jpt/12. 
     
    6348  ICE_vol_N = total(ICE_vol_N,2)/nyr 
    6449  ICE_vol_N = {arr:ICE_vol_N * 1.e-9, unit : '10^9 Km^3'} 
     50  ; 
    6551  ;ICE Area(=Surface) in SOUTH Hemisphere 
    6652  domdef, 0, jpi-1, -90, -30, /xindex  
     
    7662  ICE_vol_S = {arr:ICE_vol_S * 1.e-9, unit : '10^9 Km^3'} 
    7763  ; 
    78   title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation) + Observation (dashed)' 
     64  title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
    7965  jpt=12 
    8066  time=julday(1,15,1900)+30*lindgen(12) 
    8167  pltt, ICE_vol_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 30000. $ 
    82         , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
     68        , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
    8369; 
    84   title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation) + Observation (dashed)' 
     70  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
    8571  pltt, ICE_vol_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 11000. $ 
    86         , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
     72        , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
    8773; 
    88  
    8974  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
    9075  if KEYWORD_SET(postscript) then closeps 
     
    9681    domdef, 0, jpi-1, 30, 90, /xindex 
    9782    ;ICE Area(=Surface) in NORTH Hemisphere 
    98     ICE_N2 = rseries_ncdf(vICE2, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
    99     ICE_thick2 = rseries_ncdf(v2_Ithick, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     83    ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     84    ICE_thick2 = rseries_ncdf(v2_Ithick, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
    10085    ; Volume = Area(=Surface) * Thickness 
    10186    ICE_vol_N2 = (ICE_N2.arr < 1.e10 ) * ( ICE_thick2.arr < 1.e10) ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl 
     
    10994    ;ICE Area(=Surface) in SOUTH Hemisphere 
    11095    domdef, 0, jpi-1, -90, -30, /xindex  
    111     ICE_S2 = rseries_ncdf(vICE2, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
    112     ICE_thick2 = rseries_ncdf(v2_Ithick, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     96    ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     97    ICE_thick2 = rseries_ncdf(v2_Ithick, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
    11398    ; Volume = Area(=Surface) * Thickness 
    11499    ICE_vol_S2 = (ICE_S2.arr < 1.e10 ) * ( ICE_thick.arr < 1.e10) ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl 
     
    120105    ICE_vol_S2 = {arr:ICE_vol_S2 * 1.e-9, unit : '10^3 Km^3'} 
    121106 
    122     time = tsave   &   IF n_elements(time) NE jpt THEN stop 
     107   ; time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    123108 
    124     if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     109    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    125110 
    126     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume (North. Hemisp.)' 
    127     pltt, ICE_vol_N.arr - ICE_vol_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 
    128           , COLOR = 250, small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    129    ; 
    130     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 
    131     pltt, ICE_vol_S.arr - ICE_vol_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 
    132           , COLOR = 250, small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex     
     111    title = 'Northern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 
     112    jpt=12 
     113    time=julday(1,15,1900)+30*lindgen(12) 
     114    pltt, ICE_vol_N.arr - ICE_vol_N2.arr, 't', -125., -60., 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 
     115          , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
     116    ; 
     117    title = 'Southern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 
     118    pltt, ICE_vol_S.arr - ICE_vol_S2.arr, 't', -620., -10., 19000101, 19001231, /REMPLI, /NOERASE $ 
     119          , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex    
    133120 
    134     htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     121    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
    135122    if KEYWORD_SET(postscript) then closeps 
    136123 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro

    r2751 r3999  
    2424  ts_SSH.arr = ts_SSH.arr * 100   &   ts_SSH.unit = 'cm' 
    2525  title = prefix+' '+d1_d2+'!C'+blabla 
    26   pltt, ts_SSH, 't', -.1, .1, date1, date2, /REMPLI $ 
     26  pltt, ts_SSH, 't', -220., 220., date1, date2, /REMPLI $ 
    2727        , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 
    2828   
     
    3636 
    3737    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    38     pltt, ts_SSH.arr - ts_SSH2.arr, 't', -.1, .1, date1, date2, /REMPLI $ 
     38    pltt, ts_SSH.arr - ts_SSH2.arr, 't', -10., 10., date1, date2, /REMPLI $ 
    3939          , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 
    4040 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro

    r3610 r3999  
    4747 
    4848  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 
    49   pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $  
     49  pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
    5050        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 
    5151   
    5252  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 
    53   pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $  
     53  pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
    5454        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    5555   
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r3927 r3999  
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    7979 
    80 ; cnt = cnt+1   &   blabla = 'SNOW Volume' 
    81 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
     80;SF cnt = cnt+1   &   blabla = 'SNOW Volume' 
     81;SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    8282; 
    8383  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

    r3927 r3999  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ps_lim3ada_lim3ada 
    21 PDF_DIR=$( pwd )/pdf_lim3ada_lim3ada 
     20PS_DIR=$( pwd )/ts_500yfwb0_ps 
     21PDF_DIR=$( pwd )/ts_500yfwb0_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
     
    3434#===================== DATA ===================== 
    3535# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     36VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     37VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3838VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3939VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    4040VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
    41 VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent North Emisphere 
    42 VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent South Emisphere 
    43 VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area North Emisphere 
    44 VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area South Emisphere 
     41VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     42VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     43VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     44VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
    4545# 
    4646# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    5252READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    5353# 
    54 DATE1=20010101   ;   DATE2=20101231 
     54DATE1=20010101   ;   DATE2=21001231 
    5555#DATE1=00010101      ;   DATE2=00101231 
    5656# 
    57 VAR1_T=thetao           ;   V1T_PREF=lim3ada         ;   V1T_SUFF=_1Y_grid_T.nc 
    58 VAR1_S=so               ;   V1S_PREF=lim3ada         ;   V1S_SUFF=_1Y_grid_T.nc 
    59 VAR1_SSH=zos            ;   V1SSH_PREF=lim3ada       ;   V1SSH_SUFF=_1Y_grid_T.nc 
    60 VAR1_Q=qt               ;   V1Q_PREF=lim3ada         ;   V1Q_SUFF=_1Y_grid_T.nc 
    61 VAR1_EMP=wfo            ;   V1EMP_PREF=lim3ada       ;   V1EMP_SUFF=_1Y_grid_T.nc 
    62 VAR1_U=uocetr_eff       ;   V1U_PREF=lim3ada         ;   V1U_SUFF=_1Y_grid_U.nc 
    63 VAR1_V=vocetr_eff       ;   V1V_PREF=lim3ada         ;   V1V_SUFF=_1Y_grid_V.nc 
    64 VAR1_ICE=iiceconc       ;   V1ICE_PREF=lim3ada       ;   V1ICE_SUFF=_1M_icemod.nc 
    65 VAR1_Ithick=iicethic    ;   V1It_PREF=lim3ada        ;   V1It_SUFF=_1M_icemod.nc 
    66 VAR1_SNOW=isnowthi      ;   V1SNOW_PREF=lim3ada      ;   V1SNOW_SUFF=_1M_icemod.nc 
     57VAR1_T=votemper         ;   V1T_PREF=500yfwb0      ;   V1T_SUFF=_1Y_grid_T.nc 
     58VAR1_S=vosaline         ;   V1S_PREF=500yfwb0      ;   V1S_SUFF=_1Y_grid_T.nc 
     59VAR1_SSH=sossheig       ;   V1SSH_PREF=500yfwb0    ;   V1SSH_SUFF=_1Y_grid_T.nc 
     60VAR1_Q=sohflxdo         ;   V1Q_PREF=500yfwb0      ;   V1Q_SUFF=_1Y_grid_T.nc 
     61VAR1_EMP=sowaflup       ;   V1EMP_PREF=500yfwb0    ;   V1EMP_SUFF=_1Y_grid_T.nc 
     62VAR1_U=vozocrtx         ;   V1U_PREF=500yfwb0      ;   V1U_SUFF=_1Y_grid_U.nc 
     63VAR1_V=vomecrty         ;   V1V_PREF=500yfwb0      ;   V1V_SUFF=_1Y_grid_V.nc 
     64VAR1_ICE=iiceconc       ;   V1ICE_PREF=500yfwb0    ;   V1ICE_SUFF=_1M_icemod.nc 
     65VAR1_Ithick=iicethic    ;   V1It_PREF=500yfwb0     ;   V1It_SUFF=_1M_icemod.nc 
     66VAR1_SNOW=isnowthi      ;   V1SNOW_PREF=500yfwb0   ;   V1SNOW_SUFF=_1M_icemod.nc 
    6767# 
    6868#===================== EXP2 ===================== 
    6969# 
    70 DATE1_2=20010101   ;   DATE2_2=20101231 
     70#DATE1_2=20010101   ;   DATE2_2=20101231 
     71DATE1_2=20010101   ;   DATE2_2=21001231 
    7172# 
    72 VAR2_T=thetao           ;   V2T_PREF=lim3ada         ;   V2T_SUFF=_1Y_grid_T.nc 
    73 VAR2_S=so               ;   V2S_PREF=lim3ada         ;   V2S_SUFF=_1Y_grid_T.nc 
    74 VAR2_SSH=zos            ;   V2SSH_PREF=lim3ada       ;   V2SSH_SUFF=_1Y_grid_T.nc 
    75 VAR2_Q=qt               ;   V2Q_PREF=lim3ada         ;   V2Q_SUFF=_1Y_grid_T.nc 
    76 VAR2_EMP=wfo            ;   V2EMP_PREF=lim3ada       ;   V2EMP_SUFF=_1Y_grid_T.nc 
    77 VAR2_U=uocetr_eff       ;   V2U_PREF=lim3ada         ;   V2U_SUFF=_1Y_grid_U.nc 
    78 VAR2_V=vocetr_eff       ;   V2V_PREF=lim3ada         ;   V2V_SUFF=_1Y_grid_V.nc 
    79 VAR2_ICE=iiceconc       ;   V2ICE_PREF=lim3ada       ;   V2ICE_SUFF=_1M_icemod.nc 
    80 VAR2_Ithick=iicethic    ;   V2It_PREF=lim3ada        ;   V2It_SUFF=_1M_icemod.nc 
    81 VAR2_SNOW=isnowthi      ;   V2SNOW_PREF=lim3ada      ;   V2SNOW_SUFF=_1M_icemod.nc 
     73VAR2_T=votemper         ;   V2T_PREF=500yfwb0      ;   V2T_SUFF=_1Y_grid_T.nc 
     74VAR2_S=vosaline         ;   V2S_PREF=500yfwb0      ;   V2S_SUFF=_1Y_grid_T.nc 
     75VAR2_SSH=sossheig       ;   V2SSH_PREF=500yfwb0    ;   V2SSH_SUFF=_1Y_grid_T.nc 
     76VAR2_Q=sohflxdo         ;   V2Q_PREF=500yfwb0      ;   V2Q_SUFF=_1Y_grid_T.nc 
     77VAR2_EMP=sowaflup       ;   V2EMP_PREF=500yfwb0    ;   V2EMP_SUFF=_1Y_grid_T.nc 
     78VAR2_U=vozocrtx         ;   V2U_PREF=500yfwb0      ;   V2U_SUFF=_1Y_grid_U.nc 
     79VAR2_V=vomecrty         ;   V2V_PREF=500yfwb0      ;   V2V_SUFF=_1Y_grid_V.nc 
     80VAR2_ICE=iiceconc       ;   V2ICE_PREF=500yfwb0    ;   V2ICE_SUFF=_1M_icemod.nc 
     81VAR2_Ithick=iicethic    ;   V2It_PREF=500yfwb0     ;   V2It_SUFF=_1M_icemod.nc 
     82VAR2_SNOW=isnowthi      ;   V2SNOW_PREF=500yfwb0   ;   V2SNOW_SUFF=_1M_icemod.nc 
    8283# 
    8384######################### Export Variables ############################### 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example1

    r3928 r3999  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ps_core_tmx 
    21 PDF_DIR=$( pwd )/pdf_core_tmx 
    22 HTML_DIR=$( pwd )/html_core_tmx 
    23 SAXO_DIR=/Users/smasson/SAXO_DIR 
     20PS_DIR=$( pwd )/ps_tag33_tag331 
     21PDF_DIR=$( pwd )/pdf_tag33_tag331 
     22HTML_DIR=$( pwd )/html 
     23SAXO_DIR=/Users/sflod/SAXO_DIR 
    2424# 
    25 DIR_DATA=/Users/smasson/tmp/idlorca2/res      # path of data in NetCDF format 
    26 DIR_CLIMATO=/Users/smasson/tmp/idlorca2/clim  # path of climatological data 
    27 DIR_MASK=/Users/smasson/tmp/idlorca2/clim     # path of mask files (ex: subbasins) 
     25DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     26DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     27DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    2828# 
    2929#===================== Model GRID ===================== 
    3030# 
    31 FILE_MESH_MASK=/Users/smasson/data/vp33beta_mesh_mask_7.nc # meshmask 
     31FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 
    3232FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks 
    3333# 
    3434#===================== DATA ===================== 
    3535# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     36VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     37VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3838VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3939VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    4040VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
    41 VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent North Emisphere 
    42 VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent South Emisphere 
    43 VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area North Emisphere 
    44 VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area South Emisphere 
     41VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     42VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     43VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     44VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
    4545# 
    4646# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    5252READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    5353# 
    54 DATE1=10101   ;   DATE2=20001231 
     54DATE1=20010101   ;   DATE2=21001231 
    5555# 
    5656VAR1_T=thetao      ;   V1T_PREF=tmx_2000     ;   V1T_SUFF=_1Y_grid_T.nc 
     
    100100# 
    101101export DATE1        DATE2 
    102 export VAR1_T       V1T_PREF    V1T_SUFF 
    103 export VAR1_S       V1S_PREF    V1S_SUFF 
    104 export VAR1_SSH     V1SSH_PREF  V1SSH_SUFF 
    105 export VAR1_Q       V1Q_PREF    V1Q_SUFF 
    106 export VAR1_EMP     V1EMP_PREF  V1EMP_SUFF 
    107 export VAR1_U       V1U_PREF    V1U_SUFF 
    108 export VAR1_V       V1V_PREF    V1V_SUFF 
    109 export VAR1_ICE     V1ICE_PREF  V1ICE_SUFF 
    110 export VAR1_Ithick  V1It_PREF   V1It_SUFF 
    111 export VAR1_SNOW    V1ICE_SNOW  V1SNOW_SUFF 
     102export VAR1_T       V1T_PREF     V1T_SUFF 
     103export VAR1_S       V1S_PREF     V1S_SUFF 
     104export VAR1_SSH     V1SSH_PREF   V1SSH_SUFF 
     105export VAR1_Q       V1Q_PREF     V1Q_SUFF 
     106export VAR1_EMP     V1EMP_PREF   V1EMP_SUFF 
     107export VAR1_U       V1U_PREF     V1U_SUFF 
     108export VAR1_V       V1V_PREF     V1V_SUFF 
     109export VAR1_ICE     V1ICE_PREF   V1ICE_SUFF 
     110export VAR1_Ithick  V1It_PREF    V1It_SUFF 
     111export VAR1_SNOW    V1SNOW_PREF  V1SNOW_SUFF 
    112112 
    113113#===================== EXP2 ===================== 
     
    122122export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
    123123export VAR2_Ithick  V2It_PREF   V2It_SUFF 
    124 export VAR2_SNOW    V2ICE_SNOW V2SNOW_SUFF 
     124export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
    125125# 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example2

    r3928 r3999  
    2020PS_DIR=$( pwd )/ps_tag33_tag331 
    2121PDF_DIR=$( pwd )/pdf_tag33_tag331 
    22 HTML_DIR=$( pwd )/html_tag33_tag331 
    23 SAXO_DIR=/Users/smasson/SAXO_DIR 
     22HTML_DIR=$( pwd )/html 
     23SAXO_DIR=/Users/sflod/SAXO_DIR 
    2424# 
    25 DIR_DATA=/Users/smasson/tmp/idlorca2/res      # path of data in NetCDF format 
    26 DIR_CLIMATO=/Users/smasson/tmp/idlorca2/clim  # path of climatological data 
    27 DIR_MASK=/Users/smasson/tmp/idlorca2/clim     # path of mask files (ex: subbasins) 
     25DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     26DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     27DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    2828# 
    2929#===================== Model GRID ===================== 
    3030# 
    31 FILE_MESH_MASK=/Users/smasson/data/vp33beta_mesh_mask_7.nc # meshmask 
     31FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 
    3232FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks 
    3333# 
    3434#===================== DATA ===================== 
    3535# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     36VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     37VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3838VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3939VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    4040VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
    41 VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent North Emisphere 
    42 VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent South Emisphere 
    43 VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area North Emisphere 
    44 VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area South Emisphere 
     41VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     42VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     43VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     44VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
    4545# 
    4646# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    100100# 
    101101export DATE1        DATE2 
    102 export VAR1_T       V1T_PREF    V1T_SUFF 
    103 export VAR1_S       V1S_PREF    V1S_SUFF 
    104 export VAR1_SSH     V1SSH_PREF  V1SSH_SUFF 
    105 export VAR1_Q       V1Q_PREF    V1Q_SUFF 
    106 export VAR1_EMP     V1EMP_PREF  V1EMP_SUFF 
    107 export VAR1_U       V1U_PREF    V1U_SUFF 
    108 export VAR1_V       V1V_PREF    V1V_SUFF 
    109 export VAR1_ICE     V1ICE_PREF  V1ICE_SUFF 
    110 export VAR1_Ithick  V1It_PREF   V1It_SUFF 
    111 export VAR1_SNOW    V1ICE_SNOW  V1SNOW_SUFF 
     102export VAR1_T       V1T_PREF     V1T_SUFF 
     103export VAR1_S       V1S_PREF     V1S_SUFF 
     104export VAR1_SSH     V1SSH_PREF   V1SSH_SUFF 
     105export VAR1_Q       V1Q_PREF     V1Q_SUFF 
     106export VAR1_EMP     V1EMP_PREF   V1EMP_SUFF 
     107export VAR1_U       V1U_PREF     V1U_SUFF 
     108export VAR1_V       V1V_PREF     V1V_SUFF 
     109export VAR1_ICE     V1ICE_PREF   V1ICE_SUFF 
     110export VAR1_Ithick  V1It_PREF    V1It_SUFF 
     111export VAR1_SNOW    V1SNOW_PREF  V1SNOW_SUFF 
    112112 
    113113#===================== EXP2 ===================== 
     
    122122export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
    123123export VAR2_Ithick  V2It_PREF   V2It_SUFF 
    124 export VAR2_SNOW    V2ICE_SNOW V2SNOW_SUFF 
     124export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
    125125# 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example_ipcc

    r3927 r3999  
    3434#===================== DATA ===================== 
    3535# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     36VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     37VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3838VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3939VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    4040VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
    41 VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent North Emisphere 
    42 VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent South Emisphere 
    43 VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area North Emisphere 
    44 VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area South Emisphere 
     41VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     42VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     43VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     44VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
    4545# 
    4646# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    100100# 
    101101export DATE1        DATE2 
    102 export VAR1_T       V1T_PREF    V1T_SUFF 
    103 export VAR1_S       V1S_PREF    V1S_SUFF 
    104 export VAR1_SSH     V1SSH_PREF  V1SSH_SUFF 
    105 export VAR1_Q       V1Q_PREF    V1Q_SUFF 
    106 export VAR1_EMP     V1EMP_PREF  V1EMP_SUFF 
    107 export VAR1_U       V1U_PREF    V1U_SUFF 
    108 export VAR1_V       V1V_PREF    V1V_SUFF 
    109 export VAR1_ICE     V1ICE_PREF  V1ICE_SUFF 
    110 export VAR1_Ithick  V1It_PREF   V1It_SUFF 
    111 export VAR1_SNOW    V1ICE_SNOW  V1SNOW_SUFF 
     102export VAR1_T       V1T_PREF     V1T_SUFF 
     103export VAR1_S       V1S_PREF     V1S_SUFF 
     104export VAR1_SSH     V1SSH_PREF   V1SSH_SUFF 
     105export VAR1_Q       V1Q_PREF     V1Q_SUFF 
     106export VAR1_EMP     V1EMP_PREF   V1EMP_SUFF 
     107export VAR1_U       V1U_PREF     V1U_SUFF 
     108export VAR1_V       V1V_PREF     V1V_SUFF 
     109export VAR1_ICE     V1ICE_PREF   V1ICE_SUFF 
     110export VAR1_Ithick  V1It_PREF    V1It_SUFF 
     111export VAR1_SNOW    V1SNOW_PREF  V1SNOW_SUFF 
    112112 
    113113#===================== EXP2 ===================== 
     
    122122export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
    123123export VAR2_Ithick  V2It_PREF   V2It_SUFF 
    124 export VAR2_SNOW    V2ICE_SNOW V2SNOW_SUFF 
     124export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
    125125# 
  • branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example_old_names

    r3927 r3999  
    3434#===================== DATA ===================== 
    3535# 
    36 VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc    # Levitus_T 
    37 VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=data_1y_salinity_nomask.nc                  # Levitus_S 
     36VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     37VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
    3838VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    3939VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    4040VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
    41 VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent North Emisphere 
    42 VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index.nc                               # Ice Extent South Emisphere 
    43 VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area North Emisphere 
    44 VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index.nc                               # Ice Area South Emisphere 
     41VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     42VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     43VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     44VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
    4545# 
    4646# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
Note: See TracChangeset for help on using the changeset viewer.