Update glm_zones.F90

Fix copy_to_zone to divide by layer number (num layers above zone)
AquaticEcoDynamics · Feb 8, 2024 · b6aec13 · b6aec13
1 parent 78f84ff
commit b6aec13
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/glm_zones.F90 b/src/glm_zones.F90
@@ -309,6 +309,8 @@ SUBROUTINE copy_to_zone(x_cc, x_diag, x_diag_hz, wlev)
 
    DO zon=1,a_zones
       z_cc(zon,1:nvars) = z_cc(zon,1:nvars)/zcount(zon)
+      z_diag(zon,:)     = z_diag(zon,:)/zcount(zon)
+      z_diag_hz(zon,:)  = z_diag_hz(zon,:)/zcount(zon)
    ENDDO
 
    WHERE (zcount /= 0)