source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro @ 2751

Last change on this file since 2751 was 2751, checked in by smasson, 10 years ago

large update of IDL_scripts for Is-ENES

File size: 1.9 KB
Line 
1pro std_plot_msf, V1, V2, SUBBASIN = subbasin, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4
5@common 
6@std_common
7
8  cdti3 = string(cnt, format = '(i3.3)')
9  print, cdti3 + ') ' + blabla
10  if KEYWORD_SET(SUBBASIN) then subname = subbasin else subname = "Glo"
11  filename = cdti3 + '_MSF_'+subname+'_'+std_file1_V
12  if std_file1_V NE std_file2_V then filename = filename + '_'+std_file2_V
13  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
14;
15  if KEYWORD_SET(SUBBASIN) then begin
16    CASE subname of
17      'Atl' : var = 'atlmsk_nomed'
18      'Ind' : var = 'indmsk'
19      'IndoPac' : var = 'indpacmsk'
20      'GloNoMed' : var = 'glomsk_nomed'
21    ENDCASE
22    msk = read_ncdf( var, filename = std_file_msksub, /nostruct, _extra = ex )
23  endif else msk = tmask[*, *, 0]
24;
25  CASE subname OF
26    'GloNoMed':lat_ext = [-80, 90]
27    'Glo':lat_ext = [-80, 90]
28    'Atl':lat_ext = [-30, 90]
29    'Ind':lat_ext = [-30, 30]
30    'IndoPac':lat_ext = [-30, 70]
31  ENDCASE
32 
33  domdef, 0, 6000
34
35  mm1 = msf(V1.arr, indexboxzoom = ind, maskout = ma)
36  title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V
37  pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $
38        , small = [1, 2, 1], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so'
39
40  if std_file1_V NE std_file2_V then begin
41    mm2 = msf(V2.arr, indexboxzoom = ind, maskout = ma)
42    title = title+' - '+std_file2_V
43    pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $
44          , small = [1, 2, 2], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so'
45  endif
46 
47  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
48  if KEYWORD_SET(postscript) then closeps
49
50  return
51end
52
Note: See TracBrowser for help on using the repository browser.