c $Header: /usr/people/flittner/radtran/tomrad/pv2.2/src/RCS/iniclz1.f,v 2.22 2000/08/30 16:38:09 flittner Exp $ 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 none c commons include 'depolt.cmn' 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