[Thread Prev][Thread Next][Index]
[ferret_users] conversion of weekly data to monthly data.
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] conversion of weekly data to monthly data.
- From: STEEVEN PAUL YERRAGUNTLA <steevenpaul@xxxxxxxxx>
- Date: Thu, 28 Oct 2010 18:49:21 +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=762+a2s/09frBTgRNsmfOHQKb3bHlkjxxfPVr9wB5CI=; b=Wqe2FANXQlXdy1vZVGxrI1cY2+n65ALHvwnSiqMzZmKZ4Dy6622CZxPld0xZ7yh4Jm oR775egS4qV7Tnq3gnOQpNwpGRTFS9gVQtW3N5cBjN54E8E7i1gYYbUeVyT8wEb/qWHf Opp3T9uwW9yvwU4Xa7JIRe/G7p1mCHRxJQtiA=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:from:date:message-id:subject:to:content-type; b=Eb0bbrCiGU59sHLffr6lcXaQSPs7+mllQzXOcaQMQvY/JvniTjmrCHlXModVYvKkQt etLJS0o9VokxMYbnEizVxRaDP7iTwvhZ7SANEGILZtsa6I9DNfCnuRa31cfSnFs/HO8+ Gr0dEZZ1/6OVcD8Z069El1G6KC6g98OzwQaxo=
- Sender: owner-ferret_users@xxxxxxxx
Dear ferreters,
I wanted to convert weekly data to monthly data in a usual method in documentation and did not find any data in the resultant monthly data file and only monthly time axis can be seen. Please kindly let me know what is wrong with axis definition as i gave up after trying to sort out the problem. The jnl file lines and result as follows.
use salt_600_8697.cdf
set mem/size=400
let start_year=1986
let end_year =1997
let nyears = end_year - start_year + 1
LET indices = L[L=1:`nyears*12 + 1`]
!-> DEFINE VARIABLE indices = L[L=1:145]
!-> DEFINE VARIABLE indices = L[L=1:13]
LET month = MOD(indices-1,12)+1
LET year = start_year + INT((indices-1)/12)
DEFINE AXIS/t/UNITS=days/T0=1-jan-1900/EDGES truemonth = DAYS1900(year,month,1)
LET salt = S[gt=truemonth@ave]
save/file=salt_monthly_8697.cdf/clobber salt
LISTing to file salt_monthly_8697.cdf
yes? can var/all
yes? use salt_monthly_8697.cdf
yes? sh da
currently SET data sets:
1> ./sjava_600_8697.cdf
name title I J K L
S Salinity 1:161 1:91 1:60 1:720
2> ./salt_monthly_8697.cdf (default)
name title I J K L
SALT S[GT=TRUEMONTH@AVE] 1:161 1:91 1:60 1:144
yes? list/l=1:20/k=1 salt[x=@ave, y=@ave,d=2]
VARIABLE : S[GT=TRUEMONTH@AVE]
FILENAME : salt_monthly_8697.cdf
SUBSET : 20 points (TIME)
LONGITUDE: 39.8E to 120.2E (XY ave)
LATITUDE : 20.2S to 25.2N (XY ave)
Z (LAYERS): 1
16-JAN-1986 12 / 1:....
15-FEB-1986 00 / 2:....
16-MAR-1986 12 / 3:....
16-APR-1986 00 / 4:....
16-MAY-1986 12 / 5:....
16-JUN-1986 00 / 6:....
16-JUL-1986 12 / 7:....
16-AUG-1986 12 / 8:....
16-SEP-1986 00 / 9:....
16-OCT-1986 12 / 10:....
16-NOV-1986 00 / 11:....
16-DEC-1986 12 / 12:....
16-JAN-1987 12 / 13:....
15-FEB-1987 00 / 14:....
16-MAR-1987 12 / 15:....
16-APR-1987 00 / 16:....
16-MAY-1987 12 / 17:....
16-JUN-1987 00 / 18:....
16-JUL-1987 12 / 19:....
16-AUG-1987 12 / 20:....
-------------------------------------------------------------------------------------------------------
original weekly data as follows:
yes? list/l=1:10/k=1 s[x=@ave, y=@ave]
VARIABLE : Salinity (psu)
FILENAME : sjava_600_8697.cdf
SUBSET : 10 points (TIME)
LONGITUDE: 39.8E to 120.2E (XY ave)
LATITUDE : 20.2S to 25.2N (XY ave)
Z (LAYERS): 1
07-JUL-1957 00 / 1: 34.91
13-JUL-1957 00 / 2: 34.91
19-JUL-1957 00 / 3: 34.91
25-JUL-1957 00 / 4: 34.90
31-JUL-1957 00 / 5: 34.90
06-AUG-1957 00 / 6: 34.88
12-AUG-1957 00 / 7: 34.86
18-AUG-1957 00 / 8: 34.85
24-AUG-1957 00 / 9: 34.84
30-AUG-1957 00 / 10: 34.82
-----------------------------------------------------------------------------------
thanks in advance,
regards,
- Paul.
--
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