source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_tongue.pro @ 2711

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

add correct unit in idl plots, see ticket 724

File size: 1.8 KB
Line 
1pro xxx_Med_S_tongue, 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, '028_Med_S_tongue_'+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
14minval = min(abs(gdept -  depth),index)
15;
16varunit = S1.unit
17
18   title = sEXP1+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
19   plt, S1.arr, MININ=35., MAXIN=37., INTER=0.1,  UNIT= varunit $
20         , small=[1,2,1], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
21;
22  if sEXP1 ne sEXP2 then begin
23     
24     title = sEXP1+' - '+sEXP2+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
25     plt, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05,  UNIT = varunit $
26          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
27 
28  endif else begin
29     
30     title = sEXP1+' - Levitus Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
31     plt, S1.arr - SLev.arr, MININ=-1., MAXIN=1., INTER=0.1,  UNIT = varunit $
32          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
33     
34  endelse
35
36  title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
37  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/028_Med_S_tongue_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.png  />  ' ]
38  putfile, '/tmp/prova.txt', text
39   
40  if KEYWORD_SET(POSTSCRIPT) then begin
41      closeps
42   endif
43
44
45   return
46end
Note: See TracBrowser for help on using the repository browser.