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 branches/2014/dev_r4650_UKMO14.11_SETTE_OBSASM/NEMOGCM/CONFIG/ORCA2_LIM_OBS/IDL_scripts – NEMO

source: branches/2014/dev_r4650_UKMO14.11_SETTE_OBSASM/NEMOGCM/CONFIG/ORCA2_LIM_OBS/IDL_scripts/std_plot_bsf.pro @ 4751

Last change on this file since 4751 was 4751, checked in by djlea, 10 years ago

Changes to include an OBS test in SETTE. At the moment this uses an example profile observation.

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.