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.