[MITgcm-support] MOM_IMPLICIT_R: error when solving 3-Diag problem.

[Martin Losch]

Hi Julia,

The seaice_waterdrag value implicitly used to include the density of water, hence the factor of 1000 difference. I am sorry that this cause confusion. We decided not to add a stop but rather issue a warning ...

here’s a list of new defaults in pkg/seaice that changed a long with seaice_waterdrag defaults
https://github.com/MITgcm/MITgcm/issues/107

Martin

> On 30. Nov 2018, at 15:53, Hazel, Julia <jhazel at atmos.ucla.edu> wrote:
> 
> Yes, thanks - the model is blowing up.  It had been running fine on an older version of mitgcm, but I realize there are incompatibilities with my old data.seaice and the newer sea ice package.  (namely in the SEAICE_waterDrag value, which is now three orders of magnitude smaller than before)
> 
> 
> 
>> On Nov 29, 2018, at 9:26 PM, Matthew Mazloff <mmazloff at ucsd.edu> wrote:
>> 
>> That could do it
>> 
>> or the model is just blowing up, because I have seen a similar error when the just model crashes…
>> What is the tail of STDOUT.0000? Is your cfl number ~1?
>> 
>> -Matt
>> 
>> 
>>> On Nov 29, 2018, at 9:10 PM, Patrick Heimbach <heimbach at mit.edu> wrote:
>>> 
>>> Hello Julia,
>>> 
>>> your parameter file "data" has a typo, see line
>>> useJamartMomAdv=.TRUE.,:
>>> 
>>> That namelist error likely causes the crash.
>>> 
>>> -Patrick
>>> 
>>>> On Nov 29, 2018, at 9:51 PM, Hazel, Julia <jhazel at atmos.ucla.edu> wrote:
>>>> 
>>>> Hello- I’m trying to run on stampede2 and getting the following error after the run has executed for ~ 30 minutes.
>>>> I see the error is located in the mom_u_implicit_r.f file, but I can’t determine how to fix it.
>>>> 
>>>> TACC:  Starting up job 2532480 
>>>> TACC:  Setting up parallel environment for MVAPICH2+mpispawn. 
>>>> TACC:  Starting parallel tasks... 
>>>> MOM_IMPLICIT_R: error when solving 3-Diag problem.
>>>> [c402-031.stampede2.tacc.utexas.edu:mpispawn_1][readline] Unexpected End-Of-File on file descriptor 45. MPI process died?
>>>> [c402-031.stampede2.tacc.utexas.edu:mpispawn_1][mtpmi_processops] Error while reading PMI socket. MPI process died?
>>>> 
>>>> 
>>>> My data file is as follows— thank you!
>>>> 
>>>> # Continuous equation parameters
>>>> &PARM01
>>>> vectorInvariantMomentum=.TRUE.,
>>>> viscAr=0.00000000e+00,
>>>> viscA4=0.00000000e+00,
>>>> viscAh=0.00000000e+00,
>>>> viscA4Grid=1.00000000e-01,
>>>> viscAhGrid=0.00000000e+00,
>>>> viscA4GridMax=5.00000000e-01,
>>>> viscAhGridMax=1.00000000e+00,
>>>> useAreaViscLength=.FALSE.,
>>>> useFullLeith=.TRUE.,
>>>> viscC4leith=0.00000000e+00,
>>>> viscC4leithD=0.00000000e+00,
>>>> viscC2leith=0.00000000e+00,
>>>> viscC2leithD=0.00000000e+00,
>>>> diffKrT=0.00000000e+00,
>>>> diffKhT=0.00000000e+00,
>>>> diffK4T=0.00000000e+00,
>>>> tempAdvScheme=7,
>>>> saltAdvScheme=7,
>>>> multiDimAdvection=.TRUE.,
>>>> tempStepping=.TRUE.,
>>>> saltStepping=.TRUE.,
>>>> staggerTimeStep=.TRUE.,
>>>> eosType='JMD95Z',
>>>> no_slip_sides=.FALSE.,
>>>> no_slip_bottom=.FALSE.,
>>>> bottomDragLinear=0.00000000e+00,
>>>> bottomDragQuadratic=2.10000000e-03,
>>>> gravity=9.81000000e+00,
>>>> rhonil=1.02750000e+03,
>>>> rhoConst=1.02750000e+03,
>>>> ivdc_kappa=0.00000000e+00,
>>>> implicitDiffusion=.TRUE.,
>>>> implicitViscosity=.TRUE.,
>>>> exactConserv=.TRUE.,
>>>> useCDscheme=.FALSE.,
>>>> hFacMin=1.00000000e-01,
>>>> hFacMinDr=0.00000000e+00,
>>>> readBinaryPrec=64,
>>>> useSingleCpuIO=.TRUE.,
>>>> implicitDiffusion=.TRUE.,
>>>> implicitViscosity=.TRUE.,
>>>> debugLevel=-1,
>>>> vectorInvariantMomentum=.FALSE.,
>>>> useJamartWetPoints=.TRUE.,
>>>> useJamartMomAdv=.TRUE.,:
>>>> 
>>>> &PARM02
>>>> useSRCGSolver=.TRUE.,
>>>> cg2dMaxIters=1000,
>>>> cg2dTargetResidual=1.00000000e-12,
>>>> &
>>>> 
>>>> # Time stepping parameters
>>>> &PARM03
>>>> momDissip_In_AB=.FALSE.,
>>>> tracForcingOutAB=1,
>>>> nIter0=0,
>>>> abEps=1.00000000e-01,
>>>> chkptFreq=2.59200000e+06,
>>>> pChkptFreq=2.59200000e+06,
>>>> taveFreq=0.00000000e+00,
>>>> dumpFreq=0.00000000e+00,
>>>> monitorFreq=3.16224000e+07,
>>>> cAdjFreq=0.00000000e+00,
>>>> dumpInitAndLast=.TRUE.,
>>>> pickupStrictlyMatch=.FALSE.,
>>>> endTime=63244800,
>>>> deltaT=3.00000000e+02,
>>>> &
>>>> 
>>>> _______________________________________________
>>>> 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
>> 
>> _______________________________________________
>> 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