Soil moistures: optimized parameters

[sadiyanti]

Hi,
I have 4 layers of soils, as such 4 sets of soil moisture parameters are defined in mhm_parameter.nml:, i.e. &soilmoisture1, ... &soilmoisture4, with all attributes.
But why in the optimised FinalParameter.nml, only 1 layer is outputted?
Thanks.

[silentin]

Hello,
sorry for the late reply. I hope you could already solve your issue but if it is still unsolved I would need a bit more information about your configuration?

1. How did you specify the number of soil horizons in the &soildata section of the mhm.nml?
2. What Input data did you provide?
3. The lines from the mhm_parameter.nml and the FinalParameter.nml files.

[rohinikumar82]

You would need to change the settings of mhm_output.nml to print fluxes and states of each soil layer.

[sadiyanti]

hi Rohini, thanks.
But I can't see how to change them to represent each layer.
For example, in mHM_outputs.nml for States, I have the following:

=====================
&NLoutputResults
output_deflate_level = 6
output_double_precision = .TRUE.
timeStep_model_outputs = -1
!
!----------------
! 1. states
!----------------
! interceptional storage (L1_inter) [mm] -- case 1
outputFlxState(1)=.TRUE.

How can I change for L2_inter?
should i do the following
outputFlxState(1,1) = .TRUE. !for L1
outputFlxState(2,1) = .TRUE. !for L2

etc?

[rohinikumar82]

in mhm_outpur.nml file -- you need to make.TRUE. for outputs that print for soil layer specific outputs -- for example -- soil moisture, ET, or infiltration fluxes. Not interception, it's not a variable related to the soil layering process.

[sadiyanti]

Hi Rohini,
Of course i made it .TRUE. for the following:
The problem is that in the created FinalParam.nml only has an output for first layer (L1), whereas i have 4 Layers (L1, L2, L3, L4).

! soil water content in the single layers (L1_soilMoist) -- case 3
outputFlxState(3)=.TRUE.
!
! volumetric soil moisture in the single layers (L1_soilMoist / L1_soilMoistSat ) [mm/mm] -- case 4
outputFlxState(4)=.TRUE.
!
! mean volumetric soil moisture averaged over all soil layers (L1_soilMoist / L1_soilMoistSat ) [mm/mm] -- case 5
outputFlxState(5)=.TRUE.
!
! waterdepth in reservoir of sealed areas (L1_sealSTW) [mm] -- case 6
outputFlxState(6)=.TRUE.
!
! waterdepth in reservoir of unsat. soil zone (L1_unsatSTW) [mm] -- case 7
outputFlxState(7)=.TRUE.
!
! waterdepth in reservoir of sat. soil zone (L1_satSTW) [mm] -- case 8
! --> level of GW reservoir
outputFlxState(8)=.TRUE.

I read in mHM "known bugs": 3. If ProcessOption(3) is set to 3 and optimization is activated, the created FinalParam.nml misses the header for the namelist of the soil moisture parameters" Is this referring to the above issue?

[rohinikumar82]

FinalParam.nml has nothing to do with soil layers fluxes. It's only has the specification of global parameters valid for all soil layers.

[sadiyanti]

OFFICIAL Hi Simon, please see my response in blue. Looking forward to hear from you. Is there anything I need to change to make the outputs for each soil layers? Thanks, Dr Sri Adiyanti Senior Research Scientist – Hydrodynamic Biogeochemical Modeller | Rivers and Estuaries Science Department of Biodiversity, Conservation and Attractions 17 Dick Perry Avenue, Kensington WA 6151, AUSTRALIA P: +61 8 9278 0953 E: ***@***.******@***.***>

…

________________________________ From: Simon Lüdke ***@***.***> Sent: Wednesday, June 25, 2025 11:36 PM To: mhm-ufz/mHM ***@***.***> Cc: Sri Adiyanti ***@***.***>; Author ***@***.***> Subject: Re: [mhm-ufz/mHM] Soil moistures: optimized parameters (Discussion #124) [External Email] This email was sent from outside the department – be cautious, particularly with links and attachments. Hello, sorry for the late reply. I hope you could already solve your issue but if it is still unsolved I would need a bit more information about your configuration? 1. How did you specify the number of soil horizons in the &soildata section of the mhm.nml? iFlag_soilDB = 0 ! !> [mm] soil depth down to which organic matter is possible tillageDepth = 500 !previously 200 ! !> No. of soil horizons to be modeled nSoilHorizons_mHM = 3 ! ! soil depth information !> IF (iFlag_soilDB = 0) !> Provide below the soil_Depth() for 1,2,..,*n-1* soil horizons. Depth of the last layer(n) is taken from the soil LUT !> IF (iFlag_soilDB = 1) !> Provide below the soil_Depth() for every 1,2..n-1,*n* soil horizons. You must have soil_class-id gridded file for each layer !> Also check your tillage layer depth. It should conform with one of the below specified soil_Depth. ! !> Soil_Horizon Depth[mm] ! bottom depth of soil horizons w.r.t. ground surface (positive downwards) soil_Depth(1) = 200 soil_Depth(2) = 950 2. What Input data did you provide? Inputs: soil_class.asc and soil_classdefinition 3. The lines from the mhm_parameter.nml and the FinalParameter.nml files. I have differentiate &soilmoisture1, &soilmoisture2, &soilmoisture3 in mhm_parameter.nml, but in the FinalParameter.nml there is only 1 resulted: &soilmoisture1. In mhm_outputs.nml I set .TRUE. for outputFlState(3) and outputState(4). — Reply to this email directly, view it on GitHub<#124 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AXSZW3FZP4EZWTBWVDSFSI33FK6X3AVCNFSM6AAAAAB3VRUNCOVHI2DSMVQWIX3LMV43URDJONRXK43TNFXW4Q3PNVWWK3TUHMYTGNJXGUZDGMA>. You are receiving this because you authored the thread.

________________________________ This message is confidential and is intended for the recipient named above. If you are not the intended recipient, you must not disclose, use or copy the message or any part of it. If you received this message in error, please notify the sender immediately by replying to this message, then delete it from your system.