[MITgcm-support] (no subject)

prosperi at wicc.weizmann.ac.il prosperi at wicc.weizmann.ac.il
Sat Jun 19 04:41:21 EDT 2004


Hi guys,
I have a little problem:
I am trying to run a simplified version of the natl_box:

basically the forecing are a zonal wind, a linear meridional Qsurfaceflux 
on a 4000 meter 1000km*1000km box.
the data.pkg is

 useGMRedi=.TRUE.,
  useKPP=.TRUE.,
  useECCO=.TRUE.,

the Package.conf is

oceanic
cd_code
timeave

the data file is 

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
#
#   tRef            - Reference vertical potential temperature          
(deg C)
#   sRef            - Reference vertical salinity                         
(PSU)
#   viscAh          - Horizontal eddy viscosity coefficient             
(m^2/s)
#   viscAz          - Vertical eddy viscosity coefficient               
(m^2/s)
#   diffKhT         - Horizontal temperature diffusivity                
(m^2/s)
#   diffKzT         - Vertical temperature diffusivity                  
(m^2/s)
#   diffKhS         - Horizontal salt diffusivity                       
(m^2/s)
#   diffKzS         - Vertical salt diffusivity                         
(m^2/s)
#   f0              - Reference coriolis parameter,
#                     south edge of f on beta plane                       
(1/s)
#   beta            - df/dy                                         
(s^-1.m^-1)
#   tAlpha          - Linear EOS thermal expansion coefficient           
(1/oC)
#   sBeta           - Linear EOS haline contraction coefficient         
(1/ppt)
#   gravity         - Acceleration due to gravity                       
(m/s^2)
#   gBaro           - Accel. due to gravity used in barotropic equation 
(m/s^2)
#   rigidLid        - Set to true to use rigid lid
#   implicitFreeSurface - Set to true to use implicit free surface
#   eosType         - Flag for linear or polynomial equation of state
#   GMkbackground   - background value of GM/Redi coefficient
#   momAdvection    - On/Off flag for momentum self transport
#   momViscosity    - On/Off flag for momentum mixing
#
 &PARM01
 tRef               = 20*5.0,
 sRef               = 20*35.0,
 viscAz             = 1.E-5,
 diffKzT            = 1.E-6,
 diffKzS            = 1.E-6,
 viscAh             = 540.,
 viscA4             = 4.E8,
 diffKhT            = 20.,
 diffKhS            = 20.,
 f0                 = 1.0E-4,
 beta               = 2.E-11,
 tAlpha             = 0.28,
 sBeta              = 0.0,
 gravity            = 9.81,
 rhoConst           = 1025.0,
 rigidLid           = .FALSE.,
 implicitFreeSurface= .TRUE.,
 implicitDiffusion  = .TRUE.,
 useCDscheme        = .TRUE.,
 useNHMTerms        = .TRUE.,
 eosType            = 'LINEAR',
 momAdvection       = .TRUE.,
 momViscosity       = .TRUE.,
 implicitDiffusion  = .TRUE.,
 nonHydrostatic     = .TRUE.
 implicitViscosity  = .TRUE.,
 no_slip_bottom     = .TRUE.,
 no_slip_sides      = .FALSE.,
 readBinaryPrec     =  64.,
 &

# Elliptic solver parameters
#
#   cg2dMaxIters       - Maximum number of 2d solver iterations
#   cg2dTargetResidual - Solver target residual
#
 &PARM02
 cg2dMaxIters       = 1000,
 cg2dTargetResidual = 1.E-13,
 &

