[RegCNET] ICBC temperature Interpolation.
bixq
bixq at ictp.it
Mon Nov 23 09:03:00 CET 2009
Hi, Owanda:
You use 'OI_WK' and your version is not the latest, right?
Please do the following thing:
1. download the latest sst.wkmean.1990-present.nc file;
2. ncdump -h sst.wkmean.1990-present.nc, and get the header information
like
time = UNLIMITED ; // (1029 currently)
3. then, edit SST_1DEG.f and ICBC.f, find the COMMON /DATEWK/
In the latest version, it's like
integer WKDAY
COMMON /DATEWK/ WKDAY(427+1045)
(1045 means upto 2010010300), I guess that your's is smaller,
then increase the number of your's to a bigger number (1045, or
the number you got from step 2;
Note: you should do this kind of replace several places.
4. It should be fine now.
Of course, you can just redownload the latest code.
Regards,
Xunqiang Bi
On Sat, 21 Nov 2009, Owanda Otieno wrote:
> Dear all
> I'm trying to run ./icbc.x to interpolte sst for 2008, but I get the
> error massage shown below. When I try to interpolate sst for
> for other years prier to 2007 it works, can someone help in
> identifying the proble. Part of the interpolation is as shown;
>
> 6 ../DATA/NNRP2/2008/air.2008.nc 9.999999776482582E-003
> 477.660003662109
> 7 ../DATA/NNRP2/2008/hgt.2008.nc 1.00000000000000
> 32066.0000000000
> 8 ../DATA/NNRP2/2008/rhum.2008.nc 9.999999776482582E-003
> 302.660003662109
> 9 ../DATA/NNRP2/2008/uwnd.2008.nc 9.999999776482582E-003
> 202.660003662109
> 10 ../DATA/NNRP2/2008/vwnd.2008.nc 9.999999776482582E-003
> 202.660003662109
> READ IN fields at DATE: 2008020100
>
> SST file is not the right one
>
> Thank you
>
> --
> Owanda Otieno
> Department of Natural Resources and the Environment
> University of Connecticut
> 1376 Storrs Road, Unit 4087
> Storrs, CT 06269-4087
> Phone: (860) 486-1876
>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Xunqiang Bi email:bixq at ictp.it
Earth System Physics Group
The Abdus Salam ICTP
Strada Costiera, 11
P.O. BOX 586, 34100 Trieste, ITALY
Tel: +39-040-2240302 Fax: +39-040-2240449
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
More information about the RegCNET
mailing list