regCM version 4 error 求助

nuistxue

进行数值模拟，如下错误！


This is RegCM version 4
   SVN Revision: exported compiled at: data : Aug  3 2013  time: 20:58:23

: this run started at    : Mon Aug 19 11:07:40 2013
: it is submitted by   : xue
: it is running on     : xue
: it is using          :            1   processors
: in directory         : /home/xue/runRegCM
                       
param: RESTARTPARAM READ IN
param: TIMEPARAM READ IN
param: OUTPARAM READ IN
param: PHYSICSPARAM READ IN
param: SUBEXPARAM READ IN
param: EMANPARAM READ IN
  param: dtau =    37.500000000000000        75.000000000000000
  nrec =       289276
  param: IDATE1, IDATE2, dtmin, ktaur =   1990060100  1990070100   2.5000000000000000                0
READING HEADER FILE:/home/xue/runRegCM/input/TRS.INFO                                                                                                                                                                                                                              
param: DIMS          34          64          18
param: DOMAIN   60000.000       45.389999       13.480000       45.389999       13.480000      0.71556687   
param: PROJLAMCON
param: SIGMA   0.0000000      5.00000007E-02  0.10000000      0.16000000      0.23000000      0.31000000      0.38999999      0.47000000      0.55000001      0.63000000      0.70999998      0.77999997      0.83999997      0.88999999      0.93000001      0.95999998      0.98000002      0.99000001       1.0000000   
param: PTOP   5.0000000   
param: OUTPUT           1           1
   param: ***** mdate =   1990060100
  param: input/output parameters
   ifsave =  T  savfrq =    48.000000000000000       iftape =  T  tapfrq =    6.0000000000000000       ifprt  =  F  prtfrq =    12.000000000000000       kxout  =            6  jxsex  =           40  radisp =    6.0000000000000000       batfrq =    3.0000000000000000       nslice =          120  ifchem =  F  chemfrq =   6.0000000000000000       clmfrq =    12.000000000000000

  param: physical parameterizations
   iboudy =  5
   icup =  4
   igcc = 1
   ipptls =  1
   iocnflx =  2
   ipgf =  0
   iemiss =  0
   lakemod =  0
   Chemistry active? 0=no,1=yes  ichem = 0
   Nudge value high range   = 3.000000
   Nudge value medium range = 2.000000
   Nudge value low range    = 1.000000

  param: model parameters
    radfrq =    30.000000
    abatm  =   600.000000
    abemh  =    18.000000
    dt     =   150.000000

  ncld =  1

HT
HTSD
SATBRT
XLAT
XLONG
MSFX
MSFD
F
SNOWC
  
***************************************************
***************************************************
**** RegCM IS BEING RUN ON THE FOLLOWING GRID: ****
****     Map Projection: LAMCON                ****
****     IY=          34  JX=          64  KX=          18              ****
****     PTOP=   5.0000000      DX=   60000.000            ****
****     CLAT=    45.389999      CLON=   13.480000         ****
***************************************************
  
  Index of highest allowed pbl:  kt =  10

