Help Required in preparing SSS forcing file....

Message

jvratnam · #1 Unread post by **jvratnam** » Thu Apr 12, 2007 10:05 am

Hi

I prepared an ascii file containing the monthly Sea Surface Salinity from levitus data set using the procedure specified in the forcing.in file. I am specifying the name of the ascii data file and trying to run forcing package. I am getting the following error

GET_GRDFLD: latitude point = 39.54 is outside of input gridded data at -48.96 grid units.

I tried printing the values before the abort and I am getting

IM 89 JM 72 PLON 39.5402259826660156 PLAT 39.5402259826660156 OLON 0.500000000000000000 OLAT 89.5000000000000000

The levitus data is from 89.5N-89.5S , 0.5E-359.5E and my domain is from 20S-40N and 40E-120E.

Am I doing something wrong in generating the ascii data file ?

Thanks
Ratnam

jvratnam · #2 Unread post by **jvratnam** » Thu Apr 12, 2007 10:41 am

Hi

My problem is solved. There was a bug in the get_grid.F file. After fixing it I am able to properly generate the forcing file.

The change is to be made in line 200 of get_grid.F

get_grid.F

from:

status=nf_get_vara_FTYPE(ncgrdid,vid(idrlon),start,count,
& rlat)

to:
status=nf_get_vara_FTYPE(ncgrdid,vid(idrlat),start,count,
& rlat)

Ratnam