From b826f72507e92b4215e05e85c864368a292e0e66 Mon Sep 17 00:00:00 2001
From: Bill Sacks <sacks@ucar.edu>
Date: Thu, 4 Apr 2024 19:51:54 -0600
Subject: [PATCH] Add more prints

---
 src/biogeophys/CanopyFluxesMod.F90 | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/src/biogeophys/CanopyFluxesMod.F90 b/src/biogeophys/CanopyFluxesMod.F90
index 450ebc25b2..d04fe4b252 100644
--- a/src/biogeophys/CanopyFluxesMod.F90
+++ b/src/biogeophys/CanopyFluxesMod.F90
@@ -48,6 +48,7 @@ module CanopyFluxesMod
   use SoilWaterRetentionCurveMod, only : soil_water_retention_curve_type
   use LunaMod               , only : Update_Photosynthesis_Capacity, Acc24_Climate_LUNA,Acc240_Climate_LUNA,Clear24_Climate_LUNA
   use NumericsMod           , only : truncate_small_values
+  use spmdMod, only : iam
   !
   ! !PUBLIC TYPES:
   implicit none
@@ -643,6 +644,10 @@ subroutine CanopyFluxes(bounds,  num_exposedvegp, filter_exposedvegp,
 
       dtime = get_step_size_real()
 
+      if (iam == 20) then
+         write(iulog,*) 'WJS: Start of CanopyFluxes: ', snocan(24)
+      end if
+
       ! Make a local copy of the exposedvegp filter. With the current implementation,
       ! this is needed because the filter is modified in the iteration loop.
       !
@@ -1600,6 +1605,9 @@ subroutine CanopyFluxes(bounds,  num_exposedvegp, filter_exposedvegp,
 
          ! save before updating
          snocan_baseline(p) = snocan(p)
+         if (iam == 20 .and. p == 24) then
+            write(iulog,*) 'WJS: CanopyFluxes before: ', snocan(p)
+         end if
 
          ! Update dew accumulation (kg/m2)
          if (t_veg(p) > tfrz ) then ! above freezing, update accumulation in liqcan
@@ -1624,6 +1632,9 @@ subroutine CanopyFluxes(bounds,  num_exposedvegp, filter_exposedvegp,
             end if
          end if
 
+         if (iam == 20 .and. p == 24) then
+            write(iulog,*) 'WJS: CanopyFluxes after: ', snocan(p)
+         end if
       end do
 
       ! Remove snocan that got reduced by more than a factor of rel_epsilon
@@ -1751,6 +1762,10 @@ subroutine CanopyFluxes(bounds,  num_exposedvegp, filter_exposedvegp,
          write(iulog,*) 'energy balance in canopy ',p,', err=',err(p)
       end do
 
+      if (iam == 20) then
+         write(iulog,*) 'WJS: End of CanopyFluxes: ', snocan(24)
+      end if
+
     end associate