[MITgcm-support] adding convection counter to netcdf output file

[Ed Hill]

On Sun, 2005-11-27 at 23:49 -0500, Eli Tziperman wrote:
> Hi all, I am a beginner in this, and apologize if this is a silly
> question.  I would like to be able to plot the locations/ times where
> convection occurs.  I followed the examples in the manual in order to
> add the variable ConvectCount to ini_mnc_vars.F and write_state.F.
> However, I am getting a complaint during compilation that this
> variable does not have a type.  Is the variable only defined for
> certain convection options?  Other possible problems?  many thanks...

Hi Eli,

No, I don't think thats a silly question at all!  If you look in the
file:

  MITgcm/model/src/convective_adjustment.F

ConvectCount is a local variable.  This means that its not defined in
one of the common blocks (header files) so you don't have access to it
within the write_state() routine.  That probably explains why you're
seeing the error about a variable not having a type (its undefined).  So
one way you can fix that error is to move the ConvectCount variable into
one of the common blocks.  Or, you could create a new common block that
contains it and is only defined where you need it.

Also, looking in the convective_adjustment() routine, there appears to
be a diagnostic defined called 'CONVADJ ' that may be accumulating the
values you want.  Perhaps the easiest solution to your problem would be
to turn on the diagnostics package and then enable the 'CONVADJ '
diagnostic with a time scale appropriate for your setup.

Ed

-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS;  Rm 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3 at mit.edu                ed at eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464