[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