AUTO-CONVERSION RATE:  LAND=   0.000250                      OCEAN=   0.000250
  RELATIVE HUMIDITY THRESHOLDS:  LAND=  0.80000001192092896                                    OCEAN=  0.89999997615814209
  GULTEPE FACTORS:  LAND=  0.40000000596046448                       OCEAN=  0.40000000596046448
  MAXIMUM CLOUD COVER FOR RADIATION:   0.80000001192092896
  MAXIMUM RELATIVE HUMIDITY:    1.0099999904632568
  rh0 temperature threshold:    238.00000000000000
  Raindrop Evaporation Rate  1.00000000000000002E-003
  Raindrop Accretion Rate   3.0000000000000000


  EMANUEL (1991) CONVECTION V4.3C (20 May, 2002)
    MIN CONVECTION ORIGIN (minsig/orig):   0.94999998807907104                3
    AUTOCONVERSION THERSHOLD (elcrit):   1.09999999403953552E-003
    AUTOCONVERSION THRESHOLD TO ZERO (tlcrit):   -55.000000000000000
    ENTRAINMENT COEFFICIENT (entp):    1.5000000000000000
    FRACTIONAL AREA OF UNSATURATED DNDRAFT (sigd):   5.00000007450580597E-002
    PRECIP FRACTION OUTSIDE OF CLOUD (sigs):   0.11999999731779099
    FALL SPEED OF RAIN (omtrain):    50.000000000000000
    FALL SPEED OF SNOW (omtsnow):    5.5000000000000000
    RAIN EVAPORATION COEFFICIENT (coeffr):    1.0000000000000000
    SNOW EVAPORATION COEFFICIENT (coeffs):   0.80000001192092896
    CONVECTIVE MOMENTUM TRANSPORT COEFFICIENT (cu):   0.69999998807907104
    DOWNDRAFT VELOCITY SCALE (betae):    10.000000000000000
    MAX NEGATIVE PERTURBATION BELOW LFC (dtmax):   0.89999997615814209
    QUASI-EQUILIBRIUM APPROACH RATE (alphae):   0.20000000298023224
    QUASI-EQUILIBRIUM APPROACH RATE (damp):   0.10000000149011612



  CONVECTIVE CLOUD FRACTION/WATER
     Maximum Convective Cloud Cover
       before resolution scaling:   0.25000000000000000
     Maximum Convective Cloud Cover
       after resolution scaling:   0.25000000000000000
     Convective Cloud Water:   4.99999987368937582E-005
     the surface energy budget is used to calculate the ground temperature.   julday = 152   gmt =  0.0
          the radiation is computed every   12 time steps.

   relaxation boudnary conditions (exponential method) are used. fnudge =     0.33333E-03  gnudge =     0.48000E+06
0 k    sigma(k)     a(k)     dsigma(k)    twt(k,1)     twt(k,2)     qcon(k)

  1     0.0000     0.0250     0.0500       0.0000       0.0000       0.0000
  2     0.0500     0.0750     0.0500       0.5000       0.5000       0.5000
  3     0.1000     0.1300     0.0600       0.4545       0.5455       0.5455
  4     0.1600     0.1950     0.0700       0.4615       0.5385       0.5385
  5     0.2300     0.2700     0.0800       0.4667       0.5333       0.5333
  6     0.3100     0.3500     0.0800       0.5000       0.5000       0.5000
  7     0.3900     0.4300     0.0800       0.5000       0.5000       0.5000
  8     0.4700     0.5100     0.0800       0.5000       0.5000       0.5000
  9     0.5500     0.5900     0.0800       0.5000       0.5000       0.5000
10     0.6300     0.6700     0.0800       0.5000       0.5000       0.5000
11     0.7100     0.7450     0.0700       0.5333       0.4667       0.4667
12     0.7800     0.8100     0.0600       0.5385       0.4615       0.4615
13     0.8400     0.8650     0.0500       0.5455       0.4545       0.4545
14     0.8900     0.9100     0.0400       0.5556       0.4444       0.4444
15     0.9300     0.9450     0.0300       0.5714       0.4286       0.4286
16     0.9600     0.9700     0.0200       0.6000       0.4000       0.4000
17     0.9800     0.9850     0.0100       0.6667       0.3333       0.3333
18     0.9900     0.9950     0.0100       0.5000       0.5000       0.5000
19     1.0000


maximum time =   30.000 days.
time step =  150.00 seconds
dx =  60000. meters
grid points (x,y) = (  64,  34)
number of levels = 18
constant hor. diff. coef. =  0.36000E+05 m*m/s
maximum  hor. diff. coef. =  0.37500E+06 m*m/s
READING INITAL CONDITIONS  1990060100
U
V
TA
QV
PS
TS
            *** solar declination angle =  22.07 degrees.
ozone profiles
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.5250E-07
   0.5289E-07
   0.5775E-07
   0.6296E-07
   0.7186E-07
   0.8557E-07
   0.1074E-06
   0.1500E-06
   0.2283E-06
   0.4177E-06
   0.8028E-06
   0.1420E-05
   0.3188E-05
   0.0000E+00
  dectim =    1440.0000000000000
BCs are ready from   1990060100   to   1990060106
dt, dtau =    150.00000000000000        37.500000000000000        75.000000000000000     
  Calculating Vertical Modes
  - Linearization about standard atmosphere
0sigmaf      0.00      5.000E-02  0.100      0.160      0.230      0.310      0.390      0.470      0.550      0.630      0.710   
           0.780      0.840      0.890      0.930      0.960      0.980      0.990       1.00   
