696 |
|
|
697 |
|
!new stuff here |
698 |
|
|
699 |
< |
if (print_vac) then |
700 |
< |
vac_unit = next_unit() |
701 |
< |
vac_file = file_prefix(1:machdep_lnblnk(file_prefix))//".vac" |
702 |
< |
open (unit = vac_unit, status = 'replace', file = vac_file) |
699 |
> |
|
700 |
|
|
701 |
< |
write(vac_unit,*) n_vacancy |
702 |
< |
write(vac_unit,*) |
701 |
> |
if (print_vac) then |
702 |
> |
vac_unit = next_unit() |
703 |
> |
vac_file = file_prefix(1:machdep_lnblnk(file_prefix))//".vac" |
704 |
> |
open (unit = vac_unit, status = 'replace', file = vac_file) |
705 |
|
|
706 |
< |
endif |
706 |
> |
write(vac_unit,*) n_vacancy |
707 |
> |
write(vac_unit,*) |
708 |
|
|
709 |
+ |
endif |
710 |
+ |
|
711 |
|
do iz = start_index, nc, 1 |
712 |
|
do iy = start_index, nc, 1 |
713 |
|
do ix = start_index, nc, 1 |
759 |
|
#endif |
760 |
|
|
761 |
|
!new stuff here |
762 |
+ |
|
763 |
+ |
|
764 |
+ |
|
765 |
+ |
|
766 |
|
elseif (print_vac) then |
767 |
|
rx = box(1)*(rx - reference) |
768 |
|
ry = box(2)*(ry - reference) |