|
do i=1,521
READ(20,'(1440(I6))',end=99) (m(i,j),j=1,1440)
DO j = (1+3), (1440-3)
if(m(i,j)/=-999 .and. m(i,j+1)==-999) then
m(i,j)=2 ! 1 represents ocean
m(i,j-1)=2 ! 2 represents coastline
m(i,j-2)=2 ! -999 represents land
m(i,j-3)=2
else if(m(i,j)==-999 .and. m(i,j+1)/=-999) then
m(i,j)=2
m(i,j+1)=2
m(i,j+2)=2
m(i,j+3)=2
else if (m(i,j)<200) then
m(i,j)=2
else
m(i,j)=1
end if
END DO
write(30,'(1440(I6))') (m(i,j),j=1,1440)
end do |
|