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
| open(109,file = titre, status='unknown', form ='formatted')
open(1, file = 'Dat00210.vel',status='old',form='formatted')
read(1,*),nrtri
!------------- file 1 and vit1 ----------------------------------
write(109,*),nrtri
do i=1,nrtri
read(1,*) A,B
if ( abs(A) .gt.v1) then
A=res1
if (abs(B) .gt.v1) then
B=res1
endif
endif
write(109,*), A,B
enddo
write(109,*),'123d-4'
close(1)
close(109)
!------------file 2 and vit1 ------------------------------------
open(109,file = 'env.vel', status='old', form ='formatted'
& ,position = 'append')
open(1, file = 'Dat00220.vel',status='old',form='formatted')
read(1,*),nrtri
write(109,*),nrtri
do i=1,nrtri
read(1,*) A,B
if ( abs(A) .gt.v1) then
A=res1
if (abs(B) .gt.v1) then
B=res1
endif
endif
write(109,*), A,B
enddo
close(1)
close(109) |
Partager