Skip to content

Default soiltype namelist option #2166

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
weiwangncar merged 2 commits into from
Apr 1, 2025
Merged

Default soiltype namelist option #2166

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
1025259
FEAT: Add a default_soiltype namelist option
mathieulandreau Feb 11, 2025
b6f0053
STYLE: add default_soiltype to the README.namelist
mathieulandreau Feb 11, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2502,6 +2502,7 @@ rconfig integer ideal_xland namelist,physics 1 1
rconfig real swrad_scat namelist,physics 1 1 irh "SWRAD_SCAT" "SCATTERING FACTOR IN SWRAD" ""
rconfig integer surface_input_source namelist,physics 1 3 irh "surface_input_source" "1=static (fractional), 2=time dependent (dominant), 3=dominant cateogry from metgrid" ""
rconfig integer num_soil_layers namelist,physics 1 5 irh "num_soil_layers" "" ""
rconfig integer default_soiltype namelist,physics 1 8 rh "default_soiltype" "Soil category used to correct missing values (default 8 : silty clay loam)" ""
rconfig integer num_pft_clm namelist,physics 1 17 - "num_pft_clm" "" ""
rconfig logical input_pft namelist,physics 1 .false. h "input_pft_flag" "use input pft instead of USGS" ""
rconfig integer maxpatch namelist,physics 1 10 irh "maxpatch" "" ""
Expand Down
17 changes: 12 additions & 5 deletions dyn_em/module_initialize_real.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3114,7 +3114,7 @@ SUBROUTINE init_domain_rk ( grid &
DO i = its, MIN(ide-1,ite)
IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE
IF ( grid%landmask(i,j) .GT. 0.5 .AND. grid%isltyp(i,j) .EQ. grid%isoilwater ) THEN
grid%isltyp(i,j) = 8
grid%isltyp(i,j) = config_flags%default_soiltype
change_soilw = change_soilw + 1
iforce = iforce + 1
ELSE IF ( grid%landmask(i,j) .LT. 0.5 .AND. grid%isltyp(i,j) .NE. grid%isoilwater ) THEN
Expand All @@ -3125,10 +3125,17 @@ SUBROUTINE init_domain_rk ( grid &
END DO
END DO
IF ( change_soilw .GT. 0 .OR. change_soil .GT. 0 ) THEN
WRITE(a_message,FMT='(A,I4,A,I6)' ) &
'forcing artificial silty clay loam at ',iforce,' points, out of ',&
(MIN(ide-1,ite)-its+1)*(MIN(jde-1,jte)-jts+1)
CALL wrf_debug(0,a_message)
IF (config_flags%default_soiltype == 8) THEN
WRITE(a_message,FMT='(A,I4,A,I6)' ) &
'forcing artificial silty clay loam at ',iforce,' points, out of ',&
(MIN(ide-1,ite)-its+1)*(MIN(jde-1,jte)-jts+1)
ELSE
WRITE(a_message,FMT='(A,I4,A,I4,A,I6)' ) &
'forcing soil type to the user defined default category (',config_flags%default_soiltype,&
') at ',iforce,' points, out of ',&
(MIN(ide-1,ite)-its+1)*(MIN(jde-1,jte)-jts+1)
ENDIF
CALL wrf_debug(0,a_message)
END IF

END IF
Expand Down
1 change: 1 addition & 0 deletions run/README.namelist
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -941,6 +941,7 @@ Namelist variables for controlling the adaptive time step option:
= 10: CLM4 landsurface model
= 2: Pleim-Xu landsurface model
= 3: SSiB landsurface model
default_soiltype = 8, ! soil category applied in land cells where no soil category is defined, default is 8 = silty clay loam
num_land_cat = 21, ! number of land categories in input data.
24 - for USGS (default); 20 for MODIS
28 - for USGS if including lake category
Expand Down
13 changes: 13 additions & 0 deletions share/module_check_a_mundo.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -486,6 +486,19 @@ END FUNCTION bep_bem_ngr_u
END IF
ENDDO

!-----------------------------------------------------------------------
! Check that default_soiltype is between 1 and 19 and is not water (14)
!-----------------------------------------------------------------------
IF ((model_config_rec % default_soiltype .LT. 1) &
.OR. (model_config_rec % default_soiltype .GT. 19) &
.OR. (model_config_rec % default_soiltype .EQ. 14)) THEN
WRITE (wrf_err_message, FMT='(A,I4,A)') '--- ERROR: invalid default_soiltype category (', model_config_rec % default_soiltype, ')'
CALL wrf_message ( wrf_err_message )
WRITE (wrf_err_message, FMT='(A)') '--- ERROR: resetting default_soiltype to 8 (silty clay loam)'
CALL wrf_message ( wrf_err_message )
model_config_rec % default_soiltype = 8
END IF

!-----------------------------------------------------------------------
! Check that mosiac option cannot turn on when sf_urban_physics = 2 and 3
!-----------------------------------------------------------------------
Expand Down