cA0 "c0*(0.4201+0.2*(sa1+sa2+sb+sc))" "Initial conc. A"
cB0 "c0*0.34 - 0.1*cA0" "Initial conc. B"
ra sqrt((x-0.5)^2+(y-0.5)^2) "radius a (from center)"
rb sqrt((x-0.25)^2+(y-0.97)^2) "radius b (from 1st smaller loci)"
rc sqrt((x-0.75)^2+(y-0.97)^2) "radius c (from 2nd smaller loci)"
sa1 exp(-(20*(ra-0.13*l))^4/l^4) "inner annular front"
sa2 exp(-(20*(ra-0.3*l))^4/l^4) "outer annular front"
sb exp(-(20*(rb-0.11*l))^4/l^4) "satellite b"
sc exp(-(20*(rc-0.11*l))^4/l^4) "satellite c"