Last change
on this file since 2 was
2,
checked in by opalod, 22 years ago
|
Initial revision
|
-
Property svn:eol-style set to
native
-
Property svn:executable set to
*
-
Property svn:keywords set to
Author Date Id Revision
|
File size:
773 bytes
|
Rev | Line | |
---|
[2] | 1 | ; @initers |
---|
| 2 | ; a = a=read_ncdf('sozotaux',19920420,20000214,file='tx_ERS_19920420_20000214.nc',/nostruct) |
---|
| 3 | ; ; pb a pas de temps 222 c''est le 19960722 |
---|
| 4 | ; b=a |
---|
| 5 | ; b[*,*,222]=0. |
---|
| 6 | ; b=total(b,3)/408. |
---|
| 7 | ; mask=finite(b) |
---|
| 8 | ; |
---|
| 9 | a = read_ncdf('sozotaux',19920420,20000214,file='tx_ERS_19920420_20000214.nc',/nostruct) |
---|
| 10 | a[where(a EQ 0.)] = 1e-5 |
---|
| 11 | b = a[*,*,222] |
---|
| 12 | b[where(b GE 10)] = !values.f_nan |
---|
| 13 | a[*,*,222] = (remplit(b,nite=2,/nan)) < 300.8 |
---|
| 14 | deb = where(time GE julday(1, 1, 1993)) & deb = deb[0] |
---|
| 15 | fin = where(time lE julday(1, 1, 1997)) & fin = (reverse(fin))[0] |
---|
| 16 | print, vairdate(time[deb]), vairdate(time[fin]) |
---|
| 17 | b = a[*, *, deb:fin] |
---|
| 18 | |
---|
| 19 | b = total(b, 3)/(fin-deb+1) |
---|
| 20 | b = b[*]#replicate(1, jpt) |
---|
| 21 | a = a-temporary(b) |
---|
| 22 | |
---|
| 23 | restore, 'maskind_ers.dat' |
---|
| 24 | |
---|
| 25 | projecttx, a, mask, 2.0, FILECOEF='coeffprojtx_ers_ind.dat' |
---|
Note: See
TracBrowser
for help on using the repository browser.