[MITgcm-support] heat flux at bottom
jmc at ocean.mit.edu
Mon Aug 18 11:01:16 EDT 2008
I don't see a bottom heat flux in the present code (but this
has been used for geothermal heating, probably in a customized
version). It means you will have to modify the code.
To add a bottom heat flux, I would probably use "mypackage"
and modify/add things there:
For instance, you can change SUBROUTINE MYPACKAGE_TENDENCY_APPLY_T
(in mypackage_tendency_apply.F), to something like
IF ( k.EQ.kLowC(i,j,bi,bj) THEN
gT(i,j,k,bi,bj) = gT(i,j,k,bi,bj)
where the IF is just testing if this level is the bottom level.
And you could even store bottomHeatFlux into one of myPacakge 2-D
array (e.g., myPa_Surf1).
On Mon, Aug 18, 2008 at 05:01:06PM +0300, Hila Afargan wrote:
> I am looking for a way to impose a certain heat flux at the bottom,
> rather than at the surface.
> What do you think will be the best way to try that?
> Is there a known way that will allow me to load heat flux forcing file
> (same as surfQfile)?
> Hilla Afargan
> Department of Environmental Sciences,
> Weizmann Institute of Science,
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
More information about the MITgcm-support