[MITgcm-support] sea ice experiment
Amitabh Mitra
amitabhmitra1 at rediffmail.com
Wed Oct 31 04:54:33 EDT 2007
Sir,
I am trying to run a sea ice experiment for
the Antarctic region -70 to -40 South and
-178 to 178 East with 1x1 resolution. The
corresponding 'data' file is attached.
It runs for 10 days but when I extend it
to 20 days the model blows with the following error:
SGI_2100 11% mpirun -np 1 ./mitgcmuv
**WARNNING** CONFIG_CHECK: useMNC is TRUE and #undef
ALLOW_MNC
STOP NORMAL END
STOP
SGI_2100 12% mpirun -np 1 ./mitgcmuv
**WARNNING** CONFIG_CHECK: useMNC is TRUE and #undef
ALLOW_MNC
SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=
-2.478E+03 1.998E+03
MON_SOLUTION: STOPPING CALCULATION
S/R EEDIE: Only 0 threads have completed,
1 are expected for this configuration!
Possibly you have different setenv PARALLEL and
nThreads?
STOP MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF
SOLUTION
STOP
SGI_2100 13%
------------------------------------------------------
Please suggest suitable values for these coefficients
viscAz, viscAh, diffKhT, diffKzT, diffKhS, diffKzS
and deltaT so that the model runs efficiently.
I have read the manual but it is not helpful
in determining these values. Do you have any method
to calculate these values or it is just trial and
error method ?
Thanking you,
With Regards,
Amitabh Mitra
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20071031/8495abcf/attachment.htm>
-------------- next part --------------
# ====================
# | 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
# momAdvection - On/Off flag for momentum self transport
# momViscosity - On/Off flag for momentum mixing
#
&PARM01
tRef= 5.18, 5.13, 5.08, 5.04, 4.86,
4.62, 4.4, 4.29, 4.25, 4.3,
4.32, 4.35, 4.17, 4.0, 3.85,
sRef= 34.13, 34.14, 34.14, 34.17, 34.2,
34.23, 34.27, 34.31, 34.34, 34.4,
34.44, 34.46, 34.47, 34.48, 34.48,
no_slip_sides=.FALSE.,
no_slip_bottom=.TRUE.,
viscAz=8.33E-3,
viscAh=6.51E4,
diffKhT=39.0625,
diffKzT=7.5E-5,
diffKhS=39.0625,
diffKzS=7.5E-5,
beta=1.E-11,
tAlpha=2.E-4,
sBeta =7.4E-4,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
eosType='JMD95Z',
ivdc_kappa=0.,
implicitDiffusion=.TRUE.,
useOldFreezing=.TRUE.,
useRealFreshWaterFlux=.TRUE.,
useCDscheme=.TRUE.,
useNHMTerms=.TRUE.,
# turn on looped cells
hFacMin=.05,
hFacMindz=50.,
# set precision of data files
readBinaryPrec=32,
writeBinaryPrec=32,
saltStepping=.TRUE.,
tempStepping=.TRUE.,
momStepping=.TRUE.,
# implicitDiffusion=.TRUE.,
implicitViscosity=.TRUE.,
#globalFiles=.TRUE.,
#- not safe to use globalFiles in multi-processors runs; set instead useSingleCpuIO
useSingleCpuIO=.TRUE.,
allowFreezing=.FALSE.,
HeatCapacity_Cp = 3986.D0
gravity = 9.8156,
gBaro = 9.8156,
rhoNil = 1027.D0
# useCDscheme=.TRUE.,
# useNHMTerms=.TRUE.,
&
# Elliptic solver parameters
#
# cg2dMaxIters - Maximum number of 2d solver iterations
# cg2dTargetResidual - Solver target residual
#
&PARM02
cg2dMaxIters=500,
cg2dTargetResidual=1.D-8,
&
# 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
tauCD=172800.,
startTime=0.0,
endTime=864000.,
deltaTmom=300.0,
deltaTtracer=1800.0,
# deltaT=1200.0,
deltaTClock =1800.0,
cAdjFreq=0.,
abEps=0.1,
pChkptFreq=864000.,
chkptFreq= 864000.,
dumpFreq = 864000.,
tavefreq = 864000.,
# 6 months restoring timescale for salinity
tauSaltClimRelax = 15552000.0,
# 2 months restoring timescale for temperature
tauThetaClimRelax = 5184000.0,
monitorFreq=1.,
periodicExternalForcing=.TRUE.,
externForcingPeriod=2592000.,
externForcingCycle=31104000.,
forcing_In_AB = .FALSE.
&
# 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=.FALSE.,
usingSphericalPolarGrid=.TRUE.,
delX=356*1.E0,
delY=30*1.E0,
delZ=50., 70., 100., 140., 190.,
240., 290., 340., 390., 440.,
490., 540., 590., 640., 690.,
phimin=-70.,
thetamin=-178.,
rSphere = 6371.D3
&
# 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 = 'bathymetry.bin',
hydrogThetaFile = 'lev_monthly_ss_am_temp.bin',
hydrogSaltFile = 'lev_monthly_ss_am_salt.bin',
zonalWindFile = 'ncep_taux.bin',
meridWindFile = 'ncep_tauy.bin',
thetaClimFile = 'lev_monthly_sstemp.bin',
saltClimFile = 'lev_monthly_sssalt.bin',
surfQFile = 'ncep_qnet.bin',
## fresh water flux is turned off, uncomment next line to turn on
# (not recommend together with surface salinity restoring)
EmPmRFile = 'ncep_emp.bin',
&
More information about the MITgcm-support
mailing list