pro fastmodel, inputfile
	envistart
	maskfile='fullMasktmp'
	windFile='windMap.out'
	vegMapfile = 'meanVeg'
	
	if not file_test(maskfile) or not file_test(windFile) then return
	
	fullNDVI='fullNDVI'
;	if not file_test(fullNDVI) then $
		combunequal, inputfile, fullNDVI

	maskfile='fullMasktmp'
	windFile='windMap.out'
	vegMapfile = 'meanVeg'

;	if not file_test(vegMap) then $
		ndviToVeg, fullNDVI, vegMapfile, maskfile

	maskfile='fullMasktmp'
	windFile='windMap.out'
	vegMapfile = 'meanVeg'
	sandFile = 'SandTransportMap'
	sandmodel, vegMapfile, windFile, sandFile

	maskfile='fullMasktmp'
	windFile='windMap.out'
	vegMapfile = 'meanVeg'
		
	scaleFactor = 10
	;if not file_test('sandRS') then 
	print, scaleDown(sandFile, 'sandRS', scaleFactor)
	;if not file_test('windRS') then 
;	print, scaleDown(windFile, 'windRS', scaleFactor)
	;if not file_test('vegRS')  then 
	print, scaleDown(vegMapfile, 'vegRS', scaleFactor)
	;if not file_test('maskRS') then 
	print, scaleDown(maskfile, 'maskRS', scaleFactor)
	;if not file_test('ndviRS') then 
	print, scaleDown(fullNDVI, 'ndviRS', scaleFactor)
	
	file_delete, vegMapfile
	file_delete, sandFile
	file_delete, fullNDVI

end