Next: endstep.F
Up: The Proudman Oceanographic Laboratory
Previous: diffuse.F
Subsections
downwell.F
subroutine downwell
This subroutine updates the model temperature using heat fluxes
(hfl_in and hfl_out) calculated in (e.g.) 43.2.
These are already in units Cms
. The surface cooling is applied to the top box:

(37) 
while penetrating radiation heats down the water column according to

(38) 
The transmissivity is set to either a constant value or a simple
function of water depth (J.I. Allen , per com.)
none
 i,j,kLocal indices
 dtmp change in temperature for this time step
 rD,rLambda transmissivity
 z1,z2 depth below the sea surface at the bottom and top of the
grid box
 Rad0,Rad1,Rad2 incident radiation at the surface,
and at the bottom and top of the grid box
 surfaceout surface cooling
 first .true. first time this subroutine is called
 tmp rlambda2
 none
 set transmissivity
 calculate surface cooling
 calculate penetrating radiation
 update temperature
 none
 baroc
 vary_lambda use JIA's formula to calculate horizontally varying
transmissivity
 SCOORD Use horizontally varying vertical coordinates
What to do with the radiation reaching the sea bed is still not established.
It is lost in current version, but can be returned to the bottom box.
subroutine downwell_bio
This subroutine updates the model temperature using heat fluxes
(hfl_in and hfl_out) calculated in (e.g.) 43.2.
These are already in units Cms
. The surface cooling is applied to the top box:

(39) 
while penetrating radiation heats down the water column,
with a variable transmissivity, according to

(40) 
The transmissivity is initially set to a simple
function of water depth (J.I. Allen , per com.) put can depend (elsewhere)
on (e.g.) biological and sediment concentrations.
none
 i,j,kLocal indices
 dtmp change in temperature for this time step
 rD,rLambda transmissivity
 dz width of this grid box the
 Rad0 incident radiation at the surface
 surfaceout surface cooling
 first .true. first time this subroutine is called
 tmp rlambda3
 none
 set transmissivity
 calculate surface cooling
 calculate penetrating radiation and integrate transmissivity
 update temperature
 none
 baroc
 none
 SCOORD Use horizontally varying vertical coordinates
What to do with the radiation reaching the sea bed is still not established.
It is lost in current version, but can be returned to the bottom box.
Next: endstep.F
Up: The Proudman Oceanographic Laboratory
Previous: diffuse.F
The AMMP Project
20050420