!=================================================================
!  一次元火星大気ダストの放射計算プログラム;
!  ダスト光学パラメータの設定
!
!  2000/10/15  小高正嗣
!=================================================================
      subroutine SETDUST
     O   ( QEXT_SR, QEXT_IR1, QEXT_IR2, 
     O     QALB_SR, QALB_IR1, QALB_IR2, 
     O     QASM_SR, QASM_IR1, QASM_IR2,
     O     QTAU_ALL, QDENS  , 
     O     QREFF   , QRMOD              )

!-----------------------------------------------------------------      
      implicit real*8 ( a-h,o-z )
!-----------------------------------------------------------------      
! 消散効率(extinction efficiency)

      QEXT_SR  = 3.04D0 
!      QEXT_SR  = 2.74D0 
      QEXT_IR1 = 0.253D0 * QEXT_SR
      QEXT_IR2 = 0.166D0 * QEXT_SR

! 1 次散乱アルベド(single scattering albedo)

      QALB_SR  = 0.92D0
!      QALB_SR  = 0.86D0
      QALB_IR1 = 0.47D0
      QALB_IR2 = 0.37D0

! 非対称因子(asymmetoric factor)

      QASM_SR  = 0.55D0
!      QASM_SR  = 0.79D0
      QASM_IR1 = 0.528D0
      QASM_IR2 = 0.362D0

! 全光学的厚さ(total optical depth)

      QTAU_ALL = 0.0D0
!      QTAU_ALL = 1.0D0

! 物性パラメータ

      QDENS = 3000.0D0          ! dust density
      QREFF = 2.5D-6            ! effective radius
      QRMOD = 0.4D-6            ! mode radius

!-----------------------------------------------------------------      

      return
      end
