ratty, you highlight the proof that a specific degree in any particular area is a piece of paper that means nothing. Is there a single problem with his methods or analysis. Unlike your PHD Con Mann, the source of the data used and the programs used to create the plot are published in the article.
Technology and science ignoramus like you and Shepard have no clue as to what and how of Willis's analysis and charting. So as sniveling anonymous cowards you attack the individual. How pathetic and lame.
Fact, no IPCC models predicted the 15yr pause in warming. No IPCC model shows and of the Solar Systems cycles. IPCC predicted a false warming of the tropical upper troposphere due directly too how CO2 was going create increased radiative forcings that are all part of the
FALSIFICATION OF THE ATMOSPHERIC CO2 GREENHOUSE EFFECTS WITHIN THE FRAME OF PHYSICS
All you need do is study ratfor to understand.
watson@xen1[332]cat Is\ the\ troposphere\ worrisome.R
require("ncdf")
source('~/plotts2 function.R')
radians=function(x) x/180*pi
dsin=function(x) sin(radians(x))
# get data
mync=open.ncdf("tlt.nc")
thetempsraw=aperm(get.var.ncdf(mync,"Tlt"),c(2,1,3))
# calc variables
theheight=dim(thetempsraw)[3]
thetemps=thetempsraw[72:1,,]
thelats=-mync$dim$lat$val
thelongs=mync$dim$lon$val
latmatrix=matrix(rep(thelats,144),nrow=72)
cosmatrix=cos(latmatrix*pi/180)
latlist=list(c(67,90),c(23,67),c(-23,23),c(-67,-23),c(-90,-67))
# calculate area-weighted averages of zones
theanswer=matrix(NA,dim(thetempsraw)[3],5)
dim(theresults)
for (i in 1:length(latlist)){
j=latlist[]
mylat=which((thelats > j[1]) & (thelats <= j[2]))
mylong=which((thelongs > -180) & (thelongs < 180))
theresults=apply(thetemps[mylat,mylong,],3,function(x) weighted.mean(x,cosmatrix[mylat,mylong]))
themeans=function(x) weighted.mean(x,cosmatrix)
tempts=ts(theresults,start=c(1979,1),frequency=12)
plot(tempts,ylim=c(-3,3))
myloess=loess(as.vector(tempts)~time(tempts),span=0.6)
lines(myloess$fitted~as.vector(time(tempts)),col="blue")
theanswer[,i]=theresults
Sys.sleep(1)
}
#calculate areas as percent of globe
round((dsin(90)-dsin(67))/2*100,1)
round((dsin(67)-dsin(23))/2*100,1)
round((dsin(23)-dsin(-23))/2*100,1)
answertsraw=ts(theanswer,start=c(1979,1),frequency=12)
answerts=answertsraw-matrix(rep(colMeans(answertsraw),theheight),nrow=theheight,byrow=TRUE)
colnames(answerts)=c("N Pol (4%)","N Ext (26%)","Trop (40%)","S Ext (26%)","S Pol (4%)")
par(lwd=2)
par(col="black",cex.axis=1.3)
par(yaxp=c(-2,2,4))
plot.ts2(answerts,nc=1, ylim=c(-3,3),main="",col="salmon",yax.flip=TRUE)
title(main="UAH MSU Lower Troposphere Temperature",line=2.6)
title(main="Tropics -23° to 23°, Extratropics 23° to 67°, Poles 67° to 90°",
line=1.4,cex.main=1) |
