6,7c6
<   real(r_kind), allocatable :: szraq(:,:),raqmsoz(:,:,:,:),raqmsno2(:,:,:,:)
<   real(r_kind), allocatable :: raqmsptrop(:,:)
---
>   real(r_kind), allocatable :: szraq(:,:),raqmsoz(:,:,:,:)
12d10
<   real(r_kind) :: latob,lonob
14c12
<   character *6 nameoz,nameno2
---
>   character *6 nameoz
18,19d15
<   data nameno2/'ino2'/
<   save nameno2
37,41d32
< ! add no2 pressure for scattering weight since are a constanst
<   integer(i_kind),parameter :: nlevscatwt=35
<   real(r_kind) :: presscatwt(nlevscatwt),pintscatwt(0:nlevscatwt),pintscatwtcb(0:nlevscatwt)
<   integer(i_kind) :: no2pe=-1
<   logical :: bcastpresscatwt=.true.
58,62c49,53
< !        !do k=1,nprs
< !          !write(6,*)'top',aprprt(k)
< !          write(6,*)k,aprprs(k),'apreff',apreff(k)
< !          write(6,*)'bot',aprprb(k)
< !        end do
---
>         do k=1,nprs
>           write(6,*)'top',aprprt(k)
>           write(6,*)k,aprprs(k),'apreff',apreff(k)
>           write(6,*)'bot',aprprb(k)
>         end do
75,77c66,68
< !      do k=1,nprs
< !        write(6,*)k,aprprs(k),(apriori(9,k,ic,3),ic=1,5)
< !      end do
---
>       do k=1,nprs
>         write(6,*)k,aprprs(k),(apriori(9,k,ic,3),ic=1,5)
>       end do
146,192c137
<   end subroutine readapriori
<   subroutine sendpresscatwt
<   use kinds, only: r_kind,r_single,i_kind
<   use mpimod, only : mpi_comm_world,mpi_rtype,mpi_integer4,npe
<   use mpimod, only : mype
<   use no2info, only : ihave_no2
<   implicit none
<   integer(i_kind)           ::        nno2pe(1:npe),ierr,k
<   if(.not.ihave_no2)return
< !   write(6,*)'sendpresscatwt'
< !   call flush(6)
<   if(bcastpresscatwt)then 
< !    write(6,*)mype,'bcastpresscatwt',no2pe
<     bcastpresscatwt=.false.
<     call mpi_gather(no2pe,1,mpi_integer4,nno2pe,1,mpi_integer4,0,mpi_comm_world,ierr)
< !    write(6,*)mype,'no2pe',no2pe
< !    if(mype.eq.0)then
< !      write(6,*)'nno2pe',nno2pe
< !      call flush(6)
< !    endif
<     
<     if(mype.eq.0.and.no2pe.eq.0)then
< !      write(6,*)'zero is no2pe'
<     else
<       if(mype.eq.0)then
<         do k=1,npe
<           if(nno2pe(k)>no2pe)then
< !            write(6,*)'nno2pe not neg ',k,nno2pe(k)
<             no2pe=nno2pe(k)
<           endif
<         end do
<       endif
<     endif
<     call mpi_bcast(no2pe,1,mpi_integer4,0,mpi_comm_world,ierr)
< !    write(6,*)mype,'no2pe',no2pe
<     call mpi_bcast(presscatwt,nlevscatwt,mpi_rtype,no2pe,mpi_comm_world,ierr)
< !    call mpi_bcast(pintscatwt,nlevscatwt+1,mpi_rtype,no2pe,mpi_comm_world,ierr)
<   endif
<   pintscatwt(0)=1100.
<   pintscatwt(nlevscatwt)=0.
<   do k=1,nlevscatwt-1
<     pintscatwt(k)=.5*(presscatwt(k)+presscatwt(k+1))
<   end do
<   pintscatwtcb=pintscatwt*.1
< !  write(840+mype,*)'end of sendpressscatwt'
< 
<   end subroutine sendpresscatwt
---
>     end subroutine readapriori
