157a158
>   real(r_kind),dimension(11) :: apriori,layeff
194,195c195
< !  write(490+mype,*)'top of setupozlay'
< !  call flush(490+mype)
---
> 
345c345,350
<           call omiadjust(ozges,data(ilate,i),dpoz,dp,sfcp)
---
>           do k=1,11
>             apriori(k)=data(14+k,i)
>             layeff(k)=data(25+k,i)
> !            write(530+mype,*)'do omiadjust apr',k,apriori(k),'eff',layeff(k)
>           end do
>           call omiadjust(ozges,data(ilate,i),dpoz,dp,sfcp,apriori,layeff)
662a668
> !  if(mype.eq.0)write(6,*)'find ajl ii',ii
696d701
< !  write(530+mype,*)'bottom setupozlay'
846a852,854
> !  real pobsmax,pobsmin
> !  pobsmax=-1.
> !  pobsmin=1.e10
907a916,917
> ! ajl muse are levels that data will be assimilated while luse mean either
> ! assimilated or only diagnostic level
940a951
> !     if(pobs(kk)>10. .and. muse(i))then
1063a1075,1078
> !     if(muse(i))then
> !        write(6,*)mype,'pobs',kk,pobs(kk),' j ',j,'i',i
> !        call flush(6)
> !     endif
1194a1210,1213
> !     if(muse(i))then
> !        if(pobs(kk)<pobsmin)pobsmin=pobs(kk)
> !        if(pobs(kk)>pobsmax)pobsmax=pobs(kk)
> !     endif
1195a1215,1218
> !     if(pobs(kk)>10. .and. muse(i))then
> !       write(6,*)mype,'i',i,'pobs',kk,pobs(kk),'lat',dlat,dlon,' j ',j,'muse',muse(i),'luse',luse(i)
> !       call flush(6)
> !     endif
1230a1254,1256
> !  if(mype.eq.0)print *,'pobsmax',pobsmax,pobsmin
> !  write(6,*)mype,'pobsmax',pobsmax,pobsmin
> !  call flush(6)
1232a1259
> !   ajl this stuff only used for diagnostics files
1238c1265
< !          write(100+mype,*)mype,'pobs',kk,pobs(kk)
---
> !          write(150+mype,*)mype,'pobs',kk,pobs(kk)
