How to convert the wrfinput_d01 file into the HRLDAS_setup_2000010100_d1 format?

[OBJ780]

Dear Professor,

I am planning to use the Noah-MP model for regional simulations. I have created the wrfinput_d01 file using WPS and WRF, which contains the necessary initial state variables for the Noah-MP regional simulation at each grid point, including latitude, longitude, elevation, soil type, vegetation type, soil moisture at different depths, and soil temperature. Therefore, this file provides the required initial conditions for the regional simulation.

The current issue is how I can convert the wrfinput_d01 file into a file similar to HRLDAS_setup_2000010100_d1. Since I did not follow the standard tutorial to store the different initial variables in separate directories, I cannot use the create_forcing.exe script to generate the initial field files. As for the atmospheric forcing data, I have already prepared it using Python code.

Could you please advise on how to convert the wrfinput_d01 file into the HRLDAS_setup_2000010100_d1 format? I would greatly appreciate your guidance!

Best regards.

[cenlinhe]

You can directly use wrfinput file as the HRLDAS_setup file. Just put your wrfinput file path in the namelist.hrldas for the HRLDAS_SETUP_FILE entry.

[OBJ780]

I have encountered some difficulties:./hrldas.exe

***** ERROR: Problem reading namelist NOAHLSM_OFFLINE

At line 183 of file NoahmpReadNamelistMod.f90 (unit = 30, file = 'namelist.hrldas')
Fortran runtime error: Cannot match namelist object name .5
Here is my namelist.hrldas:
&NOAHLSM_OFFLINE

HRLDAS_SETUP_FILE = "/home/xc/IPYNB/TP_FORCING/wrfinput_d01"
INDIR = "/home/xc/IPYNB/TP_FORCING/CMFD_Forcing/"
OUTDIR = "/home/xc/IPYNB/TP_FORCING/LDASOUT_DOMAIN1/"

START_YEAR = 2016
START_MONTH = 01
START_DAY = 01
START_HOUR = 00
START_MIN = 00

! RESTART_FILENAME_REQUESTED = "./RESTART.2010010200_DOMAIN1"

KDAY = 3
SPINUP_LOOPS = 0

FORCING_NAME_T = "T2D"
FORCING_NAME_Q = "Q2D"
FORCING_NAME_U = "U2D"
FORCING_NAME_V = "V2D"
FORCING_NAME_P = "PSFC"
FORCING_NAME_LW = "LWDOWN"
FORCING_NAME_SW = "SWDOWN"
FORCING_NAME_PR = "RAINRATE"

DYNAMIC_VEG_OPTION = 4
CANOPY_STOMATAL_RESISTANCE_OPTION = 1
BTR_OPTION = 1
SURFACE_RUNOFF_OPTION = 3
SUBSURFACE_RUNOFF_OPTION = 3
DVIC_INFILTRATION_OPTION = 1
SURFACE_DRAG_OPTION = 1
FROZEN_SOIL_OPTION = 1
SUPERCOOLED_WATER_OPTION = 1
RADIATIVE_TRANSFER_OPTION = 3
SNOW_ALBEDO_OPTION = 2
PCP_PARTITION_OPTION = 1
SNOW_THERMAL_CONDUCTIVITY = 1
TBOT_OPTION = 2
TEMP_TIME_SCHEME_OPTION = 1
GLACIER_OPTION = 1
SURFACE_RESISTANCE_OPTION = 1
SOIL_DATA_OPTION = 1
PEDOTRANSFER_OPTION = 1
CROP_OPTION = 0
IRRIGATION_OPTION = 0
IRRIGATION_METHOD = 0
TILE_DRAINAGE_OPTION = 0

FORCING_TIMESTEP = 3600
NOAH_TIMESTEP = 3600
OUTPUT_TIMESTEP = 3600

SPLIT_OUTPUT_COUNT = 1
SKIP_FIRST_OUTPUT = .true.
RESTART_FREQUENCY_HOURS = 0

XSTART = 72.5
XEND = 109.5
YSTART = 25
YEND = 42

NSOIL=4
soil_thick_input(1) = 0.10
soil_thick_input(2) = 0.30
soil_thick_input(3) = 0.60
soil_thick_input(4) = 1.00

ZLVL = 10.0

SF_URBAN_PHYSICS = 0
USE_WUDAPT_LCZ = 0

/
Why is this?

[cenlinhe]

XSTART, XEND, YSTART, YEND should be the grid index (starting from 1; integer) of your subdomain, not actual lat/lon values.

[OBJ780]

