diff --git a/da/stilt/stilt.co2.simu.2015.Sep.bc.r b/da/stilt/stilt.co2.simu.2015.Sep.bc.r index c985d1c6d6a9681027d5c44ace8f8629832fae28..ea7a6d8895d2cdd118d15d5656dd877ae36cea16 100755 --- a/da/stilt/stilt.co2.simu.2015.Sep.bc.r +++ b/da/stilt/stilt.co2.simu.2015.Sep.bc.r @@ -556,12 +556,21 @@ flog.info('eventids: %s',eventids[mm], name='logger.b') #June 23 hr=substring(dd,12,13) mi=substring(dd,15,16) + flog.info('dd: %s, length dd: %s,hr: %s, mi: %s',dd,length(dd),hr,mi,name='logger.b') + n = 24-floor(as.numeric(hr) + as.numeric(mi)/60) + if (nchar(dd,type='width') ==10) + { + n=24 + flog.info('No hr and mi, set n to 24',name='logger.b') + } + 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)][1:(endtime+n)] + flog.info('n: %s, endtime: %s',n,endtime,name='logger.b') ixflux <-array(n:(endtime+n-1),dim=c(240)) flog.info('ixflux endtime: %s,%s,%s',endtime,n,as.numeric(hr)+(as.numeric(hr)/60),name='logger.b')