[MITgcm-support] balancing the heat equation

Jean-Michel Campin jmc at ocean.mit.edu
Fri Sep 7 17:25:03 EDT 2007


Hi, 

Just for this precision:
> Or are there some hidden diagnostics not written to this file ?
> david
All the diagnostics which are available in this configuration 
are listed, no exception.

Jean-Michel

On Fri, Sep 07, 2007 at 09:38:04AM -0400, David Ferreira wrote:
> Hi,
> If you turn on the diagnostics package, and run a couple of
> of time steps, the model should spit out a file called 
> "available_diagnostics.log" which, I think, has the list of all
> available dignostics for your particular configuration.
> Or are there some hidden diagnostics not written to this file ?
> david
> 
> 
> 
> Martin Losch wrote:
> >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
> >
> >
> >_______________________________________________
> >MITgcm-support mailing list
> >MITgcm-support at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support



More information about the MITgcm-support mailing list