Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

What does the soil layer variable levgrnd_class do?

A

amans

Aman Shrestha
New Member
What version of the code are you using?
ctsm5.2.023

Have you made any changes to files in the source tree?
Not yet. But looking to increase soil layer structure (either layer thickness or number of layers).

Describe your problem or question:
In src/main/initVerticalMod.F90, soil layers in each column are provided one of three classes (STANDARD, DEEP_BEDROCK, or SHALLOW_BEDROCK) in levgrnd_class. They differentiate the layers between soil and bedrock. How do they exactly impact the hydrologically active soil processes? I don't see how it is being used. Grepping for levgrnd_class in the src/ dir points to ColumnType.F90 and subgridRestMod.F90 in addition to initVerticalMod.F90, but they don't answer my question.
Apologies if I missed something obvious.
 
Sort by date Sort by votes
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I think that levgrnd_class is only used for the interpolation of initial conditions code as noted by the comment in initVerticalMod.F90:

! Possible values for levgrnd_class. The important thing is that, for a given column,
! layers that are fundamentally different (e.g., soil vs bedrock) have different
! values. This information is used in the vertical interpolation in init_interp.

See LEVGRND_CLASS in init_interp/initInterpMultilevelContainer.F90
So I don't think this is a variable used for, e.g., hydrologically active soil processes.
 
Vote Upvote 0 Downvote
You must log in or register to reply here.
Top