function LonModes5_long(DataFolder,RadarsWorking,all,p)
location0='LonFit';
location='LonMod';
lon=[-180:1:180];
LonMod.lon=lon;
if p==2
load(fullfile(DataFolder,location0,['TimeLonVecWind_Fit5',int2str(10),'_long']))
elseif p==3
load(fullfile(DataFolder,location0,['TimeLonVecWind_Fit5',int2str(10),'_long3']))
end

LonMod.Time=TimeLonVec.Time;
LonMod.S0=ones(size(TimeLonVec.Time,2),length(lon))*NaN;
LonMod.S095=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S1=ones(size(LonMod.S0,1),361)*NaN;
LonMod.S1A=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S1A95=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S1C=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S1C95=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S2=ones(size(LonMod.S0,1),361)*NaN;
LonMod.S2A=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S2A95=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S2C=ones(size(LonMod.S0,1),1)*NaN;
LonMod.S2C95=ones(size(LonMod.S0,1),1)*NaN;
LonMod.RadarsWorking=ones(size(LonMod.S0,1),1)*NaN;
MeanTest=ones(size(LonMod.S0,1),361)*NaN;

for s=10:-1:RadarsWorking

if p==2
load(fullfile(DataFolder,location0,['TimeLonVecWind_Fit5',int2str(s),'_long']))
Name=strcat('LonFit5_longAtLeast',int2str(RadarsWorking),'RadarsWorking');
elseif p==3
load(fullfile(DataFolder,location0,['TimeLonVecWind_Fit5',int2str(s),'_long3']))
Name=strcat('LonFit5_longAtLeast',int2str(RadarsWorking),'RadarsWorking_3');
end

if s==6 && all==1 && p==2
   load(fullfile(DataFolder,location0,['TimeLonVecWind_Fit5',int2str(s),'_long_all']))
   Name=[Name,'_all'];
end


ind=find(isnan(TimeLonVec.Fit.S0)==0);
for i=1:size(ind,2)
MeanTest(ind(i),:)=TimeLonVec.Fit.MeanTest(ind(i))*ones(size(lon));
LonMod.S0(ind(i),:)=TimeLonVec.Fit.S0(ind(i))*ones(size(lon));
LonMod.S095(ind(i))=TimeLonVec.Fit.S0_95(ind(i));
LonMod.S1(ind(i),:)=TimeLonVec.Fit.S1(ind(i))*sin(2*pi*(lon/360+TimeLonVec.Fit.c1(ind(i))));
LonMod.S1A(ind(i))=TimeLonVec.Fit.S1(ind(i));
LonMod.S1A95(ind(i))=TimeLonVec.Fit.S1_95(ind(i));
LonMod.S1C(ind(i))=TimeLonVec.Fit.c1(ind(i));
LonMod.S1C95(ind(i))=TimeLonVec.Fit.c1_95(ind(i));
LonMod.S2(ind(i),:)=TimeLonVec.Fit.S2(ind(i))*sin(2*pi*(lon/180+TimeLonVec.Fit.c2(ind(i))));
LonMod.S2A(ind(i))=TimeLonVec.Fit.S2(ind(i));
LonMod.S2A95(ind(i))=TimeLonVec.Fit.S2_95(ind(i));
LonMod.S2C(ind(i))=TimeLonVec.Fit.c2(ind(i));
LonMod.S2C95(ind(i))=TimeLonVec.Fit.c2_95(ind(i));
LonMod.RadarsWorking(ind(i))=s;
end


end


 
 if s==6 && all==1 && p==2
 save(fullfile(DataFolder,location,['LonModWeightedAtLeast',int2str(RadarsWorking),'RadarsWorking_long_all_Fit5'],'LonMod'))
 elseif p==2
  save(fullfile(DataFolder,location,['LonModWeightedAtLeast',int2str(RadarsWorking),'RadarsWorking_long_Fit5'],'LonMod'))
 elseif p==3
     save(fullfile(DataFolder,location,['LonModWeightedAtLeast',int2str(RadarsWorking),'RadarsWorking_long_Fit5_3']),'LonMod')
 end
