diff --git a/src/gsi/read_anowbufr.f90 b/src/gsi/read_anowbufr.f90
index e2b744eb6a..1873d0b877 100644
--- a/src/gsi/read_anowbufr.f90
+++ b/src/gsi/read_anowbufr.f90
@@ -307,6 +307,7 @@ subroutine read_anowbufr(nread,ndata,nodata,gstime,&
 
            ndata=ndata+1
            nodata=nodata+1
+           if(ndata>maxobs) exit
 
            cdata_all(iconc,ndata)  = conc                    ! pm2_5 obs     
            cdata_all(ierror,ndata) = obserror                ! pm2_5 obs error
diff --git a/src/gsi/read_dbz_nc.f90 b/src/gsi/read_dbz_nc.f90
index 89eebde8b6..147c1c246c 100644
--- a/src/gsi/read_dbz_nc.f90
+++ b/src/gsi/read_dbz_nc.f90
@@ -411,6 +411,7 @@ subroutine read_dbz_nc(nread,ndata,nodata,infile,lunout,obstype,sis,hgtl_full,no
  
  !####################       Data thinning       ###################
         icntpnt=icntpnt+1
+        if(icntpnt>maxobs) exit
 
         if(ithin > 0)then
      
diff --git a/src/gsi/read_gmi.f90 b/src/gsi/read_gmi.f90
index 0f45aa7e28..f59529662a 100644
--- a/src/gsi/read_gmi.f90
+++ b/src/gsi/read_gmi.f90
@@ -528,6 +528,7 @@ subroutine read_gmi(mype,val_gmi,ithin,rmesh,jsatid,gstime,&
         flgch = 0
 
         iobs=iobs+1
+        if(iobs>maxobs) exit
      end do read_loop
   end do read_subset
 690 continue
diff --git a/src/gsi/read_goesglm.f90 b/src/gsi/read_goesglm.f90
index e0124abbf2..bf8639c72d 100644
--- a/src/gsi/read_goesglm.f90
+++ b/src/gsi/read_goesglm.f90
@@ -276,6 +276,7 @@ subroutine read_goesglm(nread,ndata,nodata,infile,obstype,lunout,twindin,sis)
         icntpnt=icntpnt+1
 
         ndata=ndata+1
+        if(ndata>maxobs) exit
         nodata=nodata+1
         iout=ndata
         isort(icntpnt)=iout
diff --git a/src/gsi/read_radar.f90 b/src/gsi/read_radar.f90
index 8e5de5aff9..aea19586bd 100644
--- a/src/gsi/read_radar.f90
+++ b/src/gsi/read_radar.f90
@@ -2911,6 +2911,7 @@ subroutine read_radar(nread,ndata,nodata,infile,lunout,obstype,twind,sis,hgtl_fu
 !####################       Data thinning       ###################
 
                     icntpnt=icntpnt+1
+                    if(icntpnt>maxobs) exit
 
                     if(ithin > 0)then
                        if(zflag == 0)then
@@ -4031,6 +4032,7 @@ subroutine read_radar_l2rw(ndata,nodata,lunout,obstype,sis,nobs,hgtl_full)
         end if
 !####################       Data thinning       ###################
         icntpnt=icntpnt+1
+        if(icntpnt>maxobs) exit
         ithin=1 !number of obs to keep per grid box
         if(radar_no_thinning) then
           ithin=-1
diff --git a/src/gsi/read_radar_wind_ascii.f90 b/src/gsi/read_radar_wind_ascii.f90
index 2e1b06a50c..604d9d0eca 100644
--- a/src/gsi/read_radar_wind_ascii.f90
+++ b/src/gsi/read_radar_wind_ascii.f90
@@ -509,6 +509,7 @@ subroutine read_radar_wind_ascii(nread,ndata,nodata,infile,lunout,obstype,sis,hg
 !####################       Data thinning       ###################
 
            icntpnt=icntpnt+1
+           if(icntpnt>maxobs) exit
  
            if(ithin > 0)then
               if(zflag == 0)then