Dimension in soil parameter

[changliao1025]

real(r8) ,pointer :: sand3d (:,:,:) ! read in - soil texture: percent sand (needs to be a pointer for use in ncdio) allocate(sand3d(begg:endg,max_topounits,nlevsoifl))

Grid level, not pft or vegetation based.

[changliao1025]

To support the multiple hillslopes, it is also desired to keep the dimension information.

[changliao1025]

integer :: nlevsoi ! number of hydrologically active soil layers integer :: nlevsoifl ! number of soil layers on input file integer :: nlevgrnd ! number of ground layers ! (includes lower layers that are hydrologically inactive)

[changliao1025]

https://nowlin.css.msu.edu/software/triangle_form.html

[changliao1025]

Texture is : /sandy clay loam 66: Texture is : /sandy clay loam 66: Texture is : /sandy clay loam 66: Texture is : /sandy clay loam 92: Texture is : /sandy clay loam/sandy clay 92: Texture is : /clay loam 92: Texture is : /clay loam 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /sandy clay loam/sandy clay 92: Texture is : /sandy clay loam/sandy clay 92: Texture is : /clay loam 92: Texture is : /clay loam 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay 92: Texture is : /clay

[changliao1025]

Texture is : /clay 59: Texture not found for: 9.9999996E+35 9.9999996E+35 59: Texture not found for: 9.9999996E+35 9.9999996E+35 81: Texture is : /sandy clay loam/sandy clay 81: Texture is : /sandy clay loam/sandy clay 81: Texture is : /clay loam

[changliao1025]

For simplicity, we only consider one soil texture on the boundary.