Merge pull request #1376 from rgknox/trivial-zenith-defaults

fix to trivial radiation transmission

Author: glemieux

main/FatesInterfaceMod.F90

@@ -373,13 +373,13 @@ subroutine zero_bcs(fates,s)
     fates%bc_out(s)%rssun_pa(:)     = 0.0_r8
     fates%bc_out(s)%rssha_pa(:)     = 0.0_r8
 
-    fates%bc_out(s)%albd_parb(:,:) = 0.0_r8   ! zero albedo, soil absorbs all rad
-    fates%bc_out(s)%albi_parb(:,:) = 0.0_r8   ! zero albedo, soil absorbs all rad
+    fates%bc_out(s)%albd_parb(:,:) = 1.0_r8   ! zero albedo, soil absorbs all rad
+    fates%bc_out(s)%albi_parb(:,:) = 1.0_r8   ! zero albedo, soil absorbs all rad
     fates%bc_out(s)%fabd_parb(:,:) = 0.0_r8   ! no rad absorbed by veg
     fates%bc_out(s)%fabi_parb(:,:) = 0.0_r8   ! no rad absorbed by veg
-    fates%bc_out(s)%ftdd_parb(:,:) = 1.0_r8   ! rad flux to soil at bottom of veg is 100%
-    fates%bc_out(s)%ftid_parb(:,:) = 1.0_r8   ! rad flux to soil at bottom of veg is 100%
-    fates%bc_out(s)%ftii_parb(:,:) = 1.0_r8   ! rad flux to soil at bottom of veg is 100%
+    fates%bc_out(s)%ftdd_parb(:,:) = 0.0_r8   ! rad flux to soil at bottom of veg is 100%
+    fates%bc_out(s)%ftid_parb(:,:) = 0.0_r8   ! rad flux to soil at bottom of veg is 100%
+    fates%bc_out(s)%ftii_parb(:,:) = 0.0_r8   ! rad flux to soil at bottom of veg is 100%
 
     fates%bc_out(s)%elai_pa(:)   = 0.0_r8
     fates%bc_out(s)%esai_pa(:)   = 0.0_r8

radiation/FatesRadiationDriveMod.F90

@@ -107,14 +107,15 @@ subroutine FatesNormalizedCanopyRadiation(sites, bc_in, bc_out )
           if_bareground: if(currentpatch%nocomp_pft_label.ne.nocomp_bareground)then
 
              ! Initialize output boundary conditions with trivial assumption
-             ! of a black body soil and fully transmitting canopy
-             bc_out(s)%albd_parb(ifp,:)            = 0._r8
-             bc_out(s)%albi_parb(ifp,:)            = 0._r8
+             ! This matches CLM/ELM
+             ! Albedo is perfect reflector, no flux into or through canopy
+             bc_out(s)%albd_parb(ifp,:)            = 1._r8
+             bc_out(s)%albi_parb(ifp,:)            = 1._r8
              bc_out(s)%fabi_parb(ifp,:)            = 0._r8
              bc_out(s)%fabd_parb(ifp,:)            = 0._r8
-             bc_out(s)%ftdd_parb(ifp,:)            = 1._r8
-             bc_out(s)%ftid_parb(ifp,:)            = 1._r8
-             bc_out(s)%ftii_parb(ifp,:)            = 1._r8
+             bc_out(s)%ftdd_parb(ifp,:)            = 0._r8
+             bc_out(s)%ftid_parb(ifp,:)            = 0._r8
+             bc_out(s)%ftii_parb(ifp,:)            = 0._r8
 
              ! Zero diagnostics
              currentPatch%f_sun      (:,:,:) = 0._r8