Create time table
LE1202
I'm want to create a time table which interpolation in IGOR
I've try to calculate the seconds passed since 1904/01/01 00:00 by following code
But still won't work, I can't figure out what's going on
I new a table which show as
YEAR MONTH DAY HOUR MINUTE
2017 2 1 0 0
2017 2 1 1 0
2017 2 1 2 0
<br /> #pragma rtGlobals=1 // Use modern global access method.<br /> Macro test()<br /> <br /> Variable/G t1,t2,dt<br /> String sdate="20170201_0000" // Start Time: YYYYMMDD_HHMM<br /> String edate="20170201_0200" // End Time: YYYYMMDD_HHMM<br /> Variable inter_min=60 // Interpolation with N minute<br /> Variable inter_sec<br /> String ssy=sdate[0,3] , eey=edate[0,3]<br /> String ssm=sdate[4,5] , eem=edate[4,5]<br /> String ssd=sdate[6,7] , eed=edate[6,7]<br /> String ssh=sdate[9,10] , eeh=edate[9,10]<br /> String ssn=sdate[11,12] , een=edate[11,12]<br /> Variable sy=str2num(ssy) , ey=str2num(eey)<br /> Variable sm=str2num(ssm) , em=str2num(eem)<br /> Variable sd=str2num(ssd) , ed=str2num(eed)<br /> Variable sh=str2num(ssh) , eh=str2num(eeh)<br /> Variable sn=str2num(ssn) , en=str2num(een)<br /> <br /> en=en+inter_min<br /> <br /> If (inter_min==60)<br /> inter_sec=3600<br /> Endif<br /> <br /> t1=date2secs(sy,sm,sd)+sh*60*60+sn*60+0 //calculate start time seconds<br /> t2=date2secs(ey,em,ed)+eh*60*60+en*60+0 //calculate end time seconds<br /> <br /> dt=abs(t2-t1)/inter_sec // total row<br /> <br /> Make/D/N=(dt) secs,year,month,day,hour,minute<br /> Edit secs,year,month,day,hour,minute<br /> SetScale x,t1,t2,"dat",secs<br /> secs=x<br /> Variable num <br /> Variable allnum=dt<br /> <br /> Variable year_sec, month<br /> Variable totaldays,temp, temp_year, temp_month, temp_days, temp_hour, temp_min<br /> <br /> Do<br /> temp=(secs[num]/31557600) <br /> temp_year=round(temp+1904)<br /> year[num]=temp_year <br /> <br /> temp=date2secs(temp_year,01,01)+00*60*60+00*60+0<br /> temp=ceil((secs[num]-temp)/86400)+1<br /> print temp<br /> <br /> If(year_temp==2012 || year_temp==2016)<br /> totaldays=366<br /> If((temp>=1) && (temp<=31))<br /> temp_month=1<br /> temp_days=temp<br /> Endif<br /> if((temp>31) && (temp<=60))<br /> temp_month=2<br /> temp_days=temp-31<br /> endif<br /> if((temp>60) && (temp<=91))<br /> temp_month=3<br /> temp_days=temp-31-29<br /> endif<br /> Endif<br /> <br /> If((year_temp!=2012) || (year_temp!=2016))<br /> totaldays=365<br /> If((temp>=1) && (temp<=31))<br /> temp_month=1<br /> temp_days=temp<br /> endif<br /> if((temp>31) && (temp<=59))<br /> temp_month=2<br /> temp_days=temp-31<br /> endif<br /> if((temp>59) && (temp<=91))<br /> temp_month=3<br /> temp_days=temp-31-28<br /> Endif<br /> Endif<br /> <br /> month[num]=temp_month<br /> day[num]=temp_days<br /> temp=date2secs(temp_year,temp_month,temp_days)+00*60*60+00*60+0<br /> temp=(secs[num]-temp)/3600 <br /> temp_hour=floor(temp)<br /> hour[num]=temp_hour<br /> <br /> num+=1<br /> While(num<=dt)<br /> <br /> Endmacro<br /> <pre><code class="language-igor"></span>
<pre><code class="language-igor"></ccode> ... <ccode></code></pre>tags for Igor code. It's easier to read for us :-)"Won't work" is not a good description. Does it not compile? Is the result wrong? ...?
Maybe you want to have a look at
displayhelptopic "date2secs"That should do some of the work. Hours, minutes and seconds on top are easy.
Consider to use functions rather than macros.
HJ
PS: Timestamps in waves require double precision waves
April 25, 2017 at 07:36 am - Permalink