Tutorial on the interaction function details
The aim is to illustrate how to program the interaction function to have two submodels. In this example, the first grows isotropically and the second anisotropically where the growth axis is specified by a gradient of polariser.
1) Create and save a new mesh exactly as for the basic interaction function. Next, add two new morphogens, id_a and id_b.
Having saved the project (here we call it GPT_tut_interaction_20110530), click Panel: Interaction function:Edit to automatically create the new interaction function. It is shown here.
There are lots of comments shown in green. These are designed to both help us structure our thoughts on the growth model and provide reminders of useful commands and syntax without having to dive into a manual.
The comments are easily removed - but for the moment we will leave them in.
2) Add the code illustrating why we chose Matlab into the first model. This is shown in Magenta in GPT_tut_interaction_20110530 second edit. Note that it is in the space left for MODEL 1 and that MODEL 1 has been renamed 'NoPolariser'.
3) Add the second model - new code is shown in MediumBlue. Notice that the polariser is set up in the same way for both models - but in the NoPolariser (first) model the growth rates parallel and perpendicular to the polariser are identical so the polariser gradient can have no effect. The final interaction function will look like GPT_tut_interaction_20110530 second edit.
Pattern of A (blue) and B (green) morphogens at the start. Growth will be promoted by A but partially inhibited by B. |
Initial shape. |
Isotropic growth, i.e. specified growth parallel to and perpendicular to any polariser is identical (Kapar=Kbpar=Kaper=Kbper). | Anisotropic growth where specified growth perpendicular to the polariser (arrows) is 0.1 that parallel to the polariser. |