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