[MITgcm-support] Sea ice dynamics causing model to output NaN's after changing grid configuration
Mari Fjalstad Jensen
Mari.F.Jensen at uib.no
Thu Nov 21 11:58:38 EST 2019
Dear Martin and David (or others)
Did you ever solve the problem below? I am having the same problem,
NaN's after on time-step when sea ice dynamics are enabled. The
simulation runs fine when sea ice dynamics are turned off, or if there
is no wind forcing.
In addition, I recently switched from a cartesian domain with smooth
simple bathymetry (flat bottom + sloping sides) to a spherical polar-
realistic bathymetry domain. In the previous setup, the sea ice package
was working fine in the presence of wind and sea ice dynamics.
Attached are data* and *OPTIONS files
Cheers and thanks,
Mari
On 24/08/2019 08.55, Martin Losch wrote:
> Hi David,
> I would need more details to be able to help. Configuration (code directory, namelist files, version of model), and what exactly you changed from a working configuration
>
> Martin
>
>> On 24. Aug 2019, at 00:07, David Vishny <davidvish at gmail.com> wrote:
>>
>> To whom it may concern,
>>
>> I recently changed the grid configuration in my model setup, and now the model outputs all NaN’s after a single time-step if sea ice dynamics are enabled, even though there is no sea ice to begin with. If I disable sea ice dynamics, it appears that no problems occur. Thus, my advisor Malte Jansen and I believe the LSR solver is producing NaN’s.
>>
>> I know one parameter related to the LSR solver is the LSR error. For any given grid configuration, could the solver crash due to an LSR error that is either too high or too low? Are there any other parameters I should be playing around with that are related to the LSR solver?
>>
>>
>> Thanks,
>>
>> David Vishny
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
# SEAICE parameters
&SEAICE_PARM01
SEAICEwriteState = .TRUE.,
SEAICEadvScheme = 33,
SEAICE_initialHEFF = 0.0,
SEAICE_deltaTtherm = 200.,
SEAICE_deltaTdyn = 200.,
# SEAICE_deltaTevp = 20.,
SEAICE_multDim = 1,
# SEAICE_tempFrz0 = -1.96,
# SEAICE_dTempFrz_dS - 0.0,
# SEAICEuseDYNAMICS =.FALSE.,
# LSR_ERROR = 1.E-12,
# SEAICE_tave_mnc = .FALSE.,
# SEAICE_dump_mnc = .FALSE.,
# SEAICE_mon_mnc = .FALSE.,
ICE2WATR = 0.92D0,
# for backward compatibility only, in realistic runs
# this should always be turned on
SEAICEadvSnow = .TRUE.,
SEAICEuseFlooding = .TRUE.,
&
-------------- next part --------------
# Packages
&PACKAGES
usePTRACERS=.FALSE.
useEXF=.TRUE.
useSEAICE=.TRUE.
usediagnostics=.TRUE.
useKPP=.FALSE.
&
-------------- next part --------------
# $Header: /u/gcmpack/MITgcm/verification/lab_sea/input/data.exf,v 1.11 2007/04/18 19:57:21 heimbach Exp $
#
# *********************
# External Forcing Data
# *********************
&EXF_NML_01
#
useExfCheckRange = .FALSE.,
repeatPeriod = 3.1104e+7,
exf_iprec = 64,
exf_yftype = 'RL',
exf_monFreq = 2592000,
#
&
# *********************
&EXF_NML_02
#
hfluxstartdate1 = 19781216,
hfluxstartdate2 = 180000,
hfluxperiod = 3.1104e+7,
#
sfluxstartdate1 = 19781216,
sfluxstartdate2 = 180000,
sfluxperiod = 3.1104e+7,
#
ustressstartdate1 = 19781216,
ustressstartdate2 = 180000,
ustressperiod = 3.1104e+7,
#
vstressstartdate1 = 19781216,
vstressstartdate2 = 180000,
vstressperiod = 3.1104e+7,
#
atempstartdate1 = 19781216,
atempstartdate2 = 180000,
atempperiod = 3.1104e+7,
#
aqhstartdate1 = 19781216,
aqhstartdate2 = 180000,
aqhperiod = 3.1104e+7,
#
#evapstartdate1 = 19781216,
#evapstartdate2 = 180000,
#evapperiod = 2635200.0,
#
precipstartdate1 = 19781216,
precipstartdate2 = 180000,
precipperiod = 3.1104e+7,
#
uwindstartdate1 = 19781216,
uwindstartdate2 = 180000,
uwindperiod = 3.1104e+7,
#
vwindstartdate1 = 19781216,
vwindstartdate2 = 180000,
vwindperiod = 3.1104e+7,
#
swfluxstartdate1 = 19781216,
swfluxstartdate2 = 180000,
swfluxperiod = 3.1104e+7,
#
lwfluxstartdate1 = 19781216,
lwfluxstartdate2 = 180000,
lwfluxperiod = 3.1104e+7,
#
swdownstartdate1 = 19781216,
swdownstartdate2 = 180000,
swdownperiod = 3.1104e+7,
#
lwdownstartdate1 = 19781216,
lwdownstartdate2 = 180000,
lwdownperiod = 3.1104e+7,
#
climsststartdate1 = 19781216,
climsststartdate2 = 180000,
climsstperiod = 3.1104e+7,
#
climsssstartdate1 = 19781216,
climsssstartdate2 = 180000,
climsssperiod = 3.1104e+7,
#
hfluxfile = ' ',
sfluxfile = ' ',
ustressfile = ' ',
vstressfile = ' ',
atempfile = 'Tair.bin',
aqhfile = ' ',
uwindfile = 'Uwind.bin',
vwindfile = 'Vwind.bin',
# evapfile = 'Eva.bin',
precipfile = 'EmP.bin',
lwfluxfile = ' ',
swfluxfile = ' ',
lwdownfile = ' ',
swdownfile = ' ',
runoffFile = ' '
climsstfile = ' ',
climsssfile = ' ',
#
&
# *********************
&EXF_NML_03
&
# *********************
&EXF_NML_04
&
-------------- next part --------------
# Model parameters
# Continuous equation parameters
&PARM01
tRef=10.5,10.4,10.3,10.1,9.9,9.5,9.3,9.2,9.0,8.8,8.7,8.5,8.3,8.1,7.9,7.7,7.4,7.2,6.9,6.6,6.4,5.4,4.5,3.9,3.6,3.3,3.1,2.8,2.7,2.7,
sRef=30*35.0,
viscAz=1.E-4,
# viscAh=20.,
# viscA4=1.e+9,
viscC2Smag=2.5,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
# diffK4T=1.e+9,
# diffKhT=30.,
diffKzT=1.E-4,
# diffK4S=1.e+9,
# diffKhS=30.,
diffKzS=1.E-4,
tAlpha=0.2E-3,
sBeta =0.8E-3,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
staggerTimeStep=.TRUE.,
eosType='JMD95Z',
# momAdvection=.FALSE.,
readBinaryPrec=64,
hFacMin=0.1,
bottomDragQuadratic=1.e-3,
implicitDiffusion=.TRUE.,
implicitViscosity=.TRUE.,
tempAdvScheme=33,
saltAdvScheme=33,
useSingleCpuIO=.TRUE.,
# vectorInvariantMomentum=.TRUE.,
ivdc_kappa=1000.0,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-11,
&
# Time stepping parameters
&PARM03
tauCD=0.,
nIter0=0,
nTimeSteps=20412000,
# startTime=0.0,
# startTime=235872000,
# endTime=11.664e+8,
# endTime=3.1104e+7,
deltaT=200.0,
abEps=0.1,
pChkptFreq=0.0,
chkptFreq=2592000,
dumpFreq=2592000,
monitorFreq=86400,
cAdjFreq=0,
forcing_In_AB=.FALSE.,
tauThetaClimRelax=0.0,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.FALSE.,
usingSphericalPolarGrid=.TRUE.,
xgOrigin=-32.9167,
ygOrigin=58.0333,
dXspacing=.0833,
dYspacing=0.0333,
delZ=5,10,10,10,25,15*50,4*200,6*400,
&
&PARM05
checkIniTemp=.false.,
bathyFile='topog.bin',
# ThetaClimFile='SST.bin',
# zonalWindFile='taux.bin',
# EmPmRfile='EmP.bin',
&
-------------- next part --------------
A non-text attachment was scrubbed...
Name: EXF_OPTIONS.h
Type: text/x-chdr
Size: 7988 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20191121/f52bd3a5/attachment-0002.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CPP_OPTIONS.h
Type: text/x-chdr
Size: 5038 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20191121/f52bd3a5/attachment-0003.bin>
More information about the MITgcm-support
mailing list