#AMPL dat file for ECOMOD-14-746.
#Reindeer management and winter pastures in the presence of supplementary feeding and government subsidies
#May 22, 2015
#For Knitro 7.0.0 and 9.1.0
#Antti-Juhani Pekkarinen, Olli Tahvonen, Jouko Kumpula


param maxt := 100;		
param r:=0.0;			
param p:=8;		
param fn:=12;			
param mn:=8;			
param K:=1000000;   	#Change K to 10000000 when computing results without supplementary feeding with knitro 7.0.0 (otherwise use K:=1000000;)    			
param gamma:=0.48;			
param u:=0.48;							
param alpha:=0.8;   		
param sigmam:=1.1; 			
param CL:=1.14;					
param Cx:=39.54;			
param Cs:=13.35;				
param U:=0.5;
param z0:=250;	

param AQ:=0;
param Aomf:=1000;
param Aybf:=0;
param Amoh:=0;

param q:=10;
param CV:=0.4;
param supplementaryfeeding:=1;		

param wfs:=
       0 46
       1 63
       2 70
       3 76
       4 78
       5 80
       6 82
       7 82
       8 82
       9 82
       10 80
       11 80
       12 80;


param wms:=
       0 50
       1 74
       2 91
       3 108
       4 115
       5 120
       6 122
       7 122
       8 120;




param f:=                      
      0 0
      1 0
      2 0.6
      3 0.85
      4 0.9
      5 0.92
      6 0.92
      7 0.92
      8 0.93
      9 0.93
      10 0.85
      11 0.75
      12 0.6;


param fm:=                      
      0 0
      1 1
      2 7	
      3 15	
      4 20	
      5 22
      6 20
      7 18
      8 15;
   
param	xf0:=	
1	30.1862404	
2	30.178128	
3	30.178128	
4	30.1700156	
5	30.1537908	
6	30.1537908	
7	30.137566	
8	30.137566	
9	30.137566	
10	0.004913161	
11	0.002195554	
12	0.0008951	;
		
param	xm0:=	
1	3.7154792	
2	3.7154792	
3	3.7154792	
4	3.7073668	
5	3.7073668	
6	0.000123502	
7	7.86586E-05	
8	1.51499E-05	;
   

param pxf0:=
1	30.1862404	
2	30.178128	
3	30.178128	
4	30.1700156	
5	30.1537908	
6	30.1537908	
7	30.137566	
8	30.137566	
9	30.137566	
10	0.004913161	
11	0.002195554	
12	0.0008951	;

param pxm0:=
1	3.7154792	
2	3.7154792	
3	3.7154792	
4	3.7073668	
5	3.7073668	
6	0.000123502	
7	7.86586E-05	
8	1.51499E-05	;

param wc:=
      1 4.7
      2 4.7
      3 5
      4 5.4
      5 6
      6 6.2
      7 6.2
      8 6.2
      9 6.2
     10 6
     11 5.8
     12 5.6; 

# Parameters for type of the pasture rotation

param dwi:=181;					# number of the winter days when reindeer are in lichen pastures 
param dsp:=0;#31;				# number of the spring days when reindeer are in lichen pastures (0=with pasture rotation, 31=without pasture rotation) 
param dsu:=0;#92;				# number of the summer days when reindeer are in lichen pastures (0=with pasture rotation, 92=without pasture rotation) 
param dau:=0;#61;				# number of the autumn days when reindeer are in lichen pastures (0=with pasture rotation, 61=without pasture rotation) 

param wwi:=1.3;	
param wsp:=1.6;	
param wsu:=3;	
param wau:=1.6;	

param  tausp:=0.8;
param  tausu:=0.3;	
param  tauau:=0.6;