!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
!                     version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
!

      SUBROUTINE HOUGH ( KULOUT ) 1,2
#if defined (DOC)
*
***s/r HOUGH: Subroutine to compute the vertical and horizontal normal modes.
*
* Author: Luc Fillion *RPN/AES - May 93 
*
* Revision: L. Fillion - Jan 97 - Remove references to previous SEF code (inappropriate).
*Arguments
*     i : KULOUT   = Unit used for optional printing
*
#endif
      IMPLICIT NONE
*implicits
#include "comlun.cdk"
*
      INTEGER  KULOUT
C
      INTEGER  IERR
C
*MODULES
      EXTERNAL VERMOD
      EXTERNAL HPCHECK, ABORT3D
*
**
C
C* 1. COMPUTE THE NORMAL MODES
C     ------------------------
C
 100  CONTINUE
C
C*1.1   VERTICAL PROBLEM
C
 110    CONTINUE
        CALL VERMOD ( KULOUT )
C
      CALL HPCHECK(IERR)
      IF(IERR.EQ.0) THEN
      WRITE(KULOUT,FMT='(5X,10("*"),"*****************************")')
      WRITE(KULOUT,FMT='(5X,5(" "),"NO MEMORY OVERFLOW DETECTED")')
      WRITE(KULOUT,FMT='(5X,5(" "),"PROGRAM NMODES ENDED NORMALLY")')
      WRITE(KULOUT,FMT='(5X,10("*"),"****************************")')
      ELSE
      WRITE(KULOUT,FMT='(5X,10("?"),"????????????????????????")')
      WRITE(KULOUT,FMT='(5X,5(" "),"MEMORY OVERFLOW DETECTED")')
      WRITE(KULOUT,FMT='(5X,5(" "),"PROGRAM NMODES STOPS")')
      WRITE(KULOUT,FMT='(5X,10("?"),"????????????????????????")')
      CALL ABORT3D(NULOUT,'HOUGH           ')
      ENDIF
C
C
      RETURN
      END