I believe the LCZ category is an issue here. This is because the urban local climate zone (LCZ) numbers of 31-41 overlap with NLCD land categories, which may cause errors. For more details, please see
The Unix operating system usually looks at the environment variable LD_LIBRARY_PATH to search the directories in the path for the library. Try issuing these commands:
If the path where your libpnetcdf.so.4 is located is not included, then (in ksh, for example) please...
No, 'h' could be anywhere in the 8th column for a state variable.
By the way, which radiation scheme are you using? Note that not all radiation schemes calculate SWDDNI, SWDDIF, SWDDIR etc.
CAM, GODDARD, RRTMG are the schemes that can output the above variables.
By saying "point-to-point same", do you mean the wrfout files from the two runs are identical?
Can you check whether you have the file 'diffwrf' existing in in WRF/external/io_netcdf? if so, please issue the command:
WRF/external/io_netcdf/diffwrf wrfout_urban wrfout_no_urban > log
I guess you are running WRFV4.5.1 using only one processors. Please let me know if I am wrong. I ask tis because I found that
"WRF NUMBER OF TILES = 1" in your wrf451_1.0000 log file.
This is related to other issues:
(1) How did you compile WRFv4.4.2 and WRFv4.5.1?
(2) How did you run them...
Various CAMtr_volume_mixing data for different scenarios can be found in WRF/run directory.
When you run WRF under one specific scenario, you need to link the corresponding CAMtr data as "CAMtr_volume_mixing_ratio" your working directory .
For example, suppose you are running wrf..exe...
I would suggest you remove HGT from geo_em.nc files, then rerun metgrid.exe and real.exe. This is because HGT is used for vertical interpolation in REAL program. It affects both wrfinput and wrfbdy files. If you simply remove it from wrfinput, its impact still exists in wrfbdy.
I looked at your namelist.inout. Both the physics and dynamic options look appropriate.
The only concern I have is the resolution of your case. Note that cumulus scheme works fine for grid interval larger than 10km, and it should be turned off when grid interval is smaller than...
For the equation to calculate RR, we tracked it down to an OSU 1d model by Ek and Mahrt (1991).
RGAS=287 and CP=1004 works out about right. Factor becomes 8/7.
Thereby the RR calculation is correct.
SFCEVP is ACCUMULATED SURFACE EVAPORATION with positive values upward. This means that when there is evaporation from the surface to the upper air, the value is positive. If condensation occurs, then its value is negative.
I looked at the files rsl422.out and wf451. It seems that wf451.0000 is output from an initial run, while rsl422 is from a restart run.
Please let me if I am wrong.
Can you run the same case with the same settings and see whether there is large difference in the computation between them?
I have talked to our physics expert and Noah LSM expert. Unfortunately we haven't figured out yet the physical meaning of the constant " 6.48E-8" in the equation. I am sorry that I don't have an immediate answer to this question.
We will continue to pay attention to this issue and keep you...