[MITgcm-support] abnormal with extreme temperature value

李志远 oceanlizy at 163.com
Fri Mar 7 02:56:30 EST 2014


From: Li zhiyuan <oceanlizy at 163.com>
To: MITgcm <mitgcm-support at mitgcm.org>
Subject:"abnormal with extreme temperature value "


Hi,all;
     I am running a model of semi-enclosed domain with 134x165x16 grid points . the model is driven by surface wind and heat flux .when I run my mitgcm model for some steps , very large temperature values occur ,then the model blows up .
    the STDOUT file is as below :
(PID.TID 0005.0001) %MON time_tsnumber                =                   163
(PID.TID 0005.0001) %MON time_secondsf                =   7.0416000000000E+06
(PID.TID 0005.0001) %MON dynstat_eta_max              =   7.4690610491135E-01
(PID.TID 0005.0001) %MON dynstat_eta_min              =  -4.1593172249303E+00
(PID.TID 0005.0001) %MON dynstat_eta_mean             =   2.0772424301492E-02
(PID.TID 0005.0001) %MON dynstat_eta_sd               =   4.9627339137015E-01
(PID.TID 0005.0001) %MON dynstat_eta_del2             =   7.4943932050348E-04
(PID.TID 0005.0001) %MON dynstat_uvel_max             =   3.5280025421360E+00
(PID.TID 0005.0001) %MON dynstat_uvel_min             =  -1.4035655480410E+00
(PID.TID 0005.0001) %MON dynstat_uvel_mean            =  -1.8326591842827E-02
(PID.TID 0005.0001) %MON dynstat_uvel_sd              =   4.5870513051871E-02
(PID.TID 0005.0001) %MON dynstat_uvel_del2            =   9.8327354406657E-05
(PID.TID 0005.0001) %MON dynstat_vvel_max             =   2.4380767316715E+00
(PID.TID 0005.0001) %MON dynstat_vvel_min             =  -2.3087621147168E+00
(PID.TID 0005.0001) %MON dynstat_vvel_mean            =   3.4329723926579E-02
(PID.TID 0005.0001) %MON dynstat_vvel_sd              =   9.0964738317473E-02
(PID.TID 0005.0001) %MON dynstat_vvel_del2            =   1.9834779107181E-04
(PID.TID 0005.0001) %MON dynstat_wvel_max             =   2.6286782112003E-02
(PID.TID 0005.0001) %MON dynstat_wvel_min             =  -9.6061814986124E-02
(PID.TID 0005.0001) %MON dynstat_wvel_mean            =   3.6700902646170E-05
(PID.TID 0005.0001) %MON dynstat_wvel_sd              =   3.4462742242225E-03
(PID.TID 0005.0001) %MON dynstat_wvel_del2            =   7.4810854960120E-06
(PID.TID 0005.0001) %MON dynstat_theta_max            =   5.7517891007917E+01
(PID.TID 0005.0001) %MON dynstat_theta_min            =  -1.0228097331456E+03
(PID.TID 0005.0001) %MON dynstat_theta_mean           =   7.8087020365565E+00
(PID.TID 0005.0001) %MON dynstat_theta_sd             =   6.2079576234231E+00
(PID.TID 0005.0001) %MON dynstat_theta_del2           =   1.3811582622522E-02
(PID.TID 0005.0001) %MON dynstat_salt_max             =   3.5000000000152E+01
(PID.TID 0005.0001) %MON dynstat_salt_min             =   3.4999999953140E+01
(PID.TID 0005.0001) %MON dynstat_salt_mean            =   3.4999999999999E+01
 
my data file :
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 16*20.,
 sRef = 16*35.,
 viscAz=1.E-4,
 viscAh=4.E2,
 diffKhT=4.E2,
 diffKzT=3.E-5,
 diffKhS=4.E2,
 diffKzS=3.E-5,
 rhoConst=1035.,
 rhoConstFresh=1000.,
 eosType = 'LINEAR',
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.FALSE.,
 exactConserv=.TRUE.,
 hFacMin=.3,
 readBinaryPrec=32,
 &


# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.E-13,
 &


# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 200,
 deltaTmom = 30.,
 deltaTtracer = 1200.,
 deltaTClock  = 43200.,
 abEps = 0.1,
 dumpFreq=   864000.,
 monitorFreq=1.,
 &


# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delz= 10,20,25,30,50,100,100,
        200,300,400,500,500,500,
        500,500,500,
 ygOrigin=24.500,
 xgOrigin=117.525,
 delX=134*0.1,
 delY=165*0.1,
 &


# Input datasets
 &PARM05
 bathyFile=      'topo.bin',
 hydrogThetaFile='t16.bin',
# hydrogSaltFile= 'Lev_clim_salt.bin',
# zonalWindFile=  'Taux.bin',
# meridWindFile=  'Tauy.bin',
#thetaClimFile=  'SST.bin',
# saltClimFile=   'SSS.bin',
 &


I guess it  blow up because  mixing coefficients or advection scheme in my model  is incorrect , how I can avoid such problems ?


any quick reply is great appreciated !
lizhi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140307/e3cca883/attachment-0001.htm>


More information about the MITgcm-support mailing list