Project

General

Profile

« Previous | Next » 

Revision 1dc46eb9

Added by Adam Wilson almost 11 years ago

Updating figures

View differences:

climate/procedures/MOD09_CloudFigures.R
41 41 
projection(cldys)=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
42 42 

  		




  43
  43
  
    
#### Evaluate MOD35 Cloud data


  		




  44
  
  
    
mod09=brick("data/cloud_daily.nc")


  		




  
  44
  
    
mod09=brick("data/cloud_monthly.nc")


  		




  45
  45
  
    
mod09c=brick("data/cloud_ymonmean.nc",varname="CF");names(mod09c)=month.name


  		




  46
  46
  
    
mod09a=brick("data/cloud_mean.nc",varname="CF_annual")#;names(mod09c)=month.name


  		




  47
  47
  
    

  		




  ......




  76
  76
  
    

  		




  77
  77
  
    
## Figure 1: 4-panel summaries


  		




  78
  78
  
    
#- Annual average


  		




  79
  
  
    
levelplot(mod09a,col.regions=colr(100),cuts=100,at=seq(0,100,len=100),colorkey=list(space="bottom",adj=1),


  		




  
  79
  
    
levelplot(mod09a,col.regions=colr(n),cuts=100,at=seq(0,100,len=100),colorkey=list(space="bottom",adj=1),


  		




  80
  80
  
    
  margin=F,maxpixels=1e6,ylab="Latitude",xlab="Longitude",useRaster=T)+


  		




  81
  81
  
    
  layer(sp.lines(coast,col="black"),under=F)


  		




  82
  82
  
    
#- Monthly minimum


  		




  83
  83
  
    
#- Monthly maximum


  		




  84
  84
  
    
#- STDEV or Min-Max


  		




  85
  
  
    
p_mac=levelplot(mac,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),xlab="",ylab="",main=names(regs)[r],useRaster=T)


  		




  86
  
  
    
p_min=levelplot(mod09min,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  87
  
  
    
p_max=levelplot(mod09max,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  88
  
  
    
p_sd=levelplot(mod09sd,col.regions=grey(seq(0,1,len=100)),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  
  85
  
    
p_mac=levelplot(mod09a,col.regions=colr(n),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),xlab="",ylab="",main=names(regs)[r],useRaster=T)


  		




  
  86
  
    
p_min=levelplot(mod09min,col.regions=colr(n),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  
  87
  
    
p_max=levelplot(mod09max,col.regions=colr(n),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  
  88
  
    
p_sd=levelplot(mod09sd,col.regions=colr(n),cuts=99,margin=F,maxpixels=1e5,colorkey=list(space="bottom",height=.75),useRaster=T)


  		




  89
  89
  
    
p3=c("Mean Cloud Frequency (%)"=p_mac,"Max Cloud Frequency (%)"=p_max,"Min Cloud Frequency (%)"=p_min,"Cloud Frequency Variability (SD)"=p_sd,x.same=T,y.same=T,merge.legends=T,layout=c(2,2))


  		




  90
  90
  
    
print(p3)


  		




  91
  91
  
    

  		












Also available in:  Unified diff