Fixed tempunstag allocation issue

Glide runs were segfaulting because the code was trying to initialize tempstag
without allocating it. Easily fixed with 'if (associated)' logic.

Author: whlipscomb

libglide/glide_types.F90

@@ -1947,9 +1947,10 @@ subroutine glide_allocarr(model)
     ! MJH - Set temp and flwa to physically unrealistic values so we can tell later if 
     !       arrays were initialized correctly
     model%temper%temp(:,:,:) = unphys_val  ! large negative number
-    model%temper%tempunstag(:,:,:) = unphys_val
     model%temper%flwa(:,:,:) = unphys_val
     model%temper%dissip(:,:,:) = 0.d0
+    if (associated(model%temper%tempunstag)) &
+         model%temper%tempunstag(:,:,:) = unphys_val
 
     call coordsystem_allocate(model%general%ice_grid,  model%temper%bheatflx)
     call coordsystem_allocate(model%general%ice_grid,  model%temper%bwat)

libglissade/glissade_therm.F90

@@ -102,7 +102,8 @@ subroutine glissade_init_therm (temp_init,  is_restart,        &
                                   artm,                          &
                                   acab,                          &
                                   bheatflx,                      &
-                                  pmp_offset, temp, tempunstag)
+                                  pmp_offset,                    &
+                                  temp,       tempunstag)
 
     ! initialization subroutine for higher-order dycores, where temperature is defined at
     ! the midpoint of each layer plus the upper and lower surfaces