[MITgcm-support] Interior ocean heating

[Martin Losch]

Christopher,

as a quick hack you can easily modify the routine external_forcing_t.  
in external_forcing.F.
there you add you bottom heat flux (hard coded) as a tendency to the  
field gT. E.g. if you have some flux Qbot ( in W/m^2 ) you can add  
code such as
if ( klev = kLowC(i,j,bi,bj) ) then
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) + Qbot*recip_drF(klev) 
*recip_hfac(i,j,klev,bi,bj)*recip_Cp*horiVertRatio*recip_rhoConst
endif
or so. If you want to read the Qbot from a file, it's going to be a  
bit more involved.

Martin

On 7 Mar 2007, at 02:30, Christopher L. Wolfe wrote:

>
> Dear Modelers,
>
> Is there any way to add an internal heat source to the MITgcm in  
> ocean mode? What I really want to do is flux heat through the ocean  
> floor, but adding it as a source at the lowest grid point would  
> accomplish the same thing. The only way I can see to do this given  
> what I've been able to figure out so far is to fake it by relaxing  
> the model toward a hydrography which is identical to the model mean  
> state except that it's warmer in the bottom layer. However, I think  
> this would tend to suppress the formation of transient eddies,  
> which I don't really want to do.
>
> Thanks in advance,
> Christopher
>
> -----------------------------------------------------------
> Dr. Christopher L. Wolfe             		 858-534-4560
> Physical Oceanography Research Division    OAR 357
> Scripps Institution of Oceanography, UCSD  clwolfe at ucsd.edu
> -----------------------------------------------------------
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support