c	************This is chen99b.f copied from chen102b.f for the case
c	************when gama=50 and n=1, for coef for the power of 
c	when calculating variable coef below. all other variables remain the same
c       *******	Date 4/30/04 *********
c 
c	Tetas=.434, tetar=.218, alfa=.02, ks=0.367*10e-3 cm/sec, 31.6cm/d  *******
c	************* input is called chen102b.dat  *********
c	this program calculates the coef. for chen et. al. paper wat. R.R.
c	April 2001. Today's date 7,3,03
c
c	*** time in second, Ks is in cm/day, so it needs to be mul bt 24*3600 ***
c
	real tetas, tetar,ks,t100,t300,t900,t10000,alfa,kstar,ta,coef,gama
	real t5000,fourt, onefrt,an,sfourt,tpower2,t2000 
c	an is = n which is either 1 or 2. if n=1 then gama is 5
c	and when n=2 gama is = to 100. ******************
c
	open(2,file='chen99.dat')
	read(2,100)tetas,tetar,ks,t100,t900,alfa,gama,an
c	an=2.
	t100=100.00
	t900=900.00
	t2000=2000.00
	t10000=10000.00
	t40000=40000.
	t5000=5000.00
	t100000=100000.0
	t500000=500000.0
	tm=10**6.
	ttwom=2.*10**6.
100	format(6f10.5,/,2f10.5)
 	kstar=ks/(3600.*24.*(tetas-tetar))
c       tb=kstar*alfa*ttwom
c       tb=kstar*alfa*tm
c       tb=kstar*alfa*t500000
c	tb=kstar*alfa*t100000
c       tb=kstar*alfa*t40000
c	tb=kstar*alfa*t10000
c	tb=kstar*alfa*t900
c	tb=kstar*alfa*t2000
c	tb=kstar*alfa*t900
c	tb=kstar*alfa*t300
  	tb=kstar*alfa*t100
	tpower2=(4.*tb)**2.5
	fourt=4.*tb
	sfourt=(4.*tb)**.5
	onefrt=tb/4.
c	coef=tetas*gama*an*(4.*tb)**2.5
        coef=tetas*gama*an*(4.*tb)**.5
	write(6,120)tetas,tetar,ks,t100,t900,alfa,tb,kstar,coef,tpower2,
     %fourt,onefrt,sfourt
c	Def of the constants: K* is as defined in Chen's paper, tb is the dimensionless
c	t or T=(K*)**2*t/D,or T=Ks*alfa*t/(tetas-tetar) bott. of pg. 1091 . Coef 
c	 is the coef. for the second term in equ. 14. of chen.
c	 
c
120	format(1x,' tetas= ',f7.4,' tetar= ',f7.4,' ks=',f9.4,
     %' t100=',f9.2,/,' t900=',f9.2,' alfa=', f7.4,' tb=',e9.4,' kstar=',
     %e9.4,' coef=',e10.5,2x,'(4T)**2.5 = ',e9.4,/,'4T= ',e10.4,' (1/4)T=
     % ',e10.4,2x,'sqrotfourt= ', f10.6)
	stop
	end