
cruise='DYNAMO_2011';
yy=[];
y1a=load('c:\data\cwf\Dropbox\DATA\dynamo\files\DYNAMO_2011_micro_m_leg1.txt');%y1b=reshape(y1a,6,length(y1a)/6);y1a=y1b'; 
y2a=load('c:\data\cwf\Dropbox\DATA\dynamo\files\DYNAMO_2011_micro_m_leg2.txt');%y2b=reshape(y2a,6,length(y2a)/6);y2a=y2b';
y3a=load('c:\data\cwf\Dropbox\DATA\dynamo\files\DYNAMO_2011_micro_m_leg3.txt');%y3b=reshape(y3a,6,length(y3a)/6);y3a=y3b'; 
y4a=load('c:\data\cwf\Dropbox\DATA\dynamo\files\DYNAMO_2011_micro_m_leg4.txt');%read file with 10min 

yyx=[y1a' y2a' y3a' y4a']';
jdx=yyx(:,1);
ii=find(jdx<200);jdx(ii)=jdx(ii)+366;
yi=interp1(jdx,yyx,jdy);

   

   jdym=yi(:,1);%julian day at beginning of time average
   vap=yi(:,2);%column water vap (cm)
   liq=yi(:,3)*1e4;%column liq (mm), convert to g/m^2
   xnm=yi(:,4);%numbwer obs
   sigvap=yi(:,5);%
   sigliq=yi(:,6);%

ii=find(isnan(vap));vap(ii)=qa(ii)/3.3;
ii=find(sigvap>0.1 | isnan(sigvap));liq(ii)=NaN;
figure;plot(jdym,liq);


