Changeset 2751 for trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Q.pro
- 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_Q.pro
r2740 r2751 1 pro ts_mean_Q, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex1 pro std_ts_Q, 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 7 8 ; get exp1 info 9 vq = getenv('VAR1_Q') & prefix = getenv('V1Q_PREF') & suffix = getenv('V1Q_SUFF') 10 ; get exp2 info 11 vq2 = getenv('VAR2_Q') & prefix2 = getenv('V2Q_PREF') & suffix2 = getenv('V2Q_SUFF') 12 ; 13 cdti3 = string(cnt, format = '(i3.3)') 14 print, cdti3 + ') ' + blabla 15 filename = cdti3 + '_ts_Q_'+prefix 16 if prefix NE prefix2 then filename = filename + '_'+prefix2 17 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 18 19 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 20 ; 21 iodir = std_iodir_data 22 ; 23 ts_Q = rseries_ncdf(vq, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 24 title = prefix+' '+d1_d2+'!C'+blabla 25 pltt, ts_Q, 't', -4, 4, date1, date2, /REMPLI $ 26 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 6 27 7 @initorca2 8 9 if KEYWORD_SET(POSTSCRIPT) then begin 10 openps, sEXP1+'_Heat_mean.ps', portrait = 1 11 endif 12 13 ts_Heat = 0. 14 ts_Heat2 = 0. 15 ts_Time = 0. 28 IF prefix NE prefix2 THEN BEGIN 16 29 17 masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct) 30 d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 31 tsave = time 32 ts_Q2 = rseries_ncdf(vq2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 33 time = tsave & IF n_elements(time) NE jpt THEN stop 18 34 35 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 36 pltt, ts_Q.arr - ts_Q2.arr, 't', -4., 4., date1, date2, /REMPLI $ 37 , color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 19 38 20 for year=year1,year2 do begin 21 22 date1 = year * 10000L + 101 23 date2 = year * 10000L + 1231 24 filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1) 25 filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1) 39 ENDIF 26 40 27 if filename eq '' then stop 28 Heat = read_ncdf('qns+qsr',date1 ,date2 , filename = filename, /nostruct) 29 Heat2 = read_ncdf('qns+qsr',date1 ,date2 , filename = filename2, /nostruct) 41 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 42 if KEYWORD_SET(postscript) then closeps 43 44 return 45 end 30 46 31 ts_Heat = [ ts_Heat, moyenne(Heat, 'xy', mask2d = masknp)]32 ts_Heat2 = [ ts_Heat2, moyenne(Heat2, 'xy', mask2d = masknp)]33 ts_Time = [ ts_Time, Time] ; Time = date34 35 endfor36 37 ts_Heat = ts_Heat[1:*] ; remove first record of 038 ts_Heat2 = ts_Heat2[1:*] ; remove first record of 039 Time = ts_Time[1:*] ; remove first record of 040 41 jpt = n_elements(Time)42 43 title = sEXP1+' Global Annual Mean Surface Heat Flux --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)44 pltt,ts_Heat,'t', -4., 4., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $45 , small = [1,2,1], YTITLE = 'W/m2', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex46 47 if sEXP2 ne sEXP1 then begin48 49 title = sEXP1+' - '+sEXP2+' Global Annual Mean Surface Heat Flux --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)50 pltt,ts_Heat - ts_Heat2,'t', -4., 4., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $51 , COLOR = 250, small = [1,2,2], YTITLE = 'W/m2', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex52 53 endif54 55 56 if KEYWORD_SET(POSTSCRIPT) then begin57 closeps58 endif59 60 61 return62 end
Note: See TracChangeset
for help on using the changeset viewer.