Example of Parflow CLM tcl file - Wei64380/Parflow-CLM GitHub Wiki

An example of the Parflow-CLM file for the Washita catchment

  1. SCRIPT TO RUN LITTLE WASHITA DOMAIN WITH TERRAIN-FOLLOWING GRID

  2. DETAILS:

  3. Arugments are 1) runname 2) year

  4. Import the ParFlow TCL package

  5. lappend auto_path $env(PARFLOW_DIR)/bin

  6. package require parflow

  7. namespace import Parflow::*

  8. pfset FileVersion 4

  9. #-----------------------------------------------------------------------------

  10. Set Processor topology

  11. #-----------------------------------------------------------------------------

  12. pfset Process.Topology.P 2

  13. pfset Process.Topology.Q 1

  14. pfset Process.Topology.R 1

  15. #-----------------------------------------------------------------------------

  16. Computational Grid

  17. #-----------------------------------------------------------------------------

  18. pfset ComputationalGrid.Lower.X 0.0

  19. pfset ComputationalGrid.Lower.Y 0.0

  20. pfset ComputationalGrid.Lower.Z 0.0

  21. pfset ComputationalGrid.DX 1000.0

  22. pfset ComputationalGrid.DY 1000.0

  23. pfset ComputationalGrid.DZ 2.0

  24. pfset ComputationalGrid.NX 41

  25. pfset ComputationalGrid.NY 41

  26. pfset ComputationalGrid.NZ 50

  27. #-----------------------------------------------------------------------------

  28. Names of the GeomInputs

  29. #-----------------------------------------------------------------------------

  30. pfset GeomInput.Names "box_input indi_input"

  31. #-----------------------------------------------------------------------------

  32. Domain Geometry Input

  33. #-----------------------------------------------------------------------------

  34. pfset GeomInput.box_input.InputType Box

  35. pfset GeomInput.box_input.GeomName domain

  36. #-----------------------------------------------------------------------------

  37. Domain Geometry

  38. #-----------------------------------------------------------------------------

  39. pfset Geom.domain.Lower.X 0.0

  40. pfset Geom.domain.Lower.Y 0.0

  41. pfset Geom.domain.Lower.Z 0.0

  43. pfset Geom.domain.Upper.X 41000.0

  44. pfset Geom.domain.Upper.Y 41000.0

  45. pfset Geom.domain.Upper.Z 100.0

  46. pfset Geom.domain.Patches "x-lower x-upper y-lower y-upper z-lower z-upper"

  47. #-----------------------------------------------------------------------------

  48. Indicator Geometry Input

  49. #-----------------------------------------------------------------------------

  50. pfset GeomInput.indi_input.InputType IndicatorField

  51. pfset GeomInput.indi_input.GeomNames "s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12 s13 g1 g2 g3 g4 g5 g6 g7 g8"

  52. pfset Geom.indi_input.FileName "IndicatorFile_Gleeson.50z.pfb"

  53. pfset GeomInput.s1.Value 1

  54. pfset GeomInput.s2.Value 2

  55. pfset GeomInput.s3.Value 3

  56. pfset GeomInput.s4.Value 4

  57. pfset GeomInput.s5.Value 5

  58. pfset GeomInput.s6.Value 6

  59. pfset GeomInput.s7.Value 7

  60. pfset GeomInput.s8.Value 8

  61. pfset GeomInput.s9.Value 9

  62. pfset GeomInput.s10.Value 10

  63. pfset GeomInput.s11.Value 11

  64. pfset GeomInput.s12.Value 12

  65. pfset GeomInput.s13.Value 13

  66. pfset GeomInput.g1.Value 21

  67. pfset GeomInput.g2.Value 22

  68. pfset GeomInput.g3.Value 23

  69. pfset GeomInput.g4.Value 24

  70. pfset GeomInput.g5.Value 25

  71. pfset GeomInput.g6.Value 26

  72. pfset GeomInput.g7.Value 27

  73. pfset GeomInput.g8.Value 28

  74. #-----------------------------------------------------------------------------

  75. Permeability (values in m/hr)

  76. #-----------------------------------------------------------------------------

  77. pfset Geom.Perm.Names "domain s1 s2 s3 s4 s5 s6 s7 s8 s9 g2 g3 g6 g8"

  78. pfset Geom.domain.Perm.Type Constant

  79. pfset Geom.domain.Perm.Value 0.2

  80. pfset Geom.s1.Perm.Type Constant

  81. pfset Geom.s1.Perm.Value 0.269022595

  82. pfset Geom.s2.Perm.Type Constant

  83. pfset Geom.s2.Perm.Value 0.043630356

  84. pfset Geom.s3.Perm.Type Constant

  85. pfset Geom.s3.Perm.Value 0.015841225

  86. pfset Geom.s4.Perm.Type Constant

  87. pfset Geom.s4.Perm.Value 0.007582087

  88. pfset Geom.s5.Perm.Type Constant

  89. pfset Geom.s5.Perm.Value 0.01818816

  90. pfset Geom.s6.Perm.Type Constant

  91. pfset Geom.s6.Perm.Value 0.005009435

  92. pfset Geom.s7.Perm.Type Constant

  93. pfset Geom.s7.Perm.Value 0.005492736

  94. pfset Geom.s8.Perm.Type Constant

  95. pfset Geom.s8.Perm.Value 0.004675077

  96. pfset Geom.s9.Perm.Type Constant

  97. pfset Geom.s9.Perm.Value 0.003386794

  98. pfset Geom.g2.Perm.Type Constant

  99. pfset Geom.g2.Perm.Value 0.025

  100. pfset Geom.g3.Perm.Type Constant

  101. pfset Geom.g3.Perm.Value 0.059

  102. pfset Geom.g6.Perm.Type Constant

  103. pfset Geom.g6.Perm.Value 0.2

  104. pfset Geom.g8.Perm.Type Constant

  105. pfset Geom.g8.Perm.Value 0.68

  106. pfset Perm.TensorType TensorByGeom

  107. pfset Geom.Perm.TensorByGeom.Names "domain"

  108. pfset Geom.domain.Perm.TensorValX 1.0d0

  109. pfset Geom.domain.Perm.TensorValY 1.0d0

  110. pfset Geom.domain.Perm.TensorValZ 1.0d0

  111. #-----------------------------------------------------------------------------

  112. Specific Storage

  113. #-----------------------------------------------------------------------------

  114. pfset SpecificStorage.Type Constant

  115. pfset SpecificStorage.GeomNames "domain"

  116. pfset Geom.domain.SpecificStorage.Value 1.0e-5

  117. #-----------------------------------------------------------------------------

  118. Phases

  119. #-----------------------------------------------------------------------------

  120. pfset Phase.Names "water"

  121. pfset Phase.water.Density.Type Constant

  122. pfset Phase.water.Density.Value 1.0

  123. pfset Phase.water.Viscosity.Type Constant

  124. pfset Phase.water.Viscosity.Value 1.0

  125. #-----------------------------------------------------------------------------

  126. Contaminants

  127. #-----------------------------------------------------------------------------

  128. pfset Contaminants.Names ""

  129. #-----------------------------------------------------------------------------

  130. Gravity

  131. #-----------------------------------------------------------------------------

  132. pfset Gravity 1.0

  133. #-----------------------------------------------------------------------------

  134. Timing (time units is set by units of permeability)

  135. #-----------------------------------------------------------------------------

  136. pfset TimingInfo.BaseUnit 1.0

  137. pfset TimingInfo.StartCount 0.0

  138. pfset TimingInfo.StartTime 0.0

  139. pfset TimingInfo.StopTime 8760.0

  140. pfset TimingInfo.DumpInterval 24.0

  141. pfset TimeStep.Type Constant

  142. pfset TimeStep.Value 1.0

  143. #-----------------------------------------------------------------------------

  144. Porosity

  145. #-----------------------------------------------------------------------------

  146. pfset Geom.Porosity.GeomNames "domain s1 s2 s3 s4 s5 s6 s7 s8 s9"

  147. pfset Geom.domain.Porosity.Type Constant

  148. pfset Geom.domain.Porosity.Value 0.4

  149. pfset Geom.s1.Porosity.Type Constant

  150. pfset Geom.s1.Porosity.Value 0.375

  151. pfset Geom.s2.Porosity.Type Constant

  152. pfset Geom.s2.Porosity.Value 0.39

  153. pfset Geom.s3.Porosity.Type Constant

  154. pfset Geom.s3.Porosity.Value 0.387

  155. pfset Geom.s4.Porosity.Type Constant

  156. pfset Geom.s4.Porosity.Value 0.439

  157. pfset Geom.s5.Porosity.Type Constant

  158. pfset Geom.s5.Porosity.Value 0.489

  159. pfset Geom.s6.Porosity.Type Constant

  160. pfset Geom.s6.Porosity.Value 0.399

  161. pfset Geom.s7.Porosity.Type Constant

  162. pfset Geom.s7.Porosity.Value 0.384

  163. pfset Geom.s8.Porosity.Type Constant

  164. pfset Geom.s8.Porosity.Value 0.482

  165. pfset Geom.s9.Porosity.Type Constant

  166. pfset Geom.s9.Porosity.Value 0.442

  167. #-----------------------------------------------------------------------------

  168. Domain

  169. #-----------------------------------------------------------------------------

  170. pfset Domain.GeomName "domain"

  171. #----------------------------------------------------------------------------

  172. Mobility

  173. #----------------------------------------------------------------------------

  174. pfset Phase.water.Mobility.Type Constant

  175. pfset Phase.water.Mobility.Value 1.0

  176. #-----------------------------------------------------------------------------

  177. Wells

  178. #-----------------------------------------------------------------------------

  179. pfset Wells.Names ""

  180. #-----------------------------------------------------------------------------

  181. Time Cycles

  182. #-----------------------------------------------------------------------------

  183. pfset Cycle.Names "constant"

  184. pfset Cycle.constant.Names "alltime"

  185. pfset Cycle.constant.alltime.Length 1

  186. pfset Cycle.constant.Repeat -1

  187. #-----------------------------------------------------------------------------

  188. Boundary Conditions

  189. #-----------------------------------------------------------------------------

  190. pfset BCPressure.PatchNames [pfget Geom.domain.Patches]

  191. pfset Patch.x-lower.BCPressure.Type FluxConst

  192. pfset Patch.x-lower.BCPressure.Cycle "constant"

  193. pfset Patch.x-lower.BCPressure.alltime.Value 0.0

  194. pfset Patch.y-lower.BCPressure.Type FluxConst

  195. pfset Patch.y-lower.BCPressure.Cycle "constant"

  196. pfset Patch.y-lower.BCPressure.alltime.Value 0.0

  197. pfset Patch.z-lower.BCPressure.Type FluxConst

  198. pfset Patch.z-lower.BCPressure.Cycle "constant"

  199. pfset Patch.z-lower.BCPressure.alltime.Value 0.0

  200. pfset Patch.x-upper.BCPressure.Type FluxConst

  201. pfset Patch.x-upper.BCPressure.Cycle "constant"

  202. pfset Patch.x-upper.BCPressure.alltime.Value 0.0

  203. pfset Patch.y-upper.BCPressure.Type FluxConst

  204. pfset Patch.y-upper.BCPressure.Cycle "constant"

  205. pfset Patch.y-upper.BCPressure.alltime.Value 0.0

  206. pfset Patch.z-upper.BCPressure.Type OverlandFlow

  207. pfset Patch.z-upper.BCPressure.Cycle "constant"

  208. pfset Patch.z-upper.BCPressure.alltime.Value 0.0

  209. #-----------------------------------------------------------------------------

  210. Topo slopes in x-direction

  211. #-----------------------------------------------------------------------------

  212. pfset TopoSlopesX.Type "PFBFile"

  213. pfset TopoSlopesX.GeomNames "domain"

  214. pfset TopoSlopesX.FileName "LW.slopex.pfb"

  215. #-----------------------------------------------------------------------------

  216. Topo slopes in y-direction

  217. #-----------------------------------------------------------------------------

  218. pfset TopoSlopesY.Type "PFBFile"

  219. pfset TopoSlopesY.GeomNames "domain"

  220. pfset TopoSlopesY.FileName "LW.slopey.pfb"

  221. #-----------------------------------------------------------------------------

  222. Mannings coefficient

  223. #-----------------------------------------------------------------------------

  224. pfset Mannings.Type "Constant"

  225. pfset Mannings.GeomNames "domain"

  226. pfset Mannings.Geom.domain.Value 5.52e-6

  227. #-----------------------------------------------------------------------------

  228. Relative Permeability

  229. #-----------------------------------------------------------------------------

  230. pfset Phase.RelPerm.Type VanGenuchten

  231. pfset Phase.RelPerm.GeomNames "domain s1 s2 s3 s4 s5 s6 s7 s8 s9 "

  232. pfset Geom.domain.RelPerm.Alpha 3.5

  233. pfset Geom.domain.RelPerm.N 2.0

  234. pfset Geom.s1.RelPerm.Alpha 3.548

  235. pfset Geom.s1.RelPerm.N 4.162

  236. pfset Geom.s2.RelPerm.Alpha 3.467

  237. pfset Geom.s2.RelPerm.N 2.738

  238. pfset Geom.s3.RelPerm.Alpha 2.692

  239. pfset Geom.s3.RelPerm.N 2.445

  240. pfset Geom.s4.RelPerm.Alpha 0.501

  241. pfset Geom.s4.RelPerm.N 2.659

  242. pfset Geom.s5.RelPerm.Alpha 0.661

  243. pfset Geom.s5.RelPerm.N 2.659

  244. pfset Geom.s6.RelPerm.Alpha 1.122

  245. pfset Geom.s6.RelPerm.N 2.479

  246. pfset Geom.s7.RelPerm.Alpha 2.089

  247. pfset Geom.s7.RelPerm.N 2.318

  248. pfset Geom.s8.RelPerm.Alpha 0.832

  249. pfset Geom.s8.RelPerm.N 2.514

  250. pfset Geom.s9.RelPerm.Alpha 1.585

  251. pfset Geom.s9.RelPerm.N 2.413

  252. #-----------------------------------------------------------------------------

  253. Saturation

  254. #-----------------------------------------------------------------------------

  255. pfset Phase.Saturation.Type VanGenuchten

  256. pfset Phase.Saturation.GeomNames "domain s1 s2 s3 s4 s5 s6 s7 s8 s9 "

  257. pfset Geom.domain.Saturation.Alpha 3.5

  258. pfset Geom.domain.Saturation.N 2.

  259. pfset Geom.domain.Saturation.SRes 0.2

  260. pfset Geom.domain.Saturation.SSat 1.0

  261. pfset Geom.s1.Saturation.Alpha 3.548

  262. pfset Geom.s1.Saturation.N 4.162

  263. pfset Geom.s1.Saturation.SRes 0.000001

  264. pfset Geom.s1.Saturation.SSat 1.0

  265. pfset Geom.s2.Saturation.Alpha 3.467

  266. pfset Geom.s2.Saturation.N 2.738

  267. pfset Geom.s2.Saturation.SRes 0.000001

  268. pfset Geom.s2.Saturation.SSat 1.0

  269. pfset Geom.s3.Saturation.Alpha 2.692

  270. pfset Geom.s3.Saturation.N 2.445

  271. pfset Geom.s3.Saturation.SRes 0.000001

  272. pfset Geom.s3.Saturation.SSat 1.0

  273. pfset Geom.s4.Saturation.Alpha 0.501

  274. pfset Geom.s4.Saturation.N 2.659

  275. pfset Geom.s4.Saturation.SRes 0.000001

  276. pfset Geom.s4.Saturation.SSat 1.0

  277. pfset Geom.s5.Saturation.Alpha 0.661

  278. pfset Geom.s5.Saturation.N 2.659

  279. pfset Geom.s5.Saturation.SRes 0.000001

  280. pfset Geom.s5.Saturation.SSat 1.0

  281. pfset Geom.s6.Saturation.Alpha 1.122

  282. pfset Geom.s6.Saturation.N 2.479

  283. pfset Geom.s6.Saturation.SRes 0.000001

  284. pfset Geom.s6.Saturation.SSat 1.0

  285. pfset Geom.s7.Saturation.Alpha 2.089

  286. pfset Geom.s7.Saturation.N 2.318

  287. pfset Geom.s7.Saturation.SRes 0.000001

  288. pfset Geom.s7.Saturation.SSat 1.0

  289. pfset Geom.s8.Saturation.Alpha 0.832

  290. pfset Geom.s8.Saturation.N 2.514

  291. pfset Geom.s8.Saturation.SRes 0.000001

  292. pfset Geom.s8.Saturation.SSat 1.0

  293. pfset Geom.s9.Saturation.Alpha 1.585

  294. pfset Geom.s9.Saturation.N 2.413

  295. pfset Geom.s9.Saturation.SRes 0.000001

  296. pfset Geom.s9.Saturation.SSat 1.0

  297. #-----------------------------------------------------------------------------

  298. Phase sources:

  299. #-----------------------------------------------------------------------------

  300. pfset PhaseSources.water.Type "Constant"

  301. pfset PhaseSources.water.GeomNames "domain"

  302. pfset PhaseSources.water.Geom.domain.Value 0.0

  303. #----------------------------------------------------------------

  304. CLM Settings:

  305. ------------------------------------------------------------

  306. pfset Solver.LSM CLM

  307. pfset Solver.CLM.CLMFileDir "clm_output/"

  308. pfset Solver.CLM.Print1dOut False

  309. pfset Solver.BinaryOutDir False

  310. pfset Solver.CLM.CLMDumpInterval 24

  311. pfset Solver.CLM.MetForcing 1D

  312. pfset Solver.CLM.MetFileName metfile.dat

  313. pfset Solver.CLM.MetFilePath ./

  314. #pfset Solver.CLM.MetFileNT 24

  315. pfset Solver.CLM.IstepStart 1

  316. pfset Solver.CLM.EvapBeta Linear

  317. pfset Solver.CLM.VegWaterStress Saturation

  318. pfset Solver.CLM.ResSat 0.1

  319. pfset Solver.CLM.WiltingPoint 0.12

  320. pfset Solver.CLM.FieldCapacity 0.98

  321. pfset Solver.CLM.IrrigationType none

  322. #---------------------------------------------------------

  323. Initial conditions: water pressure

  324. #---------------------------------------------------------

  325. pfset ICPressure.Type PFBFile

  326. pfset ICPressure.GeomNames domain

  327. pfset Geom.domain.ICPressure.RefPatch z-upper

  328. pfset Geom.domain.ICPressure.FileName press.init.pfb

  329. #----------------------------------------------------------------

  330. Outputs

  331. ------------------------------------------------------------

  332. #Writing output (all pfb):

  333. pfset SILO.Filetype PDB

  334. pfset Solver.PrintSubsurfData False

  335. pfset Solver.PrintPressure False

  336. pfset Solver.PrintSaturation True

  337. pfset Solver.PrintMask True

  338. pfset Solver.WriteCLMBinary False

  339. pfset Solver.PrintCLM False

  340. pfset Solver.WriteSiloSpecificStorage False

  341. pfset Solver.WriteSiloMannings False

  342. pfset Solver.WriteSiloMask False

  343. pfset Solver.WriteSiloSlopes False

  344. pfset Solver.WriteSiloSubsurfData False

  345. pfset Solver.WriteSiloPressure False

  346. pfset Solver.WriteSiloSaturation True

  347. pfset Solver.WriteSiloEvapTrans False

  348. pfset Solver.WriteSiloEvapTransSum False

  349. pfset Solver.WriteSiloOverlandSum False

  350. pfset Solver.WriteSiloCLM False

  351. #-----------------------------------------------------------------------------

  352. Exact solution specification for error calculations

  353. #-----------------------------------------------------------------------------

  354. pfset KnownSolution NoKnownSolution

  355. #-----------------------------------------------------------------------------

  356. Set solver parameters

  357. #-----------------------------------------------------------------------------

  358. ParFlow Solution

  359. pfset Solver Richards

  360. pfset Solver.TerrainFollowingGrid True

  361. pfset Solver.Nonlinear.VariableDz False

  362. pfset Solver.MaxIter 25000

  363. pfset Solver.Drop 1E-20

  364. pfset Solver.AbsTol 1E-8

  365. pfset Solver.MaxConvergenceFailures 8

  366. pfset Solver.Nonlinear.MaxIter 80

  367. pfset Solver.Nonlinear.ResidualTol 1e-6

  368. new solver settings for Terrain Following Grid

  369. pfset Solver.Nonlinear.EtaChoice EtaConstant

  370. pfset Solver.Nonlinear.EtaValue 0.001

  371. pfset Solver.Nonlinear.UseJacobian True

  372. pfset Solver.Nonlinear.DerivativeEpsilon 1e-16

  373. pfset Solver.Nonlinear.StepTol 1e-30

  374. pfset Solver.Nonlinear.Globalization LineSearch

  375. pfset Solver.Linear.KrylovDimension 70

  376. pfset Solver.Linear.MaxRestarts 2

  377. pfset Solver.Linear.Preconditioner PFMG

  378. pfset Solver.Linear.Preconditioner.PCMatrixType FullJacobian

  379. #-----------------------------------------------------------------------------

  380. Distribute inputs

  381. #-----------------------------------------------------------------------------

  382. pfset ComputationalGrid.NX 41

  383. pfset ComputationalGrid.NY 41

  384. pfset ComputationalGrid.NZ 1

  385. pfdist LW.slopex.pfb

  386. pfdist LW.slopey.pfb

  387. pfset ComputationalGrid.NX 41

  388. pfset ComputationalGrid.NY 41

  389. pfset ComputationalGrid.NZ 50

  390. pfdist IndicatorFile_Gleeson.50z.pfb

  391. pfdist press.init.pfb

  392. #-----------------------------------------------------------------------------

  393. Run Simulation

  394. #-----------------------------------------------------------------------------

  395. set runname "LW"

  396. puts $runname

  397. pfrun $runname

  398. #-----------------------------------------------------------------------------

  399. Undistribute outputs

  400. #-----------------------------------------------------------------------------

  401. pfundist $runname

  402. pfundist press.init.pfb

  403. pfundist LW.slopex.pfb

  404. pfundist LW.slopey.pfb

  405. pfundist IndicatorFile_Gleeson.50z.pfb

  406. puts "ParFlow run Complete"