program Lfix
      
      real sig1,soil1,can0_1,can45_1,plant1,
     &     inco1,coho1,cohop_1,cor1,corp_1
      real sig2,soil2,can0_2,can45_2,plant2,
     &     inco2,coho2,cohop_2,cor2,corp_2
      real sig3,soil3,can0_3,can45_3,plant3,
     &     inco3,coho3,cohop_3,cor3,corp_3
C      real sig4,soil4,can0_4,can45_4,plant4,
C     &     inco4,coho4,cohop_4,cor4,corp_4

      real degrees,norm
      real sig0,sig45,cor,coho,inco,soil

      read(5,*)sig1,soil1,can0_1,can45_1,plant1,
     &     inco1,coho1,cohop_1,cor1,corp_1
      read(5,*)sig2,soil2,can0_2,can45_2,plant2,
     &     inco2,coho2,cohop_2,cor2,corp_2
      read(5,*)sig3,soil3,can0_3,can45_3,plant3,
     &     inco3,coho3,cohop_3,cor3,corp_3
C      read(5,*)sig4,soil4,can0_4,can45_4,plant4,
C     &     inco4,coho4,cohop_4,cor4,corp_4



        sig1=10.0**(sig1/10.)
        sig2=10.0**(sig2/10.)
        sig3=10.0**(sig3/10.)
C        sig4=10.0**(sig4/10.)

        soil1=10.0**(soil1/10.)
        soil2=10.0**(soil2/10.)
        soil3=10.0**(soil3/10.)
C        soil4=10.0**(soil4/10.)

        can0_1=10.0**(can0_1/10.)
        can0_2=10.0**(can0_2/10.)
        can0_3=10.0**(can0_3/10.)
C        can0_4=10.0**(can0_4/10.)

        can45_1=10.0**(can45_1/10.)
        can45_2=10.0**(can45_2/10.)
        can45_3=10.0**(can45_3/10.)
C        can45_4=10.0**(can45_4/10.)
        
        plant1=10.0**(plant1/10.)
        plant2=10.0**(plant2/10.)
        plant3=10.0**(plant3/10.)
C        plant4=10.0**(plant4/10.)

        coho1=10.0**(coho1/10.)
        coho2=10.0**(coho2/10.)
        coho3=10.0**(coho3/10.)
C        coho4=10.0**(coho4/10.)

        inco1=10.0**(inco1/10.)
        inco2=10.0**(inco2/10.)
        inco3=10.0**(inco3/10.)
C        inco4=10.0**(inco4/10.)

        cor1=10.0**(cor1/10.)
        cor2=10.0**(cor2/10.)
        cor3=10.0**(cor3/10.)
C        cor4=10.0**(cor4/10.)

        norm=10.0**(-13.69/10.)

        degrees=4.0*atan(1.0)/180.

        corp_1=corp_1*degrees
        corp_2=corp_2*degrees        
        corp_3=corp_3*degrees
C        corp_4=corp_4*degrees

        cohop_1=cohop_1*degrees
        cohop_2=cohop_2*degrees
        cohop_2=cohop_3*degrees
C        cohop_3=cohop_4*degrees

        soil=10*alog10((soil1+soil2+soil3)/3.)
        
        sig0=10*alog10((sig1+sig2+sig3)/3.)
        
        sig45=10*alog10((can45_1+can45_2+can45_3)/3.
     &       +(soil1+soil2+soil3)/3.)

        inco=10*alog10(norm*(inco1+inco2+inco3)/3.
     &       +(soil1+soil2+soil3)/3.)
                
        coho=10*alog10(norm*(coho1*cos(cohop_1)+coho2*cos(cohop_2)
     &       +coho3*cos(cohop_3))/3.
     &       +norm*(inco1+inco2+inco3)/3.
     &       +(soil1+soil2+soil3)/3.)

        cor=10*alog10(norm*(coho1*cos(cohop_1)+coho2*cos(cohop_2)
     &       +coho3*cos(cohop_3))/3.
     &       +norm*(inco1+inco2+inco3)/3.
     &       +norm*(cor*cos(corp_1)+cor2*cos(corp_2)
     &       +cor3*cos(corp_3))/3.
     &       +(soil1+soil2+soil3)/3.)

        write(6,*) sig0,sig45,cor,coho,inco,soil
        
      stop
      end