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 WRFConventions?
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.