9 parameter(maxa=16000000)
15 dimension isat(0:maxs,0:maxs)
16 real(8) arr(2,maxa),said(maxa),siid(maxa)
30 open(lunbf,file=file,form=
'unformatted')
32 CALL
openbf(lunbf,
'IN',lunbf)
34 call
ufbtab(lunbf,said,1,maxa,nret,
'SAID')
36 call
ufbtab(lunbf,siid,1,maxa,nrex,
'SIID')
44 isat(i,j) = isat(i,j)+1
50 1
format(a14,2x,a14,2x,a8)
53 if(isat(i,j).gt.0)
then
54 call satcode(i,ci,j,cj)
55 if(ci==
' ')
write(ci,
'(i4.4)')i
56 if(cj==
' ')
write(cj,
'(i4.4)')j
57 print
'(i3.3,2x,a10,2x,i10,2x,a80)',i,ci,isat(i,j),cj
67 subroutine satcode(icode,csad,jcode,csid)
69 character(16) csad,saic(1000)
70 character(80) csid,siic(1000)
72 logical first /.true./
79 read(8,*,end=2) said,saic(said)
83 read(9,4,end=3) siid,siic(siid)
subroutine openbf(LUNIT, IO, LUNDX)
This subroutine connects a new file to the BUFRLIB software for input or output operations.
subroutine ufbtab(LUNIN, TAB, I1, I2, IRET, STR)
This subroutine reads through every data subset in a BUFR file and returns one or more specified data...