|||
转自:http://bbs.keinsci.com/thread-9494-1-6.html
SUBROUTINE CONSTR_CELL_RELAX(FCELL)
USE prec
REAL(q) FCELL(3,3)
DO I=1,3
FCELL(3,I)=0
FCELL(I,3)=0
ENDDO
! just one simple example
! relaxation in x directions only
! SAVE=FCELL(1,1)
! FCELL=0 ! F90 style: set the whole array to zero
! FCELL(1,1)=SAVE
! relaxation in z direction only
! SAVE=FCELL(3,3)
! FCELL=0 ! F90 style: set the whole array to zero
! FCELL(3,3)=SAVE
RETURN
END SUBROUTINE
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-24 00:23
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社