Changeset 4320


Ignore:
Timestamp:
2013-11-28T12:41:27+01:00 (7 years ago)
Author:
flavoni
Message:

update IDL_scripts for ORCA2_LIM and ORCA2_LIM3, add plot of max_mld output 5days, see ticket: #724

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
23 added
21 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_main.sh

    r2751 r4320  
    173173       pdflist=${pdflist}' '${PDF_DIR}/${file}pdf 
    174174   done 
     175#pdfjam needed to use pdfjoin 
     176   pdfjoin ${pdflist} --outfile $PDF_DIR/${filepdf} 
    175177   texexec --pdfarrange --result=$PDF_DIR/$filepdf $pdflist 
    176178   rm -f $PDF_DIR/$( basename $filepdf .pdf ).aux $PDF_DIR/$( basename $filepdf .pdf ).log 
    177    echo ${PDF_DIR}/$filepdf done 
     179   if [ ! -e ${PDF_DIR}/$filepdf ] 
     180   then 
     181           echo "\n ATTENTION !!!!! \n" 
     182                echo "${PDF_DIR}/$filepdf do not exist \n" 
     183        else 
     184                 echo "${PDF_DIR}/$filepdf done" 
     185       fi 
    178186#commented because convert gives images of bad quality 
    179187#  convert ${PS_DIR}/*.ps /tmp/all_${PLOTTYPE}.pdf 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.pro

    r2751 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_GlobMeanTS.pro

    r2751 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r2751 r4320  
    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                         
    43    
     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 
     50 
    4451  domdef 
    4552 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro

    r2751 r4320  
    1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
    14   IF keyword_set(march) THEN var = var+'March' 
     14  IF keyword_set(april) THEN var = var+'April' 
    1515  IF keyword_set(sept) THEN var = var+'Sept' 
    1616 
    1717  filename = cdti3 + '_'+var+'_'+std_file1_I 
    18   if std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I 
    19   if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     18  IF std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I 
     19  IF KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    2020   
    2121  IF keyword_set(arc) THEN BEGIN  
    2222    domdef, 20, 380, 50, 90 
    2323    vmin = 0. 
    24     vmax = 10. 
     24    vmax = 6. 
    2525    vint = 0.5 
    26     fmt = '(i2)' 
     26    fmtd = '(f4.1)' 
    2727    vmind = -3. 
    2828    vmaxd = 3. 
    2929    vintd = 0.5 
    30     fmtd = '(i2)' 
     30    fmtd = '(f4.1)' 
     31    div = 6 
    3132    map = [90, 0, 0] 
    3233  ENDIF 
     
    4142    vintd = 0.2 
    4243    fmtd = '(f4.1)' 
     44    div = 5 
    4345    map = [-90, 0, 0] 
    4446  ENDIF 
     
    4850  title = var+'!C'+std_file1_I 
    4951  plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 
    50        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
     52       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div $ 
    5153       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 
    5254;                                ; 
     
    5456    title = title + std_file2_I 
    5557    plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    56          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
     58         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title  $ 
    5759         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
    5860  endif                         
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r2755 r4320  
    1515  PRINT, '' 
    1616; 
    17   std_iodir_data    = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format') 
    18   std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data') 
    19   std_iodir_mask    = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)') 
     17  std_iodir_data    = isadirectory(getenv('DIR_DATA'),     title = 'path of data in NetCdf format') 
     18  std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'),  title = 'path of climatological data') 
     19  std_iodir_mask    = isadirectory(getenv('DIR_MASK'),     title = 'path of mask files (ex: subbasins)') 
    2020; meshmask 
    21   std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask) 
     21  std_file_mesh = isafile(getenv('FILE_MESH_MASK'),        title = 'mesh_mask', iodir = std_iodir_mask) 
    2222  std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask) 
    2323 
    24 ; Levitus 98 
    25   std_file_Levitus_T =  isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 
    26   std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 
    27   std_file_reynolds =  isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 
    28   std_file_oaflux =  isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 
    29   std_file_mld =  isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 
     24; climatologies 
     25  std_file_Levitus_T =  isafile(getenv('FILE_TEMP_3D'),    title = 'Levitus_T', iodir = std_iodir_climato) 
     26  std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'),     title = 'Levitus_S', iodir = std_iodir_climato) 
     27  std_file_reynolds  =  isafile(getenv('FILE_SST'),        title = 'Reynolds', iodir = std_iodir_climato) 
     28  std_file_oaflux    =  isafile(getenv('FILE_FLUX'),       title = 'oaflux', iodir = std_iodir_climato) 
     29  std_file_mld       =  isafile(getenv('FILE_MLD'),        title = 'Mixed layer depth', iodir = std_iodir_climato) 
     30  std_file_ice       =  isafile(getenv('FILE_ICE'),        title = 'ICE', iodir = std_iodir_climato) 
     31  std_file_snow_arc  =  isafile(getenv('FILE_SNOW_ARC'),   title = 'SNOW_ARC', iodir = std_iodir_climato) 
     32  std_file_snow_ant  =  isafile(getenv('FILE_SNOW_ANT'),   title = 'SNOW_ANT', iodir = std_iodir_climato) 
    3033 
    3134  IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $ 
     
    155158  Ithi_1 = read_ncdf(getenv('VAR1_Ithick'), allrecords = allrec, filename = std_file1_I )  
    156159  caldat, time, mm 
    157   march = where(mm EQ 3, cnt) 
    158   Ithi_march_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(march)], 3), unit:Ithi_1.unit}  
     160  april = where(mm EQ 4, cnt) 
     161  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
    159162  sept = where(mm EQ 9, cnt) 
    160   Ithi_sept_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit}  
     163  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    161164  undefine, Ithi_1 
    162165; 
    163166  IF std_file2_I NE std_file1_I THEN BEGIN 
    164     Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )  
    165     caldat, time, mm 
    166     march = where(mm EQ 3, cnt) 
    167     Ithi_march_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(march)], 3), unit:Ithi_2.unit}  
    168     sept = where(mm EQ 9, cnt) 
    169     Ithi_sept_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(sept)], 3), unit:Ithi_2.unit}  
    170     undefine, Ithi_2 
     167     Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )  
     168     caldat, time, mm 
     169     april = where(mm EQ 4, cnt) 
     170     Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     171     sept = where(mm EQ 9, cnt) 
     172     Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     173     undefine, Ithi_2 
    171174  ENDIF ELSE BEGIN  
    172     Ithi_march_2 = {arr:-1} 
    173     Ithi_sept_2 = {arr:-1} 
     175     Ithi_april_2 = {arr:-1} 
     176     Ithi_sept_2 = {arr:-1} 
    174177  ENDELSE 
    175178; 
    176179  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
     180  help, Ifra_1.arr, jpi,jpj 
    177181  caldat, time, mm 
    178182  march = where(mm EQ 3, cnt) 
    179   Ifra_march_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(march)], 3), unit:Ifra_1.unit}  
     183  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    180184  sept = where(mm EQ 9, cnt) 
    181   Ifra_sept_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit} 
     185  Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 
    182186  undefine, Ifra_1 
    183187; 
     
    186190    caldat, time, mm 
    187191    march = where(mm EQ 3, cnt) 
    188     Ifra_march_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(march)], 3), unit:Ifra_2.unit} 
     192    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    189193    sept = where(mm EQ 9, cnt) 
    190     Ifra_sept_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit} 
     194    Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    191195    undefine, Ifra_2 
    192196  ENDIF ELSE BEGIN  
     
    194198    Ifra_sept_2 = {arr:-1} 
    195199  ENDELSE 
     200; 
     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 
    196223; 
    197224  jpt = 1 
     
    299326  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex  
    300327   
    301   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    302   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    303    
     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 
     330 
     331  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
     332  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
     333 
    304334  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    305335  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    306336   
    307   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    308   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 
     339 
     340  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
     341  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    309342 
    310343  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
     
    322355  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    323356  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
     357;; SF non funziona, dice che deve essere una struttura e non e!!!! 
     358;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
     359;;SF  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 
    324375 
    325376  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_bsf.pro

    r2751 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mld.pro

    r2751 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro

    r2760 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1

    r2753 r4320  
    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 
    21 PDF_DIR=$( pwd )/pdf 
    22 HTML_DIR=$( pwd )/html 
    23 SAXO_DIR=/Users/smasson/SAXO_DIR 
     23PS_DIR=$( pwd )/ps_plot 
     24PDF_DIR=$( pwd )/pdf_plot 
     25HTML_DIR=$( pwd )/html_plot 
     26SAXO_DIR=/Users/sflod/SAXO_DIR 
    2427# 
    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) 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     29DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     30DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    2831# 
    2932#===================== Model GRID ===================== 
    3033# 
    31 FILE_MESH_MASK=/Users/smasson/data/vp33beta_mesh_mask_7.nc # meshmask 
     34FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 
    3235FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks 
    3336# 
    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 
    38 VAR_SST=sst            ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    39 VAR_FLUX=qnet          ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    40 VAR_MLD=mld            ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     39VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     40VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
     41VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
     42VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
     43VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     44VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     45VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     46VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     47VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
     48VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic  
     49VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic 
    4150# 
    4251# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    6069VAR1_Ithick=sit 
    6170VAR1_Ifrac=sic 
     71VAR1_SNOW=isnowthi 
    6272# 
    6373#===================== EXP2 ===================== 
     
    7787VAR2_Ithick=sit 
    7888VAR2_Ifrac=sic 
     89VAR2_SNOW=isnowthi 
    7990# 
    8091######################### Export Variables ############################### 
     
    8697export FILE_MESH_MASK FILE_MASK_SUBDOMAIN 
    8798#===================== DATA ===================== 
    88 export FILE_TEMP_3D  VAR_TEMP_3D 
    89 export FILE_SAL_3D   VAR_SAL_3D 
    90 export FILE_SST      VAR_SST 
    91 export FILE_GEOHEAT  VAR_GEOHEAT 
    92 export FILE_FLUX     VAR_FLUX 
    93 export FILE_MLD      VAR_MLD 
     99export FILE_TEMP_3D   VAR_TEMP_3D 
     100export FILE_SAL_3D    VAR_SAL_3D 
     101export FILE_SST       VAR_SST 
     102export FILE_GEOHEAT   VAR_GEOHEAT 
     103export FILE_FLUX      VAR_FLUX 
     104export FILE_MLD       VAR_MLD 
     105export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 
     106export FILE_SNOW_ARC  VAR_SNOW_NH 
     107export FILE_SNOW_ANT  VAR_SNOW_SH 
    94108#===================== EXP1 ===================== 
    95109export READ_ONLY_FIRST_RECORD 
    96 export FILE1_T FILE1_U FILE1_V FILE1_I 
    97 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    98 export VAR1_U VAR1_V 
    99 export VAR1_Ithick VAR1_Ifrac 
     110export FILE1_T     FILE1_U  FILE1_V  FILE1_I 
     111export VAR1_T      VAR1_S  VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
     112export VAR1_U      VAR1_V 
     113export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
    100114#===================== EXP2 ===================== 
    101 export FILE2_T FILE2_U FILE2_V FILE2_I 
    102 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    103 export VAR2_U VAR2_V 
    104 export VAR2_Ithick VAR2_Ifrac 
     115export FILE2_T     FILE2_U    FILE2_V  FILE2_I 
     116export VAR2_T      VAR2_S    VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
     117export VAR2_U      VAR2_V 
     118export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
    105119# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2

    r2753 r4320  
    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 
    21 PDF_DIR=$( pwd )/pdf 
    22 HTML_DIR=$( pwd )/html 
    23 SAXO_DIR=/Users/smasson/SAXO_DIR 
     23PS_DIR=$( pwd )/500yfwb0_ps_plot_500y 
     24PDF_DIR=$( pwd )/500yfwb0_pdf_plot_500y 
     25HTML_DIR=$( pwd )/html_plot 
     26SAXO_DIR=/Users/sflod/SAXO_DIR 
    2427# 
    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) 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
     29DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     30DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    2831# 
    2932#===================== Model GRID ===================== 
    3033# 
    31 FILE_MESH_MASK=/Users/smasson/data/vp33beta_mesh_mask_7.nc # meshmask 
     34FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 
    3235FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks 
    3336# 
    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 
    38 VAR_SST=sst            ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    39 VAR_FLUX=qnet          ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    40 VAR_MLD=mld            ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     39VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3 
     40VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3 
     41VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
     42VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
     43VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     44VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere 
     45VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere 
     46VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere 
     47VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere 
     48VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic  
     49VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic 
    4150# 
    4251# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    6069VAR1_Ithick=sit 
    6170VAR1_Ifrac=sic 
     71VAR1_SNOW=isnowthi 
    6272# 
    6373#===================== EXP2 ===================== 
     
    7787VAR2_Ithick=sit 
    7888VAR2_Ifrac=sic 
     89VAR2_SNOW=isnowthi 
    7990# 
    8091######################### Export Variables ############################### 
     
    8697export FILE_MESH_MASK FILE_MASK_SUBDOMAIN 
    8798#===================== DATA ===================== 
    88 export FILE_TEMP_3D  VAR_TEMP_3D 
    89 export FILE_SAL_3D   VAR_SAL_3D 
    90 export FILE_SST      VAR_SST 
    91 export FILE_GEOHEAT  VAR_GEOHEAT 
    92 export FILE_FLUX     VAR_FLUX 
    93 export FILE_MLD      VAR_MLD 
     99export FILE_TEMP_3D   VAR_TEMP_3D 
     100export FILE_SAL_3D    VAR_SAL_3D 
     101export FILE_SST       VAR_SST 
     102export FILE_GEOHEAT   VAR_GEOHEAT 
     103export FILE_FLUX      VAR_FLUX 
     104export FILE_MLD       VAR_MLD 
     105export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH, VAR_ICE_area_NH, VAR_ICE_area_SH 
     106export FILE_SNOW_ARC  VAR_SNOW_NH 
     107export FILE_SNOW_ANT  VAR_SNOW_SH 
    94108#===================== EXP1 ===================== 
    95109export READ_ONLY_FIRST_RECORD 
    96 export FILE1_T FILE1_U FILE1_V FILE1_I 
    97 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    98 export VAR1_U VAR1_V 
    99 export VAR1_Ithick VAR1_Ifrac 
     110export FILE1_T     FILE1_U  FILE1_V  FILE1_I 
     111export VAR1_T      VAR1_S  VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
     112export VAR1_U      VAR1_V 
     113export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
    100114#===================== EXP2 ===================== 
    101 export FILE2_T FILE2_U FILE2_V FILE2_I 
    102 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    103 export VAR2_U VAR2_V 
    104 export VAR2_Ithick VAR2_Ifrac 
     115export FILE2_T     FILE2_U    FILE2_V  FILE2_I 
     116export VAR2_T      VAR2_S    VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
     117export VAR2_U      VAR2_V 
     118export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
    105119# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_zonal_S.pro

    r2751 r4320  
    3333; 
    3434  title = 'Salinity, '+subname+'!C'+std_file1_T 
    35   pltz, S1, MININ = 32., MAXIN = 37., INTER = .2, typein = 'yz', FORMAT = '(f4.1)' $  
     35  pltz, S1, MININ = 31., MAXIN = 36.8, INTER = .2, typein = 'yz', FORMAT = '(f4.1)' $  
    3636        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT,  _extra = ex  
    3737;  
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro

    r2751 r4320  
    1414  ts_Time = 0. 
    1515 
    16   trans = strlowcase(var_name) EQ 'vocetr_eff' 
     16  ;SF trans = strlowcase(var_name) EQ 'vocetr_eff' 
     17  trans = strlowcase(var_name) EQ 'vomeeftr' 
    1718 
     19  ; ADD definition of mask_atl_nomed to compute AMOC (with no Mediterranean Sea) 
     20  msk_atl_nomed = read_ncdf( 'atlmsk_nomed', filename = std_file_msksub, /nostruct, _extra = ex ) 
    1821  FOR i = 0, nfiles-1 DO BEGIN 
    1922    var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct) 
     
    2124; 
    2225    FOR t = 0, jpt-1 DO BEGIN 
     26      msfatl = msf(var[*, *, *, t], msk_atl_nomed, TRANSPORT = trans, /nostruct, indexboxzoom = ind) 
    2327      msftot = msf(var[*, *, *, t], TRANSPORT = trans, /nostruct, indexboxzoom = ind) 
    2428      yaxis = gphit[ind[0], ind[2]:ind[3]] 
     
    2731      indy = where(yaxis gt 40 and yaxis le 50) 
    2832      domdef, 0, 3500 
    29       t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 
     33      ;SF commented because we've to compute msftot starting from msf atlantic; i.e. without mediterranean sea. 
     34      ;SF t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 
     35      t45 = [t45, max(msfatl[indy, firstzw:lastzw], /NaN)] 
    3036       
    3137; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S 
     
    7884 
    7985; compute the MSF 
    80  
     86 
    8187  std_ts_read_AMOC, vV1, date1, date2, prefix, suffix, t45, t70, t50 
    8288  if prefix NE prefix2 then begin 
     
    114120  if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 
    115121 
    116   title = prefix+' '+d1_d2+'!C'+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 
     122  title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 
    117123  pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $ 
    118124        ,  small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Drake.pro

    r2751 r4320  
    3232  u1 = rseries_ncdf(vU1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) 
    3333  drk1 = fltarr(jpt) 
    34   trans = strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' 
     34   ;SF trans = strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' 
     35   trans = strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' 
    3536  FOR t = 0, jpt-1 DO BEGIN 
    3637    tmp = bsf(u1[*, *, *, t], TRANSPORT = trans, /nostruct) 
     
    4849    u2 = rseries_ncdf(vU2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, /nostruct) 
    4950    drk2 = fltarr(jpt) 
    50     trans = strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' 
     51    ;SF trans = strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' 
     52    trans = strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' 
    5153    FOR t = 0, jpt-1 DO BEGIN 
    5254      tmp = bsf(u2[*, *, *, t], TRANSPORT = trans, /nostruct) 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro

    r2751 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r2751 r4320  
    1 pro std_ts_ICE, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 
     1function read_arr2d, filename, varname, t1, t2 
     2;; function that read input file and return 2d array with monthly timecounter  
     3nyear = (t2-t1+1)/12 
     4arr2d = ncdf_lec(filename, VAR=varname) 
     5arr2d = arr2d[t1:t2] 
     6arr2d = reform(arr2d,12,nyear) ; put in 2D array 
     7arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     8 
     9return, arr2d 
     10end  
     11 
     12;; here start procedure that use function read_arr2d   
     13pro std_ts_ICE, masknp, s_iodir_data,  POSTSCRIPT = postscript, _extra = ex 
    214 
    315  compile_opt idl2, strictarrsubs 
     
    1022; get exp2 info 
    1123  vICE2 = getenv('VAR2_ICE')   &   prefix2 = getenv('V2ICE_PREF')   &   suffix2 = getenv('V2ICE_SUFF') 
     24; get ice climatology info 
     25  std_file_ice =  isafile(getenv('FILE_ICE'), title = 'ICE Extent Climatology', iodir = std_iodir_climato) 
    1226; 
     27  time_ice = ncdf_lec( std_file_ice, VAR='time' ) 
     28  time_ice = (time_ice - FLOOR(time_ice) ) * 12 
     29  time_ice = (round(time_ice) + 11) mod 12; round to nearest integer 
     30  t1 = where(time_ice eq 0) 
     31  t1 = t1[0] ;  jannuary 
     32  t2 = where(time_ice eq 11, count) 
     33  t2 = t2[count-1] ; last day of december 
     34 
     35  vICE_ext_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_NH'), t1, t2 ) 
     36  vICE_ext_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_SH'), t1, t2 ) 
     37; 
     38  vICE_area_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_NH'), t1, t2 ) 
     39  vICE_area_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_SH'), t1, t2 ) 
     40 
    1341  cdti3 = string(cnt, format = '(i3.3)') 
    1442  print, cdti3 + ') ' + blabla 
    15   filename = cdti3 + '_ts_AMOC_'+prefix 
     43  filename = cdti3 + '_ts_ICE_'+prefix 
    1644  if prefix NE prefix2 then filename = filename + '_'+prefix2 
    17   if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1 
     45  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    1846; 
    1947  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
    2048; 
    2149  iodir = std_iodir_data 
     50  ; ICE Area in NORTH Hemisphere 
     51  domdef, 0, jpi-1, 30, 90, /xindex 
     52  ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     53  ICE_N = grossemoyenne(ICE_N.arr, 'xy', /integration, mask2d = masknp) 
    2254 
    23   domdef, 0, jpi-1, 30, 90, /xindex 
    24   ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    25   ICE_N.arr = ICE_N.arr * 1.e-12   &   ICE_N.unit = '10^12 m^2' 
    26   domdef, 0, jpi-1, -90, -30, /xindex 
    27   ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    28   ICE_S.arr = ICE_S.arr * 1.e-12   &   ICE_S.unit = '10^12 m^2' 
     55  if jpt mod 12 ne 0 then stop 
     56  nyr=jpt/12. 
     57  ICE_N = reform(ice_n, 12,nyr) 
     58  ICE_n = total(ice_n,2)/nyr 
     59  ICE_N = {arr:ICE_N * 1.e-12, unit : '10^12 m^2'} 
    2960 
    30   title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 
    31   pltt, ICE_N, 't', 0., 15., date1, date2, /REMPLI, /PORTRAIT $ 
     61  ; ICE EXTENT (Area minus 15%) in NORTH Hemisphere 
     62  ICE_N_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     63  msk = ICE_N_15.arr gt 0.15 ; remove 0.15% for observations 
     64  ICE_N_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 
     65  if jpt mod 12 ne 0 then stop 
     66  nyr=jpt/12. 
     67  ICE_N_15 = reform(ice_n_15, 12,nyr) 
     68  ICE_n_15 = total(ice_n_15,2)/nyr 
     69  ICE_N_15 = {arr:ICE_N_15 * 1.e-12, unit : '10^12 m^2'} 
     70  ; 
     71  ;ICE Area in SOUTH Hemisphere 
     72  domdef, 0, jpi-1, -90, -30, /xindex  
     73  ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     74  ICE_S = grossemoyenne(ICE_S.arr, 'xy', /integration, mask2d = masknp) 
     75  if jpt mod 12 ne 0 then stop 
     76  nyr=jpt/12. 
     77  ICE_S = reform(ice_S, 12,nyr) 
     78  ICE_S = total(ice_S,2)/nyr 
     79  ICE_S = {arr:ICE_S * 1.e-12, unit : '10^12 m^2'} 
     80  ; ICE EXTENT (Area minus 15%) in SOUTH Hemisphere 
     81  ICE_S_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     82  msk = ICE_S_15.arr gt 0.15 ; remove 0.15% for observations 
     83  ICE_S_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 
     84  if jpt mod 12 ne 0 then stop 
     85  nyr=jpt/12. 
     86  ICE_S_15 = reform(ice_S_15, 12,nyr) 
     87  ICE_S_15 = total(ice_S_15,2)/nyr 
     88  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
     89  ; 
     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)' 
     91  jpt=12 
     92  time=julday(1,15,1900)+30*lindgen(12) 
     93  pltt, ICE_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 4., MAX = 16. $ 
    3294        , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    33    
    34   title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 
    35   pltt, ICE_S, 't', 0., 15., date1, date2, /REMPLI, /NOERASE $ 
     95  pltt, ICE_N_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; dashed lines is LINESTYLE=2  $ 
     96        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     97  pltt, vICE_area_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     98         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     99  pltt, vICE_ext_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     100        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     101; 
     102  title ='Southern 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)' 
     103  pltt, ICE_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 20. $ 
    36104        , small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    37    
    38   htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png  />  ' ] 
     105  pltt, ICE_S_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $  
     106        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     107  pltt,  vICE_area_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     108         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     109  pltt,  vICE_ext_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     110        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     111; 
     112  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
    39113  if KEYWORD_SET(postscript) then closeps 
    40114 
     
    42116 
    43117    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
    44  
    45118    tsave = time 
    46119    domdef, 0, jpi-1, 30, 90, /xindex 
    47     ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    48     ICE_N2.arr = ICE_N2.arr * 1.e-12   &   ICE_N2.unit = '10^12 m^2' 
    49     domdef, 0, jpi-1, -90, -30, /xindex 
    50     ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    51     ICE_S2.arr = ICE_S2.arr * 1.e-12   &   ICE_S2.unit = '10^12 m^2' 
    52     time = tsave   &   IF n_elements(time) NE jpt THEN stop 
     120    ;ICE Extent in NORTH Hemisphere 
     121    ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     122    ICE_N2 = grossemoyenne(ICE_N2.arr, 'xy', /integration, mask2d = masknp) 
     123    if jpt mod 12 ne 0 then stop 
     124    nyr=jpt/12. 
     125    ICE_N2 = reform(ICE_N2, 12,nyr) 
     126    ICE_N2 = total(ICE_N2,2)/nyr 
     127    ICE_N2 = {arr:ICE_N2 * 1.e-12, unit : '10^12 m^2'} 
     128    ;ICE Extent minus 15% in NORTH Hemisphere 
     129    ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     130    msk = ICE_N2_15.arr gt 0.15  ; remove 0.15% for observations 
     131    ICE_N2_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 
     132    if jpt mod 12 ne 0 then stop 
     133    nyr=jpt/12. 
     134    ICE_N2_15 = reform(ICE_N2_15, 12,nyr) 
     135    ICE_N2_15 = total(ICE_N2_15,2)/nyr 
     136    ICE_N2_15 = {arr:ICE_N2_15 * 1.e-12, unit : '10^12 m^2'} 
     137    ;ICE Extent in SOUTH Hemisphere 
     138    domdef, 0, jpi-1, -90, -30, /xindex  
     139    ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     140    ICE_S2 = grossemoyenne(ICE_S2.arr, 'xy', /integration, mask2d = masknp) 
     141    if jpt mod 12 ne 0 then stop 
     142    nyr=jpt/12. 
     143    ICE_S2 = reform(ICE_S2, 12,nyr) 
     144    ICE_S2 = total(ICE_S2,2)/nyr 
     145    ICE_S2 = {arr:ICE_S2 * 1.e-12, unit : '10^12 m^2'} 
     146    ;ICE Extent minus 15% in SOUTH Hemisphere 
     147    ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     148    msk = ICE_S2_15.arr gt 0.15  ; remove 0.15% for observations 
     149    ICE_S2_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 
     150    if jpt mod 12 ne 0 then stop 
     151    nyr=jpt/12. 
     152    ICE_S2_15 = reform(ICE_S2_15, 12,nyr) 
     153    ICE_S2_15 = total(ICE_S2_15,2)/nyr 
     154    ICE_S2_15 = {arr:ICE_S2_15 * 1.e-12, unit : '10^12 m^2'} 
     155  ; 
     156 ;   time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    53157 
    54158    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    55159 
    56     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 
    57     pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 
    58           , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    59      
    60     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 
    61     pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 
    62           , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     160 
     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; 
    63184 
    64185    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro

    r2751 r4320  
    2222; 
    2323  ts_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 
    24   ts_SSH.arr = ts_SSH.arr * 100   &   ts_SSH.unit = 'cm' 
     24  ;SF ts_SSH.arr = ts_SSH.arr * 100   &   ts_SSH.unit = 'cm' 
     25  ts_SSH.arr = ts_SSH.arr   &   ts_SSH.unit = 'm' 
    2526  title = prefix+' '+d1_d2+'!C'+blabla 
    26   pltt, ts_SSH, 't', -.1, .1, date1, date2, /REMPLI $ 
    27         , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 
     27  pltt, ts_SSH, 't', -12., 12., date1, date2, /REMPLI $ 
     28  ;SF       , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 
     29        , small = [1, 2, 1], YTITLE = 'm', TITLE = title, /PORTRAIT, _extra = ex 
    2830   
    2931  IF prefix NE prefix2 THEN BEGIN 
     
    3234    tsave = time 
    3335    ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 
    34     ts_SSH2.arr = ts_SSH2.arr * 100   &   ts_SSH2.unit = 'cm' 
     36    ;SF ts_SSH2.arr = ts_SSH2.arr * 100   &   ts_SSH2.unit = 'cm' 
     37    ts_SSH2.arr = ts_SSH2.arr   &   ts_SSH2.unit = 'm' 
    3538    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    3639 
    3740    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    38     pltt, ts_SSH.arr - ts_SSH2.arr, 't', -.1, .1, date1, date2, /REMPLI $ 
    39           , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 
     41    pltt, ts_SSH.arr - ts_SSH2.arr, 't', -10., 10., date1, date2, /REMPLI $ 
     42          ;SF , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 
     43          , color = 250, small = [1, 2, 2], YTITLE = 'm', TITLE = title, /noerase, _extra = ex 
    4044 
    4145  ENDIF  
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro

    r3610 r4320  
    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_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r2751 r4320  
    2424; reading variables 
    2525  masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill) 
     26; climatologies 
     27  std_file_Levitus_T =  isafile(getenv('FILE_TEMP_3D'),  title = 'Levitus_T', iodir = std_iodir_climato) 
     28  std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'),   title = 'Levitus_S', iodir = std_iodir_climato) 
     29  std_file_reynolds  =  isafile(getenv('FILE_SST'),      title = 'Reynolds', iodir = std_iodir_climato) 
     30  std_file_oaflux    =  isafile(getenv('FILE_FLUX'),     title = 'oaflux', iodir = std_iodir_climato) 
     31  std_file_mld       =  isafile(getenv('FILE_MLD'),      title = 'Mixed layer depth', iodir = std_iodir_climato) 
     32  std_file_ice       =  isafile(getenv('FILE_ICE'),      title = 'ICE', iodir = std_iodir_climato) 
     33  std_file_snow_arc  =  isafile(getenv('FILE_SNOW_ARC'), title = 'SNOW_ARC', iodir = std_iodir_climato) 
     34  std_file_snow_ant  =  isafile(getenv('FILE_SNOW_ANT'), title = 'SNOW_ANT', iodir = std_iodir_climato) 
     35 
    2636 
    2737  date1   = long(getenv('DATE1'))     &   date2   = long(getenv('DATE2')) 
     
    6474  cnt = cnt+1   &   blabla = 'Sea-Ice cover' 
    6575  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 
     76; 
     77  cnt = cnt+1   &   blabla = 'Sea-Ice Volume' 
     78  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    6679 
     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 
     82; 
    6783  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] 
    6884 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example1

    r2753 r4320  
    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 
    38 VAR_SST=sst            ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    39 VAR_FLUX=qnet          ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    40 VAR_MLD=mld            ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     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 
     38VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
     39VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
     40VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     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 
    4145# 
    4246# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    4852READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    4953# 
    50 DATE1=10101   ;   DATE2=20001231 
     54DATE1=20010101   ;   DATE2=21001231 
    5155# 
    52 VAR1_T=thetao     ;   V1T_PREF=tmx_2000     ;   V1T_SUFF=_1Y_grid_T.nc 
    53 VAR1_S=so         ;   V1S_PREF=tmx_2000     ;   V1S_SUFF=_1Y_grid_T.nc 
    54 VAR1_SSH=zos      ;   V1SSH_PREF=tmx_2000   ;   V1SSH_SUFF=_1Y_grid_T.nc 
    55 VAR1_Q=qns+qsr    ;   V1Q_PREF=tmx_2000     ;   V1Q_SUFF=_1Y_grid_T.nc 
    56 VAR1_EMP=wfo      ;   V1EMP_PREF=tmx_2000   ;   V1EMP_SUFF=_1Y_grid_T.nc 
    57 VAR1_U=uoce_eff   ;   V1U_PREF=tmx_2000     ;   V1U_SUFF=_1Y_grid_U.nc 
    58 VAR1_V=voce_eff   ;   V1V_PREF=tmx_2000     ;   V1V_SUFF=_1Y_grid_V.nc 
    59 VAR1_ICE=sic      ;   V1ICE_PREF=tmx_2000   ;   V1ICE_SUFF=_1Y_icemod.nc 
     56VAR1_T=thetao      ;   V1T_PREF=tmx_2000     ;   V1T_SUFF=_1Y_grid_T.nc 
     57VAR1_S=so          ;   V1S_PREF=tmx_2000     ;   V1S_SUFF=_1Y_grid_T.nc 
     58VAR1_SSH=zos       ;   V1SSH_PREF=tmx_2000   ;   V1SSH_SUFF=_1Y_grid_T.nc 
     59VAR1_Q=qns+qsr     ;   V1Q_PREF=tmx_2000     ;   V1Q_SUFF=_1Y_grid_T.nc 
     60VAR1_EMP=wfo       ;   V1EMP_PREF=tmx_2000   ;   V1EMP_SUFF=_1Y_grid_T.nc 
     61VAR1_U=uoce_eff    ;   V1U_PREF=tmx_2000     ;   V1U_SUFF=_1Y_grid_U.nc 
     62VAR1_V=voce_eff    ;   V1V_PREF=tmx_2000     ;   V1V_SUFF=_1Y_grid_V.nc 
     63VAR1_ICE=sic       ;   V1ICE_PREF=tmx_2000   ;   V1ICE_SUFF=_1M_icemod.nc 
     64VAR1_Ithick=sit    ;   V1It_PREF=tmx_2000    ;   V1It_SUFF=_1M_icemod.nc 
     65VAR1_SNOW=isnowthi ;   V1SNOW_PREF=tmx_2000  ;   V1SNOW_SUFF=_1M_icemod.nc 
    6066# 
    6167#===================== EXP2 ===================== 
     
    6369DATE1_2=10101   ;   DATE2_2=20001231 
    6470# 
    65 VAR2_T=thetao     ;   V2T_PREF=core2000     ;   V2T_SUFF=_1Y_grid_T.nc 
    66 VAR2_S=so         ;   V2S_PREF=core2000     ;   V2S_SUFF=_1Y_grid_T.nc 
    67 VAR2_SSH=zos      ;   V2SSH_PREF=core2000   ;   V2SSH_SUFF=_1Y_grid_T.nc 
    68 VAR2_Q=qns+qsr    ;   V2Q_PREF=core2000     ;   V2Q_SUFF=_1Y_grid_T.nc 
    69 VAR2_EMP=wfo      ;   V2EMP_PREF=core2000   ;   V2EMP_SUFF=_1Y_grid_T.nc 
    70 VAR2_U=uoce_eff   ;   V2U_PREF=core2000     ;   V2U_SUFF=_1Y_grid_U.nc 
    71 VAR2_V=voce_eff   ;   V2V_PREF=core2000     ;   V2V_SUFF=_1Y_grid_V.nc 
    72 VAR2_ICE=sic      ;   V2ICE_PREF=core2000   ;   V2ICE_SUFF=_1Y_icemod.nc 
     71VAR2_T=thetao      ;   V2T_PREF=core2000     ;   V2T_SUFF=_1Y_grid_T.nc 
     72VAR2_S=so          ;   V2S_PREF=core2000     ;   V2S_SUFF=_1Y_grid_T.nc 
     73VAR2_SSH=zos       ;   V2SSH_PREF=core2000   ;   V2SSH_SUFF=_1Y_grid_T.nc 
     74VAR2_Q=qns+qsr     ;   V2Q_PREF=core2000     ;   V2Q_SUFF=_1Y_grid_T.nc 
     75VAR2_EMP=wfo       ;   V2EMP_PREF=core2000   ;   V2EMP_SUFF=_1Y_grid_T.nc 
     76VAR2_U=uoce_eff    ;   V2U_PREF=core2000     ;   V2U_SUFF=_1Y_grid_U.nc 
     77VAR2_V=voce_eff    ;   V2V_PREF=core2000     ;   V2V_SUFF=_1Y_grid_V.nc 
     78VAR2_ICE=sic       ;   V2ICE_PREF=core2000   ;   V2ICE_SUFF=_1M_icemod.nc 
     79VAR2_Ithick=sit    ;   V2It_PREF=core2000    ;   V2It_SUFF=_1M_icemod.nc 
     80VAR2_SNOW=isnowthi ;   V2ICE_SNOW=core2000   ;   V2SNOW_SUFF=_1M_icemod.nc 
    7381# 
    7482######################### Export Variables ############################### 
     
    8694export FILE_FLUX     VAR_FLUX 
    8795export FILE_MLD      VAR_MLD 
     96export FILE_ICE      VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH  VAR_ICE_area_SH 
     97export FILE_SNOW     VAR_SNOW_NH     VAR_SNOW_SH    VAR_SNOW_area_NH VAR_SNOW_area_SH 
    8898#===================== EXP1 ===================== 
    8999export READ_ONLY_FIRST_RECORD 
    90100# 
    91 export DATE1 DATE2 
    92 export VAR1_T    V1T_PREF    V1T_SUFF 
    93 export VAR1_S    V1S_PREF    V1S_SUFF 
    94 export VAR1_SSH  V1SSH_PREF  V1SSH_SUFF 
    95 export VAR1_Q    V1Q_PREF    V1Q_SUFF 
    96 export VAR1_EMP  V1EMP_PREF  V1EMP_SUFF 
    97 export VAR1_U    V1U_PREF    V1U_SUFF 
    98 export VAR1_V    V1V_PREF    V1V_SUFF 
    99 export VAR1_ICE  V1ICE_PREF  V1ICE_SUFF 
     101export DATE1        DATE2 
     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 
     112 
    100113#===================== EXP2 ===================== 
    101 export DATE1_2 DATE2_2 
    102 export VAR2_T    V2T_PREF    V2T_SUFF 
    103 export VAR2_S    V2S_PREF    V2S_SUFF 
    104 export VAR2_SSH  V2SSH_PREF  V2SSH_SUFF 
    105 export VAR2_Q    V2Q_PREF    V2Q_SUFF 
    106 export VAR2_EMP  V2EMP_PREF  V2EMP_SUFF 
    107 export VAR2_U    V2U_PREF    V2U_SUFF 
    108 export VAR2_V    V2V_PREF    V2V_SUFF 
    109 export VAR2_ICE  V2ICE_PREF  V2ICE_SUFF 
     114export DATE1_2      DATE2_2 
     115export VAR2_T       V2T_PREF    V2T_SUFF 
     116export VAR2_S       V2S_PREF    V2S_SUFF 
     117export VAR2_SSH     V2SSH_PREF  V2SSH_SUFF 
     118export VAR2_Q       V2Q_PREF    V2Q_SUFF 
     119export VAR2_EMP     V2EMP_PREF  V2EMP_SUFF 
     120export VAR2_U       V2U_PREF    V2U_SUFF 
     121export VAR2_V       V2V_PREF    V2V_SUFF 
     122export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
     123export VAR2_Ithick  V2It_PREF   V2It_SUFF 
     124export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
    110125# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example2

    r2753 r4320  
    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 
    38 VAR_SST=sst            ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
    39 VAR_FLUX=qnet          ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
    40 VAR_MLD=mld            ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     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 
     38VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds 
     39VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux 
     40VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth 
     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 
    4145# 
    4246# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing 
     
    5054DATE1=20010101   ;   DATE2=21001231 
    5155# 
    52 VAR1_T=thetao     ;   V1T_PREF=tag33     ;   V1T_SUFF=_1Y_grid_T.nc 
    53 VAR1_S=so         ;   V1S_PREF=tag33     ;   V1S_SUFF=_1Y_grid_T.nc 
    54 VAR1_SSH=zos      ;   V1SSH_PREF=tag33   ;   V1SSH_SUFF=_1Y_grid_T.nc 
    55 VAR1_Q=qt         ;   V1Q_PREF=tag33     ;   V1Q_SUFF=_1Y_grid_T.nc 
    56 VAR1_EMP=wfo      ;   V1EMP_PREF=tag33   ;   V1EMP_SUFF=_1Y_grid_T.nc 
    57 VAR1_U=uocetr_eff ;   V1U_PREF=tag33     ;   V1U_SUFF=_1Y_grid_U.nc 
    58 VAR1_V=vocetr_eff ;   V1V_PREF=tag33     ;   V1V_SUFF=_1Y_grid_V.nc 
    59 VAR1_ICE=sic      ;   V1ICE_PREF=tag33   ;   V1ICE_SUFF=_1Y_icemod.nc 
     56VAR1_T=thetao      ;   V1T_PREF=tag33     ;   V1T_SUFF=_1Y_grid_T.nc 
     57VAR1_S=so          ;   V1S_PREF=tag33     ;   V1S_SUFF=_1Y_grid_T.nc 
     58VAR1_SSH=zos       ;   V1SSH_PREF=tag33   ;   V1SSH_SUFF=_1Y_grid_T.nc 
     59VAR1_Q=qt          ;   V1Q_PREF=tag33     ;   V1Q_SUFF=_1Y_grid_T.nc 
     60VAR1_EMP=wfo       ;   V1EMP_PREF=tag33   ;   V1EMP_SUFF=_1Y_grid_T.nc 
     61VAR1_U=uocetr_eff  ;   V1U_PREF=tag33     ;   V1U_SUFF=_1Y_grid_U.nc 
     62VAR1_V=vocetr_eff  ;   V1V_PREF=tag33     ;   V1V_SUFF=_1Y_grid_V.nc 
     63VAR1_ICE=sic       ;   V1ICE_PREF=tag33   ;   V1ICE_SUFF=_1M_icemod.nc 
     64VAR1_Ithick=sit    ;   V1It_PREF=tag33    ;   V1It_SUFF=_1M_icemod.nc 
     65VAR1_SNOW=isnowthi ;   V1SNOW_PREF=tag33  ;   V1SNOW_SUFF=_1M_icemod.nc 
    6066# 
    6167#===================== EXP2 ===================== 
     
    6369DATE1_2=20010101   ;   DATE2_2=21001231 
    6470# 
    65 VAR2_T=thetao     ;   V2T_PREF=trunk331     ;   V2T_SUFF=_1Y_grid_T.nc 
    66 VAR2_S=so         ;   V2S_PREF=trunk331     ;   V2S_SUFF=_1Y_grid_T.nc 
    67 VAR2_SSH=zos      ;   V2SSH_PREF=trunk331   ;   V2SSH_SUFF=_1Y_grid_T.nc 
    68 VAR2_Q=qt         ;   V2Q_PREF=trunk331     ;   V2Q_SUFF=_1Y_grid_T.nc 
    69 VAR2_EMP=wfo      ;   V2EMP_PREF=trunk331   ;   V2EMP_SUFF=_1Y_grid_T.nc 
    70 VAR2_U=uocetr_eff ;   V2U_PREF=trunk331     ;   V2U_SUFF=_1Y_grid_U.nc 
    71 VAR2_V=vocetr_eff ;   V2V_PREF=trunk331     ;   V2V_SUFF=_1Y_grid_V.nc 
    72 VAR2_ICE=sic      ;   V2ICE_PREF=trunk331   ;   V2ICE_SUFF=_1Y_icemod.nc 
     71VAR2_T=thetao      ;   V2T_PREF=trunk331     ;   V2T_SUFF=_1Y_grid_T.nc 
     72VAR2_S=so          ;   V2S_PREF=trunk331     ;   V2S_SUFF=_1Y_grid_T.nc 
     73VAR2_SSH=zos       ;   V2SSH_PREF=trunk331   ;   V2SSH_SUFF=_1Y_grid_T.nc 
     74VAR2_Q=qt          ;   V2Q_PREF=trunk331     ;   V2Q_SUFF=_1Y_grid_T.nc 
     75VAR2_EMP=wfo       ;   V2EMP_PREF=trunk331   ;   V2EMP_SUFF=_1Y_grid_T.nc 
     76VAR2_U=uocetr_eff  ;   V2U_PREF=trunk331     ;   V2U_SUFF=_1Y_grid_U.nc 
     77VAR2_V=vocetr_eff  ;   V2V_PREF=trunk331     ;   V2V_SUFF=_1Y_grid_V.nc 
     78VAR2_ICE=sic       ;   V2ICE_PREF=trunk331   ;   V2ICE_SUFF=_1M_icemod.nc 
     79VAR2_Ithick=sit    ;   V2It_PREF=trunk331    ;   V2It_SUFF=_1M_icemod.nc 
     80VAR2_SNOW=isnowthi ;   V2ICE_SNOW=trunk331   ;   V2SNOW_SUFF=_1M_icemod.nc 
    7381# 
    7482######################### Export Variables ############################### 
     
    8694export FILE_FLUX     VAR_FLUX 
    8795export FILE_MLD      VAR_MLD 
     96export FILE_ICE      VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH  VAR_ICE_area_SH 
     97export FILE_SNOW     VAR_SNOW_NH     VAR_SNOW_SH    VAR_SNOW_area_NH VAR_SNOW_area_SH 
    8898#===================== EXP1 ===================== 
    8999export READ_ONLY_FIRST_RECORD 
    90100# 
    91 export DATE1 DATE2 
    92 export VAR1_T    V1T_PREF    V1T_SUFF 
    93 export VAR1_S    V1S_PREF    V1S_SUFF 
    94 export VAR1_SSH  V1SSH_PREF  V1SSH_SUFF 
    95 export VAR1_Q    V1Q_PREF    V1Q_SUFF 
    96 export VAR1_EMP  V1EMP_PREF  V1EMP_SUFF 
    97 export VAR1_U    V1U_PREF    V1U_SUFF 
    98 export VAR1_V    V1V_PREF    V1V_SUFF 
    99 export VAR1_ICE  V1ICE_PREF  V1ICE_SUFF 
     101export DATE1        DATE2 
     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 
     112 
    100113#===================== EXP2 ===================== 
    101 export DATE1_2 DATE2_2 
    102 export VAR2_T    V2T_PREF    V2T_SUFF 
    103 export VAR2_S    V2S_PREF    V2S_SUFF 
    104 export VAR2_SSH  V2SSH_PREF  V2SSH_SUFF 
    105 export VAR2_Q    V2Q_PREF    V2Q_SUFF 
    106 export VAR2_EMP  V2EMP_PREF  V2EMP_SUFF 
    107 export VAR2_U    V2U_PREF    V2U_SUFF 
    108 export VAR2_V    V2V_PREF    V2V_SUFF 
    109 export VAR2_ICE  V2ICE_PREF  V2ICE_SUFF 
     114export DATE1_2      DATE2_2 
     115export VAR2_T       V2T_PREF    V2T_SUFF 
     116export VAR2_S       V2S_PREF    V2S_SUFF 
     117export VAR2_SSH     V2SSH_PREF  V2SSH_SUFF 
     118export VAR2_Q       V2Q_PREF    V2Q_SUFF 
     119export VAR2_EMP     V2EMP_PREF  V2EMP_SUFF 
     120export VAR2_U       V2U_PREF    V2U_SUFF 
     121export VAR2_V       V2V_PREF    V2V_SUFF 
     122export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
     123export VAR2_Ithick  V2It_PREF   V2It_SUFF 
     124export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
    110125# 
Note: See TracChangeset for help on using the changeset viewer.