function RecDur = hydrarecdur(RecSize, SampRate, RecAmp, RecComp, RecSen, RecStat, RecExt, RecParos, RecCTD, DutyCycle);
% function RecDur = hydrarecdur(RecSize, SampRate, RecAmp, RecComp, RecSen,
%                      RecStat, RecExt, RecParos, RecCTD, DutyCycle);
%
% calculates recorder duration (in days) for a Hydra deployment 
% based on the recorder capacity and deployment parameters 
% 
% Parameters:  
% RecSize, - recorder size (Mb)
% SampRate - sampling rate (Hz)
%
% the following switches take 1 as 'YES' and 0 as 'NO'
%
% RecAmp   - Record Amplitude 
% RecComp  - Record Compass
% RecSen,  - Record Sensors
% RecStat  - Record Statistics
% RecExt   - Record External Sensors
% RecParos - Record Paros
% RecCTD   - Record CTD 
% DutyCycle- Duty cycle (%)

HdrSize	 = 441;
BhdrSize = 62;
VelSize	 = 8;
AmpSize  = 4;
CompSize = 6;
SenSize	 = 4;
ExtSize  = 4;
StatSize = 36;
ParosSize= 3;
CTDSize  = 4;

DutyCycle   = DutyCycle/100;

Nbursts 	= length(RecAmp);
SampSize	= VelSize + AmpSize*RecAmp + CompSize*RecComp + SenSize*RecSen + ...
              ExtSize*RecExt + ParosSize*RecParos + CTDSize*RecCTD;

Nsamp   	= DutyCycle*3600.*SampRate;
BurstSize   = BhdrSize + StatSize*RecStat + Nsamp.*SampSize;

Nbursts     = (RecSize*1e6 - HdrSize)./BurstSize;
Blength 	= Nbursts.*DutyCycle*3600;

% total number of days 

RecDur=  Nbursts.*(SampRate>0)./24;  % convert to days