You are here

Home » Forums » Usage » Usage

wrf2geotiff.ncl is generating offsets that are way too large

18 posts / 0 new
Log in or register to post comments
Last post
October 21, 2014 - 5:53pm
#1
steve.kirby
wrf2geotiff.ncl is generating offsets that are way too large

Hi,

I found a block of NCL code online for displaying WRF data and it calls VAPOR's wrf2geotiff.ncl, but there is a problem with the offsets in the "crop" portion.

It is generating x and y offsets that seem somewhat strange. Maybe someone can identify what could be going wrong. (Note: I am using the binary version of VAPOR for Linux x86_64. Thought about building from source, but it errored out very early in the compilation so I went with binaries.)

First I create my PS workstation as required, from my WRF output, consisting of T2, U10, and V10.

a=addfile("./wrfout_d02_2010-08-20_20_00_00.nc","r");
wks = gsn_open_wks("ps","steve_plot")
plot = wrf_map_overlays(a,wks,(/contour_t,contour_p,vector/),pltres,mpres)

Looked at the resultant PS using ghostview and it looks fine.

then:
wrf2gtiff = wrf2geotiff_open(wks)
wrf2geotiff_write(wrf2gtiff,a,times(it),wks,plot,True) (admittedly, I don't know exactly what this call does yet)

then:
wrf2geotiff_close(wrf2gtiff,wks) --
This call generates this call to ImageMagick -
convert -depth 8 -rotate 0 -density 144 -crop 19171.3x19171+9611.25+-29424.9 /tmp/tiffcap3_0001.ps tiff:/tmp/tiffcap3_0001.tiff

and I get this error: "convert: geometry does not contain image `/tmp/tiffcap3_0001.ps' @ warning/transform.c/CropImage/666.
convert: TIFF: negative image positions unsupported /tmp/tiffcap3_0001.tiff @ error/tiff.c/WriteTIFFImage/3317"

Clearly, the problem is that the wrf2geotiff_close() call is generating x and y offsets of 9611.25 and -29424.9, which are both way off image, thus the error. In fact, if I use ImageMagick's "identify" on this *.ps file it tells me:

$ identify /tmp/tiffcap3_0001.ps
/tmp/tiffcap3_0001.ps PNG 612x764 612x764+0+0 16-bit sRGB 197KB 0.000u 0:00.000

Don't know why it says "PNG", but it is indicating the image is only 612x764.

If anyone can tell me why such gigantic offsets (9611.25, -29424.9) might be determined in the VAPOR NCL code, wrf2geotiff.ncl that would be excellent. Also, the height and width generated by the VAPOR NCL, (19171x19171 differ by a lot from what ImageMagick is telling me about the PS file). Maybe wrf2geotiff_write() is doing something I don't understand?

TIA for any tips.

Top
October 22, 2014 - 8:48am
steve.kirby
(Note: Forgot to format on my first submit - this is now readable :>) Hi, I found a block of NCL code online for displaying WRF data and it calls VAPOR's wrf2geotiff.ncl, but there is a problem with the offsets in the "crop" portion of the VAPOR code. It is generating x and y offsets that seem somewhat strange. Maybe someone can identify what could be going wrong. (Note: I am using the binary version of VAPOR for Linux x86_64. Thought about building from source, but it errored out very early in the compilation so I went with binaries.) First I create my PS workstation as required, from my WRF output, consisting of T2, U10, and V10. a=addfile("./wrfout_d02_2010-08-20_20_00_00.nc","r"); wks = gsn_open_wks("ps","steve_plot") plot = wrf_map_overlays(a,wks,(/contour_t,contour_p,vector/),pltres,mpres) Looked at the resultant PS using ghostview and it looks fine. then: wrf2gtiff = wrf2geotiff_open(wks) wrf2geotiff_write(wrf2gtiff,a,times(it),wks,plot,True) (admittedly, I don't know exactly what this call does yet) then: wrf2geotiff_close(wrf2gtiff,wks) -- This call generates this call to ImageMagick - convert -depth 8 -rotate 0 -density 144 -crop 19171.3x19171+9611.25+-29424.9 /tmp/tiffcap3_0001.ps tiff:/tmp/tiffcap3_0001.tiff and I get this error: "convert: geometry does not contain image `/tmp/tiffcap3_0001.ps' @ warning/transform.c/CropImage/666. convert: TIFF: negative image positions unsupported /tmp/tiffcap3_0001.tiff @ error/tiff.c/WriteTIFFImage/3317" Clearly, the problem is that the wrf2geotiff_close() call is generating x and y offsets of 9611.25 and -29424.9, which are both way off image, thus the error. In fact, if I use ImageMagick's "identify" on this *.ps file it tells me: $ identify /tmp/tiffcap3_0001.ps /tmp/tiffcap3_0001.ps PNG 612x764 612x764+0+0 16-bit sRGB 197KB 0.000u 0:00.000 Don't know why it says "PNG", but it is indicating the image is only 612x764. If anyone can tell me why such gigantic offsets (9611.25, -29424.9) might be determined in the VAPOR NCL code, wrf2geotiff.ncl that would be excellent. Also, the height and width generated by the VAPOR NCL, (19171x19171 differ by a lot from what ImageMagick is telling me about the PS file). Maybe wrf2geotiff_write() is doing something I don't understand? TIA for any tips.
Top
October 22, 2014 - 8:32pm
clyne

Hi Steve,

 

Can you post the output of the command:

 

ncdump -h wrfout_d02_2010-08-20_20_00_00.nc

 

thanks - jc

