- Timestamp:
- 2011-04-26T10:45:52+02:00 (13 years ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro
r2740 r2751 1 pro ts_mean_EMP, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex1 pro std_ts_EMP, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs 4 4 5 @common 5 @common 6 @std_common 6 7 7 @initorca2 8 9 if KEYWORD_SET(POSTSCRIPT) then begin 10 openps, sEXP1+'_EMP_mean.ps', portrait = 1 11 endif 8 ; get exp1 info 9 vemp = getenv('VAR1_EMP') & prefix = getenv('V1EMP_PREF') & suffix = getenv('V1EMP_SUFF') 10 ; get exp2 info 11 vemp2 = getenv('VAR2_EMP') & prefix2 = getenv('V2EMP_PREF') & suffix2 = getenv('V2EMP_SUFF') 12 ; 13 cdti3 = string(cnt, format = '(i3.3)') 14 print, cdti3 + ') ' + blabla 15 filename = cdti3 + '_ts_EMP_'+prefix 16 if prefix NE prefix2 then filename = filename + '_'+prefix2 17 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 12 18 13 ts_EMP = 0. 14 ts_EMP2 = 0. 15 ts_Time = 0. 19 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 20 ; 21 iodir = std_iodir_data 22 ; 23 surf_oce = e1t * e2t * tmask[*,*,0] * masknp 24 surf_oce = total(surf_oce) 25 ; 26 ts_EMP = rseries_ncdf(vemp, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 27 ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce ) & ts_EMP.unit = 'Sv' 28 title = prefix+' '+d1_d2+'!C'+blabla 29 pltt, ts_EMP, 't', -.001, .001, date1, date2, /REMPLI $ 30 , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 31 32 IF prefix NE prefix2 THEN BEGIN 16 33 17 masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct) 34 d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 35 tsave = time 36 ts_EMP2 = rseries_ncdf(vemp2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 37 ts_EMP2.arr = ts_EMP2.arr * ( 1.E-09 * surf_oce ) & ts_EMP2.unit = 'Sv' 38 time = tsave & IF n_elements(time) NE jpt THEN stop 18 39 19 for year=year1,year2 do begin 40 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', -.001, .001, date1, date2, /REMPLI $ 42 , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 20 43 21 date1 = year * 10000L + 101 22 date2 = year * 10000L + 1231 23 filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1) 24 filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1) 44 ENDIF 25 45 26 if filename eq '' then stop 27 EMP = read_ncdf('wfo',date1 ,date2 , filename = filename, /nostruct) 28 EMP2 = read_ncdf('wfo',date1 ,date2 , filename = filename2, /nostruct) 29 30 ts_EMP = [ ts_EMP, moyenne(EMP, 'xy', mask2d = masknp)] 31 ts_EMP2 = [ ts_EMP2, moyenne(EMP2, 'xy', mask2d = masknp)] 32 ts_Time = [ ts_Time, Time] ; Time = date 33 34 endfor 35 36 surf_oce = e1t * e2t * tmask[*,*,0] 37 surf_oce = total(surf_oce) 38 39 ts_EMP = ts_EMP[1:*] ; remove first record of 0 40 ts_EMP = ts_EMP * 1.E-09 * surf_oce ; to have Sv 41 42 ts_EMP2 = ts_EMP2[1:*] ; remove first record of 0 43 ts_EMP2 = ts_EMP2 * 1.E-09 * surf_oce ; to have Sv 44 45 varunit= 'Sv' 46 vol_oce = moyenne( replicate(1., jpi, jpj, jpk), /integration, 'xyz') 47 48 Time = ts_Time[1:*] ; remove first record of 0 49 jpt = n_elements(Time) 50 51 title = sEXP1+' Global Annual EMP Budget (Sv) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 52 pltt,ts_EMP,'t', -.01, 0.00, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $ 53 , small = [1,2,1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 54 55 if sEXP2 ne sEXP1 then begin 56 57 title = sEXP1+' - '+sEXP2+' Global Annual EMP Budget (Sv) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 58 pltt,ts_EMP - ts_EMP2,'t', -.002, .002, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI, /NOERASE $ 59 , COLOR = 250, small = [1,2,2], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 60 endif 61 62 if KEYWORD_SET(POSTSCRIPT) then begin 63 closeps 64 endif 65 66 46 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 47 if KEYWORD_SET(postscript) then closeps 48 67 49 return 68 50 end 51
Note: See TracChangeset
for help on using the changeset viewer.