[2751] | 1 | pro std_ts_T, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex |
---|
[2709] | 2 | |
---|
| 3 | compile_opt idl2, strictarrsubs |
---|
| 4 | |
---|
[2751] | 5 | @common |
---|
| 6 | @std_common |
---|
[2709] | 7 | |
---|
[2751] | 8 | ; get exp1 info |
---|
| 9 | vtemp = getenv('VAR1_T') & prefix = getenv('V1T_PREF') & suffix = getenv('V1T_SUFF') |
---|
| 10 | vssh = getenv('VAR1_SSH') & sshprefix = getenv('V1SSH_PREF') & sshsuffix = getenv('V1SSH_SUFF') |
---|
| 11 | ; get exp2 info |
---|
| 12 | vtemp2 = getenv('VAR2_T') & prefix2 = getenv('V2T_PREF') & suffix2 = getenv('V2T_SUFF') |
---|
| 13 | vssh2 = getenv('VAR2_SSH') & sshprefix2 = getenv('V2SSH_PREF') & sshsuffix2 = getenv('V2SSH_SUFF') |
---|
| 14 | ; |
---|
| 15 | cdti3 = string(cnt, format = '(i3.3)') |
---|
| 16 | print, cdti3 + ') ' + blabla |
---|
| 17 | filename = cdti3 + '_ts_T_'+prefix |
---|
| 18 | if prefix NE prefix2 then filename = filename + '_'+prefix2 |
---|
| 19 | if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1 |
---|
[2709] | 20 | |
---|
[2751] | 21 | d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' |
---|
[2709] | 22 | |
---|
[2751] | 23 | ; read levitus data |
---|
| 24 | std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) |
---|
| 25 | Lev = read_ncdf(getenv('VAR_TEMP_3D'), filename = std_file_Levitus_T, /nostruct ) |
---|
[2709] | 26 | Lev_xyz = moyenne(Lev, 'xyz', mask2d = masknp) |
---|
[2751] | 27 | levz = moyenne(temporary(Lev), 'xy', mask2d = masknp, /KEEPBOTTOM) |
---|
| 28 | ; |
---|
| 29 | iodir = std_iodir_data |
---|
| 30 | ; |
---|
[2709] | 31 | |
---|
[2751] | 32 | ; read exp1 data |
---|
| 33 | std_ts_read, vtemp, date1, date2, prefix, suffix, ts_Temp, ts_z $ |
---|
| 34 | , WITHSSH = vssh, SSHPREFIX = sshprefix, SSHSUFFIX = sshsuffix, LEVZ = levz |
---|
[2709] | 35 | |
---|
[2751] | 36 | title = prefix+' '+d1_d2+'!C'+blabla |
---|
| 37 | pltt, ts_Temp, 't', 1., 4.5, date1, date2, /REMPLI $ |
---|
| 38 | , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex |
---|
[2709] | 39 | |
---|
[2751] | 40 | title = prefix+' '+d1_d2+' - Levitus!C'+blabla |
---|
| 41 | pltt, ts_Temp - Lev_xyz, 't', -1., 1., date1, date2, /REMPLI $ |
---|
| 42 | , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex |
---|
[2709] | 43 | |
---|
[2751] | 44 | htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png /> ' ] |
---|
| 45 | if KEYWORD_SET(postscript) then closeps |
---|
| 46 | if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 |
---|
[2709] | 47 | |
---|
[2751] | 48 | title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla |
---|
| 49 | pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $ |
---|
| 50 | , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex |
---|
[2709] | 51 | |
---|
[2751] | 52 | title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla |
---|
| 53 | pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $ |
---|
| 54 | , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex |
---|
[2709] | 55 | |
---|
[2751] | 56 | htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png /> ' ] |
---|
| 57 | if KEYWORD_SET(postscript) then closeps |
---|
[2709] | 58 | |
---|
[2751] | 59 | IF prefix NE prefix2 THEN BEGIN |
---|
[2709] | 60 | |
---|
[2751] | 61 | d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' |
---|
| 62 | ; read exp2 data |
---|
| 63 | tsave = time |
---|
| 64 | std_ts_read, vtemp2, date1_2, date2_2, prefix2, suffix2, ts_Temp2, ts_z2 $ |
---|
| 65 | , WITHSSH = vssh2, SSHPREFIX = sshprefix2, SSHSUFFIX = sshsuffix2, LEVZ = levz |
---|
| 66 | time = tsave & IF n_elements(time) NE jpt THEN stop |
---|
[2709] | 67 | |
---|
[2751] | 68 | if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 |
---|
[2709] | 69 | |
---|
[2751] | 70 | title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla |
---|
| 71 | pltt, ts_Temp - ts_Temp2, 't', -1., 1., date1, date2, /REMPLI $ |
---|
| 72 | , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex |
---|
[2709] | 73 | |
---|
[2751] | 74 | title = prefix+' '+d1_d2+' - '+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' $ |
---|
| 76 | , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex |
---|
| 77 | |
---|
| 78 | htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_3.png /> ' ] |
---|
| 79 | if KEYWORD_SET(postscript) then closeps |
---|
| 80 | |
---|
| 81 | ENDIF |
---|
| 82 | |
---|
[2709] | 83 | return |
---|
| 84 | end |
---|