subroutine chksnow
c
c At the first physics time step after the top of each hour, check the snow
c array against the SR (snow/total precip ratio).  If SR .ge. 0.9, set this
c point to missing (so we won't do snow adjustment here).
c
      INCLUDE "parmeta"
      INCLUDE "mpif.h"
      INCLUDE "mpp.h"
C-----------------------------------------------------------------------
      INCLUDE "PPTASM.comm"
C-----------------------------------------------------------------------
      INCLUDE "CLDWTR.comm"
C-----------------------------------------------------------------------
      INCLUDE "CTLBLK.comm"
C-----------------------------------------------------------------------
      TIMES=(NTSD-1)*DT
      if (mod(times,3600.) .lt. nphs*dt) then
        ihr=(times-1.0)/3600.+1
        call loc2glb(sr,temp1)
        if (mype.eq.0) then 
           write(98,*) 'NTSD=', NTSD,' times=', times,' output SR:'
           write(90) temp1
        endif
        DO 20 J=MYJS,MYJE
          DO 10 I=MYIS,MYIE
            if (sr(i,j).ge.0.9) pptdat(i,j,ihr) = 999.
  10	  continue
  20    continue
      endif
      return
      end