Add option to write 3d soilt, soilw and soill variables to a history file

[DusanJovic-NOAA]

Description

This PR adds three 3D soil variables (soilt, soilw and soill) to a history file if corresponding entries are specified in diag_table file. The new 3d versions of these variables is identical to a set of 2d variables currently available (soilt1, soilt2, etc)

Issue(s) addressed

Link the issues to be closed with this PR, whether in this repository, or in another repository.
(Remember, issues should always be created before starting work on a PR branch!)

• fixes Add vertical dimension to certain surface fields ufs-community/ufs-weather-model#1055

Testing

How were these changes tested? Yes
What compilers / HPCs was it tested with? Intel and GNU / Hera
Are the changes covered by regression tests? Yes.
Have the ufs-weather-model regression test been run? Yes.
On what platform? Hera

• Will the code updates change regression test baseline? Yes, control_CubedSphereGrid_parallel regression test is modified to output 3d version of soilt, soilw and soill variables.

Dependencies

N/A

[FernandoAndrade-NOAA]

@jkbk2004 could you please approve this when available, thanks! Testing on #2062 has completed successfully.