Error statement failing due to mismatch in kmax_upper_shell second dimension index in Report1DError

[glemieux]

During a failure while running with hydro on, I'm finding that the error message isn't fully reporting due to crashing on line 4010:

fates/biogeophys/FatesPlantHydraulicsMod.F90

Line 4010 in 42d804b

write(fates_log(),*) ' ',1._r8/(1._r8/(csite_hydr%kmax_lower_shell(ilayer,1)*aroot_frac_plant) + 1._r8/(csite_hydr%kmax_upper_shell(ilayer,2)*aroot_frac_plant))

The issue is that the second dimension of kmax_upper_shell is being hardcoded to a value about the allocated size. I'm not sure if the write statement needs to have a reference to the second dimension variable (nshell) or if this is somehow related to the failure mode that caused the call of Report1DError in this instance.