subroutine iniclz1
c***********************************************************************
c
c        last modified by zia ahmad    9/10/93
c        purpose: to include the effect of molecular anisotropy
c
c***********************************************************************
      implicit integer*4(i-n), real*8(a-h,o-z)
c     new statement
      include "depolt.inc"
c     end of new statement
c
c     new statements
      if(ipol.eq.1)then
c         rhon=0.035d0
          gama=rhon/(2.0d0-rhon)
          q2=(1.0d0-gama)/(1.0d0+2.0d0*gama)
          q1=0.75d0*q2
          q=(2.0d0*gama)/(1.0d0-gama)
          sdp=1.0d0+q
          delp=dsqrt(2.0d0+3.0d0*q)
          q12s=q1/(2.0d0*sdp)
      endif
c     end of new statements
c
c      write(33,500)ipol
 500  format('sub iniclz1 ... ipol=',i2)
c      write(33,501)rhon,gama,q,q1,q2,sdp,delp,q12s
 501  format('sub iniclz..rhon etc.'/1p6e12.4/1p2e12.4)
      return
      end