fixing height of layers 2+ in interpolate_wind2fire_height

openwfm · Jun 12, 2020 · 0678af2 · 0678af2
1 parent 80a3806
commit 0678af2
Showing 1 changed file with 3 additions and 6 deletions.
diff --git a/phys/module_fr_sfire_atm.F b/phys/module_fr_sfire_atm.F
@@ -2072,7 +2072,7 @@ subroutine interpolate_wind2fire_height(id,       & ! to identify debugging prin
   integer::itst,itet,jtst,jtet
   integer::iftst,iftet,jftst,jftet
   real:: wjcb,wjcm,wicb,wicm,ht,i_g,loght,zr,ht_last,logwh,wh,loght_last,uk,vk,uk1,vk1,z0,logz0
-  real, dimension (its-1:ite+1,kds:kde,jts-1:jte+1):: z
+  real, dimension (its-1:ite+1,kds:kde,jts-1:jte+1)::z,zw
   character(len=128)::msg
 
 !*** executable
@@ -2099,15 +2099,12 @@ subroutine interpolate_wind2fire_height(id,       & ! to identify debugging prin
   do j = jtst,jtet
     do k=kds,kdmax+1
       do i = itst,itet       
-        z(i,k,j) = (ph(i,k,j)+phb(i,k,j))*i_g      ! altitude of the bottom w-point 
-
-        ! print *,'i,k,j=',i,k,j,'ph, phb, z=',ph(i,k,j),phb(i,k,j),z(i,k,j)
-
+        zw(i,k,j) = (ph(i,k,j)+phb(i,k,j))*i_g      ! altitude of the bottom w-point 
       enddo
     enddo
     do k=kds,kdmax
       do i = itst,itet       
-        z(i,k,j) = (z(i,k,j)+z(i,k+1,j))*0.5 - z(i,kds,j)   ! height of the cell center 
+        z(i,k,j) = (zw(i,k,j)+zw(i,k+1,j))*0.5 - zw(i,kds,j) ! height of the cell center above the bottom
       enddo
     enddo
   enddo