Project

General

Profile

« Previous | Next » 

Revision b83b17ae

Added by Adam Wilson over 11 years ago

Explored 30-day moving average

View differences:

climate/procedures/MOD35_Climatology.r
95 95 

  		




  96
  96
  
    
myear=as.integer(max(fdly$year))  #this year will be used in all dates of monthly climatologies (and day will = 15)


  		




  97
  97
  
    

  		




  
  98
  
    
## subset dates


  		




  
  99
  
    
## due to bug (?) in CDO tools, only 10 years of data can be processed at a time or strange areas of NAs appear.


  		




  
  100
  
    
datesubset="-seldate,2002-01-01,2012-01-01"


  		




  
  101
  
    

  		




  98
  102
  
    
## Monthly means


  		




  99
  103
  
    
if(verbose) print("Calculating the monthly means")


  		




  100
  
  
    
system(paste("cdo -O sorttimestamp -setyear,",myear," -setday,15 -mulc,-1 -subc,100 -ymonmean ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_ymonmean.nc",sep=""),wait=T)


  		




  
  104
  
    
system(paste("cdo -O sorttimestamp -setyear,",myear," -setday,15 -mulc,-1 -subc,100 -ymonmean  ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_ymonmean.nc",sep=""),wait=T)


  		




  
  105
  
    

  		




  
  106
  
    
#system(paste("cdo -v -O selindexbox,1,100,1100,1200 ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_dailysmall.nc",sep=""),wait=T)


  		




  
  107
  
    
#system(paste("cdo -v -O sorttimestamp -setyear,",myear," -setday,15 -mulc,-1 -subc,100 -ymonmean  ",outdir2,"/MOD35_",tile,"_dailysmall.nc ",tsdir,"/MOD35_",tile,"_ymonmean.nc",sep=""),wait=T)


  		




  
  108
  
    

  		




  
  109
  
    
system(paste("cdo -O sorttimestamp -setyear,",myear," -setday,15 -mulc,-1 -subc,100 -ymonmean  ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_ymonmean.nc",sep=""),wait=T)


  		




  
  110
  
    

  		




  
  111
  
    
system(paste("cdo -O sorttimestamp -setyear,",myear," -setday,15 -mulc,-1 -subc,100 -monmean ",datesubset," ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_monmean.nc",sep=""),wait=T)


  		




  
  112
  
    

  		




  
  113
  
    

  		




  
  114
  
    
system(paste("cdo -O sorttimestamp -setyear,",myear," -mulc,-1 -subc,100 -ydrunmean,30 ",datesubset," ",outdir2,"/MOD35_",tile,"_daily.nc ",tsdir,"/MOD35_",tile,"_ydrunmean30.nc",sep=""),wait=T)


  		




  
  115
  
    

  		




  101
  116
  
    
system(paste(ncopath,"ncrename -v PClear,PCloud ",tsdir,"/MOD35_",tile,"_ymonmean.nc",sep=""))


  		




  102
  117
  
    
system(paste(ncopath,"ncatted ",


  		




  103
  118
  
    
" -a long_name,PCloud,o,c,\"Mean Probability of Cloud\" ",


  		












Also available in:  Unified diff