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_ts_SSH.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_ts_SSH.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.8 KB
Line 
1pro std_ts_SSH, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4 
5@common
6@std_common
7
8; get exp1 info
9  vssh = getenv('VAR1_SSH')    &   prefix = getenv('V1SSH_PREF')    &   suffix = getenv('V1SSH_SUFF')
10; get exp2 info
11  vssh2 = getenv('VAR2_SSH')   &   prefix2 = getenv('V2SSH_PREF')   &   suffix2 = getenv('V2SSH_SUFF')
12;
13  cdti3 = string(cnt, format = '(i3.3)')
14  print, cdti3 + ') ' + blabla
15  filename = cdti3 + '_ts_SSH_'+prefix
16  if prefix NE prefix2 then filename = filename + '_'+prefix2
17  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
18
19  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')'
20;
21  iodir = std_iodir_data
22;
23  ts_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp)
24  ts_SSH.arr = ts_SSH.arr   &   ts_SSH.unit = 'm'
25  title = prefix+' '+d1_d2+'!C'+blabla
26  pltt, ts_SSH, 't', MIN = -12., MAX = 12., date1, date2, /REMPLI , XGRIDSTYLE = 1 $
27        , small = [1, 2, 1], YTITLE = 'm', TITLE = title, /PORTRAIT, _extra = ex
28 
29  IF prefix NE prefix2 THEN BEGIN
30
31    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')'
32    tsave = time
33    ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp)
34    ts_SSH2.arr = ts_SSH2.arr   &   ts_SSH2.unit = 'm'
35    time = tsave   &   IF n_elements(time) NE jpt THEN stop
36
37    title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla
38    pltt, ts_SSH.arr - ts_SSH2.arr, 't', MIN = -10., MAX = 10., date1, date2, /REMPLI, XGRIDSTYLE = 1 $
39          , color = 250, small = [1, 2, 2], YTITLE = 'm', TITLE = title, /noerase, _extra = ex
40
41  ENDIF
42
43  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
44  if KEYWORD_SET(postscript) then closeps
45 
46  return
47end
48
Note: See TracBrowser for help on using the repository browser.