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.
std_plot_msf.pro in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/IDL_scripts – NEMO

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/IDL_scripts/std_plot_msf.pro @ 8516

Last change on this file since 8516 was 7732, checked in by flavoni, 7 years ago

commit idl scripts deleted

File size: 2.2 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  ;formulation for variable: TRANSPORT
36  ;mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma, /transport )
37  ;
38  ;formulation for variable: VELOCITY
39  mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma )
40  title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V
41  pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $
42        , small = [1, 2, 1], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so'
43
44  if std_file1_V NE std_file2_V then begin
45    ;formulation for variable: TRANSPORT
46    ;mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma, /transport )
47    ;
48    ;formulation for variable: VELOCITY
49    mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma)
50    title = title+' - '+std_file2_V
51    pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $
52          , small = [1, 2, 2], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so'
53  endif
54 
55  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
56  if KEYWORD_SET(postscript) then closeps
57
58  return
59end
60
Note: See TracBrowser for help on using the repository browser.