Project

General

Profile

« Previous | Next » 

Revision bce95cb5

Added by Adam Wilson about 11 years ago

First global mod09 run. Updated cloud flag function

View differences:

climate/procedures/ee.MOD09.py
30 30 

  		




  31
  31
  
    
#///////////////////////////////////


  		




  32
  32
  
    
#// Function to extract cloud flags


  		




  33
  
  
    
def getmod09(img): return(img.select(['state_1km']).expression("((b(0)/1024)%2)")); 


  		




  34
  
  
    
 


  		




  
  33
  
    
def getmod09(img): return(img.select(['state_1km']).expression("((b(0)/1024)%2)>0.5")); 


  		




  
  34
  
    
# added the >0.5 because some values are coming out >1.  Need to look into this further as they should be bounded 0-1...


  		




  
  35
  
    

  		




  35
  36
  
    
#// Date ranges


  		




  36
  
  
    
yearstart=2001


  		




  37
  
  
    
yearstop=2001


  		




  
  37
  
    
yearstart=2000


  		




  
  38
  
    
yearstop=2012


  		




  38
  39
  
    
monthstart=1


  		




  39
  
  
    
monthstop=2


  		




  
  40
  
    
monthstop=12


  		




  40
  41
  
    

  		




  41
  42
  
    
#////////////////////////////////////////////////////


  		




  42
  43
  
    
# Loop through months and get monthly % missing data


  		




  43
  44
  
    

  		




  44
  45
  
    
## set a year-month if you don't want to run the loop (for testing)


  		




  45
  
  
    
year=2001


  		




  46
  
  
    
month=2


  		




  
  46
  
    
#year=2001


  		




  
  47
  
    
#month=2


  		




  47
  48
  
    

  		




  48
  49
  
    

  		




  49
  50
  
    
## define the regions to be processed


  		












Also available in:  Unified diff