c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc     Calculate q_photo                                    ccc
c      cc                   cqp.f                                ccc
c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       subroutine cqp  ( qp , rhom, lamb )  
c
        implicit real*8 ( a-h, o-z )
        real*8 lamb
        dimension qp( 1 ), rhom( 1 ), lamb( 1 )
        common /param/ ndat, n
c
        do 100 i = 1 , ndat
           qp( i ) = 0.666d0*( 1.d0 + 2.045d0*lamb( i ) )**(-2.066d0)
     &            *( 1.d0 + rhom( i )*( 1.875d8*lamb( i )
     1                                + 1.653d8*lamb( i )**2.d0 
     &                                + 8.499d8*lamb( i )**3.d0
     2                                - 1.604d8*lamb( i )**4.d0
     &                                ) **(-1.d0) )**(-1.d0)
 100    continue
        return
        end