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

[Hazel, Julia]

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