LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added some comments to the fortran code.

Browse Source
This commit is contained in:
salman 2015-09-03 10:12:12 +02:00
parent 62cf1ff6b0
commit 6815c7f089
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
trimsp/src/trimspNL.F
View File

@ -389,6 +389,7 @@ C This part reads the input file (new format)
C First line: properties of projectile C First line: properties of projectile
READ(11,*) Z1,M1,E0,Esig,ALPHA,ALPHASIG,EF,ESB,SHEATH,ERC READ(11,*) Z1,M1,E0,Esig,ALPHA,ALPHASIG,EF,ESB,SHEATH,ERC
C Second line: simulation related parameters C Second line: simulation related parameters
C Ordered as: Number of particles, seed, seed, seed, initial depth, RD, depth increment, CA, KK0, KDEE1,KDEE2,IPOT
READ(11,*) NH,RI,RI2,RI3,X0,RD,CW,CA,KK0,KK0R,KDEE1,KDEE2,IPOT READ(11,*) NH,RI,RI2,RI3,X0,RD,CW,CA,KK0,KK0R,KDEE1,KDEE2,IPOT
& ,IPOTR,IRL & ,IPOTR,IRL
C Third line: Number of layers C Third line: Number of layers
@ -523,9 +524,11 @@ C layers. It should be replaced!
C Checks wether depth interval is an integer denominator of layer thickness or not C Checks wether depth interval is an integer denominator of layer thickness or not
C If not, calculated implantation profile is not correct. C If not, calculated implantation profile is not correct.
C CW is depth increment and DX(K) is the thickness of layer K
depth_interval_flag = 1 depth_interval_flag = 1
DO K=1,L-1 DO K=1,L-1
IF(.NOT.EQUAL(DX(K)/CW-DBLE(IDINT(DX(K)/CW)),0.D0)) THEN IF(.NOT.EQUAL(DX(K)/CW-DBLE(IDINT(DX(K)/CW)),0.D0)) THEN
C If the thickness of layer K is not an integer multiples of depth increment
depth_interval_flag = 0 depth_interval_flag = 0
GO TO 44 GO TO 44
ENDIF ENDIF
@ -2293,6 +2296,7 @@ C
enddo enddo
ima = 1 ima = 1
1010 ima = MIN0(ima+2,5000) 1010 ima = MIN0(ima+2,5000)
C I am not sure what is this file for?!
open(20,file='edist') open(20,file='edist')
do ne=1,ima do ne=1,ima
write(20,1020) ne,me(ne) write(20,1020) ne,me(ne)