C*********************************************************************
C  2000/10/03  Gierasch and Goddy (1968) surface heat flux model 
C  
C*********************************************************************
C
      SUBROUTINE FXSURF
     I   ( TPOT0, CDRAGV, VSFC0,
     O     FHEAT, 
     B     DENS0, PPAI0,
     C     CPDRY, GRAV , 
     S     TSFC , 
     D     NX   , NXB  , NZ   , NZB )
C
C*********************************************************************
      IMPLICIT REAL*8 ( A-H, O-Z )
C*********************************************************************
C
      DIMENSION
     \   TPOT0( -2:NZB )

      DIMENSION
     \   CDRAGV( -2:NXB )
C
      DIMENSION
     \   FHEAT ( 2 )
C
      DIMENSION
     \   TSFC ( -2:NXB )
C
      DIMENSION
     \   PPAI0( -2:NZB ),
     \   DENS0( -2:NZB )
C
C====================================================================
C
      DIFFUS = 8.0D-4
      VISCOS = 1.0D-3
      FACT   = 0.089D0
C
C====================================================================
C
      DO 11 IX = 1, 1
         FHEAT(2) = 0.0d0
         FHEAT(1) = DENS0(0) * CPDRY * FACT * DIFFUS 
     \    * ( 2.0 * MAX( 0.0, ( TSFC(IX) - TPOT0(0) * PPAI0(0) ) ) )
     \     **(4.0D0/3.0D0)
     \    * ( GRAV / ( TPOT0(0) * PPAI0(0) * DIFFUS * VISCOS ) )
     \     **(1.0D0/3.0D0)
   11 CONTINUE
C
C====================================================================
C
      RETURN
      E N D
