[MITgcm-support] (no subject)

Patrick Heimbach heimbach at MIT.EDU
Sat Jun 19 15:50:39 EDT 2004


Hi,

I haven't looked into the natl_box/ verif. in a long time,
only the natl_box_adjoint/

Will look into it later in the coming week.

One very quick guess.
>From the std output it seems you've modified data.pkg
and hve the namelist terminator '&' in the first column.
Could you try having a blank in front?

Cheers
-Patrick



Quoting prosperi at wicc.weizmann.ac.il:

> 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
> 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> 


--------------------------------------------------------
Patrick Heimbach   Massachusetts Institute of Technology
FON: +1/617/253-5259                  EAPS, Room 54-1518
FAX: +1/617/253-4464             77 Massachusetts Avenue
mailto:heimbach at mit.edu               Cambridge MA 02139
http://www.mit.edu/~heimbach/                        USA




More information about the MITgcm-support mailing list