diff --git a/sorc/prepobs_cqcbufr.fd/cqcbufr.f b/sorc/prepobs_cqcbufr.fd/cqcbufr.f
index 0a12ac6..14ca024 100644
--- a/sorc/prepobs_cqcbufr.fd/cqcbufr.f
+++ b/sorc/prepobs_cqcbufr.fd/cqcbufr.f
@@ -927,7 +927,7 @@ SUBROUTINE AREA(IA,X,Y)
 C$$$
 
       SUBROUTINE AUXLEVS(SKIP,SAME,WIND,SINGLE,ANYS)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
@@ -1127,7 +1127,7 @@ SUBROUTINE CHANGE(ANY,SINGLE,SEQLP)
       REAL(8) BMISS
 
       SAVE IEVOLD
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
      &                 NEV, ISF(NST), NLVM, NLVW
@@ -1691,7 +1691,7 @@ SUBROUTINE CHKTMP(A,C,R,ICK)
 C$$$
 
       SUBROUTINE CKPS
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
@@ -1790,7 +1790,7 @@ SUBROUTINE CKPS
 C$$$
 
       SUBROUTINE COMPER(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -2280,7 +2280,7 @@ SUBROUTINE DHOR
       REAL(8) BMISS
 
       SAVE IEVOLD
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       COMMON /MANRES/  ZIM(21,NST),TIM(21,NST),TDIM(21,NST),QIM(21,NST),
      &                 ZHM(21,NST),THM(21,NST),TDHM(21,NST),QHM(21,NST),
      &                 ZVM(21,NST),TVM(21,NST),TDVM(21,NST),QVM(21,NST),
@@ -2408,7 +2408,7 @@ SUBROUTINE DHOR
       
       SUBROUTINE DISTR(X,MSK,XLIM,XMSG,NX,N,NDIV,DDIV,
      &  NZERO,DZERO,NS,X1,SD,SK,XK)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       INTEGER N(23), MSK(NST)
       REAL X(NST),XLIM(2)
@@ -2517,7 +2517,7 @@ SUBROUTINE DISTR(X,MSK,XLIM,XMSG,NX,N,NDIV,DDIV,
 C$$$
 
       SUBROUTINE DMA(ANY,OBS)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
      &                 NEV, ISF(NST), NLVM, NLVW
@@ -2752,7 +2752,7 @@ SUBROUTINE DRCTSL(FAALL,RAALL,DOTPRD,NDIM,MAXDIM,NXXYY,NFT,LEV,IV)
 C$$$
 
       SUBROUTINE ERR123(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -3289,7 +3289,7 @@ SUBROUTINE ERR123(L,LM)
 C$$$
 
       SUBROUTINE ERR5(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -3583,7 +3583,7 @@ SUBROUTINE ERR5(L,LM)
 C$$$
 
       SUBROUTINE ERR710(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -4211,7 +4211,7 @@ SUBROUTINE ERR710(L,LM)
 C$$$
 
       SUBROUTINE ERRTYP
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
@@ -4569,7 +4569,7 @@ SUBROUTINE EVENTW(LUNIT,EVNSTR,NLV,OBS,QMS,RCS,IND,NEVN,QCPC)
 C$$$
 
       SUBROUTINE EVPROC(I1,I2,ICALL)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       COMMON /EVENTS/  STN(2000),    SEQN(2000),  ISCAN(2000),
      &                 LEVL(2000),   PRES(2000),  LTYP(2000),
@@ -4831,7 +4831,7 @@ SUBROUTINE EVPROC(I1,I2,ICALL)
 C$$$
 
       SUBROUTINE FILALL(SAME,WIND)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -5260,7 +5260,7 @@ SUBROUTINE FILL(OBS,NLEV)
 C
 C$$$
       SUBROUTINE FULVAL(WIND)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -5613,7 +5613,7 @@ SUBROUTINE FUZZY(P,XI,XV,XH,B,PER)
 C$$$
 
       SUBROUTINE GETINC(WIND)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -5799,7 +5799,7 @@ SUBROUTINE GETINC(WIND)
 C
 
       SUBROUTINE GETLEV(WIND)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -6194,7 +6194,7 @@ SUBROUTINE GETPS(PS,ZS,Z1,T1,P1)
 C$$$
 
       SUBROUTINE HOLES
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -6287,7 +6287,7 @@ SUBROUTINE HOLES
 C$$$
 
       SUBROUTINE HORRES
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -6540,7 +6540,7 @@ SUBROUTINE HORRES
 C$$$
 
       SUBROUTINE HSC
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -6809,7 +6809,7 @@ SUBROUTINE HSC
 C$$$
 
       SUBROUTINE INCR(ITIME)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -7037,7 +7037,7 @@ SUBROUTINE INCR(ITIME)
 C$$$
 
       SUBROUTINE INCRW
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -7128,7 +7128,7 @@ SUBROUTINE INCRW
 C
 C$$$
       SUBROUTINE INIT
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -7211,7 +7211,7 @@ SUBROUTINE INIT
 
       SUBROUTINE INPUT(START,ENDIN,SKIP,SAME,WIND,ITIME,USESQN)
       SAVE
-      PARAMETER (NST=999)      ! maximum number of stations
+      PARAMETER (NST=2499)      ! maximum number of stations
       PARAMETER (MVO=5)        ! p,T,z,q,Td
       PARAMETER (MLV=255)      ! number of possible levels
       PARAMETER (MEV=13)       ! number of possible programs/events
@@ -7786,7 +7786,7 @@ SUBROUTINE INPUT(START,ENDIN,SKIP,SAME,WIND,ITIME,USESQN)
 C$$$
 
       SUBROUTINE INPUT2
-      PARAMETER(NST=999)
+      PARAMETER(NST=2499)
 
       REAL(8) BMISS,RIT_8,HDR_8(10),UPA_8(10,255)
 
@@ -7990,7 +7990,7 @@ SUBROUTINE INPUT2
 C$$$
 
       SUBROUTINE ISGOOD
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -8236,7 +8236,7 @@ SUBROUTINE ISGOOD
 C
 C$$$
       SUBROUTINE ISOLAT(ID)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
 C
 C     COLLECT LIST OF ISOLATED STATIONS.
@@ -8297,7 +8297,7 @@ SUBROUTINE ISORT(IA,INDX,N)
 C     SORT IA ACCORDING TO THE ORDER SPECIFIED BY THE
 C     INDICES IN INDX.
 
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       DIMENSION IA(*), IKSP(NST)
       INTEGER INDX(*)
       DO J=1,N
@@ -8336,7 +8336,7 @@ SUBROUTINE ISORT(IA,INDX,N)
 C$$$
 
       SUBROUTINE LAPSE
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -8458,7 +8458,7 @@ SUBROUTINE LAPSE
 C$$$
 
       SUBROUTINE LEVTYPS
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -8761,7 +8761,7 @@ FUNCTION MANLEV(P)
 C
 C$$$
       SUBROUTINE MASEVN
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -9172,7 +9172,7 @@ FUNCTION NMANLV(P)
 C$$$
 
       SUBROUTINE NOBERR
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -9397,7 +9397,7 @@ SUBROUTINE NOBERR
 C$$$
 
       SUBROUTINE OBERR
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -9989,7 +9989,7 @@ SUBROUTINE PEVENTS
 C$$$
 
       SUBROUTINE POBERR
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -10257,7 +10257,7 @@ SUBROUTINE POBERR
 C$$$
       
       SUBROUTINE PRNTOUT(SEQLP,ICALL)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -10559,7 +10559,7 @@ SUBROUTINE PRNTOUT(SEQLP,ICALL)
 C$$$
  
       SUBROUTINE PRSTNS(WIND)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
      &                 NEV, ISF(NST), NLVM, NLVW
@@ -10727,7 +10727,7 @@ SUBROUTINE PRSTNS(WIND)
 C$$$
 C-----------------------------------------------------------------------
       SUBROUTINE QCOI(LDIM,IDIM,L0,IV0,NOB1,NOB2,IDH,OINC,HRES,HSTD,WTS)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -10757,7 +10757,8 @@ SUBROUTINE QCOI(LDIM,IDIM,L0,IV0,NOB1,NOB2,IDH,OINC,HRES,HSTD,WTS)
 
       DATA MAXDIM /4/
       DATA MINDIM /2/
-      DATA MAXOBS /1000/
+C-----make MAXOBS consistent with other MAXOBS across prepobs---------      
+      DATA MAXOBS /4500/ 
       DATA NFT    /1/
 
 C-----------------------------------------------------------------------
@@ -10982,7 +10983,7 @@ SUBROUTINE QCOI(LDIM,IDIM,L0,IV0,NOB1,NOB2,IDH,OINC,HRES,HSTD,WTS)
 C$$$
 
       SUBROUTINE RBLOCKS(SKIP,WIND,ENDIN)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -11292,7 +11293,7 @@ SUBROUTINE ROUND(VAL,P,IV)
 C
 C$$$
       SUBROUTINE SEARCH(LDIM,IDIM,L0,IV0,NOB1,NOB2,IDH,OINC,OG)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -11492,7 +11493,7 @@ SUBROUTINE SEARCH(LDIM,IDIM,L0,IV0,NOB1,NOB2,IDH,OINC,OG)
 C
 C$$$
       SUBROUTINE SETTMP
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -11752,7 +11753,7 @@ SUBROUTINE SHELL(V,IV,MAX,IREV)
 C$$$
 
       SUBROUTINE SIGERR(LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -12090,7 +12091,7 @@ SUBROUTINE SORT(RA,INDX,N)
 C     SORT RA ACCORDING TO THE ORDER SPECIFIED BY THE
 C     INDICES IN INDX.
 C
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       DIMENSION RA(*), WKSP(NST)
       INTEGER INDX(*)
@@ -12139,7 +12140,7 @@ SUBROUTINE SORT(RA,INDX,N)
 C$$$
 
       SUBROUTINE STAT(ITIME)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
 C
 C     CALCULATE STATISTICS:
@@ -13069,7 +13070,7 @@ SUBROUTINE STNCNT
 C$$$
 
       SUBROUTINE STYPE
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13212,7 +13213,7 @@ SUBROUTINE STYPE
 C$$$
  
       SUBROUTINE T120(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13319,7 +13320,7 @@ SUBROUTINE T120(L,LM)
 C$$$
 
       SUBROUTINE T121(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13394,7 +13395,7 @@ SUBROUTINE T121(L,LM)
 C$$$
 
       SUBROUTINE T130(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13528,7 +13529,7 @@ SUBROUTINE T130(L,LM)
 C$$$
 
       SUBROUTINE T140(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13612,7 +13613,7 @@ SUBROUTINE T140(L,LM)
 C$$$
 
       SUBROUTINE T220(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -13770,7 +13771,7 @@ SUBROUTINE T220(L,LM)
 C$$$
 
       SUBROUTINE T240(L,LM)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -14194,7 +14195,7 @@ SUBROUTINE T240(L,LM)
 C
 C$$$
       SUBROUTINE TMPCHK
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -14326,7 +14327,7 @@ SUBROUTINE TMPCHK
 C$$$
 
       SUBROUTINE VOI(ITIME)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -14796,7 +14797,7 @@ SUBROUTINE VSOLVE (A,B,NDIM,BAD,NFT,NS,MAXDIM)
 C         
 C$$$     
       SUBROUTINE VTPEVN(DOVTMP)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -14936,7 +14937,7 @@ SUBROUTINE VTPEVN(DOVTMP)
 C$$$
 
       SUBROUTINE WBLOCKS
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       COMMON /HEADER/  SID(NST), DHR(NST), XOB(NST), YOB(NST),
      &                 ELV(NST), SQN(NST), ITP(NST), NLV,
@@ -15039,7 +15040,7 @@ SUBROUTINE WBLOCKS
 C$$$
 
       SUBROUTINE WINDATZ(SAME)
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -15279,7 +15280,7 @@ SUBROUTINE WINDATZ(SAME)
 C
 C$$$
       SUBROUTINE WTISO
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       COMMON /ISO/ IDISO(NST),NUM,ISISO
       LOGICAL ISISO
       WRITE(6,600)
@@ -15413,7 +15414,7 @@ SUBROUTINE WTSTATS
 C$$$
 
       SUBROUTINE XHORRES
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
@@ -15505,7 +15506,7 @@ SUBROUTINE XHORRES
 C$$$
 
       SUBROUTINE ZDIF
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
 
       REAL(8) BMISS
 
diff --git a/sorc/prepobs_cqcbufr.fd/radcor.f b/sorc/prepobs_cqcbufr.fd/radcor.f
index 2c88733..b5d7f04 100644
--- a/sorc/prepobs_cqcbufr.fd/radcor.f
+++ b/sorc/prepobs_cqcbufr.fd/radcor.f
@@ -82,7 +82,7 @@
 C
 C$$$
       SUBROUTINE RADEVN
-      PARAMETER (NST=999)
+      PARAMETER (NST=2499)
       PARAMETER (NRID=800) ! max # of raob ids listed per inst. type
 
       REAL(8) BMISS