copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
***s/r qntstar

      subroutine qntstar(qstr,nstr2,gotstr,ortstr,hh1,hh2,ng,kd,kf) 14

      implicit none
*
      integer ng,kd,kf
      real gotstr(ng,kf), ortstr(ng,kd:kf),
     $      nstr2(ng,kf),   qstr(ng,kd:kf),
     $        hh1(ng,kf),    hh2(ng,kd:kf)
*
#include "consdyn_8.cdk"
#include "yomdyn1.cdk"
*
      integer i,k
      real*8 p00,t00,n002,h00r,zero,pt5,one
      parameter(zero=0.d0,pt5=.5d0,one=1.d0)
*
*---------------------------------------------------------------------
***** isothermal temperature profile
***** one parameter: t00
*
      p00   = dlog(100000.0d0)
      t00   = grtstar
      n002  = grav_8*grav_8/(cpd_8*t00)
      h00r  = grav_8/(rgasd_8*t00)
*
      do k=1,kf
         do i=1,ng
            gotstr(i,k) = grav_8 / t00
            nstr2 (i,k) = n002
         end do
      end do
      do k=kd,kf
         do i=1,ng
            ortstr(i,k) = one / ( rgasd_8 * t00 )
              qstr(i,k) = p00 - hh2(i,k) * h00r
         end do
      end do
*
*---------------------------------------------------------------------
      return
      end