C	VAPOR PREESSURE CURVE FOR NORMAL DEUTERIUM
C	 Imperical equation from NBS report 9276 (Rolf Prydz)
C	M. Farkhondeh 6-7-90
c	EJB 2-99 get rid of Topdrawer stuff
c       EJB 5-00 make equivalent file for H2.
c
c
c       log_10(P) = 3.0683 - 55.256/T - 3.1282e-02*T + 6.6989e-04*T^2
c
	REAL*8 A1,A2,A3,A4,P,T,power
	real*8 power_s1,power_s2,p_s1,p_s2
C
	OPEN(UNIT=1,FILE='h2vapor.dat',TYPE='NEW')
	A1=-55.256
	A2=3.0683
	A3=-3.1282e-02
	A4=6.6989e-04
C
	do t=14.0,33.4,0.2
	   POWER = A1/T +A2 +A3*T +A4*T*T
	   power_s1 = a1/(t+1.) + a2 + a3*(t+1.) + a4*(t+1)*(t+1)
	   power_s2 = a1/(t+2.) + a2 + a3*(t+2.) + a4*(t+2)*(t+2)
	   P = 10**(POWER)*14.75 !(psi)
	   P_s1 = 10**(power_s1)*14.695 !(psia)
	   P_s2 = 10**(power_s2)*14.695 !(psia)
	   WRITE(1,100) T,P,p_s1,p_s2
 100	   FORMAT(1X,4F12.2)
	enddo
	close (unit=1)

	STOP
	END