[2751] | 1 | pro std_plot_zonal_T, T1, T2, TLev, SUBBASIN = subbasin, POSTSCRIPT = postscript, _extra = ex |
---|
| 2 | |
---|
[2709] | 3 | compile_opt idl2, strictarrsubs |
---|
| 4 | |
---|
[2751] | 5 | @common |
---|
| 6 | @std_common |
---|
[2709] | 7 | |
---|
[2751] | 8 | cdti3 = string(cnt, format = '(i3.3)') |
---|
| 9 | print, cdti3 + ') ' + blabla |
---|
| 10 | if KEYWORD_SET(SUBBASIN) then subname = subbasin else subname = "Glo" |
---|
| 11 | filename = cdti3 + '_ZonalT_'+subname+'_'+std_file1_T |
---|
| 12 | if std_file1_T NE std_file2_T then filename = filename + '_'+std_file2_T |
---|
| 13 | if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 |
---|
| 14 | ; |
---|
| 15 | ; looking for longitudinal index corresponding to the highest latitude (closest to the North Pole) |
---|
| 16 | index = where(gphit eq max(gphit)) |
---|
| 17 | nx = index[0] mod jpi |
---|
| 18 | gphi_save = gphit |
---|
| 19 | gphit[0, *] = gphit[nx, *] |
---|
| 20 | ; |
---|
| 21 | if KEYWORD_SET(SUBBASIN) then begin |
---|
| 22 | tmasksv = tmask |
---|
| 23 | CASE 1 of |
---|
| 24 | subbasin eq 'Atl' : var = 'atlmsk_nomed' |
---|
| 25 | subbasin eq 'Ind' : var = 'indmsk_nored' |
---|
| 26 | subbasin eq 'Pac' : var = 'pacmsk' |
---|
| 27 | ENDCASE |
---|
| 28 | msk = read_ncdf( var, filename = std_file_msksub, _extra = ex) |
---|
[2709] | 29 | ; from 2D array to 3D array |
---|
[2751] | 30 | msk = msk.arr[*]#replicate(1., nzt) |
---|
| 31 | tmask[firstxt:lastxt, firstyt:lastyt, firstzt:lastzt] = tmask[firstxt:lastxt, firstyt:lastyt, firstzt:lastzt] * msk |
---|
| 32 | endif |
---|
[2709] | 33 | ; |
---|
[2751] | 34 | title = 'Temperature, '+subname+'!C'+std_file1_T |
---|
| 35 | pltz, T1, MININ = -2., MAXIN = 30., INTER = 1., typein = 'yz', FORMAT = '(I2)' $ |
---|
| 36 | , small = [ 1, 2, 1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT, _extra = ex |
---|
[2709] | 37 | ; |
---|
[2751] | 38 | if std_file1_T NE std_file2_T then begin |
---|
| 39 | title = title+' - '+std_file2_T |
---|
| 40 | pltz, T1.arr-T2.arr, MIN = -2., MAX = 2., INTER = .2, typein = 'yz', FORMAT = '(f4.1)', style = 'so0so' $ |
---|
| 41 | , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex |
---|
| 42 | endif else begin |
---|
| 43 | title = title+' - Levitus' |
---|
| 44 | pltz, T1.arr-TLev.arr, MININ = -4., MAXIN = 4., INTER = .5, typein = 'yz', FORMAT = '(I2)', style = 'so0so' $ |
---|
| 45 | , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex |
---|
| 46 | endelse |
---|
[2709] | 47 | ; |
---|
[2751] | 48 | htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] |
---|
| 49 | if KEYWORD_SET(postscript) then closeps |
---|
| 50 | ; |
---|
| 51 | if KEYWORD_SET(SUBBASIN) then tmask = tmasksv |
---|
| 52 | gphit = gphi_save |
---|
| 53 | ; |
---|
| 54 | return |
---|
[2709] | 55 | end |
---|