diff --git a/fortran/trimspNL.F b/fortran/trimspNL.F
index 0f1ecba..99b49d5 100644
--- a/fortran/trimspNL.F
+++ b/fortran/trimspNL.F
@@ -387,6 +387,38 @@ C     part. refl. coeff. from Thomas et al.
       DATA  PRC/0.825D0,21.41D0,8.606D0,0.6425D0,1.907D0,1.927D0/
       DATA  number_in_layer/MAXNL*0/
 
+C=======================================================================
+C     Main run controls
+C
+C     NH       number of projectiles
+C     E0       initial projectile energy
+C     Esig     initial energy spread
+C     ALPHA    incident angle
+C     X0       initial depth/start position
+C     CW       depth-bin width
+C     L        number of target layers
+C     
+C     Projectile
+C     
+C     Z1       projectile atomic number
+C     M1       projectile mass
+C     
+C     Target layers
+C     
+C     DX(I)    thickness of layer I
+C     RHO(I)   density of layer I
+C     NJ(I)    number of elements in layer I
+C     ZT(I,J)  atomic number of element J in layer I
+C     MT(I,J)  mass of element J in layer I
+C     CO(I,J)  concentration fraction of element J in layer I
+C     
+C     Binding / displacement energies
+C     
+C     SBE(I,J) surface binding energy
+C     ED(I,J)  displacement energy
+C     BE(I,J)  bulk binding energy
+C=======================================================================
+      
       CALL getarg(1, filein)
 C     if there is no input argument take defailt file names
       if (filein.eq.'') then       
@@ -2337,7 +2369,7 @@ C
      &     ,':',A2,':',A2)
       WRITE(21,*)
       WRITE(21,10052)nh,(seconds_stop_total-seconds_start_total)
-10052 FORMAT(1x,' Simulation needed for ',I7,' muons ',I7,' seconds')
+10052 FORMAT(1x,' Simulation needed for ',I7,' projectiles ',I7,' seconds')
 
       WRITE(21,1402)innam
  1402 FORMAT(//30X,'* INPUT DATA *',5X,A12)