diff --git a/src/gsi/read_anowbufr.f90 b/src/gsi/read_anowbufr.f90 index 59ed0ffcb3..86c454a745 100644 --- a/src/gsi/read_anowbufr.f90 +++ b/src/gsi/read_anowbufr.f90 @@ -232,6 +232,8 @@ subroutine read_anowbufr(nread,ndata,nodata,gstime,& conc=indata(ncopopm) if ( iret > 0 .and. (conc < conc_missing ) .and. & (conc >= zero)) then + if(indata(nxob) > r360)cycle + if(indata(nyob) > 0.5_r_kind*r360)cycle if(indata(nxob) >= r360) indata(nxob) = indata(nxob) - r360 if(indata(nxob) < zero) indata(nxob) = indata(nxob) + r360 diff --git a/src/gsi/setuppm2_5.f90 b/src/gsi/setuppm2_5.f90 index ad940cce78..257bc44b01 100644 --- a/src/gsi/setuppm2_5.f90 +++ b/src/gsi/setuppm2_5.f90 @@ -780,8 +780,6 @@ subroutine setuppm2_5(obsLL,odiagLL,lunin,mype,nreal,nobs,isis,is,conv_diagsave) call tintrp2a11(pm25wc(:,:,:,2,nfldsig),pm25wc_ges(2),dlat,dlon,dtime,hrdifsig,& mype,nfldsig) if (pm25wc_ges(1) >= pm2_5_bg_threshold) then - pm25wc_ges(1)=1.0_r_kind - else pm25wc_ges(1)=0.0_r_kind end if if (pm25wc_ges(2) >= pm2_5_bg_threshold) then