[MITgcm-support] About Warnings

[limingjie07 at qdio.ac.cn]

Hi all,
   I meet some warnings though the model seems to go well. I am not sure that whether it will affect my results. Something may not be setted exactly in my data file. Is viscAh too large? Any help will be appreciated.                                                                         
                                                 Mingjie
 Warning reads as follows:
(PID.TID 0020.0001) ** WARNING ** CONFIG_CHECK: Implicit viscosity applies to provisional u,vVel
(PID.TID 0020.0001) ** WARNING ** => not consistent withfinal vertical shear (after appling 3-D solver solution
(PID.TID 0020.0001) ** WARNING ** CONFIG_CHECK: Implicit viscosity not implemented in CALC_GW
(PID.TID 0020.0001) ** WARNING ** CONFIG_CHECK: Explicit viscosity might become unstable if too large

 Part of my data file reads as follows:
 &PARM01
 tRef = 29.03, 28.75, 28.43, 28.04, 27.57, 27.00, 26.27, 25.54, 24.73, 23.87,
        23.11, 22.34, 21.57, 20.94, 20.28, 19.72, 19.20, 18.68, 18.15, 17.63,
        17.19, 16.79, 16.39, 16.00, 15.60, 15.26, 14.94, 14.62, 14.30, 13.98,
        13.72, 13.45, 13.15, 12.81, 12.42, 11.98, 11.47, 10.96, 10.40,  9.75,
         9.03,  8.37,  7.61,  6.94,  6.22,  5.53,  4.92,  4.30,  3.71,  3.24,
         2.83,  2.59,  2.43,  2.35,  2.29,  2.09,  1.83,  1.43,  1.21,  1.12,
 sRef = 33.86, 33.95, 34.05, 34.14, 34.24, 34.33, 34.43, 34.52, 34.57, 34.62,
        34.66, 34.69, 34.71, 34.71, 34.72, 34.73, 34.72, 34.72, 34.72, 34.71,
        34.71, 34.69, 34.68, 34.67, 34.65, 34.64, 34.63, 34.61, 34.60, 34.58,
        34.57, 34.56, 34.55, 34.53, 34.51, 34.50, 34.47, 34.45, 34.43, 34.40,
        34.38, 34.37, 34.36, 34.37, 34.39, 34.41, 34.44, 34.47, 34.51, 34.53,
        34.56, 34.60, 34.64, 34.67, 34.71, 34.72, 34.71, 34.70, 34.69, 34.71,
 
 viscAr = 1.E-5,
 viscAh =1.E3,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 bottomDragQuadratic=2.5E-3,
 diffKhT= 1.E3,
 diffKrT= 1.E-5,
 diffKhS= 1.E3,
 diffKrS= 1.E-5,
# f0=4.25e-5,
# beta=1.E-11,
 
 eosType='MDJWF',
# tAlpha=2.E-4,
# sBeta =0.E-4,
 gravity=9.81,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 implicitFreeSurface=.FALSE.,
 exactConserv=.TRUE.
 useRealFreshWaterFlux=.TRUE.,
 allowFreezing=.FALSE.,
 
 tempAdvScheme=33,
 saltAdvScheme=33,
 StaggerTimeStep=.TRUE.,
 nonHydrostatic=.TRUE.,
 hFacMin=0.2,
#implicSurfPress=0.5,
#implicDiv2DFlow=0.5,
#nonlinFreeSurf=3,
 hFacInf=0.2,
 hFacSup=1.8,
# saltStepping=.FALSE.,
 
#- not safe to use globalFiles in multi-processors runs
 globalFiles=.FALSE.,
 useSingleCpuIO=.TRUE.,
 readBinaryPrec=32,
 writeBinaryPrec=32,
 writeStatePrec=32,
 &
 
# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=400,
 cg3dTargetResidual=1.E-13,
 &
 
# Time stepping parameters
 &PARM03
 nIter0=0,
# nTimeSteps=16,
 endtime=86400.,
 deltaT=45.,
# deltaTtracer=120.,
# deltaTClock=120.,
 abEps=0.1,
 pChkptFreq=86400.,
 chkptFreq=0.,
 dumpFreq=3600.,
 monitorFreq=3600.,
 monitorSelect=2,
 
 &
# Gridding parameters
 &PARM04
# usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.TRUE.,
 ygOrigin=18.,
 dySpacing=.016667,
 dxSpacing=.016667,      
 delR=  10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        10, 11, 12, 14, 16, 18, 21, 24, 28, 33,
        38, 44, 51, 60, 70, 82, 96,113,133,156,
       184,217,256,302,356,420,495,584,620,736,
 & 
...