New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
xxx_Med_S_tongue.pro in trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts – NEMO

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

Last change on this file since 2709 was 2709, checked in by flavoni, 13 years ago

update idl scripts, 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 
16   title = sEXP1+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
17   plt, S1.arr, MININ=35., MAXIN=37., INTER=0.1,  UNIT= varunit $
18         , small=[1,2,1], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
19;
20  if sEXP1 ne sEXP2 then begin
21     
22     title = sEXP1+' - '+sEXP2+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
23     plt, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05,  UNIT = varunit $
24          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
25 
26  endif else begin
27     
28     title = sEXP1+' - Levitus Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
29     plt, S1.arr - SLev.arr, MININ=-1., MAXIN=1., INTER=0.1,  UNIT = varunit $
30          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
31     
32  endelse
33
34  title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
35  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  />  ' ]
36  putfile, '/tmp/prova.txt', text
37   
38  if KEYWORD_SET(POSTSCRIPT) then begin
39      closeps
40   endif
41
42
43   return
44end
Note: See TracBrowser for help on using the repository browser.