That's ok! Here is new questions:wrfinput_flnm: '/home/xc/IPYNB/TP_FORCING/wrfinput_d01'
wrfinput_flnm: '/home/xc/IPYNB/TP_FORCING/wrfinput_d01'
netcdf_flnm: '/home/xc/IPYNB/TP_FORCING/wrfinput_d01'
LDASIN TITLE attribute: OUTPUT FROM REAL_EM V4.1.2 PREPROCESSOR
WARNING: LDASIN file has a perverse version identifier
ldasin_version = 0
LDASIN MMNINLU attribute: MODIFIED_IGBP_MODIS_NOAH
layer_bottom(1:nsoil) = 0.1000 0.4000 1.0000 2.0000
layer_top(1:nsoil) = 0.0000 0.1000 0.4000 1.0000
Soil depth = 0.100000 0.300000 0.600000 1.000000
k, dst_centerpoint(k) = 1 5.00000007E-02
k, dst_centerpoint(k) = 2 0.250000000
k, dst_centerpoint(k) = 3 0.700000048
k, dst_centerpoint(k) = 4 1.50000000

k, src_centerpoint(k) = 1 5.00000007E-02
k, src_centerpoint(k) = 2 0.250000000
k, src_centerpoint(k) = 3 0.699999988
k, src_centerpoint(k) = 4 1.50000000
(Near) match for destination layer: Taking destination layer at 0.0500 from source layer at 0.0500
(Near) match for destination layer: Taking destination layer at 0.2500 from source layer at 0.2500
(Near) match for destination layer: Taking destination layer at 0.7000 from source layer at 0.7000
(Near) match for destination layer: Taking destination layer at 1.5000 from source layer at 1.5000
k, dst_centerpoint(k) = 1 5.00000007E-02
k, dst_centerpoint(k) = 2 0.250000000
k, dst_centerpoint(k) = 3 0.700000048
k, dst_centerpoint(k) = 4 1.50000000

k, src_centerpoint(k) = 1 5.00000007E-02
k, src_centerpoint(k) = 2 0.250000000
k, src_centerpoint(k) = 3 0.699999988
k, src_centerpoint(k) = 4 1.50000000
(Near) match for destination layer: Taking destination layer at 0.0500 from source layer at 0.0500
(Near) match for destination layer: Taking destination layer at 0.2500 from source layer at 0.2500
(Near) match for destination layer: Taking destination layer at 0.7000 from source layer at 0.7000
(Near) match for destination layer: Taking destination layer at 1.5000 from source layer at 1.5000
SNOW HEIGHT NOT FOUND - VALUE DEFINED IN LSMINIT
NTIME = 72 KHOUR= 72 dtbl = 3600.00000
A) Problem opening netcdf file: '/home/xc/IPYNB/TP_FORCING/CMFD_Forcing//2016010100.LDASIN_DOMAIN1'
Here is my forcing data format:
Variable Names and Units:
Variable Name Unit
0 time No unit
1 RAINRATE mm/s
2 U2D m s-1
3 T2D Kelvin
4 Q2D kg kg-1
5 PSFC pa
6 SWDOWN W m-2
7 LWDOWN W m-2
8 V2D m s-1
9 XLAT degrees latitude
10 XLONG degrees longitude
11 south_north No unit
12 west_east No unit
and:
<xarray.Dataset>
Dimensions: (time: 1, south_north: 340, west_east: 740)
Coordinates:

• time (time) datetime64[ns] 2018-01-01
  XLAT (south_north, west_east) float32 ...
  XLONG (south_north, west_east) float32 ...
• south_north (south_north) int64 0 1 2 3 4 5 6 ... 334 335 336 337 338 339
• west_east (west_east) int64 0 1 2 3 4 5 6 ... 733 734 735 736 737 738 739
  Data variables:
  RAINRATE (time, south_north, west_east) float32 ...
  U2D (time, south_north, west_east) float32 ...
  T2D (time, south_north, west_east) float32 ...
  Q2D (time, south_north, west_east) float32 ...
  PSFC (time, south_north, west_east) float32 ...
  SWDOWN (time, south_north, west_east) float32 ...
  LWDOWN (time, south_north, west_east) float32 ...
  V2D (time, south_north, west_east) float32 ...
  Attributes:
  title: atmospheric forcing data retrival from itpcas dataset
  history: Tue Oct 24 23:37:10 2023: ncrcat -d time,0 /home/public_fol...
  Conventions: C0ARDS
  NCO: netCDF Operators version 5.1.5 (Homepage = http://nco.sf.ne.../
  Why model can't open 2016010100.LDASIN_DOMAIN1?

[cenlinhe]

did you check if you provided the correct folder path for your forcing files?

[cenlinhe]

are you able to use ncview to open your forcing file to view those variables successfully?

[OBJ780]

When performing single-point simulations previously, there was no issue using the 30-minute atmospheric forcing data. However, when using the 1-hour forcing data, the same error occurs: A) Problem opening netcdf file.

[cenlinhe]

this is interesting. We use 1-hour forcing data a lot and did not see any problems.