  ! get co2 from a file
          open(87,file='ceq.txt')
           read(87,*) pco2a
           close(87)
	 write(6,*) pco2a
	    end
