Lines Matching refs:k
188 do k=1,ntv
361 do k=1,numImpSrfs
363 numImpt(k) = numDataImp
365 read(817,*) (ValueImpt(j,n,k),n=1,2) ! n=1 time, 2 value
385 integer nlast, nper, k, j , jstep local
388 do k =1,numImpSrfs
389 nlast=numImpt(k) ! number of time series to interpolate from
390 nper=ctime/ValueImpt(nlast,1,k)!number of periods completed to shift off
391 ptime = ctime-nper*ValueImpt(nlast,1,k) ! now time in periodic domain
395 if(ValueImpt(j,1,k).gt.ptime) then ! this is upper bound, j-1 is lower
396 wr=(ptime-ValueImpt(j-1,1,k))
397 & / ( ValueImpt(j,1,k)-ValueImpt(j-1,1,k) )
398 ValueListImp(jstep,k)= ValueImpt(j-1,2,k)*(one-wr)
399 & + ValueImpt(j,2,k)*wr
416 integer nptf, cutfreq, j, k, m, s, Filtime(nptf) local
465 do k=1,numRCRSrfs
467 numRCRt(k) = numDataRCR
469 read(818,*) ValueListRCR(j,k) ! reads Rp,C,Rd
472 read(818,*) (ValuePdist(j,n,k),n=1,2) ! n=1 time, 2 value
489 integer nlast, nper, k, j local
492 do k =1,numRCRSrfs
493 nlast=numRCRt(k) ! number of time series to interpolate from
494 nper=ctime/ValuePdist(nlast,1,k)!number of periods completed to shift off
495 ptime = ctime-nper*ValuePdist(nlast,1,k) ! now time in periodic domain
499 if(ValuePdist(j,1,k).gt.ptime) then ! this is upper bound, j-1 is lower
500 wr=(ptime-ValuePdist(j-1,1,k))
501 & / ( ValuePdist(j,1,k)-ValuePdist(j-1,1,k) )
502 Pdist(k)= ValuePdist(j-1,2,k)*(one-wr)
503 & + ValuePdist(j,2,k)*wr
526 do k=1,nSrfs
528 pressHist(k) = pressHist(k) + QHist(j,k)*betas(j,k)