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

 

 

 

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!