#include "dims.h"
      subroutine mgwinti
C-----------------------------------------------------------------------
C
C Time independent initialization for multiple gravity wave parameterization.
C
C-----------------------------------------------------------------------
      implicit none
C-----------------------------------------------------------------------
#include "pmgrid.h"
C-----------------------------------------------------------------------
#include "mgw.h"
C-----------------------------------------------------------------------
#include "comhyb.h"
C-----------------------------------------------------------------------
C
C Local variables
C
      integer k
      real
     $     oroe,                 ! efficiency factor for orographic term
     $     fac1
C
C-----------------------------------------------------------------------
C
c+
c Copy model constants
c-
      cp     = cpx
      cpvir  = cpwv / cp - 1
      g      = gx
      r      = rx
c+
c Set MGWD constants
c-
      kwv    = 6.28e-5          ! 100 km wave length
      dback  = 0.05             ! background diffusivity
      orofc2 = .5               ! critical froude number squared
      oroe   = .125             ! efficiency for orographic generation
      tauscal= .001             ! scale factor for background stress
      taubgnd= .4               ! background stress amplitude

      fracldv= 0.0              ! fraction of tau0 diverged in low level region
      zldvcon= 10.              ! constant for determining zldv
c+
c Set phase speeds 
c-
      do k = -pgwv, pgwv
         cphs(k)   = 10.* k       ! 0, +/- 10, +/- 20, ... m/s
c        cphs(k)   = 15.* k       ! 0, +/- 15, +/- 30, ... m/s walterschied
      end do
      cphs(0) = 0.
c+
c Set radiative damping times
c-
      do k = 0, plev
         alpha(k) = 1.e-6       ! about 10 days.
      end do
c+
c Min and max values to keep things reasonable
c-
      mxasym = 1.e-9           ! max factor of 10 from |tau(c)| to |tau(-c)|
      mxrange= 1.e-9           ! factor of 100 from max to min |tau(c)|
      n2min  = 1.e-8
      orohmin= 10.              ! min sdv for orographic wave drag
      orovmin=  2.              ! min wind speed for orographic wave drag
      taumin = 1.e-10
      tndmin = 1.e-10
      tndmax = 500. / 86400.    ! max permitted tendency (500 m/s/day)
      ubmc2mn= 0.01
c+
c Determine efficiency factor from gravity wave parameters
c-
      oroeko2= 0.5 * oroe * kwv
c      write(*,*)'fac1: multiply efficiency factor'
      fac1 = 1.
      efcncy = oroe * orofc2*fac1
c+
c Determine other derived constants
c-
      rog    = r/g
      efkw   = efcncy * kwv
c+
c Determine the bounds of the background and orographic stress regions
c-
      ktopbg  = 0
      kbotoro = plev
      kbotbg = plev
      ktoporo = 0
c      do k = 0, plev
c         if (hypi(k+1) .lt. 10000.) kbotbg  = k
c         if (hypi(k+1) .lt.  3000.) ktoporo = k
c      end do
c
      return
      end