From 80dc3f48052101723c394b61da68563d2b9f1aff Mon Sep 17 00:00:00 2001
From: Victor Stepanenko <stepanen@srcc.msu.ru>
Date: Tue, 9 Aug 2022 18:59:05 +0300
Subject: [PATCH] Important bug fixed in snowtemp

---
 source/model/snowtemp.f90 | 8 ++------
 1 file changed, 2 insertions(+), 6 deletions(-)

diff --git a/source/model/snowtemp.f90 b/source/model/snowtemp.f90
index 39d4c72..4d3c403 100644
--- a/source/model/snowtemp.f90
+++ b/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
 
 
-- 
GitLab