#pragma rtGlobals=3 // Use modern global access method and strict wave access. //This script calculates average concentrations of measured parameters during NAAMES 1 Cruise. //The NAAMES1 cruise times were largely divided into 4 parts based on aerosol properties, in situ sampling, and wind direction. //1. Continental: Sampled ambient air in North Atlantic was influenced by continental pollution. //2. Marine: No influence of anthropogenic or continental pollution. //3. Winds Abaft: Winds were from abaft the sampling inlet and possibly contaminated the sampled ambient air with ship plume. //4. Sea sweep: Sea sweep was deployed at stations to generate and sample insitu sea spray aerosols //This script computes the average and standard deviation for these times and creates waves in igor. //Naming of waves //1) Computed Averages: NAAMES1_(Type)_Avg(parametername) // (e.g., NAAMES1_SS_AvgOrg: Average of organic mass concentrations during Sea Sweep (SS)) //2) Computed Standard Deviation: NAAMES1_(Type)_StdDev(parametername) // (e.g., NAAMES1_SS_AvgOrg: Standard Deviation of organic mass concentrations during Sea Sweep (SS)) //3) Start Time : NAAMES1_(Type)_StartTime (e.g., NAAMES1_SS_StartTime) //4) Stop Time: NAAMES1_(Type)_StopTime (e.g., NAAMES1_SS_StopTime) //Type: a) Continental- Con // b) Marine - Mar // c)Winds Abaft - WinAb // d) Sea Sweep - SS //How it works? // 1)Load the procedure in an Igor Experiment // 2) Once you compile, a menu called NAAMES is created // 3) Click on "NAAMES" -> "NAAMES 1" // 4) A dialog box will pop up // 5) Fill the required parameters and hit "Continue" Menu "NAAMES" "NAAMES1", NAAMES1() End Function time2secs (hh, mm, ss) variable hh, mm, ss variable totalsecs totalsecs = hh*3600+mm*60+ss return totalsecs End Function NAAMES1() String StudyProject String ChemCpd String Timewave prompt StudyProject, "Category:", popup "Continental;Marine;Winds Abaft;Sea Sweep" prompt ChemCpd, "Name of Chemical Compound Wave: e.g., (Organics; Sulfate; Nitrate)" prompt Timewave, "Enter the TimeWave" doprompt "NAAMES1:", StudyProject, ChemCpd, Timewave Variable j String Cpd print StudyProject print ChemCpd print timewave for (j = 0; j < ItemsinList(ChemCpd); j +=1) Cpd = StringFromlist(j, ChemCpd) Wave AMSDistr = $Cpd Wave AMSTime= $Timewave Variable SS1start, SS1End, SS2start, SS2End, SS3start, SS3End, SS4start, SS4End, SS5start, SS5End, SS6start, SS6End, SS7start, SS7End Variable SS8start, SS8End, SS9start, SS9End, SS10start, SS10End, SS11start, SS11End, SS12start, SS12End, SS13start, SS13End, SS14start, SS14End,SS15start, SS15End, SS16start, SS16End,SS17start, SS17End, SS18start, SS18End, SS19start, SS19End String AvSS, SDSS, StrtTimeW, StpTimeW String temp Make /O/N =1 tempSS1, tempSS2, tempSS3, tempSS4, tempSS5, tempSS6, tempSS7 Make /O/N =1 tempSS8, tempSS9, tempSS10, tempSS11, tempSS12, tempSS13, tempSS14, tempSS15, tempSS16, tempSS17, tempSS18, tempSS19 Variable i ,ss1=0, ss2=0, ss3=0, ss4=0, ss5=0, ss6=0, ss7=0, ss8 =0, ss9=0, ss10=0, ss11=0, ss12=0, ss13=0, ss14=0, ss15=0, ss16=0, ss17=0, ss18=0, ss19=0 strswitch(StudyProject) Case "Sea Sweep": AvSS= "NAAMES1_SS_Avg"+Cpd SDSS= "NAAMES1_SS_StdDev"+Cpd StrtTimeW="NAAMES1_SS_StartTime" StpTimeW="NAAMES1_SS_StopTime" Make /O/N=7 $AvSS, $SDSS Make/O/D/N=7 $StrtTimeW, $StpTimeW Wave AvgSS = $AvSS Wave StdDevSS = $SDSS Wave StartTime = $StrtTimeW Wave StopTime = $StpTimeW SS1start = date2secs(2015,11,13)+ time2secs(16, 10, 00) SS1end = date2secs(2015,11,13) + time2secs(20, 45, 00) SS2start = date2secs(2015,11,14) + time2secs(10, 22, 00) SS2end = date2secs(2015,11,14) + time2secs(15, 30, 00) SS3start = date2secs(2015,11,15) + time2secs(09, 07, 00) SS3end = date2secs(2015,11,15) + time2secs(10, 12, 00) SS4start = date2secs(2015,11,18) + time2secs(09, 02, 00) SS4end = date2secs(2015,11,18) + time2secs(13, 36, 00) SS5start = date2secs(2015,11,22)+ time2secs(12, 00, 00) SS5end = date2secs(2015,11,22) + time2secs(15, 18, 00) SS6start = date2secs(2015,11,24) + time2secs(09, 04, 00) SS6end = date2secs(2015,11,24) + time2secs(11, 10, 00) SS7start = date2secs(2015,11,25) + time2secs(10, 02, 00) SS7end = date2secs(2015,11,25) + time2secs(12, 53, 00) StartTime[0] = SS1start StartTime[1] =SS2start StartTime[2] =SS3start StartTime[3] =SS4start StartTime[4] =SS5start StartTime[5] =SS6start StartTime[6] =SS7start StopTime[0] = SS1End StopTime[1] =SS2End StopTime[2] =SS3End StopTime[3] =SS4End StopTime[4] =SS5End StopTime[5] =SS6End StopTime[6] =SS7End SetScale d 0,0,"dat", StartTime SetScale d 0,0,"dat", StopTime for (i = 0; i < numpnts(AMSTime); i +=1) if (AMSTime[i] >=SS1start && AMSTime[i] < SS1end&&numtype(AMSDistr[i])== 0) tempSS1[ss1] = AMSDistr[i] ss1 +=1 insertpoints ss1, 1, tempSS1 elseif(AMSTime[i] >=SS2start && AMSTime[i] < SS2end&&numtype(AMSDistr[i])== 0) tempSS2[ss2] = AMSDistr[i] ss2 +=1 insertpoints ss2, 1, tempSS2 elseif(AMSTime[i] >=SS3start && AMSTime[i] < SS3end&&numtype(AMSDistr[i])== 0) tempSS3[ss3] = AMSDistr[i] ss3 +=1 insertpoints ss3, 1, tempSS3 elseif(AMSTime[i] >=SS4start && AMSTime[i] < SS4end&&numtype(AMSDistr[i])== 0) tempSS4[ss4] = AMSDistr[i] ss4 +=1 insertpoints ss4, 1, tempSS4 elseif(AMSTime[i] >=SS5start && AMSTime[i] < SS5end&&numtype(AMSDistr[i])== 0) tempSS5[ss5] = AMSDistr[i] ss5 +=1 insertpoints ss5, 1, tempSS5 elseif(AMSTime[i] >=SS6start && AMSTime[i] < SS6end&&numtype(AMSDistr[i])== 0) tempSS6[ss6] = AMSDistr[i] ss6 +=1 insertpoints ss6, 1, tempSS6 elseif(AMSTime[i] >=SS7start && AMSTime[i] < SS7end&&numtype(AMSDistr[i])== 0) tempSS7[ss7] = AMSDistr[i] ss7 +=1 insertpoints ss7, 1, tempSS7 endif endfor DeletePoints ss1,1, tempSS1 DeletePoints ss2,1, tempSS2 DeletePoints ss3,1, tempSS3 DeletePoints ss4,1, tempSS4 DeletePoints ss5,1, tempSS5 DeletePoints ss6,1, tempSS6 DeletePoints ss7,1, tempSS7 for (i = 1; i <8; i +=1) temp = "tempSS"+num2str(i) Wavestats /Q $temp AvgSS[i-1] = V_avg StdDevSS[i-1] = V_sdev endfor break Case "Continental": AvSS= "NAAMES1_Con_Avg"+Cpd SDSS= "NAAMES1_Con_StdDev"+Cpd StrtTimeW="NAAMES1_Con_StartTime" StpTimeW="NAAMES1_Con_StopTime" Make /O/N=5 $AvSS, $SDSS Make/O/D/N=5 $StrtTimeW, $StpTimeW Wave AvgSS = $AvSS Wave StdDevSS = $SDSS Wave StartTime = $StrtTimeW Wave StopTime = $StpTimeW SS1start = date2secs(2015,11,06)+ time2secs(00, 00, 00) SS1end = date2secs(2015,11,07) + time2secs(01, 06, 00) SS2start = date2secs(2015,11,07) + time2secs(20, 31, 00) SS2end = date2secs(2015,11,10) + time2secs(11, 57, 00) SS3start = date2secs(2015,11,10) + time2secs(23, 43, 00) SS3end = date2secs(2015,11,11) + time2secs(09, 00, 00) SS4start = date2secs(2015,11,21) + time2secs(00, 00, 00) SS4end = date2secs(2015,11,22) + time2secs(06, 00, 00) SS5start = date2secs(2015,11,29)+ time2secs(01, 16, 00) SS5end = date2secs(2015,12,01) + time2secs(00, 00, 00) StartTime[0] = SS1start StartTime[1] =SS2start StartTime[2] =SS3start StartTime[3] =SS4start StartTime[4] =SS5start StopTime[0] = SS1End StopTime[1] =SS2End StopTime[2] =SS3End StopTime[3] =SS4End StopTime[4] =SS5End SetScale d 0,0,"dat", StartTime SetScale d 0,0,"dat", StopTime for (i = 0; i < numpnts(AMSTime); i +=1) if (AMSTime[i] >=SS1start && AMSTime[i] < SS1end&&numtype(AMSDistr[i])== 0) tempSS1[ss1] = AMSDistr[i] ss1 +=1 insertpoints ss1, 1, tempSS1 elseif(AMSTime[i] >=SS2start && AMSTime[i] < SS2end&&numtype(AMSDistr[i])== 0) tempSS2[ss2] = AMSDistr[i] ss2 +=1 insertpoints ss2, 1, tempSS2 elseif(AMSTime[i] >=SS3start && AMSTime[i] < SS3end&&numtype(AMSDistr[i])== 0) tempSS3[ss3] = AMSDistr[i] ss3 +=1 insertpoints ss3, 1, tempSS3 elseif(AMSTime[i] >=SS4start && AMSTime[i] < SS4end&&numtype(AMSDistr[i])== 0) tempSS4[ss4] = AMSDistr[i] ss4 +=1 insertpoints ss4, 1, tempSS4 elseif(AMSTime[i] >=SS5start && AMSTime[i] < SS5end&&numtype(AMSDistr[i])== 0) tempSS5[ss5] = AMSDistr[i] ss5 +=1 insertpoints ss5, 1, tempSS5 endif endfor DeletePoints ss1,1, tempSS1 DeletePoints ss2,1, tempSS2 DeletePoints ss3,1, tempSS3 DeletePoints ss4,1, tempSS4 DeletePoints ss5,1, tempSS5 for (i = 1; i <6; i +=1) temp = "tempSS"+num2str(i) Wavestats /Q $temp AvgSS[i-1] = V_avg StdDevSS[i-1] = V_sdev endfor break Case "Marine": AvSS= "NAAMES1_Mar_Avg"+Cpd SDSS= "NAAMES1_Mar_StdDev"+Cpd StrtTimeW="NAAMES1_Mar_StartTime" StpTimeW="NAAMES1_Mar_StopTime" Make /O/N=14 $AvSS, $SDSS Make/O/D/N=14 $StrtTimeW, $StpTimeW Wave AvgSS = $AvSS Wave StdDevSS = $SDSS Wave StartTime = $StrtTimeW Wave StopTime = $StpTimeW SS1start = date2secs(2015,11,11)+ time2secs(09, 00, 00) SS1end = date2secs(2015,11,13) + time2secs(15, 50, 00) SS2start = date2secs(2015,11,13) + time2secs(20, 50, 00) SS2end = date2secs(2015,11,14) + time2secs(08, 40, 00) SS3start = date2secs(2015,11,14) + time2secs(15, 35, 00) SS3end = date2secs(2015,11,15) + time2secs(05, 55, 00) SS4start = date2secs(2015,11,15) + time2secs(12, 00, 00) SS4end = date2secs(2015,11,17) + time2secs(22, 35, 00) SS5start =date2secs(2015,11,18) + time2secs(14, 10, 00) SS5end = date2secs(2015,11,18) + time2secs(22, 05, 00) SS6start = date2secs(2015,11,18)+ time2secs(23, 30, 00) SS6end = date2secs(2015,11,20) + time2secs(00, 45, 00) SS7start =date2secs(2015,11,20) + time2secs(03, 00, 00) SS7end = date2secs(2015,11,20) + time2secs(15, 15, 00) SS8start = date2secs(2015,11,22) + time2secs(15, 35, 00) SS8end = date2secs(2015,11,23) + time2secs(05, 10, 00) SS9start = date2secs(2015,11,23) + time2secs(06, 45, 00) SS9end = date2secs(2015,11,24) + time2secs(07, 00, 00) SS10start = date2secs(2015,11,24) + time2secs(12, 00, 00) SS10end =date2secs(2015,11,24) + time2secs(17, 50, 00) SS11start = date2secs(2015,11,24) + time2secs(21, 45, 00) SS11end = date2secs(2015,11,25) + time2secs(08, 35, 00) SS12start = date2secs(2015,11,25) + time2secs(15, 10, 00) SS12end = date2secs(2015,11,26) + time2secs(15, 55, 00) SS13start =date2secs(2015,11,27) + time2secs(02, 25, 00) SS13end =date2secs(2015,11,27) + time2secs(09, 10, 00) SS14start =date2secs(2015,11,28) + time2secs(01, 00, 00) SS14end = date2secs(2015,11,29) + time2secs(01, 15, 00) StartTime[0] = SS1start StartTime[1] =SS2start StartTime[2] =SS3start StartTime[3] =SS4start StartTime[4] =SS5start StartTime[5] =SS6start StartTime[6] =SS7start StartTime[7] = SS8start StartTime[8] =SS9start StartTime[9] =SS10start StartTime[10] =SS11start StartTime[11] =SS12start StartTime[12] =SS13start StartTime[13] =SS14start StopTime[0] = SS1End StopTime[1] =SS2End StopTime[2] =SS3End StopTime[3] =SS4End StopTime[4] =SS5End StopTime[5] =SS6End StopTime[6] =SS7End StopTime[7] = SS8End StopTime[8] =SS9End StopTime[9] =SS10End StopTime[10] =SS11End StopTime[11] =SS12End StopTime[12] =SS13End StopTime[13] =SS14End SetScale d 0,0,"dat", StartTime SetScale d 0,0,"dat", StopTime for (i = 0; i < numpnts(AMSTime); i +=1) if (AMSTime[i] >=SS1start && AMSTime[i] < SS1end&&numtype(AMSDistr[i])== 0) tempSS1[ss1] = AMSDistr[i] ss1 +=1 insertpoints ss1, 1, tempSS1 elseif(AMSTime[i] >=SS2start && AMSTime[i] < SS2end&&numtype(AMSDistr[i])== 0) tempSS2[ss2] = AMSDistr[i] ss2 +=1 insertpoints ss2, 1, tempSS2 elseif(AMSTime[i] >=SS3start && AMSTime[i] < SS3end&&numtype(AMSDistr[i])== 0) tempSS3[ss3] = AMSDistr[i] ss3 +=1 insertpoints ss3, 1, tempSS3 elseif(AMSTime[i] >=SS4start && AMSTime[i] < SS4end&&numtype(AMSDistr[i])== 0) tempSS4[ss4] = AMSDistr[i] ss4 +=1 insertpoints ss4, 1, tempSS4 elseif(AMSTime[i] >=SS5start && AMSTime[i] < SS5end&&numtype(AMSDistr[i])== 0) tempSS5[ss5] = AMSDistr[i] ss5 +=1 insertpoints ss5, 1, tempSS5 elseif(AMSTime[i] >=SS6start && AMSTime[i] < SS6end&&numtype(AMSDistr[i])== 0) tempSS6[ss6] = AMSDistr[i] ss6 +=1 insertpoints ss6, 1, tempSS6 elseif(AMSTime[i] >=SS7start && AMSTime[i] < SS7end&&numtype(AMSDistr[i])== 0) tempSS7[ss7] = AMSDistr[i] ss7 +=1 insertpoints ss7, 1, tempSS7 elseif(AMSTime[i] >=SS8start && AMSTime[i] < SS8end&&numtype(AMSDistr[i])== 0) tempSS8[ss8] = AMSDistr[i] ss8 +=1 insertpoints ss8, 1, tempSS8 elseif(AMSTime[i] >=SS9start && AMSTime[i] < SS9end&&numtype(AMSDistr[i])== 0) tempSS9[ss9] = AMSDistr[i] ss9 +=1 insertpoints ss9, 1, tempSS9 elseif(AMSTime[i] >=SS10start && AMSTime[i] < SS10end&&numtype(AMSDistr[i])== 0) tempSS10[ss10] = AMSDistr[i] ss10 +=1 insertpoints ss10, 1, tempSS10 elseif(AMSTime[i] >=SS11start && AMSTime[i] < SS11end&&numtype(AMSDistr[i])== 0) tempSS11[ss11] = AMSDistr[i] ss11 +=1 insertpoints ss11, 1, tempSS11 elseif(AMSTime[i] >=SS12start && AMSTime[i] < SS12end&&numtype(AMSDistr[i])== 0) tempSS12[ss12] = AMSDistr[i] ss12 +=1 insertpoints ss12, 1, tempSS12 elseif(AMSTime[i] >=SS13start && AMSTime[i] < SS13end&&numtype(AMSDistr[i])== 0) tempSS13[ss13] = AMSDistr[i] ss13 +=1 insertpoints ss13, 1, tempSS13 elseif(AMSTime[i] >=SS14start && AMSTime[i] < SS14end&&numtype(AMSDistr[i])== 0) tempSS14[ss14] = AMSDistr[i] ss14 +=1 insertpoints ss14, 1, tempSS14 endif endfor DeletePoints ss1,1, tempSS1 DeletePoints ss2,1, tempSS2 DeletePoints ss3,1, tempSS3 DeletePoints ss4,1, tempSS4 DeletePoints ss5,1, tempSS5 DeletePoints ss6,1, tempSS6 DeletePoints ss7,1, tempSS7 DeletePoints ss8,1, tempSS8 DeletePoints ss9,1, tempSS9 DeletePoints ss10,1, tempSS10 DeletePoints ss11,1, tempSS11 DeletePoints ss12,1, tempSS12 DeletePoints ss13,1, tempSS13 DeletePoints ss14,1, tempSS14 for (i = 1; i <15; i +=1) temp = "tempSS"+num2str(i) Wavestats /Q $temp AvgSS[i-1] = V_avg StdDevSS[i-1] = V_sdev endfor break Case "Winds Abaft": AvSS= "NAAMES1_WinAb_Avg"+Cpd SDSS= "NAAMES1_WinAb_StdDev"+Cpd StrtTimeW="NAAMES1_WinAb_StartTime" StpTimeW="NAAMES1_WinAb_StopTime" Make /O/N=13 $AvSS, $SDSS Make/O/D/N=13 $StrtTimeW, $StpTimeW Wave AvgSS = $AvSS Wave StdDevSS = $SDSS Wave StartTime = $StrtTimeW Wave StopTime = $StpTimeW SS1start = date2secs(2015,11,07)+ time2secs(01, 07, 00) SS1end = date2secs(2015,11,07) + time2secs(20, 30, 00) SS2start = date2secs(2015,11,10) + time2secs(11, 58, 00) SS2end = date2secs(2015,11,10) + time2secs(23, 42, 00) SS3start = date2secs(2015,11,15) + time2secs(05, 58, 00) SS3end = date2secs(2015,11,15) + time2secs(07, 18, 00) SS4start = date2secs(2015,11,18) + time2secs(04, 27, 00) SS4end = date2secs(2015,11,18) + time2secs(07, 27, 00) SS5start = date2secs(2015,11,18)+ time2secs(22, 06, 00) SS5end = date2secs(2015,11,18) + time2secs(23, 27, 00) SS6start = date2secs(2015,11,20) + time2secs(00, 48, 00) SS6end = date2secs(2015,11,20) + time2secs(02, 59, 00) SS7start = date2secs(2015,11,20) + time2secs(15, 24, 00) SS7end = date2secs(2015,11,20) + time2secs(19, 07, 00) SS8start = date2secs(2015,11,20) + time2secs(20, 33, 00) SS8end = date2secs(2015,11,20) + time2secs(23, 46, 00) SS9start = date2secs(2015,11,23) + time2secs(05, 11, 00) SS9end = date2secs(2015,11,23) + time2secs(06, 42, 00) SS10start = date2secs(2015,11,24) + time2secs(17, 51, 00) SS10end = date2secs(2015,11,24) + time2secs(21, 42, 00) SS11start = date2secs(2015,11,25) + time2secs(14, 25, 00) SS11end = date2secs(2015,11,25) + time2secs(15, 10, 00) SS12start = date2secs(2015,11,26) + time2secs(15, 55, 00) SS12end = date2secs(2015,11,27) + time2secs(02, 25, 00) SS13start = date2secs(2015,11,27) + time2secs(09, 13, 00) SS13end = date2secs(2015,11,27) + time2secs(14, 50, 00) StartTime[0] = SS1start StartTime[1] =SS2start StartTime[2] =SS3start StartTime[3] =SS4start StartTime[4] =SS5start StartTime[5] =SS6start StartTime[6] =SS7start StartTime[7] = SS8start StartTime[8] =SS9start StartTime[9] =SS10start StartTime[10] =SS11start StartTime[11] =SS12start StartTime[12] =SS13start StopTime[0] = SS1End StopTime[1] =SS2End StopTime[2] =SS3End StopTime[3] =SS4End StopTime[4] =SS5End StopTime[5] =SS6End StopTime[6] =SS7End StopTime[7] = SS8End StopTime[8] =SS9End StopTime[9] =SS10End StopTime[10] =SS11End StopTime[11] =SS12End StopTime[12] =SS13End SetScale d 0,0,"dat", StartTime SetScale d 0,0,"dat", StopTime for (i = 0; i < numpnts(AMSTime); i +=1) if (AMSTime[i] >=SS1start && AMSTime[i] < SS1end&&numtype(AMSDistr[i])== 0) tempSS1[ss1] = AMSDistr[i] ss1 +=1 insertpoints ss1, 1, tempSS1 elseif(AMSTime[i] >=SS2start && AMSTime[i] < SS2end&&numtype(AMSDistr[i])== 0) tempSS2[ss2] = AMSDistr[i] ss2 +=1 insertpoints ss2, 1, tempSS2 elseif(AMSTime[i] >=SS3start && AMSTime[i] < SS3end&&numtype(AMSDistr[i])== 0) tempSS3[ss3] = AMSDistr[i] ss3 +=1 insertpoints ss3, 1, tempSS3 elseif(AMSTime[i] >=SS4start && AMSTime[i] < SS4end&&numtype(AMSDistr[i])== 0) tempSS4[ss4] = AMSDistr[i] ss4 +=1 insertpoints ss4, 1, tempSS4 elseif(AMSTime[i] >=SS5start && AMSTime[i] < SS5end&&numtype(AMSDistr[i])== 0) tempSS5[ss5] = AMSDistr[i] ss5 +=1 insertpoints ss5, 1, tempSS5 elseif(AMSTime[i] >=SS6start && AMSTime[i] < SS6end&&numtype(AMSDistr[i])== 0) tempSS6[ss6] = AMSDistr[i] ss6 +=1 insertpoints ss6, 1, tempSS6 elseif(AMSTime[i] >=SS7start && AMSTime[i] < SS7end&&numtype(AMSDistr[i])== 0) tempSS7[ss7] = AMSDistr[i] ss7 +=1 insertpoints ss7, 1, tempSS7 elseif(AMSTime[i] >=SS8start && AMSTime[i] < SS8end&&numtype(AMSDistr[i])== 0) tempSS8[ss8] = AMSDistr[i] ss8 +=1 insertpoints ss8, 1, tempSS8 elseif(AMSTime[i] >=SS9start && AMSTime[i] < SS9end&&numtype(AMSDistr[i])== 0) tempSS9[ss9] = AMSDistr[i] ss9 +=1 insertpoints ss9, 1, tempSS9 elseif(AMSTime[i] >=SS10start && AMSTime[i] < SS10end&&numtype(AMSDistr[i])== 0) tempSS10[ss10] = AMSDistr[i] ss10 +=1 insertpoints ss10, 1, tempSS10 elseif(AMSTime[i] >=SS11start && AMSTime[i] < SS11end&&numtype(AMSDistr[i])== 0) tempSS11[ss11] = AMSDistr[i] ss11 +=1 insertpoints ss11, 1, tempSS11 elseif(AMSTime[i] >=SS12start && AMSTime[i] < SS12end&&numtype(AMSDistr[i])== 0) tempSS12[ss12] = AMSDistr[i] ss12 +=1 insertpoints ss12, 1, tempSS12 elseif(AMSTime[i] >=SS13start && AMSTime[i] < SS13end&&numtype(AMSDistr[i])== 0) tempSS13[ss13] = AMSDistr[i] ss13 +=1 insertpoints ss13, 1, tempSS13 endif endfor DeletePoints ss1,1, tempSS1 DeletePoints ss2,1, tempSS2 DeletePoints ss3,1, tempSS3 DeletePoints ss4,1, tempSS4 DeletePoints ss5,1, tempSS5 DeletePoints ss6,1, tempSS6 DeletePoints ss7,1, tempSS7 DeletePoints ss8,1, tempSS8 DeletePoints ss9,1, tempSS9 DeletePoints ss10,1, tempSS10 DeletePoints ss11,1, tempSS11 DeletePoints ss12,1, tempSS12 DeletePoints ss13,1, tempSS13 for (i = 1; i <14; i +=1) temp = "tempSS"+num2str(i) Wavestats /Q $temp AvgSS[i-1] = V_avg StdDevSS[i-1] = V_sdev endfor break endswitch Killwaves tempSS1, tempSS2, tempSS3, tempSS4, tempSS5, tempSS6, tempSS7 Killwaves tempSS8, tempSS9, tempSS10, tempSS11, tempSS12, tempSS13, tempSS14, tempSS15, tempSS16, tempSS17, tempSS18, tempSS19 endfor End