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

Hazel, Julia jhazel at atmos.ucla.edu
Fri Nov 30 11:25:30 EST 2018


Great, thanks for all the help! That makes more sense since the air-ocean and air-ice drags include density.


> On Nov 30, 2018, at 7:56 AM, Martin Losch <Martin.Losch at awi.de> wrote:
> 
> 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
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support



More information about the MITgcm-support mailing list