[bas1, gmeth, fs1, clat, clong, city_radius, nn, ntotal,freq]=textread('stats.out','%s %s %s %f %f %f %f %f %f');
fs=char(fs1);
bas=char(bas1);
nt=1;
if clong < 0
   clong=360+clong;
end
%
vc=zeros(1,15);
nvd=zeros(1,15);
nvdm=zeros(1,15);
for i=1:15,
   vc(i)=30+10*i;
   nvd(i)=0;
   nvdm(i)=0;
end   
pifac=acos(-1)/180;
dfac=60.*1.852;
fs2=strcat('ascii/',fs,'/');
lg=strcat(fs2,'hurr1');
lgn=strcat(lg,'.out');
fl=load(lgn);
month=fl(1,1);
nc=0;   
stormnum=zeros(1,nn);
vstore=zeros(nn,400);
vpstore=zeros(nn,400);
pstore=zeros(nn,400);
latstore=zeros(nn,400);
longstore=zeros(nn,400);
rmstore=zeros(nn,400);
monthstore=zeros(nn,400);
daystore=zeros(nn,400);
hourstore=zeros(nn,400);
shearstore=zeros(nn,400);
for n=1:nn,
   ns=num2str(n);
   lg=strcat(fs2,'/hurr',ns);
   lgn=strcat(lg,'.out');
   fl=load(lgn);
   [i,j]=size(fl);
   m=i;
   m=min(m,399);
   mp1=m+1;
   m2=2.*m-1;
   time=[0:2:m2];
   vm=fl(:,6);
   vp=fl(:,10);
   sh=fl(:,9);
   pres=fl(:,7);
   lat=fl(:,4);
   long=fl(:,5);
   monthm=fl(:,1);
   day=fl(:,2);
   hour=fl(:,3);
   rm=fl(:,8);
         stormnum(n)=n;
         vstore(n,1:m)=vm.';
         vpstore(n,1:m)=vp.';     
         pstore(n,1:m)=pres.';     
         latstore(n,1:m)=lat.';
         longstore(n,1:m)=long.';
         rmstore(n,1:m)=rm.';
	 monthstore(n,1:m)=monthm.';
         daystore(n,1:m)=day.';
	 hourstore(n,1:m)=hour.';
         shearstore(n,1:m)=sh.';
         %
         vstore(n,mp1:400)=0;
         vpstore(n,mp1:400)=0;     
         pstore(n,mp1:400)=1000;     
         latstore(n,mp1:400)=0;
         longstore(n,mp1:400)=0;
         rmstore(n,mp1:400)=0;
	 monthstore(n,mp1:400)=0;
	 daystore(n,mp1:400)=0;
	 hourstore(n,mp1:400)=0;
	 shearstore(n,mp1:400)=0;
      end   
   %
nns=num2str(nn);
factor=1000.*freq./nn;
nvdm=nvdm.*factor;
   storestring=[fs '/' fs];
   save temp.mat freq factor stormnum clat clong city_radius vstore vpstore pstore latstore longstore rmstore monthstore daystore hourstore shearstore bas gmeth 
