Hi,
I have a
multiyear daily time series of precipitation data and
I am strugling to create composite of active (and
break) monsoon period.
I think
what I need is a mask
along the time axis that is based on
a criteria (say area averaged precipitation >1 for
3 or more successive days). I can then apply the mask
to create active phase composites for precipitation
and other variables. Let me illustrate with an
example: say daily area average precip are:
01-JUL-2007:
0.2.
02-JUL-2007: 0.8
03-JUL-2007: 1.2
04-JUL-2007: 2.1
05-JUL-2007: 2.0
06-JUL-2007: 1.5
07-JUL-2007: 1.1
08-JUL-2007: 0.4
09-JUL-2007: 0.1
10-JUL-2007: 0.6
From
above definition definition: the active phase is from 3
- 7 Jul and what I want is a mask with 1 on those dates
and undefined for other dates.
Other
wish list are: to plot active days; get a list of just
the valid days and finally save the mask for later use
on other variables.
This
is how far I have gone: For above example, I am able to
creat a mask with 1's from 3-5 Jul but not able to
extend it to 7 Jul.
let
atv=if precip GT 1 AND precip[l=@SHF:1] GT 1 AND
precip[l=@SHF:2] GT 1 then 1
list
atv
01-JUL-2007:
.....
02-JUL-2007: .....
03-JUL-2007:
1.0
04-JUL-2007: 1.0
05-JUL-2007: 1.0
06-JUL-2007: ....
07-JUL-2007: ....
08-JUL-2007: ....
09-JUL-2007: ....
10-JUL-2007: .....
What
do I modify if active period is to be valid only between
1 jul to 31 August
Many
thanks in advance.
Jagadish