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

update idl scripts for ORCA2 configurations, see ticket #1115

File:
1 edited

Legend:

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