/*
TLK Model Sample Input File: Radiosensitivity parameters correspond to an alpha/beta ratio of 3.61 Gy. DSB repair is biphasic (fast repair half-time = 0.25 h, slow repair half-time = 7 hour).  The entire cell population is actively dividing (exponential growth kinetics), i.e., GF=1.  In addition to the radiation damage, endogenous processes are presumed to create 4.33E-03 DSB/cell-hour.
*/
   MODEL: DRM=TLK CKM=QECK         !Model Specification
   CELL: DNA=5.667D+09 NC=46
   ENDO: DSB=4.3349E-03                   !Endogenous DNA damage
   IRAD: DSB=25 FCB=0.2                   !Radiation damage to DNA
   CLRP: RHT=HT4R ETA=2.00E-04 PHI=0.01 A0=AAA  GAM=0.25
   CCKM: TCC=10 TPOT=20 GF=1.0 N0=1 KAP=1.0D+38 VOL=1
   RBM: KX=6.0E+06 FMAX=1   !Analysis of PFGE data

 ! Genome instability and cell transformation
   GIMP: PGA=4.2868541E-05 PGH=0.1 NCG=1000 NCM=1
   NEO: TFS=1.1308766E-05 TFU=2.0E-03 TMAX=42 TSTP=1

!  TLK model inputs (simple and complex double strand breaks)
   HT4R = 0.25 7 /                        !Half-times for DSB repair
   AAA = 0.95 0.25 /                    !Probability DSB is rejoined correctly

!  Generate a synthetic dataset (add white noise to model-predicted responses).  The -n
!  command line option must be used to add noise to the predicted responses.
   NOISE: rerr=0.1 rngs=987654321

!  Parameters that impact on solution accuracy
   SOLVER: DBLV=0 ETOL=1.0E-09 MXSS=999999 HMAX=10

!  Parameters that control the time steps and stopping criteria
   SIMCON: FSDX=1 TSAX=5 TCUT=100 FRDL=1.0D+00 RCUT=1.0D-09

/*
   Define multiple split-dose exposure scenarios (see also the s3.inp).  In this sample
   input file, a baseline exposure scenario is defined that delivers a total
   dose of 10 Gy as two 5 Gy fractions.  The fractions are delivered using an
   instantaneous dose rate of 100 Gy/hour.  The baseline scenario is then rescaled to
   deliver 20 Gy in two 10 Gy fractions.  The second fraction is delivered starting at
   the times indicated using the RX2 data array.  NOTE: Compare to the sdrx1.inp file.
*/

   SDRX: ADR=100 TAD=RX1 TOD=RX2
   RX1= 10 20 /     !Total delivered dose (fraction size = TAD/2)
   RX2= 0 0.1 0.2 0.3 0.5 1 2 3 4 5 6 7 8 9 10 /   !Time of deliver (2nd dose)