[MITgcm-support] (no subject)
Patrick Heimbach
heimbach at MIT.EDU
Mon Jun 21 04:24:36 EDT 2004
Here's a second (a bit better) try
(the last one wasn't very good).
>From your packages.conf I take it that you have
ECCO disabled (it does not appear in your list)?
the placeholder 'oceanic' only enables
gfd gmredi kpp timeave
But in your data.pkg it is set to .TRUE.
Probably what's happening is that in packages_boot
some characters are tried to be written to an
'internal file', i.e. variable in order to generate
an error message and there is an inconsistency.
Can you just try running with useECCO = .FALSE. ?
-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