Changeset 4497 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG
- Timestamp:
- 2014-02-17T15:42:08+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro
r4495 r4497 95 95 96 96 title = prefix+' '+d1_d2+'!C'+'Max Atlantic MOC between 40N and 50N' 97 pltt, t45, 't', 0., 30., date1, date2, /REMPLI, /PORTRAIT$97 pltt, t45, 't', MIN = 0., MAX = 30., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1 $ 98 98 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, _extra = ex 99 99 IF prefix NE prefix2 then begin 100 100 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max Atlantic MOC between 40N and 50N' 101 pltt, t45, 't', 0., 30., date1, date2, /REMPLI, /NOERASE$101 pltt, t45, 't', MIN = 0., MAX = 30., date1, date2, /REMPLI , /NOERASE, XGRIDSTYLE = 1 $ 102 102 , small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 103 pltt, t45_2, 't', 0., 30.,date1, date2, /REMPLI, /NOERASE $103 pltt, t45_2, 't', date1, date2, /REMPLI, /NOERASE $ 104 104 , /ov1d, COLOR = 250, small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; RED 105 105 endif … … 110 110 111 111 title = prefix+' '+d1_d2+'!C'+'Max AntArctic Bottom Water between 80S and 65S '+d1_d2 112 pltt, -t70, 't', 0., 20., date1, date2, /REMPLI, /PORTRAIT, window = 2$112 pltt, -t70, 't', MIN = 0., MAX = 20., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2 $ 113 113 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, _extra = ex 114 114 if prefix NE prefix2 then begin 115 115 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max AntArctic Bottom Water between 80S and 65S' 116 pltt, -t70, 't', 0., 20., date1, date2, /REMPLI, /NOERASE$116 pltt, -t70, 't', MIN = 0.,MAX = 20., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 117 117 , small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 118 pltt, -t70_2, 't', 0., 20.,date1, date2, /REMPLI, /NOERASE $118 pltt, -t70_2, 't', date1, date2, /REMPLI, /NOERASE $ 119 119 , /ov1d, COLOR = 250, small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; RED 120 120 endif … … 125 125 126 126 title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 127 pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3$127 pltt, -t50, 't', MIN = 5., MAX = 30., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 3 $ 128 128 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, _extra = ex 129 129 if prefix NE prefix2 then begin 130 130 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max AntArctic Bottom Cell between 65S and 30N' 131 pltt, -t50 , 't', 5., 30., date1, date2, /REMPLI, /NOERASE$131 pltt, -t50 , 't', MIN = 5., MAX = 30., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 132 132 , small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 133 pltt, -t50_2, 't', 5., 30.,date1, date2, /REMPLI, /NOERASE $133 pltt, -t50_2, 't', date1, date2, /REMPLI, /NOERASE $ 134 134 , /ov1d, COLOR = 250, small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; RED 135 135 endif -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Drake.pro
r4495 r4497 41 41 ; 42 42 title = prefix+' '+d1_d2+'!C'+blabla 43 pltt, drk1, 't', 0., 200., date1, date2, /REMPLI $43 pltt, drk1, 't', 0., 200., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 44 44 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, /PORTRAIT, _extra = ex 45 45 … … 60 60 61 61 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla 62 pltt, drk1 , 't', 0., 200., date1, date2, /REMPLI, /NOERASE $62 pltt, drk1 , 't', 0., 200., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 63 63 , COLOR = 000, small = [1, 2, 2], TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 64 64 pltt, drk2, 't', 0., 200., date1, date2, /REMPLI, /NOERASE $ -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro
r4320 r4497 27 27 ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce ) & ts_EMP.unit = 'Sv' 28 28 title = prefix+' '+d1_d2+'!C'+blabla 29 pltt, ts_EMP, 't', -2., 2., date1, date2, /REMPLI$29 pltt, ts_EMP, 't', MIN = -2.,MAX = 2., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 30 30 , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 31 31 … … 38 38 time = tsave & IF n_elements(time) NE jpt THEN stop 39 39 40 title = prefix+' '+d1_d2+' - '+prefix2+''+d1_d2_2+'!C'+blabla41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', -2., 2., date1, date2, /REMPLI$40 title = prefix+' (BLACK) - ' +prefix2+' (RED) '+d1_d2_2+'!C'+blabla 41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', MIN = -2., MAX = 2., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 42 42 , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 43 43 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro
r4495 r4497 91 91 jpt=12 92 92 time=julday(1,15,1900)+30*lindgen(12) 93 pltt, ICE_N, 't', 0., 15., /REMPLI, /PORTRAIT, MIN = 4., MAX = 16.$93 pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 94 94 , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 95 pltt, ICE_N_15, 't', 0., 15.,/REMPLI, /PORTRAIT, LINESTYLE=2 $ ;;; dashed lines is LINESTYLE=2 $95 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ;;; dashed lines is LINESTYLE=2 $ 96 96 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 97 pltt, vICE_area_NH, 't', 0., 15., /REMPLI, /PORTRAIT$97 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT $ 98 98 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 99 pltt, vICE_ext_NH, 't', 0., 15.,/REMPLI, /PORTRAIT, LINESTYLE=2 $ ;;; dashed lines is LINESTYLE=2 $99 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ;;; dashed lines is LINESTYLE=2 $ 100 100 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 101 101 ; 102 102 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 103 pltt, ICE_S, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = 20.$103 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 104 104 ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 105 pltt, ICE_S_15, 't', 0., 15.,/REMPLI, /PORTRAIT, LINESTYLE=2 $105 pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 106 106 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 107 pltt, vICE_area_SH, 't', 0., 15., /REMPLI, /PORTRAIT$107 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT $ 108 108 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 109 pltt, vICE_ext_SH, 't', 0., 15.,/REMPLI, /PORTRAIT, LINESTYLE=2 $109 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 110 110 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 111 111 ; … … 162 162 jpt=12 163 163 time=julday(1,15,1900)+30*lindgen(12) 164 pltt, ICE_N, 't', -.2, .2, /REMPLI, /PORTRAIT, MIN = 4., MAX = 16., window = 2$164 pltt, ICE_N, 't', MIN = 4, MAX = 16, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 165 165 , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 166 pltt, ICE_N2, 't', -.2, .2, /REMPLI, /PORTRAIT$166 pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $ 167 167 , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; RED 168 pltt, ICE_N_15, 't', -.2, .2, /REMPLI, /PORTRAIT, LINESTYLE=2$ ; linee tratteggiate LINESTYLE=2 $168 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $ 169 169 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 170 pltt, ICE_N2_15, 't', -.2, .2,/REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $170 pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $ 171 171 , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 172 pltt, vICE_area_NH, 't', -2., 2., /REMPLI, /PORTRAIT$172 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT $ 173 173 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; light blue 174 pltt, vICE_ext_NH, 't', -2., 2.,/REMPLI, /PORTRAIT, LINESTYLE=2 $174 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 175 175 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; blu scuro 176 176 ; 177 177 title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (CONTINUS)'+'!C'+ 'and Extend minus 15% (DASHED)' 178 178 ; title ='Southern Hemisphere'+'!C' 179 pltt, ICE_S .arr, 't', -2., 2., /REMPLI, /NOERASE, MIN = 0., MAX = 20.$179 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 180 180 , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 181 pltt, ICE_S2 .arr, 't', -2., 2., /REMPLI, /NOERASE $181 pltt, ICE_S2, 't', /REMPLI, /NOERASE $ 182 182 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 183 pltt, ICE_S_15 .arr , 't', -2., 2.,/REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $183 pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $ 184 184 , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 185 pltt, ICE_S2_15 .arr , 't', -2., 2.,/REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $185 pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2 $ 186 186 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 187 pltt, vICE_area_SH, 't', -2., 2.,/REMPLI, /PORTRAIT $187 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT $ 188 188 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 189 pltt, vICE_ext_SH, 't', -2., 2.,/REMPLI, /PORTRAIT, LINESTYLE=2 $189 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 190 190 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 191 191 ; -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
r4495 r4497 49 49 indy2= indy1 50 50 ; 51 ;OBSERVATIONS : mean seasonal cycle/month 52 vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 53 area_obs = [0.103292, 0.0997500, 0.107625, 0.0944167, 0.0612083, 0.0262500, 0.0262500, 0.0262500, 0.0262500, 0.0843750, 0.0914583, 0.104083] 54 ; 51 55 domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices 52 56 ICE = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) … … 85 89 jpt=12 86 90 time=julday(1,15,1900)+30*lindgen(12) 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 91 ; 91 pltt, ICE_area_export, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = .5 , XGRIDSTYLE = 1 $ 92 , small = [1, 2, 1],YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 93 pltt, area_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4 $ ; light blue 94 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 95 ; 96 tot_area_expo = total(ICE_area_export.arr) 97 ; 98 xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 99 xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 100 xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 92 101 ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 93 102 title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 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 98 ; 103 pltt, ICE_vol_export, 't', /REMPLI, MIN = 0., MAX = .5 , /NOERASE, XGRIDSTYLE = 1 $ 104 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 105 pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4 $ ; light blue 106 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month',DATE_FORMAT = '%M', TITLE = title, _extra = ex 107 ; 108 tot_vol_expo = total(ICE_vol_export.arr) 109 tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 110 xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 111 xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 112 xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 113 xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 114 99 115 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 100 116 if KEYWORD_SET(postscript) then closeps … … 140 156 jpt=12 141 157 time=julday(1,15,1900)+30*lindgen(12) 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 150 ; 158 pltt, ICE_area_export, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, XGRIDSTYLE = 1, window = 2 $ 159 , small = [1, 2, 1], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 160 pltt, ICE_area_export_2 ,'t', /REMPLI, /PORTRAIT, /NOERASE $ 161 , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 162 pltt, area_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4 $ ; light blue 163 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 164 ; 165 tot_area_expo = total(ICE_area_export.arr) 166 tot_area_expo_2 = total(ICE_area_export_2.arr) 167 ; 168 xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 169 xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 170 xyouts, julday(6,15,1900), 0.44, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 171 xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 151 172 title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 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 160 ; 173 pltt, ICE_vol_export, 't', /REMPLI, MIN = 0., MAX = .5, /NOERASE, XGRIDSTYLE = 1 $ 174 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 175 pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE $ 176 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 177 pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4 $ ; light blue 178 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 179 ; 180 tot_vol_expo = total(ICE_vol_export.arr) 181 tot_vol_expo_2 = total(ICE_vol_export_2.arr) 182 tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 183 tot_vol_expo_2_Sv = tot_vol_expo_2 * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 184 xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 185 xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 186 xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 187 xyouts, julday(6,15,1900), 0.42, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 188 xyouts, julday(9,15,1900), 0.40, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 189 xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 161 190 162 191 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_ICE_Vol.pro
r4495 r4497 65 65 jpt=12 66 66 time=julday(1,15,1900)+30*lindgen(12) 67 pltt, ICE_vol_N, 't', 0., 15., /REMPLI, /PORTRAIT,MIN = 0., MAX = 30000.$67 pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 68 68 , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 69 69 ; 70 70 title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 71 pltt, ICE_vol_S, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = 11000.$71 pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 72 72 , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 73 73 ; … … 112 112 jpt=12 113 113 time=julday(1,15,1900)+30*lindgen(12) 114 pltt, ICE_vol_N, 't', -125., -60., /REMPLI, /PORTRAIT, MIN = 0., MAX = 30000., window = 2$114 pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 115 115 , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 116 pltt, ICE_vol_N2 , 't', -.2, .2,/REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $116 pltt, ICE_vol_N2 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 117 117 , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 118 118 ; 119 119 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 120 pltt, ICE_vol_S, 't', -620., -10., /REMPLI, /NOERASE, MIN = 0., MAX = 11000.$120 pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 121 121 , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 122 pltt, ICE_vol_S2, 't', -620., -10.,/REMPLI, /NOERASE $122 pltt, ICE_vol_S2, 't', /REMPLI, /NOERASE $ 123 123 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 124 124 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Q.pro
r4495 r4497 23 23 ts_Q = rseries_ncdf(vq, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 24 24 title = prefix+' '+d1_d2+'!C'+blabla 25 pltt, ts_Q, 't', -4, 4, date1, date2, /REMPLI$25 pltt, ts_Q, 't', MIN = -4, MAX = 4, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 26 26 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 27 27 … … 34 34 35 35 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla 36 pltt, ts_Q.arr , 't', -4., 4., date1, date2, /REMPLI$36 pltt, ts_Q.arr , 't', MIN = -4., MAX = 4., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 37 37 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 38 pltt, ts_Q2.arr, 't', -4., 4.,date1, date2, /REMPLI $38 pltt, ts_Q2.arr, 't', date1, date2, /REMPLI $ 39 39 , /ov1d, color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 40 40 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_S.pro
r4495 r4497 35 35 36 36 title = prefix+' '+d1_d2+'!C'+blabla 37 pltt, ts_Sal, 't', 34.54, 34.76, date1, date2, /REMPLI$37 pltt, ts_Sal, 't', MIN = 34.54, MAX = 34.76, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 38 38 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 39 39 40 40 title = prefix+' '+d1_d2+' - Levitus!C'+blabla 41 pltt, ts_Sal - Lev_xyz, 't', -.05, .05, date1, date2, /REMPLI$41 pltt, ts_Sal - Lev_xyz, 't', MIN = -.05, MAX = .05, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 42 42 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex 43 43 … … 47 47 48 48 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 49 pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so'$49 pltt, ts_z, 'zt', MIN = -.5, MAX = .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 50 50 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 51 51 52 52 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 53 pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so'$53 pltt, ts_z, 'zt', MIN = -.5, MAX = .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 54 54 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 55 55 … … 69 69 70 70 title = prefix+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 71 pltt, ts_Sal - ts_Sal2, 't', -0.1, 0.1, date1, date2, /REMPLI$71 pltt, ts_Sal - ts_Sal2, 't', MIN = -0.1, MAX = 0.1, date1, date2, /REMPLI , XGRIDSTYLE = 1 $ 72 72 , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 73 73 74 74 title = prefix+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla 75 pltt, ts_z - ts_z2, 'zt', -0.1, 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so'$75 pltt, ts_z - ts_z2, 'zt', MIN = -0.1, MAX = 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 76 76 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 77 77 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro
r4495 r4497 22 22 ; 23 23 ts_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 24 ;SF ts_SSH.arr = ts_SSH.arr * 100 & ts_SSH.unit = 'cm'25 24 ts_SSH.arr = ts_SSH.arr & ts_SSH.unit = 'm' 26 25 title = prefix+' '+d1_d2+'!C'+blabla 27 pltt, ts_SSH, 't', -12., 12., date1, date2, /REMPLI $ 28 ;SF , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 26 pltt, ts_SSH, 't', MIN = -12., MAX = 12., date1, date2, /REMPLI , XGRIDSTYLE = 1 $ 29 27 , small = [1, 2, 1], YTITLE = 'm', TITLE = title, /PORTRAIT, _extra = ex 30 28 … … 34 32 tsave = time 35 33 ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 36 ;SF ts_SSH2.arr = ts_SSH2.arr * 100 & ts_SSH2.unit = 'cm'37 34 ts_SSH2.arr = ts_SSH2.arr & ts_SSH2.unit = 'm' 38 35 time = tsave & IF n_elements(time) NE jpt THEN stop 39 36 40 37 title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla 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 38 pltt, ts_SSH.arr - ts_SSH2.arr, 't', MIN = -10., MAX = 10., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 43 39 , color = 250, small = [1, 2, 2], YTITLE = 'm', TITLE = title, /noerase, _extra = ex 44 40 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro
r4495 r4497 35 35 36 36 title = prefix+' '+d1_d2+'!C'+blabla 37 pltt, ts_Temp, 't', 1., 4.5, date1, date2, /REMPLI$37 pltt, ts_Temp, 't', MIN = 1., MAX = 4.5, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 38 38 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 39 39 40 40 title = prefix+' '+d1_d2+' - Levitus!C'+blabla 41 pltt, ts_Temp - Lev_xyz, 't', -1., 1., date1, date2, /REMPLI$41 pltt, ts_Temp - Lev_xyz, 't', MIN = -1., MAX = 1., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 42 42 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex 43 43 … … 47 47 48 48 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 49 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so'$49 pltt, ts_z, 'zt', MIN = -2.,MAX = 2., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 50 50 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 51 51 52 52 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 53 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so'$53 pltt, ts_z, 'zt', MIN = -2., MAX = 2., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 54 54 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 55 55 … … 69 69 70 70 title = prefix+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 71 pltt, ts_Temp - ts_Temp2, 't', -1., 1., date1, date2, /REMPLI$71 pltt, ts_Temp - ts_Temp2, 't', MIN = -1., MAX = 1., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 72 72 , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 73 73 74 74 title = prefix+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla 75 pltt, ts_z - ts_z2, 'zt', -1., 1., inter = .1, date1, date2, /REMPLI, style = 'so0so'$75 pltt, ts_z - ts_z2, 'zt', MIN = -1., MAX = 1., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $ 76 76 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 77 77
Note: See TracChangeset
for help on using the changeset viewer.