source: trunk/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_S.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.5 KB
Line 
1pro xxx_zonal_S, S1, S2, SLev, 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 else  subname = "Glo"
8   
9
10   if KEYWORD_SET(POSTSCRIPT) then begin
11      openps, '015_Zonal_S_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
12   endif
13
14                                ; looking for longitudinal index corresponding to the highest latitude (closest to the North Pole)
15   index=where(gphit eq max(gphit))
16   nx=index[0] mod jpi
17   gphi_save=gphit
18   gphit[0,*]=gphit[nx,*]
19
20   if KEYWORD_SET(SUBBASIN) then begin
21                                ; save tmask
22      tmasksv = tmask
23                                ;
24      CASE 1 of
25         subbasin eq 'Atl' : var = 'atlmsk_nomed'
26         subbasin eq 'Ind' : var = 'indmsk_nored'
27         subbasin eq 'Pac' : var = 'pacmsk'
28      ENDCASE
29                                ;
30      msk = read_ncdf( var, filename = mask_filename, _extra =ex)
31
32; from 2D array to 3D array
33      msk = msk.arr[*]#replicate(1., nzt)
34      tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] = tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] * msk
35   endif
36
37   varunit = S1.unit
38;
39   title=sEXP1+' '+subname+' / S year '+start_end
40   pltz, S1.arr, MININ = 32., MAXIN = 37., INTER = .2, typein = 'yz' $
41         , small = [1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT,  _extra=ex
42;
43   if sEXP1 ne sEXP2 then begin
44     
45      title = sEXP1+' - '+sEXP2+' '+subname+' / S year '+start_end
46      pltz, S1.arr-S2.arr, MIN = -1., MAX = 1., INTER = .1, typein = 'yz' $
47           , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE,  _extra=ex
48 
49   endif else begin
50     
51      title=sEXP1+' - Levitus '+subname+'/ S year'+start_end
52      pltz, S1.arr-SLev.arr, MININ = -1, MAXIN = 1., INTER = .1, typein = 'yz'  $
53            , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex
54   endelse
55;
56   title = 'Salinity differences in zonal average ' + subname +' with Levitus and '+sEXP1+'-'+sEXP2
57   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/015_Zonal_S_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
58   putfile, '/tmp/prova.txt', text
59
60   if KEYWORD_SET(POSTSCRIPT) then begin
61      closeps
62   endif
63
64
65   if KEYWORD_SET(SUBBASIN) then tmask = tmasksv
66
67   gphit=gphi_save
68
69   return
70end
Note: See TracBrowser for help on using the repository browser.