real mu_atm,m_keha
open(10,file="dim1.dat")
pi=3.1415926536
g_rask=9.814
r_keha=0.1
s_keha=pi*r_keha**2
c_keha=0.47
v_keha=4.*pi*r_keha**3/3
m_keha=5.
p_atm=100400.
t_atm=280.
mu_atm=29.
r_gaas=8314.
ro_atm_0=p_atm*mu_atm/(r_gaas*t_atm)
p1=g_rask*mu_atm/(r_gaas*t_atm)
p2=0.5*c_keha*s_keha/m_keha
p3=g_rask*v_keha/m_keha
nt=40000
dt=0.001
v=10.
x=0.
workt=0.
worka=0.

do it=1,nt
 t=dt*it
 roh=ro_atm_0*exp(-p1*x)
 at=-p2*roh*abs(v)*v
 aa=p3*roh
 a=-g_rask+aa+at
 dx=v*dt+0.5*a*dt**2
 x=x+dx
 v=v+a*dt
 ekin=0.5*m_keha*v**2
 epot=m_keha*g_rask*x
 etot=ekin+epot
 workt=workt+m_keha*at*dx
 worka=worka+m_keha*aa*dx

 if(x>0.)then
  write (10,'(15f15.5)')t,x,v,a,ekin,epot,etot,workt,worka,etot-(workt+worka)
 else
  stop
 endif
enddo

stop
end