0t mean      218.       218.       218.       218.       230.       239.       248.       255.       262.       268.       273.   
            277.       280.       283.       285.       286.       287.       287.   
0ps mean     100.   
0 vertical mode problem completed for kx= 18     0 errors detected   (should be 0)
m, fac =            4   33.333333333333336     
m, fac =            2   60.000000000000000     
Writing output files in direct access format
  
******* OPENING NEW OUTPUT FILES:  1990060100
OPENING NEW OUT FILE: /home/xue/runRegCM/output//ATM.1990060100
OPENING NEW SRF FILE: /home/xue/runRegCM/output//SRF.1990060100
OPENING NEW RAD FILE: /home/xue/runRegCM/output//RAD.1990060100
OUT-history written date =    1990060100.0000000     
LAND variables written at   1990060100   2.5000000000000000     
Writing rad fields at ktau =            1  1990060100
     at day =    0.0851, ktau =         50 :  1st, 2nd time deriv of ps =  0.26715E-04 0.98959E-07,  no. of points w/convection =     148
LAND variables written at   1990060103   180.00000000000000     
     at day =    0.1719, ktau =        100 :  1st, 2nd time deriv of ps =  0.14142E-04 0.29969E-07,  no. of points w/convection =     180
            *** solar declination angle =  22.11 degrees.
  dectim =    2880.0000000000000     
OUT-history written date =    1990060106.0000000     
LAND variables written at   1990060106   0.0000000000000000     
Writing rad fields at ktau =          144  1990060106
BCs are ready from   1990060106   to   1990060112
     at day =    0.2587, ktau =        150 :  1st, 2nd time deriv of ps =  0.98340E-05 0.19447E-07,  no. of points w/convection =     286
     at day =    0.3455, ktau =        200 :  1st, 2nd time deriv of ps =  0.80761E-05 0.13971E-07,  no. of points w/convection =     684
LAND variables written at   1990060109   180.00000000000000     
     at day =    0.4323, ktau =        250 :  1st, 2nd time deriv of ps =  0.79871E-05 0.11274E-07,  no. of points w/convection =     807
            *** solar declination angle =  22.14 degrees.
  dectim =    4320.0000000000000     
OUT-history written date =    1990060112.0000000     
LAND variables written at   1990060112   0.0000000000000000     
Writing rad fields at ktau =          288  1990060112
BCs are ready from   1990060112   to   1990060118
     at day =    0.5191, ktau =        300 :  1st, 2nd time deriv of ps =  0.11038E-04 0.13815E-07,  no. of points w/convection =     866
     at day =    0.6059, ktau =        350 :  1st, 2nd time deriv of ps =  0.99502E-05 0.13795E-07,  no. of points w/convection =     855
LAND variables written at   1990060115   180.00000000000000     
     at day =    0.6927, ktau =        400 :  1st, 2nd time deriv of ps =  0.84321E-05 0.10493E-07,  no. of points w/convection =     665
            *** solar declination angle =  22.17 degrees.
  dectim =    5760.0000000000000     

Program received signal 8 (SIGFPE): Floating-point exception.

Backtrace for this error:
  + /lib64/libc.so.6() [0x3f86432920]
  + function outtap_ (0x50E0BD)
    at line 214 of file outtap.F90
  + function output_ (0x5571D3)
    at line 2042 of file output.F90
  + in the main program
    at line 342 of file regcm.F90
  + /lib64/libc.so.6(__libc_start_main+0xfd) [0x3f8641ecdd]
** Warning: a core dump was requested, but the core sizelimit
**          is currently zero.

Quit

bjj948

计算不稳定。
可以减少时间步长试试。

nuistxue

bjj948 发表于 2013-8-20 08:10

                               
登录/注册后可看大图

计算不稳定。
可以减少时间步长试试。

   多谢老师，确与步长相关后采用这样一组参数（  
     radfrq =    30.000000
    abatm  =   600.000000
    abemh  =    18.000000
    dt     =   100.000000）并设置
  truelatl = 10.0,
truelath = 30 （ LAMCON  投影， 替换原来的30.0和60）
    running 无错。
补充问题：计算不稳定，主要是我PC计算资源限制造成还是算法复杂度问题，目前近会去尝试不同的参数，其中原理不甚明白。