alltiles=system("ls -r MODTILES/ | grep tif$ | cut -c1-6 | sort | uniq - ",intern=T)

## or run all tiles

## subset to tiles in global region (not outside global boundary in sinusoidal projection)

tiles=tiles[tiles%in%alltiles]

queue="normal" #"devel"

queue="long" #"devel"

nodes=120

### write qsub script to disk

cat(paste("

queue="devel"

walltime=2

### qsub script