revised

da/stilt/stilt.co2.simu.2015.Sep.bc.r

@@ -520,20 +520,22 @@ for(mm in 1:length(pfbfns))
     #print(dateall)  #reversed order
     
     # note: ixflux should not be from 1:ntimes, should be more accurate in hours
-    dd=as.character(min(yy))
+    #dd=as.character(min(yy))
+    dd=as.character(xx)
     #yrlist=substring(dd,1,4)
     #monlist=substring(dd,6,7)
     #daylist=substring(dd,9,10)
     hr=substring(dd,12,13)
     mi=substring(dd,15,16)
 
-    n = ceiling(as.numeric(hr) + as.numeric(mi)/60)
+    n = 24-floor(as.numeric(hr) + as.numeric(mi)/60)
     print("XXXn") 
     print(n)
     datevalid=dateall[1]+0.5*3600-(0:(ntimes-1))*3600
     #ixflux=(1:ntimes)[(datevalid<=xx)][1:endtime]   #time backward
     ixflux=(1:ntimes)[(datevalid<=xx)][n:(endtime+n-1)]
-   
+    print(ixflux)
+    
     if(bioflux_flag == "SiB3")   #footprint is backward
     {
         gpptmp=foot*gppflux[,,ixflux]