Example 4.4 (Quarantine as independent concern): Model 4 SEIRS MultiSpecies Spatial Model with quarantine strategy

Refers to the model:

KendrickModel SEIRS
    attribute: #(status -> S E I R);
    parameters: #(beta lambda gamma mu1 mu2 sigma nu);
    lambda: #(beta*I/N);
    transitions: #(
        S -- lambda --> E.
        E -- sigma --> I.
        I -- gamma --> R.
        R -- nu --> S.
        status -- mu2 --> Empty.
        Empty -- mu1 --> S.).

KendrickModel Quarantine
	attribute: #(qState -> F Q);
	parameters: #(delta epsilon);
	transitions: #(
		F -- delta --> Q.
		Q -- epsilon --> F.).
		
KendrickModel MultiSpecies
    attribute: #(species -> human bird).

Map IndoChina
    for: #(country -> Thailand Laos Vietnam Cambodia MyanmarBurma);
    borders: #(
        #(0 1 0 1 1)
        #(1 0 1 1 1)
        #(0 1 0 1 0)
        #(1 1 1 0 0)
        #(1 1 0 0 0)).

KendrickModel Spatial
    maps: 'IndoChina';
    transitionRate: #rho.

Composition Influenza
    model: 'SEIRS';
	 model: 'Quarantine';
    model: 'MultiSpecies';
    model: 'Spatial'.

Scenario Scr1
    on: 'Influenza';
	 populationSize: 27500;
	 N: #(country_species);
	 beta_species_qState: #(
	 	#(#(0 0.21) #(0 0))
		#(#(0 0.42) #(0 0)));
	 lambda: #(beta*I_F_country/(N-I_Q_country) sum);
	 gamma_species_qState: #(
	 	#(0.25 0.096)
		#(0.233 0.082));
    sigma_species_qState: #(
		#(0.5 0)
		#(0.67 0));
	 mu1_species_qState: #(
		#(0.0000365 0)
		#(0.00137 0));
	 mu2_species: #(0.0000365 0.00137);
	 nu_qState: #(0.00274 0);
	 delta_species_status: #(
		#(0 0 0.068 0)
		#(0 0 0.055 0));
	 epsilon_status: #(0 0 0 1);
	 rho_species_qState: #(
		#(0.03 0) 
		#(0.1 0));
	 S_species_country_F: #(
		#(500 500 500 500 500) 
		#(5000 5000 5000 4990 5000));
	 I_species_country_F: #(
	 	#(0 0 0 0 0) 
		#(0 0 0 10 0)).

Simulation InfluenzaSim rungeKutta
    scenario: 'Scr1';
    from: 0;
    to: 500;
    step: 1.

Visualization InfluenzaVizQ diagram
    for: 'InfluenzaSim';
    legends: #('humans' 'birds');
    data: #(I_F_species);
    legendTitle: 'Infectious';
    xLabel: 'Time (days)';
    exportToPng.