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_bsf.pro in NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_PISCES/IDL_scripts – NEMO

source: NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_PISCES/IDL_scripts/std_plot_bsf.pro @ 12808

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

commit idl scripts deleted

File size: 1.5 KB
Line 
1pro std_plot_bsf, U1, U2, 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  filename = cdti3 + '_BSF_'+std_file1_U
11  if std_file1_U NE std_file2_U then filename = filename + '_'+std_file2_U
12  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
13 
14  domdef, 0, 6000
15
16  ;formulation for variable: TRANSPORT
17  ;bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0], /transport )
18  ;
19  ;formulation for variable: VELOCITY
20  bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0])
21  title = 'Barotropic Stream Function!C'+std_file1_U
22  plt, bb1, min = -200., max = 200., int = 10., /portrait, FORMAT = '(I4)', STYLE = 'so0so' $
23       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title,  _extra = ex
24
25  if std_file1_U NE std_file2_U then BEGIN
26    ;formulation for variable: TRANSPORT
27    ;bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0], /transport )
28    ;formulation for variable: VELOCITY
29    bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0] )
30    title = 'Barotropic Stream Function!C'+std_file1_U+' - '+std_file2_U
31    plt, bb1.arr - bb2.arr, min = -20., max = 20., int = 2.,  FORMAT = '(I2)', STYLE = 'so0so' $
32         , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, /NOERASE, _extra = ex
33  endif
34 
35  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
36  if KEYWORD_SET(postscript) then closeps
37
38  return
39end
40
Note: See TracBrowser for help on using the repository browser.