program felszinlegkor
c     Keszitette: Lazar Dora /ELTE/
c     ldora1989@gmail.com
c     2012.04.25.
      implicit none
      integer i
      integer zr
      real theta1,theta2,theta3
      real thetas1,thetas2,thetas3
      real c1,c2,c3
      real z0,k,rb
      real lamda,gamma,ro,row,cp
      real KS1,KS2,KS3
      real B1,B2,B3
      real K1,K2,K3
      real D1,D2,D3
      real Kef1,Kef2
      real pszi1,pszi2,pszi3
      real pszis1,pszis2,pszis3
      real Q1,Q2,Q3
      real Rn(574),G(574),S(574),Ur(574)
      real d(574),es(574),T(574)
      real ucs(574),ra(574),Ae(574),E0(574)
      integer ev(574),h(574),n1(574),n2(574)
      integer x1(574),o(574),p(574)
      real l(574),x2(574)
      real r(574),e(574),di(574),cs(574),v(574)
      real par(574),belso(574),oszto(574)
c40    format(6i4,i2,f8.1,f6.1,i6,3i5,i7,f6.1)

      open(10,file='M0210_BS.928')
      open(25,file='Qk.txt')
      open(30,file='parolgas.txt')

      do i=1,574
      read(10,*)ev(i),h(i),n1(i),n2(i),
     + x1(i),o(i),p(i),l(i),x2(i),
     + T(i),r(i),e(i),Ur(i),
     + di(i),S(i),cs(i),v(i)

      theta1=0.10
      theta2=0.08
      theta3=0.02

      c1=30.
      c2=3.5
      c3=2.3

      zr=20.
      z0=0.01
      k=0.378

      lamda=2.5*(10.**6.)
      gamma=0.65
      ro=1.2
      cp=1004.

      Rn(i)=-122+0.808*S(i)
      G(i)=0.15*Rn(i)

      d(i)=(lamda/461.)*(es(i)/(T(i)**2.))
      es(i)=6.11*exp((17.4*T(i))/(T(i)+239.16))

      ucs(i)=(k*Ur(i))/(log(zr/z0))
      ra(i)=(0.74*log(zr/z0))/(k*ucs(i))

      rb=c1+c2*(thetas1/theta1)**c3
      Ae(i)=Rn(i)-G(i)
      belso(i)=ro*cp*(d(i)/ra(i))
      oszto(i)=d(i)+gamma*(1.+rb/ra(i))
      E0(i)=(d(i)*Ae(i)+belso(i))/oszto(i)
      par(i)=E0(i)*lamda

      KS1=176.*(10.**(-6.))
      KS2=7.*(10.**(-6.))
      KS3=7.*(10.**(-6.))
      thetas1=0.395
      thetas2=0.451
      thetas3=0.451
      B1=4.05
      B2=5.39
      B3=5.39

      K1=KS1*(theta1/thetas1)**(2.*B1+3.)
      K2=KS2*(theta2/thetas2)**(2.*B2+3.)
      K3=KS3*(theta3/thetas3)**(2.*B3+3.)
      D1=0.1
      D2=0.2
      D3=0.3
      Kef1=(K1*D1+K2*D2)/(D1+D2)
      Kef2=(K2*D2+K3*D3)/(D2+D3)

      pszis1=-0.121
      pszis2=-0.478
      pszis3=-0.478
      pszi1=pszis1*((theta1/thetas1)**((-1.)*B1))
      pszi2=pszis2*((theta2/thetas2)**((-1.)*B2))
      pszi3=pszis3*((theta3/thetas3)**((-1.)*B3))

      row=1000.

      Q1=row*Kef1*(2.*((pszi1-pszi2)/(D1+D2))+1.)
      Q2=row*Kef2*(2.*((pszi2-pszi3)/(D2+D3))+1.)
      Q3=row*KS3*((theta3/thetas3)**(2.*B3+3.))

      write(*,*)'Q1: ',Q1,'Q2: ',Q2,'Q3: ',Q3
      write(25,*)Q1,' ',Q2,' ',Q3
      write(30,*)par(i)
      enddo
      close(10)
      close(25)
      close(30)
      stop