shr_map_getWts error for a single point simulation

[koichi_sakaguchi@pnnl_gov]

I am trying to run CLM4 for a single point with my own surface and atmospheric forcing data.
I have made the required surface forcing data sets with "getregional_datasets.pl"
For atmospheric forcing, I downloaded an example monthly forcing data for Vancouver, Canada and changed its lat, lon, and grid edges to be the same as the surface forcing data for testing. The coordinates look like:

forcing files: atmo, surface
EDGEN: -1.894737, -1.89473684210535
EDGES: -3.789474, -3.78947368421061
ESGEE: 306.25, 306.25
EDGEW: 303.75, 303.75
LATIXY: -2.842105, -2.84210526315798
LONGXY: 305, 305


But I get this error message at runtime.

('shr_map_getWts') ERROR: yd outside bounds -2.84210526315798
(shr_sys_abort) ERROR: ('shr_map_getWts') ERROR yd outside 90 degree bounds

Somehow the subroutine shr_map_getWts complains that my latitude coordinate is beyond + or -90 degree bounds, but they are obviously not.

If anyone knows how to solve this problem, please let me know.

Thanks,

Koichi

 

[slevis]

Dear Koichi,

I haven't seen this before, and I don't have experience with clm4's methods of running single point BUT...

Is it possible that the error comes from your surface's north edge being further north than your atmosphere's north edge?

Sam

 

[erik]

Koichi

This sounds like a known problem with ccsm4 and cesm1.0, that should be fixed in the latest cesm1.0.1 release. I would update to the cesm1.0.1 release and first try the case outlined for Vancouver from the User's Guide.

http://www.cesm.ucar.edu/models/cesm1.0/clm/models/lnd/clm/doc/UsersGuide/x3279.html#sp_atm_forcing

make sure that works and then move onto your case.


Erik

 

[koichi_sakaguchi@pnnl_gov]

Hi Sam and Erik, thank you both for your replies.

I updated the model to CESM1.0.1 and ran the example for Vancouver given in the CLM4 user guide, and it worked! Before updating the model, I got the same error message I posted for Vancouver example. So hopefully my own simulation for the Amazon single point also work. Probably I need to modify the lat and lon as Sam pointed out.

I will post the result when I'm done with the Amazon.

Again, thank you so much for your help,

Koichi

 

[jlchen@lzb_ac_cn]

Hi, Koichi.Sakaguchi, how did you solved this problem?