namespace mom_geothermal

Overview

Implemented geothermal heating at the ocean bottom. More…

namespace mom_geothermal {

// global functions

subroutine, public geothermal(
    h h,
    tv tv,
    dt dt,
    ea ea,
    eb eb,
    G G,
    GV GV,
    CS CS,
    halo halo
    );

subroutine, public geothermal_init(Time Time, G G, param_file param_file, diag diag, CS CS);
subroutine, public geothermal_end(CS CS);

} // namespace mom_geothermal

Detailed Documentation

Implemented geothermal heating at the ocean bottom.

Geothermal heating can be added either in a layered isopycnal mode, in which the heating raises the density of the layer to the target density of the layer above, and then moves the water into that layer, or in a simple Eulerian mode, in which the bottommost GEOTHERMAL_THICKNESS are heated. Geothermal heating will also provide a buoyant source of bottom TKE that can be used to further mix the near-bottom water. In cold fresh water lakes where heating increases density, water should be moved into deeper layers, but this is not implemented yet.

Global Functions

subroutine, public geothermal(
    h h,
    tv tv,
    dt dt,
    ea ea,
    eb eb,
    G G,
    GV GV,
    CS CS,
    halo halo
    )

Applies geothermal heating, including the movement of water between isopycnal layers to match the target densities. The heating is applied to the bottommost layers that occur within ### of the bottom. If the partial derivative of the coordinate density with temperature is positive or very small, the layers are simply heated in place. Any heat that can not be applied to the ocean is returned (WHERE)?

Parameters:

g	The ocean’s grid structure.
gv	The ocean’s vertical grid structure.
h	Layer thicknesses [H ~> m or kg m-2]
tv	A structure containing pointers to any available thermodynamic fields. Absent fields have NULL ptrs.
dt	Time increment [s].
ea	The amount of fluid moved downward into a layer; this should be increased due to mixed layer detrainment [H ~> m or kg m-2]
eb	The amount of fluid moved upward into a layer; this should be increased due to mixed layer entrainment [H ~> m or kg m-2].
cs	The control structure returned by a previous call to geothermal_init.
halo	Halo width over which to work

subroutine, public geothermal_init(
    Time Time,
    G G,
    param_file param_file,
    diag diag,
    CS CS
    )

Initialize parameters and allocate memory associated with the geothermal heating module.

Parameters:

time	Current model time.
g	The ocean’s grid structure.
param_file	A structure to parse for run-time parameters.
diag	Structure used to regulate diagnostic output.
cs	Pointer pointing to the module control structure.

subroutine, public geothermal_end(CS CS)

Clean up and deallocate memory associated with the geothermal heating module.

Parameters:

cs	Geothermal heating control structure that will be deallocated in this subroutine.