[MITgcm-support] heat flux at bottom

[Jean-Michel Campin]

Hi Hilla,

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)
                   + bottomHeatFlux(i,j,bi,bj)*recip_Cp*mass2rUnit
                    *recip_drF(k)*recip_hFacC(i,j,k,bi,bj)
                    *recip_rA(i,j,bi,bj)
  ENDIF
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).
   
Jean-Michel

On Mon, Aug 18, 2008 at 05:01:06PM +0300, Hila Afargan wrote:
> Hi,
>
> 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)?
>
> Thanks,
> Hilla
>
>
> -- 
> Hilla Afargan
> Department of Environmental Sciences,
> Weizmann Institute of Science,
> Rehovot,76100
> Israel
> Tel:+972-8-9344924
> Fax:+972-8-9344124
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support