Changeset 4495 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
- Timestamp:
- 2014-02-10T16:51:42+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
r4395 r4495 71 71 ICE_area_export = reform(ICE_area_export, 12, nyr) 72 72 ICE_area_export = total(ICE_area_export,2)/nyr 73 ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} 73 ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} ; annual mean 74 ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'} ; monthly mean 75 74 76 ; 75 77 ICE_vol_export = reform(ICE_vol_export, 12, nyr) 76 78 ICE_vol_export = total(ICE_vol_export,2)/nyr 77 ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} 79 ; ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 80 ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 78 81 79 82 ; … … 82 85 jpt=12 83 86 time=julday(1,15,1900)+30*lindgen(12) 84 pltt, ICE_area_export, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 2. $ 85 , small = [1, 2, 1],YTITLE = '10^6 Km^2/year', TITLE = title, _extra = ex 87 ;; pltt, ICE_area_export, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = .5 $ 88 pltt, ICE_area_export, 't', 0., 15., /REMPLI, /PORTRAIT,MIN = 0., MAX = .5 $ 89 ;; , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 90 , small = [1, 2, 1],YTITLE = '10^6 Km^2/month', TITLE = title, _extra = ex 86 91 ; 87 92 ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 88 93 title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 89 pltt, ICE_vol_export, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 6. $ 90 , small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 94 ;; pltt, ICE_vol_export, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 95 pltt, ICE_vol_export, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 96 ;; , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 97 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 91 98 ; 92 99 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] … … 119 126 ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 120 127 ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 121 ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} 128 ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} ; annual mean 129 ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'} ; monthly mean 122 130 ; 123 131 ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 124 132 ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 125 ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} 133 ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 134 ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 126 135 ; 127 136 ; … … 131 140 jpt=12 132 141 time=julday(1,15,1900)+30*lindgen(12) 133 pltt, ICE_area_export.arr, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = 2., window = 2 $ 134 , small = [1, 2, 1],YTITLE = '10^6 Km^2/year', TITLE = title, _extra = ex 135 pltt, ICE_area_export_2.arr ,'t', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ 136 , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/year', TITLE = title, _extra = ex 137 ; 142 ;; pltt, ICE_area_export.arr, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, window = 2 $ 143 pltt, ICE_area_export, 't', 0., 15., /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, window = 2 $ 144 ;; , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 145 , small = [1, 2, 1], YTITLE = '10^6 Km^2/month', TITLE = title, _extra = ex 146 ;; pltt, ICE_area_export_2.arr ,'t', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ 147 pltt, ICE_area_export_2 ,'t', 0., 15., /REMPLI, /PORTRAIT $ 148 ;; , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 149 , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/month', TITLE = title, _extra = ex 138 150 ; 139 151 title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 140 pltt, ICE_vol_export.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 6. $ 141 , small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 142 pltt, ICE_vol_export_2.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE $ 143 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 152 ;; pltt, ICE_vol_export.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 153 pltt, ICE_vol_export, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 154 ;; , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 155 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 156 ;; pltt, ICE_vol_export_2.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE $ 157 pltt, ICE_vol_export_2, 't', 0., 15., /REMPLI, /NOERASE $ 158 ;; , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 159 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 144 160 ; 145 161
Note: See TracChangeset
for help on using the changeset viewer.