[MITgcm-support] Layers package : nothing but zeros for output?

[Dan Jones]

Hi all,

I am interested in implementing the layers package in an idealized model.
 I can get it to compile, run, and produce output files such as Hs, Hw, UH,
and VH.  However, these output files contain nothing but zeros!  I'm very
likely making some trivial error somewhere - please take a look at my setup
below:

-------------------- /code/LAYERS_SIZE.h
-------------------------------------------------------------------------------

C $Header: /u/gcmpack/MITgcm/pkg/layers/LAYERS_SIZE.h,v 1.3 2013/01/08
21:52:34 jmc Exp $
C $Name: checkpoint64t $
C * Compiled-in size options for the LAYERS package *
C
C  - Just as you have to define Nr in SIZE.h, you must define the number
C    of vertical layers for isopycnal averaging so that the proper array
C    sizes can be declared in the LAYERS.h header file.
C
C  - Variables -
C      NLayers        :: the number of isopycnal layers (must match
data.layers)
C      FineGridFact   :: how many fine-grid cells per dF cell
C      FineGridMax    :: the number of points in the finer vertical grid
C                        used for interpolation
C      layers_maxNum  :: max number of tracer fields used for layer
averaging
      INTEGER    Nlayers, FineGridFact, FineGridMax, layers_maxNum
      PARAMETER( Nlayers = 10 )
      PARAMETER( FineGridFact = 10 )
      PARAMETER( FineGridMax = Nr * FineGridFact )
      PARAMETER( layers_maxNum = 1 )

-------------------- /input/data.layers
-------------------------------------------------------------------------------

 &LAYERS_PARM01
# averaging
  layers_taveFreq=45000.,
# no instantaneous output
  layers_diagFreq=0.,
  layers_G=-0.2,0.6,1.5,2.3,3.2,4.0,4.8,5.7,6.5,7.4,8.2,
 &
# Note: Some systems use & as the
# namelist terminator. Other systems
# use a / character.

As you can see in the header file, I'm running MITgcm c64t.  Based on the
thread linked below, I assume that I'm using the default layer coordinate
(i.e. theta), so as long as I'm picking reasonable values of theta I should
have well-defined layers.  Is that correct?

http://mitgcm.org/pipermail/mitgcm-support/2011-December/007569.html

Thanks in advance for taking a look!
Dan

---------------------------------------------------------

Dr Dan Jones
British Antarctic Survey, NERC
Cambridge, UK

Phone: +44 (0)1223 221505
Fax: +44 (0)1223 362616
Skype:  dcjones.work

---------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140611/7604f85e/attachment.htm>