| 1 | pro std_plot_zonal_S, S1, S2, SLev, SUBBASIN = subbasin, POSTSCRIPT = postscript, _extra = ex |
|---|
| 2 | |
|---|
| 3 | compile_opt idl2, strictarrsubs |
|---|
| 4 | |
|---|
| 5 | @common |
|---|
| 6 | @std_common |
|---|
| 7 | |
|---|
| 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 + '_ZonalS_'+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) |
|---|
| 29 | ; from 2D array to 3D array |
|---|
| 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 |
|---|
| 33 | ; |
|---|
| 34 | title = 'Salinity, '+subname+'!C'+std_file1_T |
|---|
| 35 | pltz, S1, MININ = 31., MAXIN = 36.8, INTER = .2, typein = 'yz', FORMAT = '(f4.1)' $ |
|---|
| 36 | , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT, _extra = ex |
|---|
| 37 | ; |
|---|
| 38 | if std_file1_T NE std_file2_T then begin |
|---|
| 39 | title = title+' - '+std_file2_T |
|---|
| 40 | pltz, S1.arr-S2.arr, MIN = -1., MAX = 1., INTER = .1, 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, S1.arr-SLev.arr, MININ = -1, MAXIN = 1., INTER = .1, typein = 'yz', FORMAT = '(f4.1)', style = 'so0so' $ |
|---|
| 45 | , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex |
|---|
| 46 | endelse |
|---|
| 47 | ; |
|---|
| 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 |
|---|
| 55 | end |
|---|
