[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