a = csvread('C:/Users/Mohammed/Desktop/SAM trials/SPT + WT/Design TES range.csv', {'skip'=1, 'numeric'=true});
b = csvread('C:/Users/Mohammed/Desktop/SAM trials/SPT + WT/Design SM range.csv', {'skip'=1, 'numeric'=true});
c = csvread('C:/Users/Mohammed/Desktop/SAM trials/SPT + WT/cut in speed.csv', {'skip'=1, 'numeric'=true});
d = csvread('C:/Users/Mohammed/Desktop/SAM trials/SPT + WT/cut out speed.csv', {'skip'=1, 'numeric'=true});

active_case = ('SPT');

for (w=0; w<2; w++)

	 {        for (x=0; x<2; x++)
	 
		 {    R = a[w][0];
              set ('tshours', R);
     
	          S = b[x][0];
              set ('solarm', S);
     
              simulate ();   
            
              SPT_output_array = get('system_pre_curtailment_kwac'); 
			  outln (SPT_output_array);
										  }
			                                    } 
active_case = ('WT');

for (y=0; y<2; y++)

	 {        for (z=0; z<2; z++)
	 
		 {    T = c[y][0];
              set ('wind_turbine_cutin', T);
              
	          U = d[z][0];
              set ('wind_turbine_cutout', U); 
     
			  simulate ();   
            
              WT_output_array = get ('system_pre_curtailment_kwac'); 
              outln (WT_output_array);
                                         }    
											 }            
                   
for ( i=0; i<4; i++ ) 
	
	  {       {total_output = SPT_output_array[i] + WT_output_array[i];     } 
      
			  active_case('WT + SPT'); 

			  set( 'spec_mode', 1 ); // sets mode to import hourly or subhourly profile from file

              set('energy_output_array', total_output);

              simulate();                               
            
              energy_output_array = get ('total_output');

              outln (energy_output_array);
              
                                              }