On Jan 24, 2013, at 7:57 AM, Laurence wrote: > SOLVED!!!
Good! > I realised I had to create the constants that control my reactions > into a CellVariable which are therefore localised to the reaction area. Yes. If the reaction constant is both localized and constant, then the .setValue() approach works, e.g., kRx = CellVariable(..., value=0.) kRx.setValue(kRx0, where=(X > XA) & (Y < YB) ...) If, instead, you wanted a coefficient that depended on a threshold value of your concentration variable, you could write kRx = kRx0 * (C > C0) or you could combine that threshold with your spatial localization: kRx = kRx0 * ((C > C0) & (X > XA) & (Y < YB) ...) > Just need to work out what it all means :S That, unfortunately, lies outside FiPy's design specifications. _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
