Error: "The volume can not be fuond"

 

Function SV()
	Variable/G experimenttime
	Variable/G recordingtime
	Variable/G duration
	Variable/G rate
	Variable/G savedrate
	Variable/G plotresolution
	// variable comment
	Variable/G com
	string/G wct
end

 //Principal
Function master(experimenttime, recordingtime, duration, rate, savedrate, plotresolution)
variable experimenttime, recordingtime, duration, rate, savedrate, plotresolution
variable psdflag=1 // first time, make the periodogram
NewPath/O Path1, "C:Users:Carlo:SD" //defines where to save the files

// rate= sampling rate for data acquiring, like 40000 Hz
// savedrate= rate for high time-resolution data recording, like 200 Hz
// plotresolution= interval between each point on strip graph
// recordingtime= interval for saving on disk, like 10s or 30s
// experimenttime= global experiment time, after which acquisition stops, like 3600 s (break stops before and saves)
// savepath = forlder where save the files, like "C:Users:Carlo:SD"


string timenow=date()+" "+time() // date and time now
string filename
variable t0 = stopMSTimer(-2)/10^6 // consults the microseconds timer, from the moment of PC startup
variable t1=0, dt1=0, dt2=0 // time init
if (plotresolutionduration
make /O/N=0 variancestrip //where to store the variance every plotresolution>duration
setscale/P x, t1, (plotresolution), "s", timestrip, deflstrip, variancestrip // sets the spacing of the waves, just to have "s" as unit
setscale/P y, 0,1,"V" deflstrip, variancestrip 
setscale/P y,0,1,"s",timestrip
//display
display/W=(220,1,930,150) deflstrip vs timestrip
AppendToGraph/L=L2 variancestrip vs timestrip
ModifyGraph axisEnab(left)={0,0.45}
ModifyGraph axisEnab(L2)={0.55,1}
ModifyGraph freePos(L2)=0
ModifyGraph lblPosMode(L2)=2
ModifyGraph rgb(deflstrip)=(0,0,0)


make /O/N=0 storeaverage, storevariance, storetime, waveaverage, wavevariance, wavetime // to store average and variance for the recordingtime
display/W=(220,220, 550,450) storeaverage vs storetime
variable acq_cycle=0 // counter for acquisition cycle

do //big loop for acquisition and plot
	t1=stopMSTimer(-2)/10^6-t0 // evaluate the experimental time
	//print "t1=",t1
	waveacquire(t1, duration, rate, savedrate, psdflag) // DATA ACQUISITION
	dt1+=stopMSTimer(-2)/10^6-t0-t1 // increment the delta time for plotting
	dt2+=stopMSTimer(-2)/10^6-t0-t1 // increment the delta time for saving
	
	insertpoints numpnts(storeaverage), numpnts(waveaverage), storeaverage // extends storeaverage
	storeaverage[numpnts(storeaverage)-numpnts(waveaverage), numpnts(storeaverage)-1]=waveaverage[p-numpnts(storeaverage)+numpnts(waveaverage)-1] // stores the averaged deflection points
	insertpoints numpnts(storevariance), numpnts(wavevariance), storevariance // extends storevariance
	storevariance[numpnts(storevariance)-numpnts(wavevariance), numpnts(storevariance)-1]=wavevariance[p-numpnts(storeaverage)+numpnts(waveaverage)-1] // stores the variance time points too
	insertpoints numpnts(storetime), numpnts(wavetime), storetime // extends storetime
	storetime[numpnts(storetime)-numpnts(wavetime), numpnts(storetime)-1]=wavetime[p-numpnts(storeaverage)+numpnts(waveaverage)-1] // and the times
	
	acq_cycle=+1 //increment the counter for acquisition cycle
	if (psdflag==1)
		psdflag=0
	endif // switches off psd
	//print dt1," ",dt2
	if (dt1>=plotresolution) // after exiting this loops, it adds a plot point ... 
		// insertpoints numpnts(deflstrip), 1, deflstrip, variancestrip // makes room
		timestrip[numpnts(timestrip)]={stopMSTimer(-2)/10^6-t0} // to store the time too
		deflstrip[numpnts(deflstrip)]={mean(storeaverage, pnt2x(storeaverage, numpnts(storeaverage)-acq_cycle*rate/savedrate), pnt2x(storeaverage, numpnts(storeaverage)))}
		variancestrip[numpnts(variancestrip)]={mean(storevariance, pnt2x(storevariance, numpnts(storevariance)-acq_cycle*rate/savedrate), pnt2x(storevariance, numpnts(storevariance)))/sqrt(acq_cycle*rate/savedrate)}
		doupdate // explicitly ask to update the plots (in hope it has not done it before at insertpoint)
		// print "plotplot"
		dt1=0 //resets delta time for plotting
		acq_cycle=0 // resets acquisition cycle, so it later knows what to use for plotting
	endif
	if  (dt2>=recordingtime) // after exiting this loop, it records on file and resets storeaverage and storevariance
		// saves storeaverage and storevariance on file with automatic names 
		// empty storeaverage and storevariance for next acquisitions
		 timenow=date()+" "+time() // updates time
		 filename="DeflectionStrip "+num2str(savedrate)+"Hz on "+timenow+"-"+num2istr(dt2)+".txt" // for saving
		save /G/P=Path1 storetime storeaverage storevariance as filename		// print "I want to save storeaverage and storevariance"
		redimension/N=0 storetime, storeaverage, storevariance //clean up
		dt2=0 //resets delta time for recording
		psdflag=1 // nextg cycle, record the PSD (if there is no next cycle, it records it at the end)
	endif

while (t18192*4)
	segmentsize=8192// segment size for Periodogram calculation
else
	segmentsize=(duration*rate)/5
endif
//print segmentsize
make/O/N=(rate*duration) wave0
make/O/N=(rate*duration/blocksize) wavevariance // wave variance decimated
make/O/N=(rate*duration/blocksize) waveaverage // wave after decimation
make/O/N=(rate*duration/blocksize) wavetime // to store the time points
make/O/N=0 media //to store the recordings
make/O/N=0 timepoints //to store the time of recordings
make/O/N=(segmentsize/2+1) periodogramdomain


setscale/P x, t1, (1/rate), "s", wave0 // sets the spacing of the wave, so the sampling rate

DAQmx_scan/DEV="dev3" WAVES="wave0, 2;" //acquires from DEV device on channel 2 and puts in wave0

if (psdflag==1)
	DSPPeriodogram/Q/DLSG/NOR=(segmentsize*segmentsize/2)/SEGN={(segmentsize),(segmentsize/2)}/Win=hanning wave0
	timenow=date()+" "+time() // updates time
	filename="PSD "+num2str(savedrate)+"Hz on "+timenow+".txt" // for saving
	periodogramdomain=pnt2x(W_periodogram, p)
	save /G/P=Path1 periodogramdomain W_periodogram as filename		// print "I want to save storeaverage and storevariance"

endif

wavevariance=variance(wave0, pnt2x(wave0, blocksize*p), pnt2x(wave0, (blocksize*(p+1))-1))/blocksize
setscale/P x, t1, blocksize*(1/rate), "s", wavevariance // sets the spacing of the wave, so the sampling rate
waveaverage=sum(wave0, pnt2x(wave0, blocksize*p), pnt2x(wave0, (blocksize*(p+1))-1))/blocksize
setscale/P x, t1, blocksize*(1/rate), "s", waveaverage // sets the spacing of the wave, so the sampling rate
wavetime=pnt2x(wavevariance, p) // timings from the scale (?)
variable stoptimer=stopMSTimer(-2) // stops the microseconds timer

end