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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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)' 
Note: See TracChangeset for help on using the changeset viewer.