long_photo_rate.f90

module long_photo_rate
	
  use precision, only: dp
  use input, only: level_py, use_which_source
  use array, only: long_source_HI_density_array, long_source_HeI_density_array, long_source_HeII_density_array, &
                   long_source_HI_column_density_in, long_source_HeI_column_density_in, long_source_HeII_column_density_in, &
                   long_source_HI_column_density_out, long_source_HeI_column_density_out, long_source_HeII_column_density_out, &
                   long_source_HI_photoionization_rate_array, long_source_HeI_photoionization_rate_array, &			   
                   long_source_HeII_photoionization_rate_array, long_source_photoheating_rate_array, &				   
                   long_source_shell_volume
  use radiation, only: photoion_shell, photrates

contains
	
  subroutine long_photo_rate_calculation()
		
    type(photrates) :: phi	
    integer :: i,j,k
		
    do i = -level_py,level_py
      do j = -level_py,level_py
        do k = -level_py,level_py
				
			call photoion_shell (phi, &
			                     long_source_HI_column_density_in(i,j,k),long_source_HI_column_density_out(i,j,k),&
			                     long_source_HeI_column_density_in(i,j,k),long_source_HeI_column_density_out(i,j,k),&
			                     long_source_HeII_column_density_in(i,j,k),long_source_HeII_column_density_out(i,j,k),&
			                     long_source_shell_volume(i,j,k),use_which_source)
long_source_HI_photoionization_rate_array(i,j,k) = phi%photo_cell_HI/long_source_HI_density_array(i,j,k)
long_source_HeI_photoionization_rate_array(i,j,k) = phi%photo_cell_HeI/long_source_HeI_density_array(i,j,k)		  
long_source_HeII_photoionization_rate_array(i,j,k) = phi%photo_cell_HeII/long_source_HeII_density_array(i,j,k)		  
long_source_photoheating_rate_array(i,j,k) = phi%heat	 	
		  							  
        enddo
      enddo
    enddo

  end subroutine long_photo_rate_calculation
  	
end module long_photo_rate