source: trunk/CONFIG/ORCA2_LIM/IDL_scripts/plt_msf.pro @ 2150

Last change on this file since 2150 was 2150, checked in by flavoni, 10 years ago

commit idl tools for post-processing, for ORCA2_LIM configuration, see ticket #724

File size: 2.0 KB
Line 
1pro plt_msf, V1, V2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename,_extra = ex
2
3  compile_opt idl2, strictarrsubs
4
5  @common
6
7  if KEYWORD_SET(SUBBASIN) then subname = subbasin
8
9
10  if KEYWORD_SET(POSTSCRIPT) then begin
11     openps, '019_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps',  portrait = 1
12  endif
13
14  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]
26;
27
28 
29  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     
35  ENDCASE
36
37 if subbasin eq 'Pac' then subname = "IndoPac"
38 
39  domdef, 0, 6000
40
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'
46
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     
53  endif
54 
55  title = + subname+ 'Meridional Stream Function'
56  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/019_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
63
64  return
65end
66
Note: See TracBrowser for help on using the repository browser.