mirror of
https://github.com/CHN-beta/nixos.git
synced 2024-10-24 20:18:56 +08:00
23 lines
441 B
Fortran
23 lines
441 B
Fortran
SUBROUTINE CONSTR_CELL_RELAX(FCELL)
|
|
USE prec
|
|
REAL(q) FCELL(3,3)
|
|
|
|
LOGICAL FILFLG
|
|
INTEGER ICELL(3,3)
|
|
INQUIRE(FILE='OPTCELL',EXIST=FILFLG)
|
|
IF (FILFLG) THEN
|
|
OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD')
|
|
DO J=1,3
|
|
READ(67,"(3I1)") (ICELL(I,J),I=1,3)
|
|
ENDDO
|
|
CLOSE(67)
|
|
DO J=1,3
|
|
DO I=1,3
|
|
IF (ICELL(I,J)==0) FCELL(I,J)=0.0
|
|
ENDDO
|
|
ENDDO
|
|
ENDIF
|
|
|
|
RETURN
|
|
END SUBROUTINE
|