!=================================================================
!  一次元火星大気ダストの放射計算プログラム;
!  δ-Eddington 近似
!
!  2000/10/17  小高正嗣
!=================================================================
      subroutine EDDIN
     I   ( QASM  , QALB  , COSZ,
     O     gamma1, gamma2, gamma3 )

!-----------------------------------------------------------------      
      implicit real*8 ( a-h,o-z )
!-----------------------------------------------------------------      

      gamma1 = 0.25D0 * ( 7.0D0 - ( 4.0D0 + 3.0D0 * QASM ) * QALB )
      gamma2 = - 0.25D0 * ( 1.0D0 - ( 4.0D0 - 3.0D0 * QASM ) * QALB )
      gamma3 = 0.25D0 * ( 2.0D0 - 3.0D0 * QASM * COSZ ) 


      return
      end
