tshf

 TSELEV reads elevation time series from TSEPIC.CDF style files

 USAGE:  [e,jd]=tselev(cdf);

0001 function [e,jd]=tselev(cdf);
0002 % TSELEV reads elevation time series from TSEPIC.CDF style files
0003 %
0004 % USAGE:  [e,jd]=tselev(cdf);
0005 
0006 [t]=mcvgt(cdf,'time');
0007 nt=length(t);
0008 %
0009 [stations]=mcvgt(cdf,'stations');
0010 nsta=length(stations);
0011 corner=[0 0 0 0];
0012 edges=[nt, 1, 1, nsta];
0013 e=mcvgt(cdf,'heat_flux',corner,edges);
0014 %
0015 base_date=zeros(1,6);
0016 base_date(1:3)=mcagt(cdf,'global','base_date');
0017 jd0=julian(base_date);
0018 jd=jd0+t/(3600*24);