adding temporary diagnostics

Author: glemieux

biogeochem/FatesSoilBGCFluxMod.F90

@@ -673,6 +673,8 @@ subroutine FluxIntoLitterPools(csite)
     currentPatch => csite%oldest_patch
     flux_patch_loop: do while (associated(currentPatch))
 
+      write(*,*) 'Entering FluxIntoLitterPools for patch ', currentPatch%patchno
+
       associate( &
         bc_out => csite%bc_out(currentPatch%patchno), &
         bc_in => csite%bc_in(currentPatch%patchno) &
@@ -835,6 +837,12 @@ subroutine FluxIntoLitterPools(csite)
           flux_all_si = sum(flux_cel_si(:) * bc_in%dz_decomp_sisl(:)) + &
                         sum(flux_lig_si(:) * bc_in%dz_decomp_sisl(:)) + &
                         sum(flux_lab_si(:) * bc_in%dz_decomp_sisl(:))
+         !  if (flux_all_si > 1e35_r8) then
+             write(fates_log(), *) 'flux_all: ', flux_all_si
+             write(fates_log(), *) 'sum cel, lig, lab: ', sum(flux_cel_si), sum(flux_lig_si), sum(flux_lab_si)
+             write(fates_log(), *) 'sum dz: ', sum(bc_in%dz_decomp_sisl(:))
+            !  call endrun(msg=errMsg(__FILE__, __LINE__))
+         !  end if
 
       end do flux_elem_loop
 

main/FatesInterfaceTypesMod.F90

@@ -1115,6 +1115,7 @@ subroutine CheckInterfaceVariables(this)
       ! Check that the rank matches
       if (this%hlm_vars(i)%data_rank /= this%fates_vars(i)%data_rank) then
         write(*,*) "Rank mismatch for variable: ", this%key(i)
+        write(*,*) "Rank: hlm, fates: ", this%hlm_vars(i)%data_rank, this%fates_vars(i)%data_rank
         call endrun(msg=errMsg(__FILE__, __LINE__))
       end if
 
@@ -1689,20 +1690,43 @@ subroutine UpdateLitterFluxes(this, dtime)
       j = this%filter_litter_flux(i)
 
       ! Update the hlm variables with the fates variables
+      ! write(fates_log(),*) 'Updating litterfluxes: key:', this%GetRegistryVariableKey(j)
       call this%hlm_vars(j)%Update(this%fates_vars(j), scalar=dtime)
 
+      ! Dump HLM variable
+      ! write(fates_log(),*) 'Dumping HLM variable'
+      ! call this%hlm_vars(j)%Dump()
+      ! write(fates_log(),*) 'Dumping FATES variable'
+      ! call this%fates_vars(j)%Dump()
+
     end do
 
     ! Update the HLM variable with the total litterfall
+    ! write(fates_log(),*) 'Updating total litterfall'
     j = this%GetRegistryVariableIndex(hlm_fates_litter_carbon_total)
     call this%hlm_vars(j)%Update(this%fates_vars(j))
 
     if (hlm_parteh_mode == prt_cnp_flex_allom_hyp) then
       j = this%GetRegistryVariableIndex(hlm_fates_litter_phosphorus_total)
       call this%hlm_vars(j)%Update(this%fates_vars(j))
 
+      ! Dump HLM variable
+      ! write(fates_log(),*) 'Dumping HLM variable'
+      ! call this%hlm_vars(j)%Dump()
+      ! write(fates_log(),*) 'Dumping FATES variable'
+      ! call this%fates_vars(j)%Dump()
+
       j = this%GetRegistryVariableIndex(hlm_fates_litter_nitrogen_total)
       call this%hlm_vars(j)%Update(this%fates_vars(j))
+
+      ! Dump HLM variable
+      ! write(fates_log(),*) 'Dumping HLM variable'
+      ! call this%hlm_vars(j)%Dump()
+      ! write(fates_log(),*) 'Dumping FATES variable'
+      ! call this%fates_vars(j)%Dump()
+
+      ! write(fates_log(),*) 'Finished updating litter fluxes'
+
     end if