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_read.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_ts_read.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_ts_read, var_name, dt1, dt2, prefix, suffix, ts, ts_z, masknp $
2                 , WITHSSH = withssh, SSHPREFIX = sshprefix, SSHSUFFIX = sshsuffix, LEVZ = levz
3
4  compile_opt idl2, strictarrsubs
5 
6@common
7@std_common
8
9
10  list = rseries_ncdf(var_name, dt1, dt2, prefix, suffix, /fileslist)
11  nfiles = n_elements(list)
12  IF keyword_set(withssh) THEN BEGIN
13    sshlist = rseries_ncdf(withssh, dt1, dt2, sshprefix, sshsuffix, /fileslist)
14    IF nfiles NE n_elements(sshlist) THEN stop
15  ENDIF
16
17  ts = 0.
18  ts_Time = 0.
19  ts_z = fltarr(jpk)
20
21  FOR i = 0, nfiles-1 DO BEGIN
22    IF keyword_set(withssh) THEN ssh = read_ncdf(withssh, allrecords = allrec, filename = sshlist[i], /nostruct)
23    var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct)
24    ts_Time = [ ts_Time, Time]
25;
26    IF jpt EQ 1 THEN txyz = moyenne(var, 'xyz', mask2d = masknp, ssh = ssh) $
27    ELSE txyz = grossemoyenne(var, 'xyz', mask2d = masknp, ssh = ssh)
28    ts = [ ts, txyz ]
29;
30    IF jpt EQ 1 THEN tz = moyenne(temporary(var), 'xy', mask2d = masknp, ssh = ssh, /KEEPBOTTOM) $
31    ELSE tz = grossemoyenne(temporary(var), 'xy', mask2d = masknp, ssh = ssh, /KEEPBOTTOM)
32    IF keyword_set(levz) THEN ts_z = [ [ts_z], [tz - (levz[*] # replicate(1., jpt))] ] $
33    ELSE ts_z = [ [ts_z], [tz] ]
34  ENDFOR
35
36  time = ts_Time[1:*]           ; remove first record of 0
37  jpt = n_elements(time)
38 
39  ts = ts[1:*]        ; remove first record of 0
40  ts_z = ts_z[*, 1:*]           ; remove first record of 0
41
42
43return
44end
Note: See TracBrowser for help on using the repository browser.