[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: WRFConvention.java problem.

Hi Todd:

I am on vacation and cant look at this till I get back around July 8. Generally you should send these problems to address@hidden or address@hidden In this case it will probably wait until i get back anyway.


Todd Plessel wrote:
Hi John,

We are displaying WRF polar stereographic gridded files
in VERDI - a Java application that uses
netcdf-4.1 Java which contains WRFConvention.java.


Apparently the map is not being projected correctly.
The attached zip file contains a sample WRF NetCDF
data file wrfinput_d01 and screen.gif which shows
correct CMAQ (M3IO NetCDF) data and map on the left and
WRF land-use data and mis-matched map on the right.
The grids are almost identical and so should display
almost the same. Note that South America and
Africa are clearly mis-registered with respect to the
land-use data where (most/all) red is ocean.

Some relevant attributes are:

west_east = 199 ;
south_north = 199 ;
float LU_INDEX(Time, south_north, west_east) ;
:DX = 108000.f ;
:DY = 108000.f ; :GRIDTYPE = "C" ;
:MAP_PROJ = 2 ;
:CEN_LAT = 90.f ;
:CEN_LON = -8.f ;
:TRUELAT1 = 45.f ;
:TRUELAT2 = 45.f ;
:MOAD_CEN_LAT = 90.f ;
:STAND_LON = -98.f ;

CEN_LON looks suspect. Perhaps it should be -98?
In any case, which attributes are actually used by

I don't see any information about the ellipsoid in the file.
CMAQ usually uses a sphere of radius 6,370,000 meters.
Perhaps WRF does too.
It would seem that the constructors for all projectors should
have ellipsoid input parameters e.g.,
double major_semiaxis,  double minor_semiaxis.

Is there a fix available?



NOTE: All email exchanges with Unidata User Support are recorded in the Unidata inquiry tracking system and then made publicly available through the web. If you do not want to have your interactions made available in this way, you must let us know in each email you send to us.