|
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
|
| Line | |
|---|
| 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.
