Floating-point divide by zero exception in ssnow%smp_hys computation

[SeanBryan51]

Hacking a temporary fix for #395 and running CABLE-MPI offline (main branch - commit 95b9b5e) using the crujra_accessN96_1h configuration results in the following divide by zero exception:

[gadi-cpu-clx-2663:671237:0:671237] Caught signal 8 (Floating point exception: floating-point divide by zero)
==== backtrace (tid: 671237) ====
 0 0x0000000000012d20 __funlockfile()  :0
 1 0x00000000006c2338 cable_param_module_mp_derived_parameters_()  /home/189/sb8430/cable/src/offline/cable_parameters.F90:2320
 2 0x00000000006228a2 cable_input_module_mp_load_parameters_()  /home/189/sb8430/cable/src/offline/cable_input.F90:2829
 3 0x0000000000416cd2 cable_mpimaster_mp_mpidrv_master_()  /home/189/sb8430/cable/src/offline/cable_mpimaster.F90:601
 4 0x000000000040e5fc MAIN__()  /home/189/sb8430/cable/src/offline/cable_mpidrv.F90:54
 5 0x000000000040da22 main()  ???:0
 6 0x000000000003a7e5 __libc_start_main()  ???:0
 7 0x000000000040d92e _start()  ???:0
=================================
forrtl: error (75): floating point exception
Image              PC                Routine            Line        Source             
cable-mpi          0000000000C7E474  Unknown               Unknown  Unknown
libpthread-2.28.s  0000155548397D20  Unknown               Unknown  Unknown
cable-mpi          00000000006C2338  cable_param_modul        2320  cable_parameters.F90
cable-mpi          00000000006228A2  cable_input_modul        2829  cable_input.F90
cable-mpi          0000000000416CD2  cable_mpimaster_m         601  cable_mpimaster.F90
cable-mpi          000000000040E5FC  MAIN__                     54  cable_mpidrv.F90
cable-mpi          000000000040DA22  Unknown               Unknown  Unknown
libc-2.28.so       0000155547C677E5  __libc_start_main     Unknown  Unknown
cable-mpi          000000000040D92E  Unknown               Unknown  Unknown

The exception occurs on this line of the code:

CABLE/src/offline/cable_parameters.F90

Line 2320 in 95b9b5e

(ssnow%ssat_hys(i,k)-ssnow%watr_hys(i,k)) )**&

It looks like ssnow%ssat_hys(i,k) and ssnow%watr_hys(i,k) are both uninitialised and contain the same garbage value, causing the subtraction of the two values to result in divide by zero.

Steps to reproduce (Gadi)

Apply the following patch to fix the error described in #395 (WARNING - this patch is untested and should not be used for work other than reproducing this issue):

diff --git a/src/offline/cable_parameters.F90 b/src/offline/cable_parameters.F90
index b6133f6..c741eaf 100644
--- a/src/offline/cable_parameters.F90
+++ b/src/offline/cable_parameters.F90
@@ -3340,11 +3340,11 @@ CONTAINS
     totdepth = 0.0
     DO is = 1, ms-1
        totdepth = totdepth + soil_zse(is) * 100.0  ! unit in centimetres
-       veg%froot(:, is) = MIN( 1.0, 1.0-veg%rootbeta(:)**totdepth )
+       veg%froot(ifmp:fmp, is) = MIN( 1.0, 1.0-veg%rootbeta(ifmp:fmp)**totdepth )
     END DO
-    veg%froot(:, ms) = 1.0 - veg%froot(:, ms-1)
+    veg%froot(ifmp:fmp, ms) = 1.0 - veg%froot(ifmp:fmp, ms-1)
     DO is = ms-1, 2, -1
-       veg%froot(:, is) = veg%froot(:, is)-veg%froot(:,is-1)
+    veg%froot(ifmp:fmp, is) = veg%froot(ifmp:fmp, is)-veg%froot(ifmp:fmp,is-1)
     END DO

   END SUBROUTINE init_veg_from_vegin

The steps to reproduce the error are the same as that described in #395.