Project

General

Profile

« Previous | Next » 

Revision d86b0a4a

Added by Adam Wilson almost 11 years ago

updates to MOD09 visualizations

View differences:

climate/procedures/ee_compile.R
76 76 

  		




  77
  77
  
    

  		




  78
  78
  
    
### generate the monthly mean and sd


  		




  79
  
  
    
system(paste("cdo -O merge -ymonmean data/mod09.nc -chname,CF,CF_sd -ymonstd data/mod09.nc data/mod09_clim.nc"))


  		




  80
  
  
    
system(paste("cdo -O -ymonmean data/mod09.nc data/mod09_clim2.nc"))


  		




  
  79
  
    
#system(paste("cdo -P 10 -O merge -ymonmean data/mod09.nc -chname,CF,CF_sd -ymonstd data/mod09.nc data/mod09_clim.nc"))


  		




  
  80
  
    
system(paste("cdo  -O -ymonmean data/mod09.nc data/mod09_clim_mean.nc"))


  		




  
  81
  
    
system(paste("cdo  -O -chname,CF,CF_sd -ymonstd data/mod09.nc data/mod09_clim_sd.nc"))


  		




  81
  82
  
    

  		




  82
  83
  
    
#  Overall mean


  		




  83
  
  
    
system(paste("cdo -O  -chname,CF,CF_annual -timmean data/mod09.nc  data/mod09_clim2.nc"))


  		




  
  84
  
    
system(paste("cdo -O  -chname,CF,CF_annual -timmean data/mod09.nc  data/mod09_clim_mac.nc"))


  		




  84
  85
  
    

  		




  85
  86
  
    
### Long term summaries


  		




  86
  87
  
    
seasconc <- function(x,return.Pc=T,return.thetat=F) {


  		




  ......




  114
  115
  
    

  		




  115
  116
  
    

  		




  116
  117
  
    
## read in monthly dataset


  		




  117
  
  
    
mod09=brick("data/mod09_clim.nc",varname="CF")


  		




  
  118
  
    
mod09=brick("data/mod09_clim_mean.nc",varname="CF")


  		




  118
  119
  
    
plot(mod09[1])


  		




  119
  120
  
    

  		




  120
  121
  
    
mod09_seas=calc(mod09,seasconc,return.Pc=T,return.thetat=F,overwrite=T,filename="data/mod09_seas.nc",NAflag=255,datatype="INT1U")


  		




  
  122
  
    
mod09_seas2=calc(mod09,seasconc,return.Pc=F,return.thetat=T,overwrite=T,filename="data/mod09_seas_theta.nc",datatype="INT1U")


  		




  121
  123
  
    

  		




  122
  124
  
    
plot(mod09_seas)


  		












Also available in:  Unified diff