Question about simulation

28 views
Skip to first unread message

Francis Oussou

unread,
Mar 26, 2025, 9:12:58 PMMar 26
to ParFlow
Dear all,
I believe you are doing well.

Can you please check my simulation outputs (see attached)? I want to know why the fluxes do not reflect the land use spatial variability. See the solver and CLM settings below.
I tried multiple times but could not solve the problem. Any suggestions?

Best regards,
Francis


#==
# CLM Settings:
# ------------------------------------------------------------
nalohou_run.Solver.LSM = 'CLM'
nalohou_run.Solver.CLM.CLMFileDir  =    "clm_output/"
nalohou_run.Solver.CLM.Print1dOut  = True
#nalohou_run.Solver.BinaryOutDir  =   True
#nalohou_run.Solver.WriteCLMBinary = True
nalohou_run.Solver.CLM.DailyRST = False
nalohou_run.Solver.CLM.SingleFile  = True
nalohou_run.Solver.CLM.CLMDumpInterval = 1
#nalohou_run.Solver.CLM.UseSlopeAspect = True
#nalohou_run.Solver.SlopeAccountingCLM = True

nalohou_run.Solver.CLM.MetForcing  = '3D'#'NC' # 3D
nalohou_run.Solver.CLM.MetFileName  = 'ERA5' #"NLDAS"
nalohou_run.Solver.CLM.MetFilePath  ="/bg/home/ousso-f/others/processing_@123/FORCING/nalohou_pfb" # already distributed
nalohou_run.Solver.CLM.MetFileNT  =  24
nalohou_run.Solver.CLM.IstepStart  = 1

nalohou_run.Solver.CLM.EvapBeta  =   'Linear'
nalohou_run.Solver.CLM.VegWaterStress  = 'Saturation'
nalohou_run.Solver.CLM.ResSat  =     0.1
nalohou_run.Solver.CLM.WiltingPoint  =  0.12
nalohou_run.Solver.CLM.FieldCapacity  = 0.98
nalohou_run.Solver.CLM.IrrigationType  = 'none'

nalohou_run.Solver.CLM.RootZoneNZ  = 4#10
nalohou_run.Solver.CLM.SoiLayer  =   4


nalohou_run.NetCDF.CLMNumStepsPerFile = 720#24
nalohou_run.NetCDF.WriteCLM = True

# Prof Reed lines
nalohou_run.Solver.PrintLSMSink        = False

#nalohou_run.Solver.CLM.CLMDumpInterval = 1
##nalohou_run.Solver.CLM.CLMFileDir      = 'output/'

nalohou_run.Solver.CLM.BinaryOutDir    = False

nalohou_run.Solver.CLM.IstepStart      = 1

nalohou_run.Solver.WriteCLMBinary      = False

nalohou_run.Solver.WriteSiloCLM        = False

nalohou_run.Solver.CLM.WriteLogs       = False

nalohou_run.Solver.CLM.WriteLastRST    = True

#nalohou_run.Solver.CLM.DailyRST        = False

#nalohou_run.Solver.CLM.SingleFile      = True


#-----------------------------------------------------------------------------
# Set solver parameters
#-----------------------------------------------------------------------------
#Run with only these solver settings, comment the rest out: Prof Reed lines

nalohou_run.Solver.TerrainFollowingGrid = True

nalohou_run.Solver = 'Richards'

nalohou_run.Solver.Nonlinear.MaxIter = 100#700

nalohou_run.Solver.Nonlinear.ResidualTol = 1e-5#1e-2

nalohou_run.Solver.Nonlinear.EtaChoice = 'EtaConstant'

nalohou_run.Solver.Nonlinear.EtaValue = 0.001

nalohou_run.Solver.Nonlinear.UseJacobian = True

nalohou_run.Solver.Linear.Preconditioner = 'PFMG'

#nalohou_run.Solver.PrintInitalConditions = False


latent_heat.png
evap_tot.png
sensible_heat.png
lclu.png
ground_heat.png
wttd.png

Reed M. Maxwell

unread,
Mar 26, 2025, 9:32:43 PMMar 26
to ParFlow

It looks like one of your inputs is mismatched.  That is, ParFlow is expecting something like NX, NY, NZ and you have either the order off e.g. NZ, NY, NX or something is off by a value e.g. NX-1, NY, NZ.

 

 

 

--
You received this message because you are subscribed to the Google Groups "ParFlow" group.
To unsubscribe from this group and stop receiving emails from it, send an email to parflow+u...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/parflow/CAE5G0iqA8QxUekgoTh7gXJXQTsS-ygOGsdCSV1ZFDBXF7WxbQA%40mail.gmail.com.

Reply all
Reply to author
Forward
0 new messages