Changeset 481 for experiments

Timestamp:

03/26/16 00:26:52 (9 years ago)

Author:

oriona

Message:

Food seasonality added.

File:

• experiments/frams/foraminifera/data/scripts/foraminifera.expdef (modified) (9 diffs)

experiments/frams/foraminifera/data/scripts/foraminifera.expdef

@@ -31 +31 @@
 global o;
 global reprocounter;
+global changePeriod;
+global phase;
 
 @include "foraminifera.inc"
@@ -91 +93 @@
 
         //reproduction
-        ExpParams.foramPop = 4; 
+        ExpParams.foramPop = 20;        
         ExpParams.crossprob = 0;
         ExpParams.mutationprob = 0;
@@ -110 +112 @@
         ExpParams.chamber_difference_haplo = 0.0;
         ExpParams.chamber_proculus_diplo = micronsToFrams(20);
-        ExpParams.chamber_difference_diplo = 0.2;
+        ExpParams.chamber_difference_diplo = micronsToFrams(8);
 
         max_chamber_volume = [Vector.new(), Vector.new()];
@@ -129 +131 @@
         ExpParams.energy_move = 0.0000005;
 
-        ExpParams.energies0_haplo = energyFromVolume(ExpParams.chamber_proculus_haplo,1);
+        ExpParams.energies0_haplo = energyFromVolume(micronsToFrams(20),1);
         ExpParams.energies0_diplo = energyFromVolume(micronsToFrams(1.25),1);
         ExpParams.feedtrans = 0.001;
@@ -139 +141 @@
 
         //nutrients
-        ExpParams.foodperiod = 480;
+        changePeriod = 0;
+        phase = "high";
+        ExpParams.foodperiod = 21600;
+        ExpParams.foodPeriodChange = 0;
         ExpParams.nutrientradius = micronsToFrams(10);
         ExpParams.energy_nut = energyFromVolume(ExpParams.nutrientradius,1);
@@ -226 +231 @@
         {
                 addInitialForam(0, i);  
-                addInitialForam(1, i);
+                //addInitialForam(1, i);
         }
         o = Populations[0][0].getMechPart(0).orient.clone();
@@ -652 +657 @@
 function nutrientGrowth()
 {
+        if (ExpParams.foodPeriodChange > 0)
+        {
+                        changePeriod += 1;
+                        if (phase=="low" && (changePeriod*ExpParams.secPerStep) >= 23328000) //9 months
+                        {
+                                ExpParams.foodperiod = ExpParams.foodperiod/ExpParams.foodPeriodChange;
+                                phase = "high";
+                                changePeriod = 0;
+                        }
+                
+                        else if (phase == "high" && (changePeriod*ExpParams.secPerStep) >= 7776000) //3 months
+                        {
+                                ExpParams.foodperiod = ExpParams.foodperiod*ExpParams.foodPeriodChange;
+                                phase = "low";
+                                changePeriod = 0;
+                        }
+        }
         nutrientenergywaiting = nutrientenergywaiting + 1;
         if (nutrientenergywaiting*ExpParams.secPerStep >= ExpParams.foodperiod)
         {
-                for (var i = 0; i < ExpParams.feedrate*(ExpParams.secPerStep/ExpParams.foodperiod); i++)
+                for (var i = 0; i < ExpParams.feedrate; i++)
                 {    
                         addNutrient();
@@ -976 +998 @@
 
 prop:
+id:foodPeriodChange
+name:Set variable feed rate
+type:f 0
+group:Energy
+
+prop:
 id:ingestion
 name:Ingestion rate
@@ -996 +1024 @@
 id:foodperiod
 name:Time between food occurrences
-type:f 0 100000
+type:f 0 1000000
 group:Energy