1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
|
integer,parameter :: nzonal = 43200
integer,parameter :: nmerid = 21600
integer,parameter :: nland = 2
character(len=264) :: fnin
character(len=264) :: fnout
! inout file stuff
integer,parameter :: iunit = 10
integer,parameter :: junit = 20
integer :: ios,iopen
! misc
integer,dimension(nland) :: land
integer(kind=1), dimension(nzonal) :: indictot
fnin="file.txt"
fnout="out.txt"
open (iunit,file=trim(fnin),iostat=ios,status='old',action='read')
if (ios /= 0) then
stop 'Problem opening input file'
end if
open (junit,file=trim(fnout))
do i=1,nmerid
if(mod(i,100).eq.0) print *,int(100.*float(i)/float(nmerid)),"%"
indictot=0
read(iunit,*) indictot(1:nzonal)
do j=1,nzonal
if(indictot(j).ne.1) indictot(j) = 0
enddo
write(junit,*) (indictot(j),j=1,nzonal)
enddo
close(junit)
close(iunit) |
Partager