%% Starvation study for Sardine %Clear workspace clc clear all % Deleting the previous file %delete('D:/UCT/PhD/Model output/Sardinestarve.csv') dbstop if error globalvariables_Sardine; Varmat_Sardine; SO_Dev(1) = 0; % Final Oocyte development time %Make it between 11 and 18 SO_DevTime1 = 18; SO_DevTime2 = 11; Starve(1) = 0; SpM3(1) = 0; %Test variables mydata_Sardinops_sagax %Parameters_Sardine %Pre-assign matrix Sam = 10; %number of years the fish live %Froese et al. 2010 j = 1; %number of times fish feed M = 365; %number of days in a year N = (M*Sam*j); %number of days * number of years * number of feedings fn = 8; %number of values for f Vn = 38; %number of variables in matrix n = (N-1)*fn; %Max number of rows in matrix SNC = zeros (Vn,n); %Figure_Hold % Setting up the Matrix for saving Sardine_label; fname = 'Sardinestarve0.8_20.csv'; txt=sprintf('%s,', hdrS{:}); txt(end)=''; dlmwrite(['D:\UCT\PhD\Model output\',fname],txt,''); Parameters_Sardine Startvalue_Sardine for spawndays = [1887 2040 2162 2252 2405 2527 2618 2771 2893 2983 3136 3258 3348 3501 3623]; for d = 2:1:N; %2191 if they eat once a day T1D(d) = 20; %Reference Temperature TK(d) = T1D(d)+273; Sparameter_calc % Peak spawning between August and March % Atresia Sept & March % Starving start of peak (8th month) % Starving middle of peak (12th month) % Starving end of peak (3rd month) % Temp = 26deg First spawn 290 days % 1st % middle 336 days % Temp = 24deg First spawn 360 days % end 426 days % Temp = 22deg first spawn 448 days % 2nd % f = 1.0 & Temp = 20deg First spawn 559 days % start 579 days % Temp = 18 deg first spawn 699 days % middle 701 days % f = 0.9 first spawn 744 days % end 791 days % Temp = 16 deg first spawn 879 days % 3rd % start 944 days % middle 1066 days % Temp = 14 deg first spawn 1108 days % end 1157 days % 4th % start 1310 days % Temp = 12 deg first spawn 1401 days % middle 1432 days % end 1522 days % 5th % start 1675 days % middle 1797 % f = 0.8 first spawn 1843 days % end 1887 days % 6th % start 2040 days % middle 2162 % end 2252 days % 7th % start 2405 days % middle 2527 days % end 2618 days % 8th % start 2771 days % middle 2893 days % end 2983 days % 9th % start 3136 days % middle 3258 days % end 3348 days % 10th % start 3501 days % middle 3623 days if spawndays == 354; Spawn = 1.25; elseif spawndays == 426; Spawn = 1.5; elseif spawndays == 579; Spawn = 2; elseif spawndays == 701; Spawn = 2.25; elseif spawndays == 791; Spawn = 2.5; elseif spawndays == 944; Spawn = 3; elseif spawndays == 1066; Spawn = 3.25; elseif spawndays == 1157; Spawn = 3.5; elseif spawndays == 1310; Spawn = 4; elseif spawndays == 1432; Spawn = 4.25; elseif spawndays == 1522; Spawn = 4.5; elseif spawndays == 1675; Spawn = 5; elseif spawndays == 1797; Spawn = 5.25; elseif spawndays == 1887; Spawn = 5.5; elseif spawndays == 2040; Spawn = 6; elseif spawndays == 2162; Spawn = 6.25; elseif spawndays == 2252; Spawn = 6.5; elseif spawndays == 2405; Spawn = 7; elseif spawndays == 2527; Spawn = 7.25; elseif spawndays == 2618; Spawn = 7.5; elseif spawndays == 2771; Spawn = 8; elseif spawndays == 2893; Spawn = 8.25; elseif spawndays == 2983; Spawn = 8.5; elseif spawndays == 3136; Spawn = 9; elseif spawndays == 3258; Spawn = 9.25; elseif spawndays == 3348; Spawn = 9.5; elseif spawndays == 3501; Spawn = 10; else Spawn = 10.25; end if d >= spawndays; f = 0; Starve(d) = Starve(d-1) + 1; Sardine SO_Dev(d) = SO_DevTime2-1; if SER(d) <= 0; disp('The fish is dead') %Stop the starve and reset for new spawn d = 2; % Break out of inner loop and go to next spawn break end else f = 0.8; Starve(d) = Starve(d-1) + 0; Sardine end Sardine_Save fprintf('d = %f\nf = %f\nStarve(d) = %f\nER = %f\nLength = %f\nWet Weight = %f\n\n',d,f,Starve(d),SER(d),SL(d),SW_ww(d)); end end %hold off