145,146d144
<   use gridmod, only : raqmsadjust
<   use mpimod, only : mype
168d165
<   real(r_kind),dimension(nsig):: val_lay_raqms
337,355d333
<        if(raqmsadjust)then
<          val1= zero
<          do kk=nsig,1,-1
<            val_lay_raqms(kk)= zero_quad
<            w1=ozptr%wij(1,kk)
<            w2=ozptr%wij(2,kk)
<            w3=ozptr%wij(3,kk)
<            w4=ozptr%wij(4,kk)
<            val_lay_raqms(kk)=val_lay_raqms(kk) + &
<                 w1* soz(j1,kk)+ &
<                 w2* soz(j2,kk)+ &
<                 w3* soz(j3,kk)+ &
<                 w4* soz(j4,kk)
< !           if(mype.eq.0)then
< !             write(6,*)'kk intoz',kk,'val_layraqnms',val_lay_raqms(kk),'eff',ozptr%efficiency(kk)
< !           endif
<            val1=val1+val_lay_raqms(kk)*ozptr%efficiency(kk)
<           end do
<        else
383,385d360
< !           if(mype.eq.0)then
< !               write(6,*)'vallyagsichem',j,val_lay(j),'eff',ozptr%efficiency(j)
< !           endif
388d362
<       endif
453,481d426
<          if(raqmsadjust)then
<            do kk=nsig,1,-1
<               w1=ozptr%wij(1,kk)
<               w2=ozptr%wij(2,kk)
<               w3=ozptr%wij(3,kk)
<               w4=ozptr%wij(4,kk)
<               roz(j1,kk)  = roz(j1,kk) + valx*w1*ozptr%efficiency(kk)
<               roz(j2,kk)  = roz(j2,kk) + valx*w2*ozptr%efficiency(kk)
<               roz(j3,kk)  = roz(j3,kk) + valx*w3*ozptr%efficiency(kk)
<               roz(j4,kk)  = roz(j4,kk) + valx*w4*ozptr%efficiency(kk)
<            enddo
<            if (l_foto) then
<               do kk=nsig,1,-1
<                  w1=ozptr%wij(1,kk)
<                  w2=ozptr%wij(2,kk)
<                  w3=ozptr%wij(3,kk)
<                  w4=ozptr%wij(4,kk)
<                  j1x=w1+(kk-1)*lat2*lon2
<                  j2x=w2+(kk-1)*lat2*lon2
<                  j3x=w3+(kk-1)*lat2*lon2
<                  j4x=w4+(kk-1)*lat2*lon2
<                  dhat_dt_oz(j1x) =dhat_dt_oz(j1x) + valx*w1*time_oz*ozptr%efficiency(kk)
<                  dhat_dt_oz(j2x) =dhat_dt_oz(j2x) + valx*w2*time_oz*ozptr%efficiency(kk)
<                  dhat_dt_oz(j3x) =dhat_dt_oz(j3x) + valx*w3*time_oz*ozptr%efficiency(kk)
<                  dhat_dt_oz(j4x) =dhat_dt_oz(j4x) + valx*w4*time_oz*ozptr%efficiency(kk)
<               enddo
<            endif
<          else
<            dz1=nsig+1
486,487d430
< ! ajl bug fix 12/16/2015 need to redfine dz1 since was used in loop above
< ! else the do kl=1 loop gets skiped since dz1 will be like 1
493d435
< !              write(700+mype,*)'kl ',kl,' eff ',ozptr%efficiency(kl), 'iz1 ',iz1,iz2
502d443
< !                 write(700+mype,*)'kk ',kk,' add ',val_lay(kl)*w1*delz,' delz ',delz
510d450
<          endif
