Open LluisFB opened 1 day ago
Hi Lluís,
where did you get the required units from? in the FPS (and in CORDEX in general), they seem to be kg m-2
The REFSMC parameter leads to a single soil column value, depending on the soil type.
But, then, there are places in the code where the field capacity (read into SMCREF) is considered by soil layer, as a 3D variable: https://github.com/NCAR/noahmp/blob/c232448ef191de0ab0357ebd249636f504bb7a24/src/module_sf_noahmplsm.F#L8089
Maybe a NoahMP developer can help us here to understand where we can get the field capacity actually used by the model @cenlinhe @barlage
The soil properties in the parameters structure have dimension NSOIL, regardless of iopt_soil setting
Depending on the iopt_soil setting, it will either use dominant for all layers (option 1) or use classes that can be different in each layer (option 2)
I believe the default in WRF is to use option 1
Thank you for the swift reply!
we are using the default opt_soil = 1, yes (one never ends learning about parameter settings in WRF) https://github.com/CORDEX-WRF-community/WRF/blob/v4.5.1.5/run/README.namelist#L1346
@LluisFB, then I think you are OK with your computation, except for the units. I'm no expert in soil processes, but I'd say that to convert from volumetric fraction:
mrsofc(i,j) = REFSMC(i,j) * tot_dz * 1000. [kg m-2]
Also, not sure why this variable was requested in the FPS, since soil moisture itself is not requested.
Regarding the CMORization of the WRF output, I created a python/Fortran based post-processor: python-cordex.
Accordingly to the requested CORDEX variables, we need to provide
mrsofc
being `Capacity of Soil to Store Water (Field Capacity).WRF does not provide this variable, thus, what I implemented is the following:
ISLTYP
SOILPARM.TBL
ISLTYP
with the correspondent entries inSOILPARM.TBL
REFSMC
(whicih accroding to this NoahMP- web, corresponds to Reference soil moisture (field capacity), where transpiration begins to stress [volumetric fraction])m3m-2
by integrating to all soil depth (multiplying by totDZ) mrsofc(i,j) = REFSMC(i,j)*tot_dzPlease if anyone from the soil/land community can revise it, will be highly appreciated
Hope it make sense.
Lluís