Important bug fixed in snowtemp

source/model/snowtemp.f90

@@ -39,6 +39,8 @@ real(kind=ireals)   , intent(in) :: hour
 real(kind=ireals)   , intent(in) :: phi
 !real(kind=ireals)   , intent(in) :: extice
 real(kind=ireals)   , intent(in) :: fetch
+real(kind=ireals), intent(inout) :: a(1:nveclen), b(1:nveclen), &
+&                                   c(1:nveclen), d(1:nveclen), Temp(1:nveclen)
 
 
 real(kind=ireals) :: xx
@@ -52,7 +54,6 @@ real(kind=ireals) :: snowmass_init
 real(kind=ireals) :: CCT(ML)
 real(kind=ireals) :: dt
 
-real(kind=ireals), allocatable :: a(:), b(:), c(:), d(:), Temp(:)
 
 integer(kind=iintegers) :: i
 integer(kind=iintegers) :: iyear
@@ -61,9 +62,6 @@ integer(kind=iintegers) :: iday
 
 SAVE
 
-allocate(a(1:nveclen), b(1:nveclen), c(1:nveclen), &
-&        d(1:nveclen), Temp(1:nveclen))
-
 do i = itop,ms
   T(i) = Tsn(i)    
 enddo
@@ -91,8 +89,6 @@ totalmelts = totalmelts + snmelt*dt
 
 !if (hs1 == 0) hs1=0.00001
 
-deallocate(a, b, c, d, Temp)
-
 END SUBROUTINE SNOWTEMP