Skip to content

Commit

Permalink
Progress towards conservative boundary flux tendencies
Browse files Browse the repository at this point in the history
RE: Issue NOAA-GFDL/MOM6-examples#161
Without the v_extensive flag, the boundary_forcing_heat_tendency and
boundary_forcing_salt_tendency terms were incorrectly being output onto
the diagnostic grid via the ALE remapping core instead of the
reintegration routines. Additionally, the diagnostic vertical grid
needed to be built after applyBoundaryFluxesInOut. Despite these two fixes,
there remains  ~10% error between the native and remapped output.
native and remapped outputs.
  • Loading branch information
Andrew Shao committed Sep 20, 2017
1 parent 3d60b91 commit 2987394
Showing 1 changed file with 5 additions and 2 deletions.
7 changes: 5 additions & 2 deletions src/parameterizations/vertical/MOM_diabatic_driver.F90
Original file line number Diff line number Diff line change
Expand Up @@ -809,6 +809,7 @@ subroutine diabatic(u, v, h, tv, Hml, fluxes, visc, ADp, CDp, dt, G, GV, CS)

! diagnose the tendencies due to boundary forcing
if(CS%boundary_forcing_tendency_diag) then
call diag_update_remap_grids(CS%diag)
call diagnose_boundary_forcing_tendency(tv, h, temp_diag, saln_diag, h_diag, dt, G, GV, CS)
endif

Expand Down Expand Up @@ -2183,14 +2184,16 @@ subroutine diabatic_driver_init(Time, G, GV, param_file, useALEalgorithm, diag,

CS%id_boundary_forcing_heat_tend = register_diag_field('ocean_model',&
'boundary_forcing_heat_tendency', diag%axesTL, Time, &
'Boundary forcing heat tendency','Watts/m2')
'Boundary forcing heat tendency','Watts/m2', &
v_extensive = .true.)
if (CS%id_boundary_forcing_heat_tend > 0) then
CS%boundary_forcing_tendency_diag = .true.
endif

CS%id_boundary_forcing_salt_tend = register_diag_field('ocean_model',&
'boundary_forcing_salt_tendency', diag%axesTL, Time, &
'Boundary forcing salt tendency','kg m-2 s-1')
'Boundary forcing salt tendency','kg m-2 s-1', &
v_extensive = .true.)
if (CS%id_boundary_forcing_salt_tend > 0) then
CS%boundary_forcing_tendency_diag = .true.
endif
Expand Down

0 comments on commit 2987394

Please sign in to comment.