better readable
This commit is contained in:
parent
f415a8cad2
commit
2fc66cff5b
|
@ -747,13 +747,13 @@ pure function lattice_symmetrizeC66(struct,C66)
|
||||||
case (LATTICE_iso_ID)
|
case (LATTICE_iso_ID)
|
||||||
do k=1,3
|
do k=1,3
|
||||||
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
|
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
|
||||||
lattice_symmetrizeC66(k,k) = C66(1,1)
|
lattice_symmetrizeC66(k,k) = C66(1,1)
|
||||||
lattice_symmetrizeC66(k+3,k+3) = 0.5_pReal*(C66(1,1)-C66(1,2))
|
lattice_symmetrizeC66(k+3,k+3) = 0.5_pReal*(C66(1,1)-C66(1,2))
|
||||||
enddo
|
enddo
|
||||||
case (LATTICE_fcc_ID,LATTICE_bcc_ID)
|
case (LATTICE_fcc_ID,LATTICE_bcc_ID)
|
||||||
do k=1,3
|
do k=1,3
|
||||||
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
|
forall(j=1:3) lattice_symmetrizeC66(k,j) = C66(1,2)
|
||||||
lattice_symmetrizeC66(k,k) = C66(1,1)
|
lattice_symmetrizeC66(k,k) = C66(1,1)
|
||||||
lattice_symmetrizeC66(k+3,k+3) = C66(4,4)
|
lattice_symmetrizeC66(k+3,k+3) = C66(4,4)
|
||||||
enddo
|
enddo
|
||||||
case (LATTICE_hex_ID)
|
case (LATTICE_hex_ID)
|
||||||
|
|
Loading…
Reference in New Issue