Methane production refinement at the talik depth

source/model/methane_mod.f90

@@ -399,7 +399,7 @@ forall (i = 1:gs%ns, z(i) > rnroot) forg(i) = veg * exp((rnroot - z(i))*10.)
 !     diff(i) = DIFF_WATER_METHANE(Tsoil(i)) * &
 !     & tortuosity_coef*fcoarse
 !   end if                                             ! soil layers
-!   diff(i) = 0.25*( por(i) + por(i-1) )*tortuosity_coef * &
+!   diff(i) = 0.25/\cf*( por(i) + por(i-1) )*tortuosity_coef * &
 !   & ( (water_porvolrat(i) + water_porvolrat(i-1)) * &
 !   & DIFF_WATER_METHANE( 0.5*(Tsoil(i) + Tsoil(i-1)) ) + &
 !   & (air_porvolrat(i) + air_porvolrat(i-1)) * &
@@ -416,14 +416,14 @@ forall (i = 1:gs%ns, z(i) > rnroot) forg(i) = veg * exp((rnroot - z(i))*10.)
      ft1 = 1. ! Note, methane generation through old organics decomposition
               ! is zero below talik automatically, see the code below
      if (Tsoil(i) .gt. MELTINGPOINT(Sals(i)/wl(i),pressoil(i),tricemethhydr%par)) then
-       if (i /= i_talik) then
+       if (i /= i_talik .or. i == 1 .or. i == gs%ns) then
          ft = 1.
        else
          ft = max((h_talik - 0.5*(z(i-1) + z(i)) ) / &
          & (0.5*(z(i+1) - z(i-1))),0.e0_ireals)
        endif
      else
-       if (i /= i_talik) then
+       if (i /= i_talik .or. i == 1 .or. i == gs%ns) then
          ft = 0.
        else
          if (i == gs%ns) then