; @initers
; a = a=read_ncdf('sozotaux',19920420,20000214,file='tx_ERS_19920420_20000214.nc',/nostruct)
; ; pb a pas de temps 222 c''est le 19960722
; b=a
; b[*,*,222]=0.
; b=total(b,3)/408.
; mask=finite(b)
;
a = read_ncdf('sozotaux',19920420,20000214,file='tx_ERS_19920420_20000214.nc',/nostruct)
a[where(a EQ 0.)] = 1e-5
b = a[*,*,222]
b[where(b GE 10)] = !values.f_nan
a[*,*,222] = (remplit(b,nite=2,/nan)) < 300.8
deb = where(time GE julday(1, 1, 1993)) & deb = deb[0]
fin = where(time lE julday(1, 1, 1997)) & fin = (reverse(fin))[0]
print, vairdate(time[deb]), vairdate(time[fin])
b = a[*, *, deb:fin]

b = total(b, 3)/(fin-deb+1)
b = b[*]#replicate(1, jpt)
a = a-temporary(b)

restore, 'maskind_ers.dat'

projecttx, a, mask, 2.0, FILECOEF='coeffprojtx_ers_ind.dat'