fix restart problem with nitrogen flux history fields

modified: src/physics/cam/cam_diagnostics.F90 modified: src/physics/cam/physpkg.F90
fvitt · Dec 26, 2024 · 82a3c9c · 82a3c9c
1 parent 24e012d
commit 82a3c9c
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 7 deletions.
diff --git a/src/physics/cam/cam_diagnostics.F90 b/src/physics/cam/cam_diagnostics.F90
@@ -1906,13 +1906,6 @@ subroutine diag_export(cam_out)
       call outfld('a2x_DSTDRY4',  cam_out%dstdry4,  pcols, lchnk)
     end if
 
-    if (associated(cam_out%nhx_nitrogen_flx)) then
-       call outfld('a2x_NHXDEP', cam_out%nhx_nitrogen_flx, pcols, lchnk)
-    end if
-    if (associated(cam_out%noy_nitrogen_flx)) then
-       call outfld('a2x_NOYDEP', cam_out%noy_nitrogen_flx, pcols, lchnk)
-    end if
-
   end subroutine diag_export
 
 !#######################################################################

diff --git a/src/physics/cam/physpkg.F90 b/src/physics/cam/physpkg.F90
@@ -1988,6 +1988,14 @@ subroutine tphysac (ztodt,   cam_in,  &
 
     call clybry_fam_set( ncol, lchnk, map2chm, state%q, pbuf )
 
+    ! output these here -- after updates by chem_timestep_tend or export_fields within the current time step
+    if (associated(cam_out%nhx_nitrogen_flx)) then
+       call outfld('a2x_NHXDEP', cam_out%nhx_nitrogen_flx, pcols, lchnk)
+    end if
+    if (associated(cam_out%noy_nitrogen_flx)) then
+       call outfld('a2x_NOYDEP', cam_out%noy_nitrogen_flx, pcols, lchnk)
+    end if
+
   end subroutine tphysac
 
   subroutine tphysbc (ztodt, state,  &