Changeset 3999 for branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro
- Timestamp:
- 2013-08-06T13:39:29+02:00 (11 years ago)
- 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 26 26 v2_Ithick = getenv('VAR2_Ithick') & prefix2 = getenv('V2It_PREF') & suffix2 = getenv('V2It_SUFF') 27 27 28 ; get ice VOLUME climatology info29 ;;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) ) * 1234 ; time_ice = round(time_ice) ; round to nearest integer35 ; t1 = where(time_ice eq 0)36 ; t1 = t1[0] ; jannuary37 ; t2 = where(time_ice eq 11, count)38 ; t2 = t2[count-1] ; last day of december39 ; nyear = (t2-t1+1)/1240 ; 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 ;46 28 cdti3 = string(cnt, format = '(i3.3)') 47 29 print, cdti3 + ') ' + blabla 48 filename = cdti3 + '_ts_ICE_ '+prefix30 filename = cdti3 + '_ts_ICE_Vol'+prefix 49 31 if prefix NE prefix2 then filename = filename + '_'+prefix2 50 32 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 33 ; 51 34 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 35 ; 52 36 iodir = std_iodir_data 53 37 ; ICE Area(=Surface) in NORTH Hemisphere … … 58 42 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 59 43 ICE_vol_N = grossemoyenne(ICE_vol_N, 'xy', /integration, mask2d = masknp) 44 ; 60 45 if jpt mod 12 ne 0 then stop 61 46 nyr=jpt/12. … … 63 48 ICE_vol_N = total(ICE_vol_N,2)/nyr 64 49 ICE_vol_N = {arr:ICE_vol_N * 1.e-9, unit : '10^9 Km^3'} 50 ; 65 51 ;ICE Area(=Surface) in SOUTH Hemisphere 66 52 domdef, 0, jpi-1, -90, -30, /xindex … … 76 62 ICE_vol_S = {arr:ICE_vol_S * 1.e-9, unit : '10^9 Km^3'} 77 63 ; 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)' 79 65 jpt=12 80 66 time=julday(1,15,1900)+30*lindgen(12) 81 67 pltt, ICE_vol_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 30000. $ 82 , small = [1, 2, 1], YTITLE = ' 10^9Km^3 ', TITLE = title, _extra = ex68 , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 83 69 ; 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)' 85 71 pltt, ICE_vol_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 11000. $ 86 , small = [1, 2, 2], YTITLE = ' 10^9Km^3 ', TITLE = title, _extra = ex72 , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 87 73 ; 88 89 74 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 90 75 if KEYWORD_SET(postscript) then closeps … … 96 81 domdef, 0, jpi-1, 30, 90, /xindex 97 82 ;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) 100 85 ; Volume = Area(=Surface) * Thickness 101 86 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 … … 109 94 ;ICE Area(=Surface) in SOUTH Hemisphere 110 95 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) 113 98 ; Volume = Area(=Surface) * Thickness 114 99 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 … … 120 105 ICE_vol_S2 = {arr:ICE_vol_S2 * 1.e-9, unit : '10^3 Km^3'} 121 106 122 time = tsave & IF n_elements(time) NE jpt THEN stop107 ; time = tsave & IF n_elements(time) NE jpt THEN stop 123 108 124 if KEYWORD_SET(postscript) then openps, filename+' .ps', portrait = 1109 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 125 110 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 133 120 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 /> ' ] 135 122 if KEYWORD_SET(postscript) then closeps 136 123
Note: See TracChangeset
for help on using the changeset viewer.