Garbled output from ARWPost

[Oscar van der Velde]

I am now trying ARWPost 2.0 to get output for GrADS for my WRF 3.0 run.
The program runs successfully though the files and writes one output file and the ctl. I ran GrADS 1.8, opened the ctl file and tried to display the parameters (which I have requested for an array of pressure levels only). The output is very odd.

For example:
open my.ctl
set t 7
d dbz
output: "Constant field. Value = 0"

set z 2
d dbz
output: a plot of dBZ values.

set z 6
d dbz
output: a plot with okay contours but with very low values up to 0.45!

equally, for another time where I specify the same level I got output for, it shows "constant field" again.
Not only dbz but any other parameter: u and v are at different z levels, or "constant field".
For your reference here is the ctl file:

Code: Select all

dset ^arwpost_d01.dat
undef 1.e30
options byteswapped
title  OUTPUT FROM WRF V3.0 MODEL
pdef   89  75 lcc  42.094    1.049   45.000   38.000  42.09400  42.09400    1.04900   8000.000   8000.000
xdef  266 linear   -3.76067   0.03603604
ydef  169 linear   38.99683   0.03603604
zdef   25 levels 
1000.00000
 975.00000
 950.00000
 925.00000
 900.00000
 875.00000
 850.00000
 825.00000
 800.00000
 775.00000
 750.00000
 725.00000
 700.00000
 650.00000
 600.00000
 550.00000
 500.00000
 450.00000
 400.00000
 350.00000
 300.00000
 250.00000
 200.00000
 150.00000
 100.00000
tdef   14 linear 06Z18APR2008      60MN     
VARS   49
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             25  0  x-wind component (m s-1)
V             25  0  y-wind component (m s-1)
W             25  0  z-wind component (m s-1)
PH            25  0  perturbation geopotential (m2 s-2)
PHB           25  0  base-state geopotential (m2 s-2)
T             25  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
P             25  0  perturbation pressure (Pa)
PB            25  0  BASE STATE PRESSURE (Pa)
T2             1  0  TEMP at 2 M (K)
TH2            1  0  POT TEMP at 2 M (K)
PSFC           1  0  SFC PRESSURE (Pa)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
PBLH           1  0  PBL HEIGHT (m)
pressure      25  0  Model pressure (hPa)
height        25  0  Model height (KM)
tk            25  0  Temperature (K)
theta         25  0  Potential Temperature (K)
rh            25  0  Relative Humidity (%)
clflo          1  0  Low Cloud Fraction (%)
clfmi          1  0  Mid Cloud Fraction (%)
clfhi          1  0  High Cloud Fraction (%)
rh2            1  0  Relative Humidity at 2m (%)
umet          25  0  Rotated wind component (m s-1)
vmet          25  0  Rotated wind component (m s-1)
u10m           1  0  Rotated wind component (m s-1)
v10m           1  0  Rotated wind component (m s-1)
slp            1  0  Sea Levelp Pressure (hPa)
dbz           25  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)
cape          25  0  CAPE (J/kg)
cin           25  0  CIN (J/kg)
mcape          1  0  MCAPE (J/kg)
mcin           1  0  MCIN (J/kg)
lcl            1  0  LCL (meters AGL)
lfc            1  0  LFC (meters AGL)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.0 MODEL
@ global String comment SIMULATION_START_DATE = 2008-04-18_06:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION =    90
@ global String comment SOUTH-NORTH_GRID_DIMENSION =    76
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    50
@ global String comment MAP_PROJ =     1
@ global String comment DX =      8000.00
@ global String comment DY =      8000.00
@ global String comment CEN_LAT =        42.09
@ global String comment CEN_LON =         1.05
@ global String comment TRUELAT1 =        42.09
@ global String comment TRUELAT2 =        42.09
@ global String comment MOAD_CEN_LAT =        42.09
@ global String comment STAND_LON =         1.05
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =    10
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     1
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     1
@ global String comment CU_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment UCMCALL =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment W_DAMPING =     0
@ global String comment PD_MOIST =     1
@ global String comment PD_SCALAR =     1
@ global String comment PD_TKE =     0
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment GRID_ID =     1
@ global String comment PARENT_ID =     0
@ global String comment I_PARENT_START =     1
@ global String comment J_PARENT_START =     1
@ global String comment PARENT_GRID_RATIO =     1
@ global String comment DT =        48.00
@ global String comment ISWATER =    16
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14


