source: trunk/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_depth.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: 1.8 KB
Line 
1pro xxx_Med_S_depth, S1, S2, SLev, sEXP1, sEXP2, start_end, DEPTH=depth, LAT=lat, POSTSCRIPT = postscript,  _extra = ex
2 
3  compile_opt idl2, strictarrsubs
4
5   @common
6
7
8   if KEYWORD_SET(POSTSCRIPT) then begin
9      openps, '024_Med_S_depth_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.ps'   
10   endif
11;domdef, 0, 6000
12
13;domdef,270,365,5,70
14   minval = min(abs(gdept -  depth),index)
15   
16   title = sEXP1+'  Mediterranean water at ' +strtrim(lat,1)+'°N / year '+start_end
17   pltz, S1.arr, MININ = 35., MAXIN = 37., INTER = 0.1,  UNIT = varunit $
18         , small = [1,2,1], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., /PORTRAIT, _extra = ex
19;
20
21   if sEXP1 ne sEXP2 then begin
22     
23      title = sEXP1+' - '+sEXP2+' Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
24      pltz, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05, UNIT = varunit $
25            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000.,  _extra = ex
26     
27   endif else begin
28   
29      title = sEXP1+' - Levitus Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
30      pltz, S1.arr-SLev.arr, MININ = -1., MAXIN = 1., INTER=0.1, UNIT = varunit $
31            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., _extra = ex
32   endelse
33
34
35   title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
36   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/024_Med_Levitus_S_depth_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.png  />  ' ]
37   putfile, '/tmp/prova.txt', text
38
39   if KEYWORD_SET(POSTSCRIPT) then begin
40      closeps
41   endif
42
43
44   return
45end
Note: See TracBrowser for help on using the repository browser.