Project

General

Profile

« Previous | Next » 

Revision d91f0663

Added by Adam Wilson about 11 years ago

nearly finished first global run of MOD35 cloud frequency

View differences:

climate/procedures/MOD35_Climatology.r
34 34 
 fdly=data.frame(path=list.files(outdir,pattern="nc$",full=T),stringsAsFactors=F)
35 35 
  fdly$file=basename(fdly$path)
36 36 
  fdly$dateid=substr(fdly$file,14,21)
37 
  fdly$date=as.Date(substr(fdly$file,14,21),"%Y%m%d")
37 
  fdly$date=as.Date(substr(fdly$file,14,21),"%Y%m%d") #23,30 for alltiles
38 38 
  fdly$month=format(fdly$date,"%m")
39 39 
  fdly$year=format(fdly$date,"%Y")
40 40 
nrow(fdly)
......
55 55 

  		




  56
  56
  
    
## merge all daily files to create a single file with all dates


  		




  57
  57
  
    
system(paste(ncopath,"ncrcat -O ",outdir,"/*nc ",outdir2,"/MOD35_",tile,"_daily.nc",sep=""))


  		




  
  58
  
    
#system(paste(ncopath,"ncrcat -O ",outdir,"/*alltests*nc ",outdir2,"/MOD35_",tile,"_daily_alltiles.nc",sep=""))


  		




  58
  59
  
    
system(paste("ncdump -h ",outdir2,"/MOD35_",tile,"_daily.nc",sep=""))


  		




  59
  60
  
    
 


  		




  60
  61
  
    
## Update attributes


  		












Also available in:  Unified diff