source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_depth.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_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, '029_'+sEXP1+'_Med_Levitus_S_depth'+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   varunit = S1.unit
17
18   title = sEXP1+'  Mediterranean water at ' +strtrim(lat,1)+'°N / year '+start_end
19   pltz, S1.arr, MININ = 35., MAXIN = 37., INTER = 0.1,  UNIT = varunit $
20         , small = [1,2,1], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., /PORTRAIT, _extra = ex
21;
22
23   if sEXP1 ne sEXP2 then begin
24     
25      title = sEXP1+' - '+sEXP2+' Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
26      pltz, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05, UNIT = varunit $
27            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000.,  _extra = ex
28     
29   endif else begin
30   
31      title = sEXP1+' - Levitus Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
32      pltz, S1.arr-SLev.arr, MININ = -1., MAXIN = 1., INTER=0.1, UNIT = varunit $
33            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., _extra = ex
34   endelse
35
36
37   title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
38   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/029_'+sEXP1+'_Med_Levitus_S_depth'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.png  />  ' ]
39   putfile, '/tmp/prova.txt', text
40
41   if KEYWORD_SET(POSTSCRIPT) then begin
42      closeps
43   endif
44
45
46   return
47end
Note: See TracBrowser for help on using the repository browser.