|
1
|
## script to visualize cloud frequency data
|
|
2
|
|
|
3
|
setwd("~/acrobates/adamw/projects/cloud/")
|
|
4
|
|
|
5
|
library(rasterVis)
|
|
6
|
|
|
7
|
|
|
8
|
#### Evaluate MOD35 Cloud data
|
|
9
|
mod09=brick("~/acrobates/adamw/projects/cloud/data/mod09.nc")
|
|
10
|
NAvalue(mod09)=-1
|
|
11
|
|
|
12
|
cols=colorRampPalette(c("black","blue","red"))
|
|
13
|
|
|
14
|
r=mod09[[10]]
|
|
15
|
levelplot(r,col.regions=cols(100),at=seq(0,100,len=100),margin=F,maxpixels=1e6)
|
|
16
|
|
|
17
|
|
|
18
|
|
|
19
|
## climatologies
|
|
20
|
mod09c=brick("~/acrobates/adamw/projects/cloud/data/mod09_clim.nc",varname="CF")
|
|
21
|
levelplot(mod09c,col.regions=cols(100),at=seq(0,100,len=100),margin=F,max.pixels=1e7)
|
|
22
|
|
|
23
|
|
|
24
|
|
|
25
|
|
