<div dir="ltr"><div>Hi All,</div><div><br></div><div>I have a model domain at 7km resolution in north Indian ocean with a southern open boundary. I get an error when the simulation begins.<br></div><div>Error: 'MOM_IMPLICIT_R: error when solving 3-Diag problem.'</div><div>Upon further investigation, the error message changes to: 'STOP ABNORMAL END: S/R MON_SOLUTION, stops due to EXTREME Pot.Temp' ( I did check and found extreme +ve and -ve values for temperature and salinity at deeper depths)<br></div><div><br></div><div>From previous support resolutions, I know it has to do with CFL criteria, and I have tried changing deltaT, viscAh/viscAhgrid, viscAz, initial conditions files and even delZ levels. Somehow I am not able to figure it out.</div><div><br></div><div>I have a second question. In<strong> </strong>calc_oce_mxlayer.F, if I am not mistaken, model will calc 'hMixCriteria<span></span>' first and based on its value, will decide to go for method 1 or 2. If so, how can I set the 2nd method for MLD calculation?</div><div><br></div><div><br></div><div># ====================<br># | Model parameters |<br># ====================<br>#<br># Continuous equation parameters<br> &PARM01<br><br>## viscAhgrid=9.8E-7,<br> viscAh=10.E3,<br> viscAz=1.E-4,<br>## viscAz=1.E-7,<br>## viscA4=0.,<br> diffKhT=1.,<br> diffKzT=1.E-5,<br>## diffK4T=0.,<br> diffKhS=1.,<br> diffKzS=1.E-5,<br>## diffK4S=0.,<br><br>## useFullLeith=.TRUE.,<br>## This uses the Strain-Tension form for calculating viscous terms<br>## useStrainTensionVisc=.TRUE.,<br><br>## This is the standard Leith Viscosity coefficient, <br>## visc=(visc2leith/pi)**3 grad(curl v_h)<br>## viscC2leith=2.0, <br><br>## This is an additional Leith Viscosity, prop to grad(div.v_h)<br>## It traps instabilities that have no signal in grad(curl v_h)<br>## viscC2leithD=2.0,<br><br>## This is the max Viscosity used, be it Smag or Leith, and it scales <br>## with grid size and timestep<br>## viscAhGridMax=1.E-2,<br>## viscA4GridMax=1.E-9,<br>##<br> no_slip_sides=.FALSE.,<br> no_slip_bottom=.TRUE.,<br> tempAdvScheme=33,<br> saltAdvScheme=33,<br> staggerTimeStep=.TRUE.,<br>#<br> tAlpha=2.E-4,<br> sBeta =7.E-4,<br> gravity=9.81,<br> f0=0.0000101456,<br> beta=2.2209E-11,<br> rigidLid=.FALSE.,<br> implicitFreeSurface=.TRUE.,<br> implicitDiffusion=.TRUE.,<br> implicitViscosity=.TRUE.,<br> nonHydrostatic=.FALSE.,<br> bottomDragLinear=0.0002,<br> eosType='MDJWF',<br> useRealFreshWaterFlux=.TRUE.,<br> exactConserv=.TRUE.,<br>## nonlinFreeSurf=4,<br>## select_rStar=2,<br> hFacMin=0.5,<br># hFacMinDz=5.,<br> hFacInf=0.2,<br> hFacSup=2.0,<br> useSingleCpuIO=.TRUE<br> readBinaryPrec=64,<br> writeBinaryPrec=64,<br> &<br><br># Elliptic solver parameters<br> &PARM02<br> cg2dMaxIters=1000,<br> cg2dTargetResidual=1.E-13,<br> cg3dMaxIters=400,<br> cg3dTargetResidual=1.E-13,<br>## cg2dUseMinResSol=1,<br> &<br><br># Time stepping parameters<br> &PARM03<br> startTime= 0.,<br> endTime= 86400.,<br> deltaTmom= 180.0,<br> deltaTtracer   = 3600.0,<br> deltaTClock    = 3600.0,<br> deltaTfreesurf = 3600.0,<br> pChkptFreq= 86400.,<br> dumpFreq= 7200.,<br># monitorFreq= 180.0,<br> &<br><br># Gridding parameters<br> &PARM04<br> usingCurvilinearGrid=.FALSE.,<br> usingCartesianGrid=.TRUE.,<br> usingSphericalPolarGrid=.FALSE.,<br> dxSpacing=7695.06,<br> dySpacing=7045.008,<br> delZ=    5.034 , 10.067 , 10.118 , 10.139 , 10.218 , 10.277 , 10.409 , 10.541 , 10.774 , 11.046,<br>       11.473 , 12.010 , 12.802 , 13.838 , 15.307 , 17.259 , 19.951 , 23.514 , 28.291 , 34.483,<br>      42.457 , 52.357 , 64.361 , 78.256 , 93.721 , 109.99 , 126.30 , 141.72 , 155.66 , 167.62,<br>      177.57 , 185.50 , 191.73 , 196.46 , 200.07 , 202.72 , 204.71 , 206.14 , 207.22 , 207.97,<br>      208.55 , 208.93 , 209.25 , 209.44 , 209.61 , 209.70 , 209.80 , 209.84 , 209.9,<br><br> radius_fromHorizGrid=6370.E3,<br> &<br><br># Input datasets<br> &PARM05<br> hydrogThetaFile='Twoa_dec_init_UPD_7km.bin',<br> hydrogSaltFile='Swoa_dec_init_UPD_7km.bin',<br> uVelInitFile='Usoda_dec_init_UPD_7km.bin',<br> vVelInitFile='Vsoda_dec_init_UPD_7km.bin',<br> bathyFile='bathy_NIO_7km_UPD.bin',<br> checkIniTemp=.false.,<br> checkIniSalt=.false.,<br> &</div><div><br></div><div>Kindly advise. Thanks</div><div>Kunal<br></div><div><br></div><div><br></div></div>