cleaning up some print statements on zenith processing improvements f…

…or FATES
rgknox · Jan 29, 2025 · 9940262 · 9940262
1 parent 2c751e9
commit 9940262
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 6 deletions.
diff --git a/src/biogeophys/BalanceCheckMod.F90 b/src/biogeophys/BalanceCheckMod.F90
@@ -956,6 +956,7 @@ subroutine BalanceCheck( bounds, &
 
            indexp = maxloc( abs(errsol(bounds%begp:bounds%endp)), 1 , mask = (errsol(bounds%begp:bounds%endp) /= spval) ) + bounds%begp -1
            indexg = patch%gridcell(indexp)
+           indexc = patch%column(indexp)
            write(iulog,*)'WARNING:: BalanceCheck, solar radiation balance error (W/m2)'
            write(iulog,*)'nstep         = ',nstep
            write(iulog,*)'errsol        = ',errsol(indexp)

diff --git a/src/biogeophys/SurfaceAlbedoMod.F90 b/src/biogeophys/SurfaceAlbedoMod.F90
@@ -267,8 +267,6 @@ subroutine UpdateZenithAngles(bounds, surfalb_inst, nextsw_cday, declinp1)
          else
             coszen_col(c) = coszen_grc(g)
          endif
-
-         print*,"COSZEN_COL:",c,coszen_col(c)
 
       end do
 
@@ -481,14 +479,11 @@ subroutine SurfaceAlbedo(bounds,nc,  &
           fabi_sha_z    =>    surfalb_inst%fabi_sha_z_patch         & ! Output:  [real(r8) (:,:) ]  absorbed shaded leaf diffuse PAR (per unit lai+sai) for each canopy layer
           )
 
-
-     print*,"PEFORMING ALBEDO"
-
     ! Apply column level zenith angles to the patch level
     do fp = 1,num_nourbanp
        p = filter_nourbanp(fp)
        c = patch%column(p)
-       if(col%is_fates(p))then
+       if(patch%is_fates(p))then
           coszen_patch(p) = spval
        else
           coszen_patch(p) = coszen_col(c)