doc/plots.rwg.R

In wendellopes/rvswf: RVSWF: Expansion in Vector Spherical Wave Functions written in R

#-------------------------------------------------------------------------------
# NAMES
#-------------------------------------------------------------------------------
md<-ifelse(TM,"tm","te")
nm.pwe<-ifelse(lmax<10, paste("rwg.vswf.",md,".0",lmax,sep=""),
                        paste("rwg.vswf.",md,"." ,lmax,sep=""))
nm.vwf<-paste("rwg.vwfd.",md,".00",sep="")
nm.vwf.i<-paste(nm.vwf,".png",sep="")
nm.vwf.d<-paste(nm.vwf,".Rdata",sep="")
nm.pwe.i<-paste(nm.pwe,".png",sep="")
nm.pwe.d<-paste(nm.pwe,".Rdata",sep="")
#-------------------------------------------------------------------------------
# IMAGE
#-------------------------------------------------------------------------------
source("Image.R")
if(TM){
   zl<-range(Re(tez.RWG))
   #1
   if(!file.exists(nm.vwf.i)){
      png(nm.vwf.i)
      Image((y+yo)/lambda,(x+xo)/lambda,z=Re(tez.RWG),nlevels=256,axes=TRUE,color.palette=cm.colors,#zlim=zl,
         plot.axes={axis(1);axis(2);abline(h=xo/lambda,v=yo/lambda,col='green')},
         xlab=expression(y/lambda),ylab=expression(x/lambda))
      dev.off()
   }
   #2
   png(nm.pwe.i)
   Image((y+yo)/lambda,(x+xo)/lambda,z=Re(tez.PWE),nlevels=256,axes=TRUE,color.palette=cm.colors,#zlim=zl,
      plot.axes={axis(1);axis(2);abline(h=xo/lambda,v=yo/lambda,col='green')},
      xlab=expression(y/lambda),ylab=expression(x/lambda))
   dev.off()
}else{
   zl<-range(Re(thz.RWG))
   #1
   if(!file.exists(nm.vwf.i)){
      png(nm.vwf.i)
      Image((y+yo)/lambda,(x+xo)/lambda,z=Re(thz.RWG),nlevels=256,axes=TRUE,color.palette=cm.colors,#zlim=zl,
         plot.axes={axis(1);axis(2);abline(h=xo/lambda,v=yo/lambda,col='green')},
         xlab=expression(y/lambda),ylab=expression(x/lambda))
      dev.off()
   }
   #2
   png(nm.pwe.i)
   Image((y+yo)/lambda,(x+xo)/lambda,z=Re(thz.PWE),nlevels=256,axes=TRUE,color.palette=cm.colors,#zlim=zl,
      plot.axes={axis(1);axis(2);abline(h=xo/lambda,v=yo/lambda,col='green')},
      xlab=expression(y/lambda),ylab=expression(x/lambda))
   dev.off()
}
#-------------------------------------------------------------------------------
# DATASETS
#-------------------------------------------------------------------------------
if(TM){
   save(RWG, file=nm.vwf.d)
   save(PWE, file=nm.pwe.d)
}else{
   save(RWG, file=nm.vwf.d)
   save(PWE, file=nm.pwe.d)
}