copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
***s/r cxxpar -- Computes parameters c01 to c07
*

      subroutine cxxpar 1
      implicit none
*
#include "consdyn_8.cdk"
#include "yomdyn.cdk"
#include "yomdyn1.cdk"
#include "dynpar.cdk"
      integer gnelast
*
      real*8 pt5,one
      parameter(pt5=0.5d0,one=1.0d0)
*---------------------------------------------------------------
*
* hydrostatic approximation: gnnohyd = 0
* Boussinesq  approximation: gnbousq = 0
* homogeneous fluid approx.: gnhomog = 0
*
* - modify these two at your own risk:
*     initial conditions must be prepared consistently !
* - you may not set gnhomog=0 and gnnohyd=0 simultaneously !
*
      gnbousq = 1
      gnhomog = 1
      gnelast = 1
*
      if (theocase.eq.'BROOK') gnbousq = 0
      if (theocase.eq.'BROOK') gnhomog = 0
      if (theocase.eq.'BROOK') gnelast = 0
*
      c01  =  one                       *dble(gnelast)
      c02  =  one                       *dble(gnbousq)
      c03  =  one                                     *dble(gnnohyd)
      c04  =  one / grav_8              *dble(gnbousq)
      c05  =  pt5*c02
      c06  =  grav_8                                  *dble(gnhomog)
      c07  =  cappa_8/(one-cappa_8)     *dble(gnbousq)
*
      gama_star = grav_8 / ( cpd_8 * grtstar ) * c02
      c2r_star = (one-cappa_8) / ( rgasd_8 * grtstar )
*---------------------------------------------------------------
      return
      end