My coordinate system was viewed as Polar coordinates, but it is now

[qhgogogo]

Dear NCAR scientist,

I recently created my own forcing and domain data for running CLM45 which I got the following error:

'ERROR: ('shr_map_getWts') ERROR yd outside 90 degree bounds'.

in shr_abour_mod.f90.

I found my coordinate units were treated as radiance but they are not, they are lat/lon degrees.

I am not sure what caused this error in my domain/forcing files. Could you help me resolve this?

Thanks in advance!

Han

 

[erik]

I assume you mean "radians" rather than "radiance". But, if your units are degrees make sure the "unit" attribute labels it as such on your lat/lon coordinates. There's also some information output besides what you show above that might help you to figure it out. Depending on what version of CLM you are using you should be able to look at the file shr_map_mod.F90 to see more about the error message and what's going wrong.

Is your forcing data global or just over a given region? If the units aren't provided running a regional grid can cause it to make an error in detecting what the units are.

 

[qhgogogo]

I assume you mean "radians" rather than "radiance". But, if your units are degrees make sure the "unit" attribute labels it as such on your lat/lon coordinates. There's also some information output besides what you show above that might help you to figure it out. Depending on what version of CLM you are using you should be able to look at the file shr_map_mod.F90 to see more about the error message and what's going wrong.

Is your forcing data global or just over a given region? If the units aren't provided running a regional grid can cause it to make an error in detecting what the units are.

Thanks a lot, Erik! The issue is resolved after debugging into shr_map_modF90. My forcing is for a very small region such that the function 'shr_map_checkRad' in shr_map_mod.F90 detected it as radians. I modified this and it is now correct. Thank you!

Best,

Han