[Thread Prev][Thread Next][Index]

[ferret_users] conversion of weekly data to monthly data.



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