/*
RMR Model Sample Input File: Radiosensitivity parameters correspond to an alpha/beta ratio of 11.43 Gy. The half-time for DSB repair is 1.0 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=RMR CKM=QECK !Model Specification
CELL: DNA=5.667D+09 NC=46
ENDO: DSB=4.3349E-03 !Endogenous DNA damage
IRAD: DSB=25 !Radiation damage to DNA
CLRP: RHT=1.0 ETA=1.50E-04 PHI= 0.01 A0=0.25 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
! 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=2 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
/*
In this example, the dose rate decreases exponentially with a half-life
of 1443.36 h (= half-life for 125I). The dose rate function is truncated
at 2886.72 h (2X the value of RHL). The baseline 1 Gy scenario is
rescaled so that the dose delivered in 2886.72 h equals 140,
145, or 150 Gy. A relatively small step tolerance (STOL) is used
so that the shape of the dose rate function is very similar to a
true exponential. Values for STOL in the range from about 0.05 to
1E-03 should be adequate for most biological simulations.
NOTE: See the table of selected radioisotope half-lives for other
radioisotopes.
*/
DECAY: TAD=1 RHL=1443.36 TCUT=2886.72 STOL=1.0D-02 SAD=RX1
RX1= 140 145 150 / !Array of delivered doses (Gy)
