Back to home page

darwin3

 
 

    


Warning, /verification/global_ocean.cs32x15/input.viscA4/data is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit b15f6f1e on 2026-01-19 15:35:47 UTC
a0691689e0 Jean*0001 # ====================
                0002 # | Model parameters |
                0003 # ====================
                0004 #
                0005 # Continuous equation parameters
                0006  &PARM01
                0007  tRef=15*20.,
                0008  sRef=15*35.,
e6f4278565 Jean*0009  viscAh =1.E5,
05a00b43fa Jean*0010 #- biharmonic Viscosity: 3.e15 is close to the stability limit with deltaTMom=20mn
e6f4278565 Jean*0011  viscA4 =3.E14,
a0691689e0 Jean*0012  viscAr =1.E-3,
d68faeaed2 Jean*0013  diffKhT=0.,
                0014  diffK4T=0.,
a0691689e0 Jean*0015  diffKrT=3.E-5,
d68faeaed2 Jean*0016  diffKhS=0.,
                0017  diffK4S=0.,
a0691689e0 Jean*0018  diffKrS=3.E-5,
                0019  ivdc_kappa=10.,
                0020  implicitDiffusion=.TRUE.,
                0021  gravity=9.81,
05a00b43fa Jean*0022  rhoConst=1035.,
a0691689e0 Jean*0023  rhoConstFresh=1000.,
                0024  eosType='JMD95Z',
                0025  staggerTimeStep=.TRUE.,
                0026  vectorInvariantMomentum=.TRUE.,
                0027  implicitFreeSurface=.TRUE.,
                0028  exactConserv=.TRUE.,
79b5d5775c Jean*0029 #-- just to test (not physically justified):
d958cb0c95 Jean*0030  selectP_inEOS_Zc=3,
e6f4278565 Jean*0031  nonHydrostatic=.TRUE.,
d958cb0c95 Jean*0032  selectNHfreeSurf=1,
e6f4278565 Jean*0033  useNHMTerms=.TRUE.,
79b5d5775c Jean*0034  implicDiv2DFlow=0.5,
                0035  selectAddFluid=1,
                0036 #---
a0691689e0 Jean*0037  useRealFreshWaterFlux=.TRUE.,
                0038  allowFreezing=.TRUE.,
                0039  hFacMin=.1,
                0040  hFacMinDr=20.,
                0041  readBinaryPrec=64,
d958cb0c95 Jean*0042 #writeBinaryPrec=64,
cc30a7a996 Jean*0043  useSingleCpuIO=.TRUE.,
246faa4842 Jean*0044  plotLevel=0,
a0691689e0 Jean*0045  &
                0046 
                0047 # Elliptic solver parameters
                0048  &PARM02
                0049  cg2dMaxIters=200,
e6e223b277 Jean*0050  cg3dMaxIters=60,
a0691689e0 Jean*0051 #cg2dTargetResidual=1.E-9,
e6e223b277 Jean*0052 #cg3dTargetResidual=1.E-9,
b15f6f1e9f Jean*0053  cg2dTargetResWunit=7.E-13,
                0054  cg3dTargetResWunit=4.E-13,
a0691689e0 Jean*0055  &
                0056 
                0057 # Time stepping parameters
                0058  &PARM03
e6f4278565 Jean*0059  nIter0=86400,
a0691689e0 Jean*0060  nTimeSteps=10,
e6f4278565 Jean*0061  deltaTClock =3600.,
                0062  abEps = 0.12,
                0063  forcing_In_AB=.FALSE.,
f4b54c2014 Jean*0064  pickupStrictlyMatch=.FALSE.,
                0065  pChkptFreq  =31104000.,
e6f4278565 Jean*0066  dumpFreq    =31104000.,
                0067  monitorFreq =86400.,
79b5d5775c Jean*0068 #-- just to test: switch off periodic forcing --> constant forcing
                0069 # periodicExternalForcing=.TRUE.,
                0070 # externForcingPeriod=2592000.,
                0071 # externForcingCycle=31104000.,
                0072 #--
a0691689e0 Jean*0073 # 2 months restoring timescale for temperature
                0074  tauThetaClimRelax = 5184000.,
                0075 # 2yrs restoring timescale for salinity
                0076  tauSaltClimRelax = 62208000.,
                0077  monitorFreq =1.,
                0078  &
                0079 
                0080 # Gridding parameters
                0081  &PARM04
                0082  usingCurvilinearGrid=.TRUE.,
e6f4278565 Jean*0083  horizGridFile='grid_cs32',
4e59369e7e Jean*0084  radius_fromHorizGrid=6370.E3,
a0691689e0 Jean*0085  delR= 50., 70., 100., 140., 190.,
                0086        240., 290., 340., 390., 440.,
                0087        490., 540., 590., 640., 690.,
4e59369e7e Jean*0088  rSphere=6370.E3,
a0691689e0 Jean*0089  &
                0090 
                0091 # Input datasets
                0092  &PARM05
                0093  bathyFile      ='bathy_Hmin50.bin',
                0094  hydrogThetaFile='lev_T_cs_15k.bin',
                0095  hydrogSaltFile ='lev_S_cs_15k.bin',
                0096  zonalWindFile  ='trenberth_taux.bin',
                0097  meridWindFile  ='trenberth_tauy.bin',
                0098  thetaClimFile  ='lev_surfT_cs_12m.bin',
                0099  saltClimFile   ='lev_surfS_cs_12m.bin',
7130bed1c3 Jean*0100  surfQnetFile   ='shiQnet_cs32.bin',
79b5d5775c Jean*0101 #EmPmRFile      ='shiEmPR_cs32.bin',
                0102 #- just to test using addMassFile:
                0103  addMassFile    ='addMass_k1.bin',
                0104 #EmPmRFile      ='large_Precip.bin',
a0691689e0 Jean*0105  &