[3331] | 1 | PRO std_ts_read, var_name, dt1, dt2, prefix, suffix, ts, ts_z $ |
---|
| 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 | |
---|
| 43 | return |
---|
| 44 | end |
---|