         do k=1,navgkernel
           ! JDE believe pressavgkernel is bottom pressures, see trantropomi.no2.f90 and read_no2.f90
           pressavgkernel(k)=data(ipressavgkernel+k,i)*.01 ! make mb
           presscbavgkernel(k)=data(ipressavgkernel+k,i)*.001 ! make cb
           avgkernel(k)=data(iavgkernel+k,i)*data(iamftot,i) ! make avgkernel to use to with slant column
         pressavgkernel(navgkernel+1)=0.0 ! top pressure, pressavgkernel is now interfaces
         presscbavgkernel(navgkernel+1)=0.0
           do k=2,min(navgkernel,lvltrop) ! JDE loop to set bottom pres level used 
             if(pressavgkernel(k-1)*.1>=prsitmp(1).and.pressavgkernel(k)*.1<prsitmp(1))then
           ptropobs=pressavgkernel(lvltrop+1)
             if (presscbavgkernel(kbot+nz-2)*psi < one )then 
               no2p_omi(nz)=presscbavgkernel(kbot+nz-2)
             no2ges(1)=no2ges(1)+no2ges1(k)*avgkernel(ktop-k+1)
                       no2tail(ibin)%head%avgkernel(nlayers), &
                  no2tail(ibin)%head%avgkernel(k)=avgkernel(ktop-k+1)