By the way, this is for d01, which was set to run microphysics option 5. The metadata apparently take my d02 value.
The GrADS output suggests the control file is not matching the binary output correctly.
I appreciate your help.

Oscar

[wrfhelp]

What happens if you only compute a small number of fields, like tc,height,XLAT,XLONG.

The meta data comes directly from the wrfout file, so it cannot mix information from d02.

wrfhelp

[Oscar van der Velde]

For these parameters I get the following ctl file (for my nest and for splitted files this time, fewer levels, but doesn't matter for the results):
I have also tried omitting the PDEF line and putting the grid size directly in XDEF and YDEF, but does not fix it.

Code: Select all

dset ^arwpost_d02_%y4-%m2-%d2_%h2:%n2.dat
options  byteswapped template
undef 1.e30
title  OUTPUT FROM WRF V3.0 MODEL
pdef   87  75 lcc  41.443    1.913   44.000   38.000  42.09400  42.09400    1.04900   2666.667   2666.667
xdef  252 linear    0.40986   0.01201201
ydef  169 linear   40.41742   0.01201201
zdef    3 levels 
 900.00000
 700.00000
 500.00000
tdef    6 linear 13Z18APR2008      12MN     
VARS    4
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
height         3  0  Model height (KM)
tc             3  0  Temperature (C)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.0 MODEL
@ global String comment SIMULATION_START_DATE = 2008-04-18_06:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION =    88
@ global String comment SOUTH-NORTH_GRID_DIMENSION =    76
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    50
@ global String comment MAP_PROJ =     1
@ global String comment DX =      2666.67
@ global String comment DY =      2666.67
@ global String comment CEN_LAT =        41.44
@ global String comment CEN_LON =         1.91
@ global String comment TRUELAT1 =        42.09
@ global String comment TRUELAT2 =        42.09
@ global String comment MOAD_CEN_LAT =        42.09
@ global String comment STAND_LON =         1.05
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =    10
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     1
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     1
@ global String comment CU_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment UCMCALL =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment W_DAMPING =     0
@ global String comment PD_MOIST =     1
@ global String comment PD_SCALAR =     1
@ global String comment PD_TKE =     0
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment GRID_ID =     2
@ global String comment PARENT_ID =     1
@ global String comment I_PARENT_START =    40
@ global String comment J_PARENT_START =    17
@ global String comment PARENT_GRID_RATIO =     3
@ global String comment DT =        16.00
@ global String comment ISWATER =    16
@ global String comment ISICE =     0
@ global String comment ISURBAN =     0
@ global String comment ISOILWATER =     0


For this I get in GrADS:

ga-> open arwpost_d02
Scanning description file: arwpost_d02
Data file arwpost_d02_%y4-%m2-%d2_%h2:%n2.dat is open as file 1
LON set to 0.40986 3.42487
LAT set to 40.4174 42.4354
LEV set to 900 900
Time values set: 2008:4:18:13 2008:4:18:13
Notice: Implied interpolation for file arwpost_d02
Interpolation will be performed on any data displayed from this file
d xlat
(draws a correct plot of latitude)
d xlong
(Constant field. Value = 0)
d height
(Constant field. Value = 0)
d tc
(Constant field. Value = 0)

set lev 700
d xlat
ga-> d xlat
Notice: Automatic Grid Interpolation Taking Place
Data Request Warning: Request beyond file limits
Entire grid contents are set to missing data
Dimension ranges are: X = 1 252 Y = 0 169 Z = 2 1 T = 1 6
Warning issued for variable = xlat
Cannot contour grid - all undefined values

looking a bit further, at lev = 500 'd height' gives exactly the field of longitude.

After some testing it appears to be working for every first parameter in the file, but only for the lowest level (if lcl,tc,height are requested, height is the first one, gives only valid output at 900 hPa in my example, but for all output times).

Oscar

[Oscar van der Velde]

Thanks to wrfhelp it now works properly. Do this if you encounter the problem:

Could you please remove the recl from the configure.arwp file.
clean (NOT clean -a), and compile again.
This should solve the problem. Let me know if it does not.

wrfhelp

best wishes
Oscar

[cioran]

Hi!

We get the same
Cannot contour grid - all undefined values message,
removing the recl option in configure.arwp brought up an error in arwpost Fortran runtime error: Write exceeds length of DIRECT access record

Any ideas? Thank you!

[maldmm5]

How did you manage to fix this problem

[camilath]

Well, I made this, but I had the same problem of Cioran.. "Fortran runtime error: Write exceeds length of DIRECT access record"

Someone could resolve this?

Thanks!