From 6815c7f0892981ebd485caccf92a01bddd6ec66f Mon Sep 17 00:00:00 2001 From: Zaher Salman Date: Thu, 3 Sep 2015 10:12:12 +0200 Subject: [PATCH] Added some comments to the fortran code. --- trimsp/src/trimspNL.F | 4 ++++ 1 file changed, 4 insertions(+) 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)