[Thread Prev][Thread Next][Index]

[ferret_users] Extract gridded SST data on space-time scattered buoy SST data



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


--
--------------------------------------
Dr. Theodoros Yemenis
Center for Marine and Atmospheric Research
University of Hamburg
Institute of Oceanography
Bundesstrasse 53
D-20146 Hamburg
Germany
---------------------------------------

phone: +49 40 42838 6104
fax: +49 40 42838 7063
+49 40 36 801 801
email: yemenis@ifm.zmaw.de ---------------------------------------




[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement