[Thread Prev][Thread Next][Index]
Re: [ferret_users] Extract gridded SST data on space-time scatteredbuoy SST data
Hi Theodor,
In your definition,
LET SST_AMSRE_T =
SAMPLET_DATE(SST_AMSRE_XY,2005,MONTH[D=3],DAY[D=3],12,0,0)
The arguments year, month, day, hour, minute, and second should all have
the same length. Set up variables year, hour, minute, day with the same
number of points as MONTH and DAY:
let year = 2005 + 0*month[d=3]
let hour = 12 + 0*month[d=3]
let minut = 0*month[d=3]
let secnd = 0*month[d=3]
and then
LET SST_AMSRE_T =
SAMPLET_DATE(SST_AMSRE_XY,year,MONTH[D=3],DAY[D=3],hour,minut,secnd)
I realize that the example in the Users Guide probably led you to call
the function with the single 0's in some of the arguments. I have
changed the Users Guide to reflect what the function does.
Ansley
Dr. Theodor Yemenis wrote:
Hi Ferreters,
I want to compare BUOY-SST data with satellite SST data.
The buoy data are in ascii format. I read them in FERRET with:
!
-------------------------------------------------------------------------
!
-------------------------------------------------------------------------
!
-------------------------------------------------------------------------
! READ BUOY VARIABLES ( SST )
set data/ez/skip=10/var="No, Buoy, Deziday, Year, Month, Day, Hour,
Minute, Se, P, TEMP, SST, Lat, Lon, UW, VW"
argos_lofzy05_000($BUOY_NR).dat
SET VARIABLE/BAD=8192.0/UNITS="grad C" SST[d=3]
!
-------------------------------------------------------------------------
sh da
3> ./argos_lofzy05_00009.dat
name title I J
K L
NO NO 1:4967 ...
... ...
BUOY BUOY 1:4967 ...
... ...
DEZIDAY DEZIDAY 1:4967 ...
... ...
YEAR YEAR 1:4967 ...
... ...
MONTH MONTH 1:4967 ...
... ...
DAY DAY 1:4967 ...
... ...
HOUR HOUR 1:4967 ...
... ...
MINUTE MINUTE 1:4967 ...
... ...
SE SE 1:4967 ...
... ...
P P 1:4967 ...
... ...
TEMP TEMP 1:4967 ...
... ...
SST SST 1:4967 ...
... ...
LAT LAT 1:4967 ...
... ...
LON LON 1:4967 ...
... ...
UW UW 1:4967 ...
... ...
VW VW 1:4967 ...
... ...
The satellite data are defined in a regular .25 degree grid.
My region is the norwegian sea and the time interval FEB-AUG 2005:
SET REGION/X=10W:26E/Y=60N:76N/T=1-FEB-2005:10-AUG-2005
!
-------------------------------------------------------------------------
! READ THE SST TMI/AMSR-E DATA
USE sst.tmi_amsre.LOFZY05.long.cdf
!
-------------------------------------------------------------------------
sh da
2> ./sst.tmi_amsre.LOFZY05.long.cdf
name title I J
K L
SST TMI/AMSR-E SST 1:144 1:64
... 1:1180
Now I want to extract the satellite data to the irregular space and
time positions of the buoy data.
I tried:
LET SST_AMSRE_XY = SAMPLEXY(SST[D=2],LON[D=3],LAT[D=3])
LET SST_AMSRE_T =
SAMPLET_DATE(SST_AMSRE_XY,2005,MONTH[D=3],DAY[D=3],12,0,0)
but it doesn't work.
The variable SST_AMSRE_XY is a two dimensional array (time,points),
the second variable SST_AMSRE_T is a one-dimensional array, but the
extracted values belong all to the first time-point
MONTH[D=3,i=1],DAY[D=3,i=1].
I would very much appreciate if someone can help me to extract the
satellite data correctly in space and time?
Theodor
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement