[RegCNET] regcm NaN ds=45 dt=120

[xmengliu]

 Dear all: 
     When I running RegCM from Dec23 2006 to Feb1 2007, I meet a problm:
 
 
                         ...........
 
 BATS variables written at   2006122300   2.00000000000000     
Writing rad fields at ktau =            1  2006122300 
    at day =    0.0681, ktau =         50 :  1st, 2nd time deriv of ps =  0.74663E-04 0.27378E-06,  no. of points w/convection =     123 
BATS variables written at   2006122303   180.000000000000     
    at day =    0.1375, ktau =        100 :  1st, 2nd time deriv of ps =  0.61828E-04 0.25774E-06,  no. of points w/convection =     148 
    at day =    0.2069, ktau =        150 :  1st, 2nd time deriv of ps =  0.83809E-04 0.41155E-06,  no. of points w/convection =     110 
          *** solar declination angle = -23.41 degrees. 
 dectim =    2880.00000000000     
OUT-history written date =    2006122306.00000     
BATS variables written at   2006122306  0.000000000000000E+000 
Writing rad fields at ktau =          180  2006122306 
BCs are ready from   2006122306   to   2006122312 
    at day =    0.2764, ktau =        200 :  1st, 2nd time deriv of ps =  0.17990E-03 0.10383E-05,  no. of points w/convection =     126 
    at day =    0.3458, ktau =        250 :  1st, 2nd time deriv of ps =  0.40411E-03 0.25749E-05,  no. of points w/convection =     110 
BATS variables written at   2006122309   180.000000000000     
    at day =    0.4153, ktau =        300 :  1st, 2nd time deriv of ps =  0.86163E-03 0.56394E-05,  no. of points w/convection =     117 
    at day =    0.4847, ktau =        350 :  1st, 2nd time deriv of ps =  NaN         NaN        ,  no. of points w/convection =      53 
          *** solar declination angle = -23.41 degrees. 
 dectim =    4320.00000000000     
OUT-history written date =    2006122312.00000     
BATS variables written at   2006122312  0.000000000000000E+000 
Writing rad fields at ktau =          360  2006122312 
BCs are ready from   2006122312   to   2006122318 
    at day =    0.5542, ktau =        400 :  1st, 2nd time deriv of ps =  NaN         NaN        ,  no. of points w/convection =       0 
    at day =    0.6236, ktau =        450 :  1st, 2nd time deriv of ps =  NaN         NaN        ,  no. of points w/convection =       0 
BATS variables written at   2006122315   180.000000000000     
    at day =    0.6931, ktau =        500 :  1st, 2nd time deriv of ps =  NaN         NaN        ,  no. of points w/convection =       0 

 
 
 
I wonder what the problem is ? When I running from Jun25 2006 to Aug1 2006(in summer). No such problem.
 
in getdat.f
 
calday = dble(julday)+(NNNNNN-NSTRT0)/4.+(xtime/60.+gmt)/24. 
 
 
in file  regcm.in 

$restartparam 
ifrest  = .false. , 
idate0  = 2006122300, 
idate1  = 2006122300, 
idate2  = 2007010100, 
nslice  =        120, 
$ 
$timeparam 
radfrq  =    30., 
abemh   =    18., 
abatm   =   240., 
dt      =   120., 
ibdyfrq =          6, 
$ 

I have ever run by time step is 90s ,the NaN is remain.
 
in file domain.param
 
 parameter(iproj='LAMCON') ! Map projection (LAMCON,NORMER,POLSTR,ROTMER)
      parameter(iy=100)          ! number of grid points in y direction (i)
      parameter(jx=120)          ! number of grid points in x direction (j)
      parameter(kz=18)          ! number of vertical levels (14, 18 or 23)
      parameter(nsg=1)          ! number of subgrids in 1 grid (one direction)
      parameter(ds=45.0)        ! grid point separation in km
      parameter(ptop= 10.0)      ! pressure of model top in cb.
      parameter(clat= 25.00)    ! Central latitude (North positive)
      parameter(clon= 110.00)    ! Central longitude (East positive)
      parameter(plat=clat)      ! Pole latitude  (used for ROTMER)
      parameter(plon=clon)      ! Pole longitude (used for ROTMER)
      parameter(truelatL= 30.)  ! LAMCON true latitude (low  latitude side)
      parameter(truelatH= 60.)  ! LAMCON true latitude (high latitude side)
      parameter(ntypec=10)      ! Global topo/landuse data resolution (min)
      parameter(ntypec_s=10)    ! Global topo/landuse data resolution (min)
      parameter(h2opct=75.)     ! Minimum H2O percent to be considered water
      parameter(ifanal =.true.) ! Type of interpolation for topography
      parameter(smthbdy=.true.)! Smoothing of boundaries
      parameter(lakadj=.false.) ! Adjust Great Lakes levels (for US domain)
      parameter(igrads=1)       ! Output GrADS control file?
      parameter(ibigend=1)      ! 1=big-endian; (always 1)
      parameter(ibyte=1)        ! For direct access open statements (1 or 4)

 
 
 Please help me! 
 
 
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20080703/52a0abc9/attachment-0002.html>