Top
October 23, 2014 - 3:03pm
(Reply to #3)
steve.kirby

Hi,

Here it is --

netcdf wrfout_d02_2010-08-20_20_00_00 {
dimensions:
Time = UNLIMITED ; // (1 currently)
DateStrLen = 19 ;
west_east = 120 ;
south_north = 120 ;
bottom_top = 89 ;
bottom_top_stag = 90 ;
soil_layers_stag = 4 ;
west_east_stag = 121 ;
south_north_stag = 121 ;
variables:
char Times(Time, DateStrLen) ;
float LU_INDEX(Time, south_north, west_east) ;
LU_INDEX:FieldType = 104 ;
LU_INDEX:MemoryOrder = "XY " ;
LU_INDEX:description = "LAND USE CATEGORY" ;
LU_INDEX:units = "" ;
LU_INDEX:stagger = "" ;
LU_INDEX:coordinates = "XLONG XLAT" ;
float ZNU(Time, bottom_top) ;
ZNU:FieldType = 104 ;
ZNU:MemoryOrder = "Z " ;
ZNU:description = "eta values on half (mass) levels" ;
ZNU:units = "" ;
ZNU:stagger = "" ;
float ZNW(Time, bottom_top_stag) ;
ZNW:FieldType = 104 ;
ZNW:MemoryOrder = "Z " ;
ZNW:description = "eta values on full (w) levels" ;
ZNW:units = "" ;
ZNW:stagger = "Z" ;
float ZS(Time, soil_layers_stag) ;
ZS:FieldType = 104 ;
ZS:MemoryOrder = "Z " ;
ZS:description = "DEPTHS OF CENTERS OF SOIL LAYERS" ;
ZS:units = "m" ;
ZS:stagger = "Z" ;
float DZS(Time, soil_layers_stag) ;
DZS:FieldType = 104 ;
DZS:MemoryOrder = "Z " ;
DZS:description = "THICKNESSES OF SOIL LAYERS" ;
DZS:units = "m" ;
DZS:stagger = "Z" ;
float U(Time, bottom_top, south_north, west_east_stag) ;
U:FieldType = 104 ;
U:MemoryOrder = "XYZ" ;
U:description = "x-wind component" ;
U:units = "m s-1" ;
U:stagger = "X" ;
U:coordinates = "XLONG_U XLAT_U" ;
float V(Time, bottom_top, south_north_stag, west_east) ;
V:FieldType = 104 ;
V:MemoryOrder = "XYZ" ;
V:description = "y-wind component" ;
V:units = "m s-1" ;
V:stagger = "Y" ;
V:coordinates = "XLONG_V XLAT_V" ;
float W(Time, bottom_top_stag, south_north, west_east) ;
W:FieldType = 104 ;
W:MemoryOrder = "XYZ" ;
W:description = "z-wind component" ;
W:units = "m s-1" ;
W:stagger = "Z" ;
W:coordinates = "XLONG XLAT" ;
float PH(Time, bottom_top_stag, south_north, west_east) ;
PH:FieldType = 104 ;
PH:MemoryOrder = "XYZ" ;
PH:description = "perturbation geopotential" ;
PH:units = "m2 s-2" ;
PH:stagger = "Z" ;
PH:coordinates = "XLONG XLAT" ;
float PHB(Time, bottom_top_stag, south_north, west_east) ;
PHB:FieldType = 104 ;
PHB:MemoryOrder = "XYZ" ;
PHB:description = "base-state geopotential" ;
PHB:units = "m2 s-2" ;
PHB:stagger = "Z" ;
PHB:coordinates = "XLONG XLAT" ;
float T(Time, bottom_top, south_north, west_east) ;
T:FieldType = 104 ;
T:MemoryOrder = "XYZ" ;
T:description = "perturbation potential temperature (theta-t0)" ;
T:units = "K" ;
T:stagger = "" ;
T:coordinates = "XLONG XLAT" ;
float MU(Time, south_north, west_east) ;
MU:FieldType = 104 ;
MU:MemoryOrder = "XY " ;
MU:description = "perturbation dry air mass in column" ;
MU:units = "Pa" ;
MU:stagger = "" ;
MU:coordinates = "XLONG XLAT" ;
float MUB(Time, south_north, west_east) ;
MUB:FieldType = 104 ;
MUB:MemoryOrder = "XY " ;
MUB:description = "base state dry air mass in column" ;
MUB:units = "Pa" ;
MUB:stagger = "" ;
MUB:coordinates = "XLONG XLAT" ;
float NEST_POS(Time, south_north, west_east) ;
NEST_POS:FieldType = 104 ;
NEST_POS:MemoryOrder = "XY " ;
NEST_POS:description = "-" ;
NEST_POS:units = "-" ;
NEST_POS:stagger = "" ;
NEST_POS:coordinates = "XLONG XLAT" ;
float P(Time, bottom_top, south_north, west_east) ;
P:FieldType = 104 ;
P:MemoryOrder = "XYZ" ;
P:description = "perturbation pressure" ;
P:units = "Pa" ;
P:stagger = "" ;
P:coordinates = "XLONG XLAT" ;
float PB(Time, bottom_top, south_north, west_east) ;
PB:FieldType = 104 ;
PB:MemoryOrder = "XYZ" ;
PB:description = "BASE STATE PRESSURE" ;
PB:units = "Pa" ;
PB:stagger = "" ;
PB:coordinates = "XLONG XLAT" ;
float FNM(Time, bottom_top) ;
FNM:FieldType = 104 ;
FNM:MemoryOrder = "Z " ;
FNM:description = "upper weight for vertical stretching" ;
FNM:units = "" ;
FNM:stagger = "" ;
float FNP(Time, bottom_top) ;
FNP:FieldType = 104 ;
FNP:MemoryOrder = "Z " ;
FNP:description = "lower weight for vertical stretching" ;
FNP:units = "" ;
FNP:stagger = "" ;
float RDNW(Time, bottom_top) ;
RDNW:FieldType = 104 ;
RDNW:MemoryOrder = "Z " ;
RDNW:description = "inverse d(eta) values between full (w) levels" ;
RDNW:units = "" ;
RDNW:stagger = "" ;
float RDN(Time, bottom_top) ;
RDN:FieldType = 104 ;
RDN:MemoryOrder = "Z " ;
RDN:description = "inverse d(eta) values between half (mass) levels" ;
RDN:units = "" ;
RDN:stagger = "" ;
float DNW(Time, bottom_top) ;
DNW:FieldType = 104 ;
DNW:MemoryOrder = "Z " ;
DNW:description = "d(eta) values between full (w) levels" ;
DNW:units = "" ;
DNW:stagger = "" ;
float DN(Time, bottom_top) ;
DN:FieldType = 104 ;
DN:MemoryOrder = "Z " ;
DN:description = "d(eta) values between half (mass) levels" ;
DN:units = "" ;
DN:stagger = "" ;
float CFN(Time) ;
CFN:FieldType = 104 ;
CFN:MemoryOrder = "0 " ;
CFN:description = "extrapolation constant" ;
CFN:units = "" ;
CFN:stagger = "" ;
float CFN1(Time) ;
CFN1:FieldType = 104 ;
CFN1:MemoryOrder = "0 " ;
CFN1:description = "extrapolation constant" ;
CFN1:units = "" ;
CFN1:stagger = "" ;
float P_HYD(Time, bottom_top, south_north, west_east) ;
P_HYD:FieldType = 104 ;
P_HYD:MemoryOrder = "XYZ" ;
P_HYD:description = "hydrostatic pressure" ;
P_HYD:units = "Pa" ;
P_HYD:stagger = "" ;
P_HYD:coordinates = "XLONG XLAT" ;
float Q2(Time, south_north, west_east) ;
Q2:FieldType = 104 ;
Q2:MemoryOrder = "XY " ;
Q2:description = "QV at 2 M" ;
Q2:units = "kg kg-1" ;
Q2:stagger = "" ;
Q2:coordinates = "XLONG XLAT" ;
float T2(Time, south_north, west_east) ;
T2:FieldType = 104 ;
T2:MemoryOrder = "XY " ;
T2:description = "TEMP at 2 M" ;
T2:units = "K" ;
T2:stagger = "" ;
T2:coordinates = "XLONG XLAT" ;
float TH2(Time, south_north, west_east) ;
TH2:FieldType = 104 ;
TH2:MemoryOrder = "XY " ;
TH2:description = "POT TEMP at 2 M" ;
TH2:units = "K" ;
TH2:stagger = "" ;
TH2:coordinates = "XLONG XLAT" ;
float PSFC(Time, south_north, west_east) ;
PSFC:FieldType = 104 ;
PSFC:MemoryOrder = "XY " ;
PSFC:description = "SFC PRESSURE" ;
PSFC:units = "Pa" ;
PSFC:stagger = "" ;
PSFC:coordinates = "XLONG XLAT" ;
float U10(Time, south_north, west_east) ;
U10:FieldType = 104 ;
U10:MemoryOrder = "XY " ;
U10:description = "U at 10 M" ;
U10:units = "m s-1" ;
U10:stagger = "" ;
U10:coordinates = "XLONG XLAT" ;
float V10(Time, south_north, west_east) ;
V10:FieldType = 104 ;
V10:MemoryOrder = "XY " ;
V10:description = "V at 10 M" ;
V10:units = "m s-1" ;
V10:stagger = "" ;
V10:coordinates = "XLONG XLAT" ;
float RDX(Time) ;
RDX:FieldType = 104 ;
RDX:MemoryOrder = "0 " ;
RDX:description = "INVERSE X GRID LENGTH" ;
RDX:units = "" ;
RDX:stagger = "" ;
float RDY(Time) ;
RDY:FieldType = 104 ;
RDY:MemoryOrder = "0 " ;
RDY:description = "INVERSE Y GRID LENGTH" ;
RDY:units = "" ;
RDY:stagger = "" ;
float RESM(Time) ;
RESM:FieldType = 104 ;
RESM:MemoryOrder = "0 " ;
RESM:description = "TIME WEIGHT CONSTANT FOR SMALL STEPS" ;
RESM:units = "" ;
RESM:stagger = "" ;
float ZETATOP(Time) ;
ZETATOP:FieldType = 104 ;
ZETATOP:MemoryOrder = "0 " ;
ZETATOP:description = "ZETA AT MODEL TOP" ;
ZETATOP:units = "" ;
ZETATOP:stagger = "" ;
float CF1(Time) ;
CF1:FieldType = 104 ;
CF1:MemoryOrder = "0 " ;
CF1:description = "2nd order extrapolation constant" ;
CF1:units = "" ;
CF1:stagger = "" ;
float CF2(Time) ;
CF2:FieldType = 104 ;
CF2:MemoryOrder = "0 " ;
CF2:description = "2nd order extrapolation constant" ;
CF2:units = "" ;
CF2:stagger = "" ;
float CF3(Time) ;
CF3:FieldType = 104 ;
CF3:MemoryOrder = "0 " ;
CF3:description = "2nd order extrapolation constant" ;
CF3:units = "" ;
CF3:stagger = "" ;
int ITIMESTEP(Time) ;
ITIMESTEP:FieldType = 106 ;
ITIMESTEP:MemoryOrder = "0 " ;
ITIMESTEP:description = "" ;
ITIMESTEP:units = "" ;
ITIMESTEP:stagger = "" ;
float XTIME(Time) ;
XTIME:FieldType = 104 ;
XTIME:MemoryOrder = "0 " ;
XTIME:description = "minutes since simulation start" ;
XTIME:units = "" ;
XTIME:stagger = "" ;
float QVAPOR(Time, bottom_top, south_north, west_east) ;
QVAPOR:FieldType = 104 ;
QVAPOR:MemoryOrder = "XYZ" ;
QVAPOR:description = "Water vapor mixing ratio" ;
QVAPOR:units = "kg kg-1" ;
QVAPOR:stagger = "" ;
QVAPOR:coordinates = "XLONG XLAT" ;
float QCLOUD(Time, bottom_top, south_north, west_east) ;
QCLOUD:FieldType = 104 ;
QCLOUD:MemoryOrder = "XYZ" ;
QCLOUD:description = "Cloud water mixing ratio" ;
QCLOUD:units = "kg kg-1" ;
QCLOUD:stagger = "" ;
QCLOUD:coordinates = "XLONG XLAT" ;
float QRAIN(Time, bottom_top, south_north, west_east) ;
QRAIN:FieldType = 104 ;
QRAIN:MemoryOrder = "XYZ" ;
QRAIN:description = "Rain water mixing ratio" ;
QRAIN:units = "kg kg-1" ;
QRAIN:stagger = "" ;
QRAIN:coordinates = "XLONG XLAT" ;
float QICE(Time, bottom_top, south_north, west_east) ;
QICE:FieldType = 104 ;
QICE:MemoryOrder = "XYZ" ;
QICE:description = "Ice mixing ratio" ;
QICE:units = "kg kg-1" ;
QICE:stagger = "" ;
QICE:coordinates = "XLONG XLAT" ;
float QSNOW(Time, bottom_top, south_north, west_east) ;
QSNOW:FieldType = 104 ;
QSNOW:MemoryOrder = "XYZ" ;
QSNOW:description = "Snow mixing ratio" ;
QSNOW:units = "kg kg-1" ;
QSNOW:stagger = "" ;
QSNOW:coordinates = "XLONG XLAT" ;
float QGRAUP(Time, bottom_top, south_north, west_east) ;
QGRAUP:FieldType = 104 ;
QGRAUP:MemoryOrder = "XYZ" ;
QGRAUP:description = "Graupel mixing ratio" ;
QGRAUP:units = "kg kg-1" ;
QGRAUP:stagger = "" ;
QGRAUP:coordinates = "XLONG XLAT" ;
float QNICE(Time, bottom_top, south_north, west_east) ;
QNICE:FieldType = 104 ;
QNICE:MemoryOrder = "XYZ" ;
QNICE:description = "Ice Number concentration" ;
QNICE:units = " kg-1" ;
QNICE:stagger = "" ;
QNICE:coordinates = "XLONG XLAT" ;
float QNRAIN(Time, bottom_top, south_north, west_east) ;
QNRAIN:FieldType = 104 ;
QNRAIN:MemoryOrder = "XYZ" ;
QNRAIN:description = "Rain Number concentration" ;
QNRAIN:units = " kg(-1)" ;
QNRAIN:stagger = "" ;
QNRAIN:coordinates = "XLONG XLAT" ;
float LANDMASK(Time, south_north, west_east) ;
LANDMASK:FieldType = 104 ;
LANDMASK:MemoryOrder = "XY " ;
LANDMASK:description = "LAND MASK (1 FOR LAND, 0 FOR WATER)" ;
LANDMASK:units = "" ;
LANDMASK:stagger = "" ;
LANDMASK:coordinates = "XLONG XLAT" ;
float TSLB(Time, soil_layers_stag, south_north, west_east) ;
TSLB:FieldType = 104 ;
TSLB:MemoryOrder = "XYZ" ;
TSLB:description = "SOIL TEMPERATURE" ;
TSLB:units = "K" ;
TSLB:stagger = "Z" ;
TSLB:coordinates = "XLONG XLAT" ;
float SMOIS(Time, soil_layers_stag, south_north, west_east) ;
SMOIS:FieldType = 104 ;
SMOIS:MemoryOrder = "XYZ" ;
SMOIS:description = "SOIL MOISTURE" ;
SMOIS:units = "m3 m-3" ;
SMOIS:stagger = "Z" ;
SMOIS:coordinates = "XLONG XLAT" ;
float SH2O(Time, soil_layers_stag, south_north, west_east) ;
SH2O:FieldType = 104 ;
SH2O:MemoryOrder = "XYZ" ;
SH2O:description = "SOIL LIQUID WATER" ;
SH2O:units = "m3 m-3" ;
SH2O:stagger = "Z" ;
SH2O:coordinates = "XLONG XLAT" ;
float SMCREL(Time, soil_layers_stag, south_north, west_east) ;
SMCREL:FieldType = 104 ;
SMCREL:MemoryOrder = "XYZ" ;
SMCREL:description = "RELATIVE SOIL MOISTURE" ;
SMCREL:units = "" ;
SMCREL:stagger = "Z" ;
SMCREL:coordinates = "XLONG XLAT" ;
float SEAICE(Time, south_north, west_east) ;
SEAICE:FieldType = 104 ;
SEAICE:MemoryOrder = "XY " ;
SEAICE:description = "SEA ICE FLAG" ;
SEAICE:units = "" ;
SEAICE:stagger = "" ;
SEAICE:coordinates = "XLONG XLAT" ;
float XICEM(Time, south_north, west_east) ;
XICEM:FieldType = 104 ;
XICEM:MemoryOrder = "XY " ;
XICEM:description = "SEA ICE FLAG (PREVIOUS STEP)" ;
XICEM:units = "" ;
XICEM:stagger = "" ;
XICEM:coordinates = "XLONG XLAT" ;
float SFROFF(Time, south_north, west_east) ;
SFROFF:FieldType = 104 ;
SFROFF:MemoryOrder = "XY " ;
SFROFF:description = "SURFACE RUNOFF" ;
SFROFF:units = "mm" ;
SFROFF:stagger = "" ;
SFROFF:coordinates = "XLONG XLAT" ;
float UDROFF(Time, south_north, west_east) ;
UDROFF:FieldType = 104 ;
UDROFF:MemoryOrder = "XY " ;
UDROFF:description = "UNDERGROUND RUNOFF" ;
UDROFF:units = "mm" ;
UDROFF:stagger = "" ;
UDROFF:coordinates = "XLONG XLAT" ;
int IVGTYP(Time, south_north, west_east) ;
IVGTYP:FieldType = 106 ;
IVGTYP:MemoryOrder = "XY " ;
IVGTYP:description = "DOMINANT VEGETATION CATEGORY" ;
IVGTYP:units = "" ;
IVGTYP:stagger = "" ;
IVGTYP:coordinates = "XLONG XLAT" ;
int ISLTYP(Time, south_north, west_east) ;
ISLTYP:FieldType = 106 ;
ISLTYP:MemoryOrder = "XY " ;
ISLTYP:description = "DOMINANT SOIL CATEGORY" ;
ISLTYP:units = "" ;
ISLTYP:stagger = "" ;
ISLTYP:coordinates = "XLONG XLAT" ;
float VEGFRA(Time, south_north, west_east) ;
VEGFRA:FieldType = 104 ;
VEGFRA:MemoryOrder = "XY " ;
VEGFRA:description = "VEGETATION FRACTION" ;
VEGFRA:units = "" ;
VEGFRA:stagger = "" ;
VEGFRA:coordinates = "XLONG XLAT" ;
float GRDFLX(Time, south_north, west_east) ;
GRDFLX:FieldType = 104 ;
GRDFLX:MemoryOrder = "XY " ;
GRDFLX:description = "GROUND HEAT FLUX" ;
GRDFLX:units = "W m-2" ;
GRDFLX:stagger = "" ;
GRDFLX:coordinates = "XLONG XLAT" ;
float ACGRDFLX(Time, south_north, west_east) ;
ACGRDFLX:FieldType = 104 ;
ACGRDFLX:MemoryOrder = "XY " ;
ACGRDFLX:description = "ACCUMULATED GROUND HEAT FLUX" ;
ACGRDFLX:units = "J m-2" ;
ACGRDFLX:stagger = "" ;
ACGRDFLX:coordinates = "XLONG XLAT" ;
float SNOW(Time, south_north, west_east) ;
SNOW:FieldType = 104 ;
SNOW:MemoryOrder = "XY " ;
SNOW:description = "SNOW WATER EQUIVALENT" ;
SNOW:units = "kg m-2" ;
SNOW:stagger = "" ;
SNOW:coordinates = "XLONG XLAT" ;
float SNOWH(Time, south_north, west_east) ;
SNOWH:FieldType = 104 ;
SNOWH:MemoryOrder = "XY " ;
SNOWH:description = "PHYSICAL SNOW DEPTH" ;
SNOWH:units = "m" ;
SNOWH:stagger = "" ;
SNOWH:coordinates = "XLONG XLAT" ;
float CANWAT(Time, south_north, west_east) ;
CANWAT:FieldType = 104 ;
CANWAT:MemoryOrder = "XY " ;
CANWAT:description = "CANOPY WATER" ;
CANWAT:units = "kg m-2" ;
CANWAT:stagger = "" ;
CANWAT:coordinates = "XLONG XLAT" ;
float SST(Time, south_north, west_east) ;
SST:FieldType = 104 ;
SST:MemoryOrder = "XY " ;
SST:description = "SEA SURFACE TEMPERATURE" ;
SST:units = "K" ;
SST:stagger = "" ;
SST:coordinates = "XLONG XLAT" ;
float SSTSK(Time, south_north, west_east) ;
SSTSK:FieldType = 104 ;
SSTSK:MemoryOrder = "XY " ;
SSTSK:description = "SKIN SEA SURFACE TEMPERATURE" ;
SSTSK:units = "K" ;
SSTSK:stagger = "" ;
SSTSK:coordinates = "XLONG XLAT" ;
float LAI(Time, south_north, west_east) ;
LAI:FieldType = 104 ;
LAI:MemoryOrder = "XY " ;
LAI:description = "Leaf area index" ;
LAI:units = "area/area" ;
LAI:stagger = "" ;
LAI:coordinates = "XLONG XLAT" ;
float TKE_MYJ(Time, bottom_top, south_north, west_east) ;
TKE_MYJ:FieldType = 104 ;
TKE_MYJ:MemoryOrder = "XYZ" ;
TKE_MYJ:description = "TKE FROM MELLOR-YAMADA-JANJIC" ;
TKE_MYJ:units = "m2 s-2" ;
TKE_MYJ:stagger = "" ;
TKE_MYJ:coordinates = "XLONG XLAT" ;
float EL_MYJ(Time, bottom_top, south_north, west_east) ;
EL_MYJ:FieldType = 104 ;
EL_MYJ:MemoryOrder = "XYZ" ;
EL_MYJ:description = "MIXING LENGTH FROM MELLOR-YAMADA-JANJIC" ;
EL_MYJ:units = "m" ;
EL_MYJ:stagger = "" ;
EL_MYJ:coordinates = "XLONG XLAT" ;
float MAPFAC_M(Time, south_north, west_east) ;
MAPFAC_M:FieldType = 104 ;
MAPFAC_M:MemoryOrder = "XY " ;
MAPFAC_M:description = "Map scale factor on mass grid" ;
MAPFAC_M:units = "" ;
MAPFAC_M:stagger = "" ;
MAPFAC_M:coordinates = "XLONG XLAT" ;
float MAPFAC_U(Time, south_north, west_east_stag) ;
MAPFAC_U:FieldType = 104 ;
MAPFAC_U:MemoryOrder = "XY " ;
MAPFAC_U:description = "Map scale factor on u-grid" ;
MAPFAC_U:units = "" ;
MAPFAC_U:stagger = "X" ;
MAPFAC_U:coordinates = "XLONG_U XLAT_U" ;
float MAPFAC_V(Time, south_north_stag, west_east) ;
MAPFAC_V:FieldType = 104 ;
MAPFAC_V:MemoryOrder = "XY " ;
MAPFAC_V:description = "Map scale factor on v-grid" ;
MAPFAC_V:units = "" ;
MAPFAC_V:stagger = "Y" ;
MAPFAC_V:coordinates = "XLONG_V XLAT_V" ;
float MAPFAC_MX(Time, south_north, west_east) ;
MAPFAC_MX:FieldType = 104 ;
MAPFAC_MX:MemoryOrder = "XY " ;
MAPFAC_MX:description = "Map scale factor on mass grid, x direction" ;
MAPFAC_MX:units = "" ;
MAPFAC_MX:stagger = "" ;
MAPFAC_MX:coordinates = "XLONG XLAT" ;
float MAPFAC_MY(Time, south_north, west_east) ;
MAPFAC_MY:FieldType = 104 ;
MAPFAC_MY:MemoryOrder = "XY " ;
MAPFAC_MY:description = "Map scale factor on mass grid, y direction" ;
MAPFAC_MY:units = "" ;
MAPFAC_MY:stagger = "" ;
MAPFAC_MY:coordinates = "XLONG XLAT" ;
float MAPFAC_UX(Time, south_north, west_east_stag) ;
MAPFAC_UX:FieldType = 104 ;
MAPFAC_UX:MemoryOrder = "XY " ;
MAPFAC_UX:description = "Map scale factor on u-grid, x direction" ;
MAPFAC_UX:units = "" ;
MAPFAC_UX:stagger = "X" ;
MAPFAC_UX:coordinates = "XLONG_U XLAT_U" ;
float MAPFAC_UY(Time, south_north, west_east_stag) ;
MAPFAC_UY:FieldType = 104 ;
MAPFAC_UY:MemoryOrder = "XY " ;
MAPFAC_UY:description = "Map scale factor on u-grid, y direction" ;
MAPFAC_UY:units = "" ;
MAPFAC_UY:stagger = "X" ;
MAPFAC_UY:coordinates = "XLONG_U XLAT_U" ;
float MAPFAC_VX(Time, south_north_stag, west_east) ;
MAPFAC_VX:FieldType = 104 ;
MAPFAC_VX:MemoryOrder = "XY " ;
MAPFAC_VX:description = "Map scale factor on v-grid, x direction" ;
MAPFAC_VX:units = "" ;
MAPFAC_VX:stagger = "Y" ;
MAPFAC_VX:coordinates = "XLONG_V XLAT_V" ;
float MF_VX_INV(Time, south_north_stag, west_east) ;
MF_VX_INV:FieldType = 104 ;
MF_VX_INV:MemoryOrder = "XY " ;
MF_VX_INV:description = "Inverse map scale factor on v-grid, x direction" ;
MF_VX_INV:units = "" ;
MF_VX_INV:stagger = "Y" ;
MF_VX_INV:coordinates = "XLONG_V XLAT_V" ;
float MAPFAC_VY(Time, south_north_stag, west_east) ;
MAPFAC_VY:FieldType = 104 ;
MAPFAC_VY:MemoryOrder = "XY " ;
MAPFAC_VY:description = "Map scale factor on v-grid, y direction" ;
MAPFAC_VY:units = "" ;
MAPFAC_VY:stagger = "Y" ;
MAPFAC_VY:coordinates = "XLONG_V XLAT_V" ;
float F(Time, south_north, west_east) ;
F:FieldType = 104 ;
F:MemoryOrder = "XY " ;
F:description = "Coriolis sine latitude term" ;
F:units = "s-1" ;
F:stagger = "" ;
F:coordinates = "XLONG XLAT" ;
float E(Time, south_north, west_east) ;
E:FieldType = 104 ;
E:MemoryOrder = "XY " ;
E:description = "Coriolis cosine latitude term" ;
E:units = "s-1" ;
E:stagger = "" ;
E:coordinates = "XLONG XLAT" ;
float SINALPHA(Time, south_north, west_east) ;
SINALPHA:FieldType = 104 ;
SINALPHA:MemoryOrder = "XY " ;
SINALPHA:description = "Local sine of map rotation" ;
SINALPHA:units = "" ;
SINALPHA:stagger = "" ;
SINALPHA:coordinates = "XLONG XLAT" ;
float COSALPHA(Time, south_north, west_east) ;
COSALPHA:FieldType = 104 ;
COSALPHA:MemoryOrder = "XY " ;
COSALPHA:description = "Local cosine of map rotation" ;
COSALPHA:units = "" ;
COSALPHA:stagger = "" ;
COSALPHA:coordinates = "XLONG XLAT" ;
float HGT(Time, south_north, west_east) ;
HGT:FieldType = 104 ;
HGT:MemoryOrder = "XY " ;
HGT:description = "Terrain Height" ;
HGT:units = "m" ;
HGT:stagger = "" ;
HGT:coordinates = "XLONG XLAT" ;
float TSK(Time, south_north, west_east) ;
TSK:FieldType = 104 ;
TSK:MemoryOrder = "XY " ;
TSK:description = "SURFACE SKIN TEMPERATURE" ;
TSK:units = "K" ;
TSK:stagger = "" ;
TSK:coordinates = "XLONG XLAT" ;
float P_TOP(Time) ;
P_TOP:FieldType = 104 ;
P_TOP:MemoryOrder = "0 " ;
P_TOP:description = "PRESSURE TOP OF THE MODEL" ;
P_TOP:units = "Pa" ;
P_TOP:stagger = "" ;
float T00(Time) ;
T00:FieldType = 104 ;
T00:MemoryOrder = "0 " ;
T00:description = "BASE STATE TEMPERATURE" ;
T00:units = "K" ;
T00:stagger = "" ;
float P00(Time) ;
P00:FieldType = 104 ;
P00:MemoryOrder = "0 " ;
P00:description = "BASE STATE PRESURE" ;
P00:units = "Pa" ;
P00:stagger = "" ;
float TLP(Time) ;
TLP:FieldType = 104 ;
TLP:MemoryOrder = "0 " ;
TLP:description = "BASE STATE LAPSE RATE" ;
TLP:units = "" ;
TLP:stagger = "" ;
float TISO(Time) ;
TISO:FieldType = 104 ;
TISO:MemoryOrder = "0 " ;
TISO:description = "TEMP AT WHICH THE BASE T TURNS CONST" ;
TISO:units = "K" ;
TISO:stagger = "" ;
float MAX_MSTFX(Time) ;
MAX_MSTFX:FieldType = 104 ;
MAX_MSTFX:MemoryOrder = "0 " ;
MAX_MSTFX:description = "Max map factor in domain" ;
MAX_MSTFX:units = "" ;
MAX_MSTFX:stagger = "" ;
float MAX_MSTFY(Time) ;
MAX_MSTFY:FieldType = 104 ;
MAX_MSTFY:MemoryOrder = "0 " ;
MAX_MSTFY:description = "Max map factor in domain" ;
MAX_MSTFY:units = "" ;
MAX_MSTFY:stagger = "" ;
float RAINC(Time, south_north, west_east) ;
RAINC:FieldType = 104 ;
RAINC:MemoryOrder = "XY " ;
RAINC:description = "ACCUMULATED TOTAL CUMULUS PRECIPITATION" ;
RAINC:units = "mm" ;
RAINC:stagger = "" ;
RAINC:coordinates = "XLONG XLAT" ;
float RAINNC(Time, south_north, west_east) ;
RAINNC:FieldType = 104 ;
RAINNC:MemoryOrder = "XY " ;
RAINNC:description = "ACCUMULATED TOTAL GRID SCALE PRECIPITATION" ;
RAINNC:units = "mm" ;
RAINNC:stagger = "" ;
RAINNC:coordinates = "XLONG XLAT" ;
float SNOWNC(Time, south_north, west_east) ;
SNOWNC:FieldType = 104 ;
SNOWNC:MemoryOrder = "XY " ;
SNOWNC:description = "ACCUMULATED TOTAL GRID SCALE SNOW AND ICE" ;
SNOWNC:units = "mm" ;
SNOWNC:stagger = "" ;
SNOWNC:coordinates = "XLONG XLAT" ;
float GRAUPELNC(Time, south_north, west_east) ;
GRAUPELNC:FieldType = 104 ;
GRAUPELNC:MemoryOrder = "XY " ;
GRAUPELNC:description = "ACCUMULATED TOTAL GRID SCALE GRAUPEL" ;
GRAUPELNC:units = "mm" ;
GRAUPELNC:stagger = "" ;
GRAUPELNC:coordinates = "XLONG XLAT" ;
float HAILNC(Time, south_north, west_east) ;
HAILNC:FieldType = 104 ;
HAILNC:MemoryOrder = "XY " ;
HAILNC:description = "ACCUMULATED TOTAL GRID SCALE HAIL" ;
HAILNC:units = "mm" ;
HAILNC:stagger = "" ;
HAILNC:coordinates = "XLONG XLAT" ;
float CLDFRA(Time, bottom_top, south_north, west_east) ;
CLDFRA:FieldType = 104 ;
CLDFRA:MemoryOrder = "XYZ" ;
CLDFRA:description = "CLOUD FRACTION" ;
CLDFRA:units = "" ;
CLDFRA:stagger = "" ;
CLDFRA:coordinates = "XLONG XLAT" ;
float SWDOWN(Time, south_north, west_east) ;
SWDOWN:FieldType = 104 ;
SWDOWN:MemoryOrder = "XY " ;
SWDOWN:description = "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE" ;
SWDOWN:units = "W m-2" ;
SWDOWN:stagger = "" ;
SWDOWN:coordinates = "XLONG XLAT" ;
float GLW(Time, south_north, west_east) ;
GLW:FieldType = 104 ;
GLW:MemoryOrder = "XY " ;
GLW:description = "DOWNWARD LONG WAVE FLUX AT GROUND SURFACE" ;
GLW:units = "W m-2" ;
GLW:stagger = "" ;
GLW:coordinates = "XLONG XLAT" ;
float SWNORM(Time, south_north, west_east) ;
SWNORM:FieldType = 104 ;
SWNORM:MemoryOrder = "XY " ;
SWNORM:description = "NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT)" ;
SWNORM:units = "W m-2" ;
SWNORM:stagger = "" ;
SWNORM:coordinates = "XLONG XLAT" ;
float OLR(Time, south_north, west_east) ;
OLR:FieldType = 104 ;
OLR:MemoryOrder = "XY " ;
OLR:description = "TOA OUTGOING LONG WAVE" ;
OLR:units = "W m-2" ;
OLR:stagger = "" ;
OLR:coordinates = "XLONG XLAT" ;
float XLAT(Time, south_north, west_east) ;
XLAT:FieldType = 104 ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degree_north" ;
XLAT:stagger = "" ;
float XLONG(Time, south_north, west_east) ;
XLONG:FieldType = 104 ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degree_east" ;
XLONG:stagger = "" ;
float XLAT_U(Time, south_north, west_east_stag) ;
XLAT_U:FieldType = 104 ;
XLAT_U:MemoryOrder = "XY " ;
XLAT_U:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT_U:units = "degree_north" ;
XLAT_U:stagger = "X" ;
XLAT_U:coordinates = "XLONG_U XLAT_U" ;
float XLONG_U(Time, south_north, west_east_stag) ;
XLONG_U:FieldType = 104 ;
XLONG_U:MemoryOrder = "XY " ;
XLONG_U:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG_U:units = "degree_east" ;
XLONG_U:stagger = "X" ;
XLONG_U:coordinates = "XLONG_U XLAT_U" ;
float XLAT_V(Time, south_north_stag, west_east) ;
XLAT_V:FieldType = 104 ;
XLAT_V:MemoryOrder = "XY " ;
XLAT_V:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT_V:units = "degree_north" ;
XLAT_V:stagger = "Y" ;
XLAT_V:coordinates = "XLONG_V XLAT_V" ;
float XLONG_V(Time, south_north_stag, west_east) ;
XLONG_V:FieldType = 104 ;
XLONG_V:MemoryOrder = "XY " ;
XLONG_V:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG_V:units = "degree_east" ;
XLONG_V:stagger = "Y" ;
XLONG_V:coordinates = "XLONG_V XLAT_V" ;
float ALBEDO(Time, south_north, west_east) ;
ALBEDO:FieldType = 104 ;
ALBEDO:MemoryOrder = "XY " ;
ALBEDO:description = "ALBEDO" ;
ALBEDO:units = "-" ;
ALBEDO:stagger = "" ;
ALBEDO:coordinates = "XLONG XLAT" ;
float ALBBCK(Time, south_north, west_east) ;
ALBBCK:FieldType = 104 ;
ALBBCK:MemoryOrder = "XY " ;
ALBBCK:description = "BACKGROUND ALBEDO" ;
ALBBCK:units = "" ;
ALBBCK:stagger = "" ;
ALBBCK:coordinates = "XLONG XLAT" ;
float EMISS(Time, south_north, west_east) ;
EMISS:FieldType = 104 ;
EMISS:MemoryOrder = "XY " ;
EMISS:description = "SURFACE EMISSIVITY" ;
EMISS:units = "" ;
EMISS:stagger = "" ;
EMISS:coordinates = "XLONG XLAT" ;
float NOAHRES(Time, south_north, west_east) ;
NOAHRES:FieldType = 104 ;
NOAHRES:MemoryOrder = "XY " ;
NOAHRES:description = "RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET" ;
NOAHRES:units = "W m{-2}" ;
NOAHRES:stagger = "" ;
NOAHRES:coordinates = "XLONG XLAT" ;
float TMN(Time, south_north, west_east) ;
TMN:FieldType = 104 ;
TMN:MemoryOrder = "XY " ;
TMN:description = "SOIL TEMPERATURE AT LOWER BOUNDARY" ;
TMN:units = "K" ;
TMN:stagger = "" ;
TMN:coordinates = "XLONG XLAT" ;
float XLAND(Time, south_north, west_east) ;
XLAND:FieldType = 104 ;
XLAND:MemoryOrder = "XY " ;
XLAND:description = "LAND MASK (1 FOR LAND, 2 FOR WATER)" ;
XLAND:units = "" ;
XLAND:stagger = "" ;
XLAND:coordinates = "XLONG XLAT" ;
float UST(Time, south_north, west_east) ;
UST:FieldType = 104 ;
UST:MemoryOrder = "XY " ;
UST:description = "U* IN SIMILARITY THEORY" ;
UST:units = "m s-1" ;
UST:stagger = "" ;
UST:coordinates = "XLONG XLAT" ;
float PBLH(Time, south_north, west_east) ;
PBLH:FieldType = 104 ;
PBLH:MemoryOrder = "XY " ;
PBLH:description = "PBL HEIGHT" ;
PBLH:units = "m" ;
PBLH:stagger = "" ;
PBLH:coordinates = "XLONG XLAT" ;
float HFX(Time, south_north, west_east) ;
HFX:FieldType = 104 ;
HFX:MemoryOrder = "XY " ;
HFX:description = "UPWARD HEAT FLUX AT THE SURFACE" ;
HFX:units = "W m-2" ;
HFX:stagger = "" ;
HFX:coordinates = "XLONG XLAT" ;
float QFX(Time, south_north, west_east) ;
QFX:FieldType = 104 ;
QFX:MemoryOrder = "XY " ;
QFX:description = "UPWARD MOISTURE FLUX AT THE SURFACE" ;
QFX:units = "kg m-2 s-1" ;
QFX:stagger = "" ;
QFX:coordinates = "XLONG XLAT" ;
float LH(Time, south_north, west_east) ;
LH:FieldType = 104 ;
LH:MemoryOrder = "XY " ;
LH:description = "LATENT HEAT FLUX AT THE SURFACE" ;
LH:units = "W m-2" ;
LH:stagger = "" ;
LH:coordinates = "XLONG XLAT" ;
float ACHFX(Time, south_north, west_east) ;
ACHFX:FieldType = 104 ;
ACHFX:MemoryOrder = "XY " ;
ACHFX:description = "ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE" ;
ACHFX:units = "J m-2" ;
ACHFX:stagger = "" ;
ACHFX:coordinates = "XLONG XLAT" ;
float ACLHF(Time, south_north, west_east) ;
ACLHF:FieldType = 104 ;
ACLHF:MemoryOrder = "XY " ;
ACLHF:description = "ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE" ;
ACLHF:units = "J m-2" ;
ACLHF:stagger = "" ;
ACLHF:coordinates = "XLONG XLAT" ;
float SNOWC(Time, south_north, west_east) ;
SNOWC:FieldType = 104 ;
SNOWC:MemoryOrder = "XY " ;
SNOWC:description = "FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)" ;
SNOWC:units = "" ;
SNOWC:stagger = "" ;
SNOWC:coordinates = "XLONG XLAT" ;
float SR(Time, south_north, west_east) ;
SR:FieldType = 104 ;
SR:MemoryOrder = "XY " ;
SR:description = "fraction of frozen precipitation" ;
SR:units = "-" ;
SR:stagger = "" ;
SR:coordinates = "XLONG XLAT" ;
float POTEVP(Time, south_north, west_east) ;
POTEVP:FieldType = 104 ;
POTEVP:MemoryOrder = "XY " ;
POTEVP:description = "accumulated potential evaporation" ;
POTEVP:units = "W m-2" ;
POTEVP:stagger = "" ;
POTEVP:coordinates = "XLONG XLAT" ;
float SNOPCX(Time, south_north, west_east) ;
SNOPCX:FieldType = 104 ;
SNOPCX:MemoryOrder = "XY " ;
SNOPCX:description = "snow phase change heat flux" ;
SNOPCX:units = "W m-2" ;
SNOPCX:stagger = "" ;
SNOPCX:coordinates = "XLONG XLAT" ;
float SOILTB(Time, south_north, west_east) ;
SOILTB:FieldType = 104 ;
SOILTB:MemoryOrder = "XY " ;
SOILTB:description = "bottom soil temperature" ;
SOILTB:units = "K" ;
SOILTB:stagger = "" ;
SOILTB:coordinates = "XLONG XLAT" ;
int SAVE_TOPO_FROM_REAL(Time) ;
SAVE_TOPO_FROM_REAL:FieldType = 106 ;
SAVE_TOPO_FROM_REAL:MemoryOrder = "0 " ;
SAVE_TOPO_FROM_REAL:description = "1=original topo from real/0=topo modified by WRF" ;
SAVE_TOPO_FROM_REAL:units = "flag" ;
SAVE_TOPO_FROM_REAL:stagger = "" ;

// global attributes:
:TITLE = " OUTPUT FROM WRF V3.2.1 MODEL" ;
:START_DATE = "2010-08-20_13:00:00" ;
:SIMULATION_START_DATE = "2010-08-20_13:00:00" ;
:WEST-EAST_GRID_DIMENSION = 121 ;
:SOUTH-NORTH_GRID_DIMENSION = 121 ;
:BOTTOM-TOP_GRID_DIMENSION = 90 ;
:DX = 1000.f ;
:DY = 1000.f ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 1 ;
:KM_OPT = 4 ;
:DAMP_OPT = 3 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 8 ;
:RA_LW_PHYSICS = 1 ;
:RA_SW_PHYSICS = 1 ;
:SF_SFCLAY_PHYSICS = 4 ;
:SF_SURFACE_PHYSICS = 2 ;
:BL_PBL_PHYSICS = 4 ;
:CU_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 1 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 0 ;
:GFDDA_INTERVAL_M = 0 ;
:GFDDA_END_H = 0 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:FEEDBACK = 0 ;
:SMOOTH_OPTION = 0 ;
:SWRAD_SCAT = 1.f ;
:W_DAMPING = 1 ;
:MOIST_ADV_OPT = 1 ;
:SCALAR_ADV_OPT = 1 ;
:TKE_ADV_OPT = 1 ;
:DIFF_6TH_OPT = 2 ;
:DIFF_6TH_FACTOR = 0.12f ;
:OBS_NUDGE_OPT = 1 ;
:FDDA_END = 240.f ;
:OBS_NUDGE_WIND = 1 ;
:OBS_COEF_WIND = 0.001f ;
:OBS_NUDGE_TEMP = 1 ;
:OBS_COEF_TEMP = 0.001f ;
:OBS_NUDGE_MOIS = 1 ;
:OBS_COEF_MOIS = 0.001f ;
:OBS_NUDGE_PSTR = 0 ;
:OBS_COEF_PSTR = 0.f ;
:OBS_IONF = 1 ;
:OBS_IDYNIN = 1 ;
:OBS_DTRAMP = 40.f ;
:BUCKET_MM = -1.f ;
:BUCKET_J = -1.f ;
:PREC_ACC_DT = 0.f ;
:OMLCALL = 0 ;
:ISFTCFLX = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 120 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 121 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 120 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 121 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 89 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 90 ;
:GRID_ID = 2 ;
:PARENT_ID = 1 ;
:I_PARENT_START = 61 ;
:J_PARENT_START = 61 ;
:PARENT_GRID_RATIO = 3 ;
:DT = 1.f ;
:CEN_LAT = 32.217f ;
:CEN_LON = -107.167f ;
:TRUELAT1 = 32.217f ;
:TRUELAT2 = 32.217f ;
:MOAD_CEN_LAT = 32.217f ;
:STAND_LON = -107.167f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 13.f ;
:JULYR = 2010 ;
:JULDAY = 232 ;
:MAP_PROJ = 1 ;
:MMINLU = "USGS" ;
:NUM_LAND_CAT = 24 ;
:ISWATER = 16 ;
:ISLAKE = -1 ;
:ISICE = 24 ;
:ISURBAN = 1 ;
:ISOILWATER = 14 ;
}

Top
October 23, 2014 - 5:12pm
clyne

Steve,

 

This is a bug. A bug report has been opened, and it describes a partial workaround for the problem:

 

https://sourceforge.net/p/vapor/bugs/1041/

 

We will try to have this fixed in the 2.4 release of VAPOR, which should be available by the end of the year.

Top
October 27, 2014 - 9:47am
(Reply to #5)
steve.kirby
Hi, Thanks for letting me know. Looking forward to a code fix for this. In the interim, as a temporary workaround, I turned off cropping by setting the last arg of wrf2geotiff_write to FALSE. However, now my final TIFF file is not being created. (The "_temp.tiff" TIFF file is ok). This call is the culprit: tiff2geotiff -4 " +proj=lcc +ellps=WGS84 +lon_0=-107.167 +lat_0=32.217 +lat_1=32.217 +lat_2=32.217" -m ./steve_plot.dat ./steve_plot_temp.tiff ./steve_plot.tiff It is giving me a segmentation fault. Could the fact that all three lat values being the same be an issue? steve_plot.dat shows these values: 2010-08-20_20:00:00 -107.796 31.6802 -106.531 32.7506 7.85233 7.47339 23.5151 19.5763 Let me know if you see anything that could be a problem. TIA, Steve
Top
October 27, 2014 - 4:47pm
(Reply to #6)
clyne

I can't reproduce this. Can you post your ./steve_plot_temp.tiff file somehow?

Top
October 27, 2014 - 6:27pm
(Reply to #7)
steve.kirby
Hi, Is there a way to attach files in this forum? If not, I could email it to you. thanks, Steve
Top
October 28, 2014 - 9:28am
(Reply to #8)
steve.kirby
Hi, I clicked on your email link in this forum, but didn't see a way to attach a file. I'm surprised there is no way to attach files, either in the forum itself, or via the email links. Could you provide an email? Or is there a better way?
Top
October 28, 2014 - 9:35am
(Reply to #9)
clyne

You can mail it to vapor@ucar.edu

Top
October 28, 2014 - 4:14pm
(Reply to #10)
steve.kirby
Just tried to send, but it was rejected. Here's what I got: The message you are sending is too large. All messages must be less than 5.00 MBytes. The max size was set by the owners of this list (vapor+owner@ucar.edu). If you have questions related to this or any other Google Group, visit the Help Center at http://support.google.com/a/ucar.edu/bin/topic.py?topic=25838 The file size is about 7MB.
Top
October 28, 2014 - 4:25pm
(Reply to #11)
clyne

Can you put it on DropBox or GoogleDrive?

Top
October 30, 2014 - 12:52pm
(Reply to #12)
steve.kirby

I got the file sent to vapor@ucar.edu

Top
November 3, 2014 - 8:40am
(Reply to #13)
clyne

Steve, I'venot received anything.

Top
November 3, 2014 - 10:01am
(Reply to #14)
steve.kirby

Hi,

I just uploaded it again. Can you check please. There should be a password for you to access.

--Steve

Top
November 3, 2014 - 10:34am
(Reply to #15)
clyne

Uploaded it where? I thought you were trying to email the image. I've not received any email.

Top
November 4, 2014 - 9:26am
(Reply to #16)
clyne

This is a bug in tiff-4.0.0. It has been fixed in tiff-4.0.3. The next release of VAPOR (version 2.4.0) which will be released in December, will be built against tiff-4.0.3.

A possible workaround for vapor2.3.0 pre-compiled binaries may be to replace the tiff shared library in $VAPOR_HOME/lib (libtiff.so) with the tiff-4.0.3 libtiff.so

Top
August 27, 2015 - 10:09am
alan

We now have a new version of wrf2geotiff.ncl that fixes the image alignment problem.  This will be available in the next vapor release, later this year.  If you need it sooner, send an email to vapor@ucar.edu and I'll send you a copy.

-Alan

Top

VAPOR Sponsors

SPONSORS

VAPOR sponsors