# Time stepping parameters
#
#   startTime         - Integration starting time                (s)
#   endTime           - Integration ending time                  (s)
#   tauCD             - CD scheme coupling timescale             (s)
#   deltaTMom         - Timestep for momemtum equations          (s)
#   deltaTtracer      - Tracer timestep                          (s)
#   deltaTClock       - Timestep used as model "clock"           (s)
#   abEps             - Adams-Bashforth stabilising factor
#   pChkPtFreq        - Frequency of permanent check pointing    (s)
#   chkPtFreq         - Frequency of rolling check pointing      (s)
#   dumpFreq          - Frequency at which model state is stored (s)
#   tauThetaClimRelax - Relaxation to climatology time scale     (s)
#   tauSaltClimRelax  - Relaxation to climatology time scale     (s)
#
 &PARM03
 startTime  = 0.0,
 endTime    = 2808000.,
 deltaTmom  = 2400.0,
 TauCD      = 172800.
 cAdjFreq   = 0.,
 abEps      = 0.1,
 pChkptFreq = 0.0,
 chkptFreq  = 0.0,
 dumpFreq   = 864000.0,
 monitorFreq=1.,
 &
=
# Gridding parameters
#
#   usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
#   usingCartesianGrid      - On/Off flag for selecting cartesian 
coordinates
#   delX                    - Zonal grid spacing         (degrees)
#   delY                    - Meridional grid spacing    (degrees)
#   delZ                    - Vertical grid spacing      (m)
#   phiMin                  - Southern boundary latitude (degrees)
#
 &PARM04
 usingCartesianGrid      = .TRUE.,
 usingSphericalPolarGrid = .FALSE.,
 delX   = 20*50.E3,
 delY   = 20*50.E3,
 delZ   = 30*133
 &

# Input datasets
#
#   bathyFile       - File containing bathymetry
#   hydrogThetaFile - File containing initial potential temperature data
#   hydrogSaltFile  - File containing initial salinity data
#   zonalWindFile   - File containing zonal wind data
#   meridWindFile   - File containing meridional wind data
#   thetaClimFile   - File containing theta climatology used for 
relaxation
#   saltClimFile    - File containing salt climatology used for relaxation
#
 &PARM05
 bathyFile        = 'topog.box',
 zonalWindFile    = 'windx.sin_y',
 surfQFile        = 'Qsurf.bin',

 &
 and the SIZE.h

  &           sNx =  20,
     &           sNy =  20,
     &           OLx =   3,
     &           OLy =   3,
     &           nSx =   1,
     &           nSy =   1,
     &           nPx =   1,
     &           nPy =   1,
     &           Nx  = sNx*nSx*nPx,
     &           Ny  = sNy*nSy*nPy,
     &           Nr  =  30)

the 
CPP_EEOPTIONS.h
CPP_OPTIONS.h 
are the defaults ones.

In the compilation everything seems ok but I got thi error message

forrtl: severe (24): end-of-file during read, unit 11, file 
/tmp/fortEAAPcOYNM
   0: __FINI_00_remove_gp_range [0x3ff81a6de38]
   1: __FINI_00_remove_gp_range [0x3ff81a6e45c]
   2: __FINI_00_remove_gp_range [0x3ff81aa6e14]
   3: packages_boot_ [packages_boot.f: 1409, 0x120142e00]
   4: initialise_fixed_ [initialise_fixed.f: 1485, 0x1201383d4]
   5: the_model_main_ [the_model_main.f: 1989, 0x12014cb44]
   6: main_ [main.f: 3749, 0x1200d6bf0]
   7: main [for_main.c: 203, 0x12015249c]
   8: __start [0x12001abb8]

the output.txt finish in this way

 PACKAGES_BOOT: opening data.pkg
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.pkg
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) // Parameter file "data.pkg"
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) ># Packages
(PID.TID 0000.0001) > &PACKAGES
(PID.TID 0000.0001) >  useGMRedi=.TRUE.,
(PID.TID 0000.0001) >  useKPP=.TRUE.,
(PID.TID 0000.0001) >  useECCO=.TRUE.,
(PID.TID 0000.0001) >&

so seems to me that I have a problem with the packages dependence.

How can I figure out what's going on?

Thank you very much in advance 

Davide






More information about the MITgcm-support mailing list