Project

General

Profile

« Previous | Next » 

Revision 1dc46eb9

Added by Adam Wilson almost 11 years ago

Updating figures

View differences:

climate/procedures/MOD09_Visualize.R
45 45 

  		




  46
  46
  
    

  		




  47
  47
  
    
## climatologies


  		




  48
  
  
    
mac=brick("~/acrobates/adamw/projects/cloud/data/mod09_clim_mac.nc",varname="CF_annual")


  		




  
  48
  
    
mac=brick("~/acrobates/adamw/projects/cloud/data/cloud_mean.nc",varname="CF_annual")


  		




  49
  49
  
    

  		




  50
  50
  
    
pdf("output/mod09_climatology.pdf",width=11,height=8.5)


  		




  51
  51
  
    
levelplot(mac,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e6)+


  		




  ......




  82
  82
  
    
## reduced resolution


  		




  83
  83
  
    

  		




  84
  84
  
    
## read in GEWEX 1-degree data


  		




  85
  
  
    
gewex=raster("data/gewex/CA_PATMOSX_NOAA.nc")


  		




  
  85
  
    
gewex=mean(brick("data/gewex/CA_PATMOSX_NOAA.nc",varname="a_CA"))


  		




  86
  86
  
    

  		




  87
  87
  
    
mod09_8km=aggregate(mod09_mac,8)


  		




  88
  88
  
    

  		




  89
  89
  
    
pdf("output/mod09_resolution.pdf",width=11,height=8.5)


  		




  90
  90
  
    
p1=levelplot(mod09_mac,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  91
  
  
    
p2=levelplot(mod09_8km,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  92
  
  
    
p3=levelplot(mod09_1deg,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  93
  
  
    
print(c(p1,p2,p3,x.same=T,y.same=T,merge.legends=F))


  		




  
  91
  
    
#p2=levelplot(mod09_8km,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  92
  
    
p3=levelplot(gewex,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  93
  
    
print(c(p1,p3,x.same=T,y.same=T,merge.legends=F))


  		




  94
  94
  
    

  		




  95
  
  
    

  		




  96
  
  
    
p1=levelplot(crop(mod09_mac,reg2),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  97
  
  
    
p2=levelplot(crop(mod09_8km,reg2),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  98
  
  
    
p3=levelplot(crop(mod09_1deg,reg2),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  99
  
  
    
print(c(p1,p2,p3,x.same=T,y.same=T,merge.legends=F))


  		




  
  95
  
    
p1=levelplot(crop(mac,regs[["Venezuela"]]),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  96
  
    
#p2=levelplot(crop(mod09_8km,reg2),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  97
  
    
p3=levelplot(crop(gewex,regs[["Venezuela"]]),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  98
  
    
print(c(MOD09=p1,GEWEX=p3,x.same=T,y.same=T,merge.legends=F))


  		




  100
  99
  
    

  		




  101
  100
  
    
p1=levelplot(crop(mod09_mac,reg3),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  102
  
  
    
p2=levelplot(crop(mod09_8km,reg3),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  
  101
  
    
#p2=levelplot(crop(mod09_8km,reg3),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  103
  102
  
    
p3=levelplot(crop(mod09_1deg,reg3),col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,max.pixels=1e5)


  		




  104
  
  
    
print(c(p1,p2,p3,x.same=T,y.same=T,merge.legends=F))


  		




  
  103
  
    
print(c(p1,p3,x.same=T,y.same=T,merge.legends=F))


  		




  105
  104
  
    

  		




  106
  105
  
    
dev.off()


  		




  107
  106
  
    

  		












Also available in:  Unified diff