- 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_plot_msf.pro
r2733 r2751 1 pro plt_msf, V1, V2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename,_extra = ex1 pro std_plot_msf, V1, V2, SUBBASIN = subbasin, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs 4 4 5 @common 5 @common 6 @std_common 6 7 7 if KEYWORD_SET(SUBBASIN) then subname = subbasin8 9 10 if KEYWORD_SET(POSTSCRIPT) then begin11 openps, '024_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 112 endif13 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 + '_MSF_'+subname+'_'+std_file1_V 12 if std_file1_V NE std_file2_V then filename = filename + '_'+std_file2_V 13 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 14 ; 14 15 if KEYWORD_SET(SUBBASIN) then begin 15 16 CASE 1 of 17 subbasin eq 'Atl' : var = 'atlmsk_nomed' 18 subbasin eq 'Ind' : var = 'indmsk' 19 subbasin eq 'Pac' : var = 'indpacmsk' 20 subbasin eq 'Glo' : var = 'glomsk_nomed' 21 ENDCASE 22 ; 23 msk = read_ncdf( var, filename = mask_filename, /nostruct, _extra = ex ) 24 25 endif else msk = tmask[*,*,0] 16 CASE subname of 17 'Atl' : var = 'atlmsk_nomed' 18 'Ind' : var = 'indmsk' 19 'IndoPac' : var = 'indpacmsk' 20 'GloNoMed' : var = 'glomsk_nomed' 21 ENDCASE 22 msk = read_ncdf( var, filename = std_file_msksub, /nostruct, _extra = ex ) 23 endif else msk = tmask[*, *, 0] 26 24 ; 27 28 29 25 CASE subname OF 30 'Glo':lat_ext=[-80,90]31 'Atl':lat_ext=[-30,90]32 'Ind':lat_ext=[-30,30]33 'Pac':lat_ext=[-30,70]34 26 'GloNoMed':lat_ext = [-80, 90] 27 'Glo':lat_ext = [-80, 90] 28 'Atl':lat_ext = [-30, 90] 29 'Ind':lat_ext = [-30, 30] 30 'IndoPac':lat_ext = [-30, 70] 35 31 ENDCASE 36 37 if subbasin eq 'Pac' then subname = "IndoPac"38 32 39 33 domdef, 0, 6000 40 34 41 mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma) 42 mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma) 43 title=sEXP1+' Meridional Stream Function '+subname+'/ year '+start_end 44 pltz, mm1, -20., 20., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex , FORMAT = '(I3)' $ 45 , small=[1,2,1], COAST_THICK=2 , zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so' 35 mm1 = msf(V1.arr, indexboxzoom = ind, maskout = ma) 36 title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V 37 pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $ 38 , small = [1, 2, 1], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so' 46 39 47 if sEXP1 ne sEXP2 then begin 48 49 title=sEXP1+' - '+sEXP2+subname+'/ year '+start_end 50 pltz, mm1.arr-mm2.arr, -10., 10., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex, FORMAT = '(I3)' $ 51 , small=[1,2,2], COAST_THICK=2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so' 52 40 if std_file1_V NE std_file2_V then begin 41 mm2 = msf(V2.arr, indexboxzoom = ind, maskout = ma) 42 title = title+' - '+std_file2_V 43 pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)' $ 44 , small = [1, 2, 2], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so' 53 45 endif 54 46 55 title = + subname+ 'Meridional Stream Function' 56 text = [ title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/024_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png /> ' ] 57 putfile, '/tmp/prova.txt', text 58 59 if KEYWORD_SET(POSTSCRIPT) then begin 60 closeps 61 endif 62 47 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 48 if KEYWORD_SET(postscript) then closeps 63 49 64 50 return
Note: See TracChangeset
for help on using the changeset viewer.