[MITgcm-support] balancing the heat equation
Martin Losch
Martin.Losch at awi.de
Fri Sep 7 03:13:10 EDT 2007
Hi Nicolas/Allison,
the safest way to find out what diagnostics are available is to have
a look in the respective "diagnostics_init" files. The are usually:
pkg/diagnostics/diagsnostics_init_main.F (for general stuff, e.g.
TOTTTEND)
pkg/$pkg/$pkg_diagnostics_init.F
so for all the advective and diffusive stuff you'd have to look in
pkg/generic_advdiff/gad_diagnostics_init.F (which is already a
counter example to the about scheme, oh well).
Anyway, for checking the balance offline you'll find everything in
these two files except for the non-local transport of KPP. Since you
have all the other terms, you could assume that the residual is the
non-local transport of KPP. If you don't trust the MITgcm (shame on
you (-:) then you'll have to do the modifications suggested earlier
in this thread.
Martin
On 6 Sep 2007, at 22:59, Nicolas Wienders wrote:
>
>> Hi Allison,
>> without checking details in the code I would suggest that you use
>> the diagnostics package for diagnosing the different
>> contributions. As far as I can see, the diagnostics are all
>> available:
>> dT/dt:
>> diagName = 'TOTTTEND'
>> diagTitle = 'Tendency of Potential Temperature'
>> uTx/vTy/wTz: 'ADVx/y/r_TH '
>> KTxx/yy/zz: 'DFx/y/rE_TH ' (and 'DFrI_TH ' for the implicit part
>> which is relevant for KPP)
>
>
> I didn't see the variable 'TOTTTEND' in the manual (pages 376-385)
> but could find it in the code.
>
> - Are some available variables missing in the package documentation
> (manual) ?
>
> - Has anyone done a complete Heat equation balance and in that case
> what diagnostics variables
> where you outputting?
>
> Thank you.
>
>
>
>
>
>
> --
>
> Dr. Nicolas Wienders (850) 644-1987 (office)
> Florida State University 294-3542 (cell)
> Oceanography - OSB415 644-2581 (fax)
> P.O. Box 3064320 117 N. Woodward Avenue
> Tallahassee, FL 32306-4320 wienders at ocean.fsu.edu
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list