Question about running 20th Century Reanalysis V3 data

[Zachary Esch]

Hello, I've been having issues with trying to get 20CRV3 input data to run in WRF v4.7.1. Since 20th century reanalysis V3 is missing specific humidity data for levels above 100 mb I'm using 1000-100 mb pressure levels. Yet when running real.exe I keep getting this error and I can't figure out what is causing it. I attached my namelist file.

*** -> Check your wrfinput_d01 file, you might have input pressure levels too close together (********* Pa and ********* Pa) at ( 38, 103) for variable type Z
*** -> Check your wrfinput_d01 file, you might have input pressure levels too close together (********* Pa and ********* Pa) at ( 39, 103) for variable type Z
*** -> Check your wrfinput_d01 file, you might have input pressure levels too close together (********* Pa and ********* Pa) at ( 40, 103) for variable type Z
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 6321
troubles, could not find trapping x locations
-------------------------------------------

i,j = 40 103
target pressure and value = 59.0746689 -4.83910707E+26
column of pressure and value = 59.0799675 39957.9453
column of pressure and value = 59.0799675 40177.0977
column of pressure and value = 59.0799675 40401.1992
column of pressure and value = 59.0799675 40630.3242
column of pressure and value = 59.0871162 6.36811194E+26
column of pressure and value = 59.0871162 3.29212637E+27
column of pressure and value = 59.0871162 6.10528790E+27
column of pressure and value = 59.0871162 9.09714672E+27
column of pressure and value = 59.0912895 1.06955698E+28
column of pressure and value = 59.0912895 1.06955698E+28
column of pressure and value = 59.0912895 1.06955698E+28
column of pressure and value = 59.0912895 1.06955698E+28
column of pressure and value = 59.0799675 0.00000000
column of pressure and value = NaN -1.38638918E+30
column of pressure and value = NaN -4.47504655E+30
column of pressure and value = NaN -6.19864746E+30
column of pressure and value = NaN -6.19247105E+30
column of pressure and value = 5.70378256 -6.18890593E+30
column of pressure and value = 5.52146101 -6.18890593E+30
column of pressure and value = 5.29831743 -6.18890593E+30
column of pressure and value = 5.01063538 -6.18890593E+30
column of pressure and value = 4.60517025 -6.18890593E+30

 

[Ming Chen]

If you set p_top_requested = 15000, can you run real.exe successfully?

 

[Zachary Esch]

15000

I tried that and got the exact same error.

 

[Ming Chen]

In this case, I am suspicous that your input data is either not correct or sufficient. Please check whether all the variables below are included in your intermediate files:

Field name in intermediate file	Units​	Description​	Notes​
TT	K	3-d air temperature
RH	%	3-d relative humidity	Not needed if SPECHUMD is available
SPECHUMD	kg kg-1	3-d specific humidity	Not needed if RH is available
UU	m s-1	3-d wind u-component
VV	m s-1	3-d wind v-component
GHT	m	3-d geopotential height
PRESSURE	Pa	3-d pressure	Only needed for non-isobaric datasets
PSFC	Pa	Surface pressure
PMSL	Pa	Mean sea-level pressure
SKINTEMP	K	Skin temperature
SOILHGT	m	Soil height
TT	K	2-meter air temperature
RH	%	2-meter relative humidity	Not needed if SPECHUMD is available
SPECHUMD	kg kg-1	2-meter specific humidity	Not needed if RH is available
UU	m s-1	2-meter wind u-component
VV	m s-1	2-meter wind v-component
LANDSEA	fraction	Land-sea mask (0=water, 1=land)
SMtttbbb	m3 m-3	Soil moisture	'ttt' is the layer top depth in cm, and 'bbb' is the layer bottom depth in cm
STtttbbb	K	Soil temperature
SOILMmmm	kg m-3	Soil moisture	'mmm' is the level depth in cm, not needed if SMtttbbb available
SOILTmmm	K	Soil temperature

 

[Zachary Esch]

In this case, I am suspicous that your input data is either not correct or sufficient. Please check whether all the variables below are included in your intermediate files:

Field name in intermediate file	Units​	Description​	Notes​
TT	K	3-d air temperature
RH	%	3-d relative humidity	Not needed if SPECHUMD is available
SPECHUMD	kg kg-1	3-d specific humidity	Not needed if RH is available
UU	m s-1	3-d wind u-component
VV	m s-1	3-d wind v-component
GHT	m	3-d geopotential height
PRESSURE	Pa	3-d pressure	Only needed for non-isobaric datasets
PSFC	Pa	Surface pressure
PMSL	Pa	Mean sea-level pressure
SKINTEMP	K	Skin temperature
SOILHGT	m	Soil height
TT	K	2-meter air temperature
RH	%	2-meter relative humidity	Not needed if SPECHUMD is available
SPECHUMD	kg kg-1	2-meter specific humidity	Not needed if RH is available
UU	m s-1	2-meter wind u-component
VV	m s-1	2-meter wind v-component
LANDSEA	fraction	Land-sea mask (0=water, 1=land)
SMtttbbb	m3 m-3	Soil moisture	'ttt' is the layer top depth in cm, and 'bbb' is the layer bottom depth in cm
STtttbbb	K	Soil temperature
SOILMmmm	kg m-3	Soil moisture	'mmm' is the level depth in cm, not needed if SMtttbbb available
SOILTmmm	K	Soil temperature

Hi, I figured it out the error and managed to fix it. It was the way I was creating the WRF intermediate files that was the issue.

 

[Ming Chen]

Thank you for the update. Would you please clarify what issue you had previously? I suppose such kind of information will be helpful for other users who may also want to use 20CRV3 as input to WRF. Thanks in advance.