<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><font size="4">Hi, </font><div dir="auto" style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"></div><div><font size="4"><br></font></div><div><font size="4">I’m trying to simulate the magma ocean in tidally locked lava planet.</font></div><div><font size="4">To speed up the simulation, the initial temperature is set to decrease from 3000 K at the substellar point (0,0) to 50 K at the side boundaries. </font></div><div><font size="4">Meanwhile, the initial temperature is vertically uniform.</font></div><div><font size="4"><br></font></div><div><font size="4">However, there is strong instability<span class="Apple-tab-span" style="white-space:pre"> </span> when I start this simulation. I tried to decrease the maximum initial temperature to 1700 K, and the horizontal contrast becomes to (1700-50) K, but it is still unstable.</font></div><div><font size="4"><br></font></div><div><font size="4">Is there any way I can simulate stably with this initial temperature setup?</font></div><div><font size="4">By the way, the tilmestep deltaT is 20 s.</font></div><div><font size="4"><br></font></div><div><font size="4">Thanks very much!!</font></div><div><font size="4">Yanhong</font></div><div><font size="4"><br></font></div><div><font size="4">Below is my data file:</font></div><div><div># ====================</div><div># | Model parameters |</div><div># ====================</div><div>#</div><div># Continuous equation parameters</div><div> &PARM01</div><div># tRef= 18.8, 16.3, 13.1, 10.4, 8.1, 6.0, 4.5,</div><div># 3.4, 2.7, 2.2, 1.8, 1.5, 1.2, 1.0, 0.9,</div><div> tRefFile='T0solid_diffusion_160x72.bin',</div><div># tRef=44*1600.,</div><div> sRef= 44*34.7,</div><div> cosPower=1.,</div><div> viscAr=1.E-3,</div><div> viscAh=12.E5,</div><div> no_slip_sides=.FALSE.,</div><div> no_slip_bottom=.TRUE.,</div><div> bottomDragLinear=1.E-3,</div><div>#diffK4T=2.E12,</div><div> diffKhT=1.E3,</div><div> diffKrNrT=44*1.E-4,</div><div> ivdc_kappa=100.,</div><div> implicitDiffusion=.TRUE.,</div><div># eosType='JMD95Z',</div><div> eosType='LINEAR',</div><div> tAlpha=4.E-5,</div><div> sBeta =0.E-4,</div><div> saltStepping=.FALSE.,</div><div> gravity=22,</div><div> rhonil=3300.,</div><div> rhoConst=2600.,</div><div> HeatCapacity_Cp=1800.,</div><div> implicitFreeSurface=.TRUE.,</div><div> nonlinFreeSurf=4,</div><div> select_rStar=2,</div><div> hFacInf=0.2,</div><div> hFacSup=2.0,</div><div> exactConserv=.TRUE.,</div><div> staggerTimeStep=.TRUE.,</div><div> useCDscheme=.TRUE.,</div><div> readBinaryPrec=64,</div><div> tempAdvScheme=77,</div><div># usePickupBeforeC54=.TRUE.,</div><div> rotationPeriod=69120.,</div><div># useCoriolis=.FALSE.,</div><div> Tsolidus=1700.,</div><div> Tliquidus=2000.,</div><div> &</div><div><br></div><div># Elliptic solver parameters</div><div> &PARM02</div><div> cg2dMaxIters=1000,</div><div>#cg2dTargetResidual=1.E-13,</div><div> cg2dTargetResWunit=1.E-17,</div><div> &</div><div><br></div><div># Time stepping parameters</div><div> &PARM03</div><div> niter0=0,</div><div> nTimeSteps=20000000,</div><div>#endTime=3110400000.,</div><div> deltaT=10.0,</div><div># deltaTmom=100.0,</div><div># deltaTfreesurf=2000.0,</div><div># deltaTtracer=600.0,</div><div># deltaTClock =600.0,</div><div> forcing_In_AB=.FALSE.,</div><div> abEps=0.1,</div><div> tauCD=321428.,</div><div> tauThetaClimRelax=2000.0,</div><div># pickupStrictlyMatch=.FALSE.,</div><div> chkptFreq=0.0,</div><div> pChkptFreq=5000000.0,</div><div> taveFreq=0.0,</div><div> dumpFreq=5000000.0,</div><div> diagFreq=31104000000.0,</div><div> monitorFreq=31104000000.0,</div><div># monitorFreq=1.0,</div><div> monitorSelect=2,</div><div> &</div><div><br></div><div># Gridding parameters</div><div> &PARM04</div><div> usingSphericalPolarGrid=.TRUE.,</div><div># delRc= 20., 60., 85., 120., 165., 215., 265., 315.,</div><div># 365., 415., 465., 515., 565., 615., 665., 350.,</div><div># delRcFile='delRc.bin',</div><div> delR= 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 20., 20., 20., 20., 20., 20.,</div><div> 40., 40., 40., 40., 40., 60., 60., 60., 60., 80., 80., 80., 100., 100., 100., 100., </div><div> 100., 150., 150., 200., 200., 200., 200., 200., 300., 300., 300., 400.,</div><div> rSphere=9E6,</div><div> ygOrigin=-81.,</div><div> xgOrigin=-180,</div><div> delY=72*2.25,</div><div> delX=160*2.25,</div><div> &</div><div><br></div><div># Input datasets</div><div> &PARM05</div><div> bathyFile='topo3D_160x72_5400m_4walls.bin',</div><div># zonalWindFile='taux_Az1E-5mod_160x72.bin',</div><div># meridWindFile='tauy_Az1E-5mod_160x72.bin',</div><div> thetaClimFile='sst3D_mod_160x72.bin',</div><div> &</div></div><div><br></div><div><br></div></body></html>