program count_grid
c
c=================================================
c
c read the country-code, and
c count the number of the grid
c
c country and country-code are from the table
c=================================================
c..........................................
c country_code.asc
c country_code
c
c country_code.list
c..........................................
character*28 country(10)
parameter( infile = 11, iofile = 12, nc = 720, nr= 360 )
parameter(infile2 = 13)
integer data(2000,4000),count(10),code(10),region(10)
c
c ------ read from country_code.asc file ------
open(infile, file= 'region_code3.asc', status='old')
c
do j=1,nr
read(infile,*) (data(i,j),i=1,nc)
end do
c
c ------ read from country-code-table file ------
open(infile2,file='region_code3', status='old')
c read(infile2,*)
do n =1, 10
read(infile2,200) code(n), country(n), region(n)
end do
200 format(i4,a28,i8)
close(infile2)
c
c ------ count grid ------
do n=1,10
count(n) =0
end do
do n=1,10
do j =1 ,nr
do i = 1, nc
if(data(i,j) .eq. code(n)) then
count(n) = count(n) +1
end if
end do
end do
end do
c
close(infile)
c
open(iofile, file='region_code3.txt', status='new')
do n =1, 10
write(iofile,300) code(n),country(n),region(n),count(n)
end do
300 format(i4,a28,i8,i8)
close(iofile)
end