forall is deprecated in Fortran 2018

This commit is contained in:
Martin Diehl 2019-09-19 14:24:13 -07:00
parent b41d6182a5
commit 08d2cdcb25
1 changed files with 6 additions and 2 deletions

View File

@ -736,13 +736,17 @@ pure function lattice_symmetrizeC66(struct,C66)
select case(struct)
case (LATTICE_iso_ID)
do k=1,3
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
do j=1,3
lattice_symmetrizeC66(k,j) = C66(1,2)
enddo
lattice_symmetrizeC66(k,k) = C66(1,1)
lattice_symmetrizeC66(k+3,k+3) = 0.5_pReal*(C66(1,1)-C66(1,2))
enddo
case (LATTICE_fcc_ID,LATTICE_bcc_ID)
do k=1,3
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
do j=1,3
lattice_symmetrizeC66(k,j) = C66(1,2)
enddo
lattice_symmetrizeC66(k,k) = C66(1,1)
lattice_symmetrizeC66(k+3,k+3) = C66(4,4)
enddo