diff --git a/trimsp/src/trimspNL.F b/trimsp/src/trimspNL.F index 83fcec1..13073fc 100644 --- a/trimsp/src/trimspNL.F +++ b/trimsp/src/trimspNL.F @@ -389,6 +389,7 @@ C This part reads the input file (new format) C First line: properties of projectile READ(11,*) Z1,M1,E0,Esig,ALPHA,ALPHASIG,EF,ESB,SHEATH,ERC 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 & ,IPOTR,IRL 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 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 DO K=1,L-1 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 GO TO 44 ENDIF @@ -2293,6 +2296,7 @@ C enddo ima = 1 1010 ima = MIN0(ima+2,5000) +C I am not sure what is this file for?! open(20,file='edist') do ne=1,ima write(20,1020) ne,me(ne)