fix missing czil1 in vege_flux

wzzheng90 · Mar 7, 2022 · f093f77 · f093f77
1 parent 53c0c7a
commit f093f77
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/physics/module_sf_noahmplsm.f90 b/physics/module_sf_noahmplsm.f90
@@ -3887,6 +3887,7 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
   real (kind=kind_phys), intent(in   ) :: shdfac      ! greeness vegetation fraction (-)
   real (kind=kind_phys), intent(inout) :: ustarx      ! friction velocity
   real (kind=kind_phys), intent(  out) :: csigmaf1    !
+  real (kind=kind_phys)                :: czil1       ! canopy based czil
 
   real (kind=kind_phys) :: t, tdc       !kelvin to degree celsius with limit -50 to +50