mom_geothermal mom_geothermal::geothermal_cs subroutine, public subroutine, public mom_geothermal::geothermal_entraining (h, tv, dt, ea, eb, G, GV, US, CS, halo) geothermal_entraining h h tv tv dt dt ea ea eb eb G G GV GV US US 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 GEOTHERMAL_THICKNESS 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)? 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 [T ~> 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]. us A dimensional unit scaling type cs The control structure returned by a previous call to geothermal_init. halo Halo width over which to work mom_error_handler::mom_error subroutine, public subroutine, public mom_geothermal::geothermal_in_place (h, tv, dt, G, GV, US, CS, halo) geothermal_in_place h h tv tv dt dt G G GV GV US US CS CS halo halo Applies geothermal heating to the bottommost layers that occur within GEOTHERMAL_THICKNESS of the bottom, by simply heating the water in place. Any heat that can not be applied to the ocean is returned (WHERE)? 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 [T ~> s]. us A dimensional unit scaling type cs The control structure returned by a previous call to geothermal_init. halo Halo width over which to work subroutine, public subroutine, public mom_geothermal::geothermal_init (Time, G, GV, US, param_file, diag, CS, useALEalgorithm) geothermal_init Time Time G G GV GV US US param_file param_file diag diag CS CS useALEalgorithm useALEalgorithm Initialize parameters and allocate memory associated with the geothermal heating module. time Current model time. g The ocean's grid structure. gv The ocean's vertical grid structure. us A dimensional unit scaling type 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. usealealgorithm logical for whether to use ALE remapping mom_diabatic_driver::diabatic_driver_init subroutine, public subroutine, public mom_geothermal::geothermal_end (CS) geothermal_end CS CS Clean up and deallocate memory associated with the geothermal heating module. cs Geothermal heating control structure that will be deallocated in this subroutine. 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.