@@ -151,8 +151,8 @@ plotsurf <- function(freqs=0,
151151 # make log plot to expand lower values. Plot 0 as 0.1
152152 badcw [zcw ,nfreq ] <- 0.1
153153 }
154- plot( badcw [, nfreq ], type = " b " , xlim = c( t1 , t2 ),
155- log = ifelse(logbad ," y" ," " ))
154+
155+ plot( badcw [, nfreq ], type = " b " , xlim = c( t1 , t2 ), log = ifelse(logbad ," y" ," " ))
156156 if (! do_dat ) {
157157 timeaxis(3 , labels = TRUE , time.zone = badcw @ time.zone , date.too = FALSE ,
158158 xlab = FALSE )
@@ -165,30 +165,30 @@ plotsurf <- function(freqs=0,
165165 pcu [zcw ,] <- 0
166166 colnames(pcu ) <- " UncorrCW"
167167 units(pcu ) <- " %"
168+
168169 plot(pcu ,type = " b" ,xlim = c(t1 ,t2 ))
169170 if (! do_dat ) {
170- timeaxis(3 , labels = FALSE , time.zone = badcw @ time.zone , date.too = FALSE ,
171+ timeaxis(3 , labels = FALSE , time.zone = pcu @ time.zone , date.too = FALSE ,
171172 xlab = FALSE )
172173 axis(4 )
173174 }
174175
175176 plot(snr [,nfreq ],type = " b" ,xlim = c(t1 ,t2 ))
176177 if (! do_dat ) {
177- timeaxis(3 , labels = FALSE , time.zone = badcw @ time.zone , date.too = FALSE ,
178+ timeaxis(3 , labels = FALSE , time.zone = snr @ time.zone , date.too = FALSE ,
178179 xlab = FALSE )
179180 axis(4 )
180181 }
181182
182183 if (! do_dat ) {
183184 plot(pow [,nfreq ],type = " b" ,xlim = c(t1 ,t2 ))
184- timeaxis(3 , labels = FALSE , time.zone = badcw @ time.zone , date.too = FALSE ,
185+ timeaxis(3 , labels = FALSE , time.zone = pow @ time.zone , date.too = FALSE ,
185186 xlab = FALSE )
186187 axis(4 )
187188 title(main = titlestr , line = - par(" cex.main" ), outer = TRUE )
188189 }
189190 else plot(pow [,nfreq ],type = " b" ,xlim = c(t1 ,t2 ), title = titlestr )
190191
191-
192192 par(ask = ask )
193193 }
194194
@@ -279,8 +279,8 @@ plotsurf <- function(freqs=0,
279279 zcw <- ! is.na(badcw [,1 ]) & badcw [,1 ] == 0
280280 badcw [zcw ,1 ] <- 0.1
281281 }
282- plot( badcw [, 1 ], type = " b " , xlim = c( t1 , t2 ),
283- log = ifelse(logbad ," y" ," " ))
282+
283+ plot( badcw [, 1 ], type = " b " , xlim = c( t1 , t2 ), log = ifelse(logbad ," y" ," " ))
284284 if (! do_dat ) {
285285 timeaxis(3 , time.zone = badcw @ time.zone , date.too = FALSE , xlab = FALSE )
286286 axis(side = 4 )
@@ -292,6 +292,7 @@ plotsurf <- function(freqs=0,
292292 pcu [zcw ,] <- 0
293293 colnames(pcu ) <- " UncorrCW"
294294 units(pcu ) <- " %"
295+
295296 plot(pcu , type = " b" ,xlim = c(t1 ,t2 ))
296297 if (! do_dat ) {
297298 timeaxis(3 , labels = FALSE , time.zone = badcw @ time.zone , date.too = FALSE )
@@ -309,6 +310,7 @@ plotsurf <- function(freqs=0,
309310 pow @ data [,1 ] <- apply(pow , 1 , function (x ) { mean(x , na.rm = T ) })
310311
311312 titlestr <- " All frequencies"
313+
312314 plot(pow [,1 ], type = " b" , xlim = c(t1 ,t2 ))
313315 if (! do_dat ) {
314316 timeaxis(3 , labels = FALSE , time.zone = badcw @ time.zone , date.too = FALSE )
0 commit comments