- 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_SSH.pro
r2740 r2751 1 pro ts_mean_SSH, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex1 pro std_ts_SSH, 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 vssh = getenv('VAR1_SSH') & prefix = getenv('V1SSH_PREF') & suffix = getenv('V1SSH_SUFF') 10 ; get exp2 info 11 vssh2 = getenv('VAR2_SSH') & prefix2 = getenv('V2SSH_PREF') & suffix2 = getenv('V2SSH_SUFF') 12 ; 13 cdti3 = string(cnt, format = '(i3.3)') 14 print, cdti3 + ') ' + blabla 15 filename = cdti3 + '_ts_SSH_'+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_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 24 ts_SSH.arr = ts_SSH.arr * 100 & ts_SSH.unit = 'cm' 25 title = prefix+' '+d1_d2+'!C'+blabla 26 pltt, ts_SSH, 't', -.1, .1, date1, date2, /REMPLI $ 27 , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 6 28 7 @initorca2 29 IF prefix NE prefix2 THEN BEGIN 30 31 d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 32 tsave = time 33 ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 34 ts_SSH2.arr = ts_SSH2.arr * 100 & ts_SSH2.unit = 'cm' 35 time = tsave & IF n_elements(time) NE jpt THEN stop 36 37 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 38 pltt, ts_SSH.arr - ts_SSH2.arr, 't', -.1, .1, date1, date2, /REMPLI $ 39 , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 40 41 ENDIF 42 43 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 44 if KEYWORD_SET(postscript) then closeps 8 45 9 if KEYWORD_SET(POSTSCRIPT) then begin10 openps, sEXP1+'_SSH_mean.ps', portrait = 111 endif12 13 ts_SSH = 0.14 ts_SSH2 = 0.15 ts_Time = 0.16 17 masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct)18 19 for year=year1,year2 do begin20 21 date1 = year * 10000L + 10122 date2 = year * 10000L + 123123 filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)24 filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)25 26 SSH = read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)27 SSH2 = read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)28 29 ts_SSH = [ ts_SSH, moyenne(SSH, 'xy', mask2d = masknp)]30 ts_SSH2 = [ ts_SSH2, moyenne(SSH2, 'xy', mask2d = masknp)]31 ts_Time = [ ts_Time, Time] ; Time = date32 33 endfor34 35 ts_SSH = ts_SSH[1:*] ; remove first record of 036 ts_SSH2 = ts_SSH2[1:*] ; remove first record of 037 Time = ts_Time[1:*] ; remove first record of 038 39 jpt = n_elements(Time)40 41 title = sEXP1+' Global Annual Mean SSH --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)42 pltt,ts_SSH,'t', -.001, .001, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $43 , small = [1,2,1], YTITLE = 'm', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex44 45 if sEXP2 ne sEXP1 then begin46 47 title = sEXP1+' - '+sEXP2+' Global Annual Mean SSH --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)48 pltt, ts_SSH - ts_SSH2,'t', -.0001, .0001, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $49 , COLOR = 250, small = [1,2,2], YTITLE = 'm', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex50 51 endif52 53 54 55 if KEYWORD_SET(POSTSCRIPT) then begin56 closeps57 endif58 59 60 46 return 61 47 end 48
Note: See TracChangeset
for help on using the changeset viewer.