[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