[Thread Prev][Thread Next][Index]
[ferret_users] How to extract DJF average time series from a monthly time series?
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] How to extract DJF average time series from a monthly time series?
- From: steeven paul yerraguntla <steevenpaul@xxxxxxxxx>
- Date: Mon, 11 Oct 2010 18:42:10 +0530
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:mime-version:received:from:date :message-id:subject:to:content-type; bh=4ASTAG/jLKPEMqyxocd1cNjZNGbP18XC2Md8TFKDASg=; b=soNjwZnWVw8qOQkIWDPHw0eWICEg4pFSL7mFCGh5dR6mIG/MVm+PpHXjN6NahJ9/98 LWAy1NGDTu2slHyhvjTkPkI+QqGsHoBqKvX7PtUd6+67iqzbhkBRdS0Mw+uL+gkjd0L3 K8fOzs7pm29x9kxGAGOAiC2oatHSQyzEyPrz4=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:from:date:message-id:subject:to:content-type; b=B3NTTwL+4heJIAh5cr1NPs1eQAJOh31xrcelHskpzWfoOrE8joE0dXbKZ9umpXpOxQ 8w768fdzV7RvL9GWVIbnBmWjqPNzStAGYkwy1YfG8naOgMIyL2HMnwwMO69ZmtbY6LoG 7XymG1FLcHV+gNfq4sEEx/IHQ3bAkP9arxDEw=
- Sender: owner-ferret_users@xxxxxxxx
Dear ferreters,
������ I want to extract the DJF months time series in the following format.
������ 15-dec-1980
������ 15-jan-1981
������ 15-feb-1981
������ 15-dec-1981
������ 15-jan-1982
������ 15-feb-1982
������������ ----
�������� and so on
I tried to follow the examples in archives, provided by Jaison Kurian and Ansley which are for extracting JJA months and i could not get time series for DJF months though making some changes to the scripts. Hope someone will help me in this connection.��� Ansley's� Scripts as follows
--------------������ -----------------����������� ---------------------------------------���������� ------------------------------------------
use monthly_navy_winds
let nmon = `uwnd,return=Lend`
let tt = t[t=1:`nmon`]
let june = if mod(tt,12) ge 6 then tt
let june_aug = if mod(june,12) le 8 then june
list june_aug
! The number of summer months.
let nja = `june_aug[t=@ngd]`
let just_summer = compressl(june_aug)
list just_summer[L=1:`nja`]
! Now sample the Lth point of uwnd for the summer months
let uwnd_jja = samplel(uwnd,just_summer[l=1:`nja`])
! Define a calendar axis containing these months.
def axis/t0="`uwnd,return=t0`"/units="months" tsummer =
just_summer[l=1:`nja`]
! Put the sampled summer data on the calendar axis
let uwnd_summer = uwnd_jja[gt=tsummer@asn]
list/x=180/y=0 uwnd_summer
--------------------------------------------------------------- end of script --------------------------------------------------------------------------------------------------
Thanks in advance,
�
--
Steeven Paul Y
Research Fellow
NIO RC Vizag
176 Lawsons Bay Colony
Visakhapatnam - 530017
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement