pro doitquick
	if file_test('fullNDVI') then begin
		NDVIfile='fullNDVI'
		print, NDVIfile
	endif
	if file_test('fullMasktmp') then begin
		maskfile='fullMasktmp'
		print, maskfile
	endif
	if file_test('classNDVItmp') then begin
		classfile = 'classNDVItmp'
		envistart
		envi_open_file, classfile, r_fid=c_fid
		envi_file_query, c_fid, num_classes=nclasses
		envi_file_mng, id=c_fid, /remove
		print, classfile, nclasses
	endif
	if file_test('spatCent.out') then begin
		spatCfile='spatCent.out'
		print, spatCfile
	endif
	if file_test('windMap.out') then begin
		windFile='windMap.out'
		print, windFile
	endif
	if file_test('NvCcorrelation') then begin
		corrFile='NvCcorrelation'
		print, corrFile
	endif
	if file_test('FutureVegMapCorrected') then begin
		FVcorrected='FutureVegMapCorrected'
		print, FVcorrected
	endif
	HighPlainsDuneModel, 'ModelInput', /nocent, fullMask=maskfile, $
		fullNDVI=NDVIfile, n_Classes=nclasses,classFile=classfile, $
		spatCfile=spatCfile, windFile=windFile, corrFile=corrFile, $
		FVCmap=FVcorrected
end		

pro doitREALLYFAST
	if not file_test('combmeta') then return
	fastmodel, 'combmeta'
end


;; currently set up to test the effects of 6-41 maximum classes in increments of 5
pro doitmoderate
	for max_classes=6,41,5 do begin
	if file_test('fullNDVI') then begin
		NDVIfile='fullNDVI'
		print, NDVIfile
	endif
	if file_test('fullMasktmp') then begin
		maskfile='fullMasktmp'
		print, maskfile
	endif
;	if file_test('classNDVItmp') then begin
;		classfile = 'classNDVItmp'
;		envistart
;		envi_open_file, classfile, r_fid=c_fid
;		envi_file_query, c_fid, num_classes=nclasses
;		envi_file_mng, id=c_fid, /remove
;		print, classfile, nclasses
;	endif
	if file_test('spatCent.out') then begin
		spatCfile='spatCent.out'
		print, spatCfile
	endif
	if file_test('windMap.out') then begin
		windFile='windMap.out'
		print, windFile
	endif

	mc=max_classes
	HighPlainsDuneModel, 'ModelInput', /nocent, fullMask=maskfile, $
		fullNDVI=NDVIfile, $; n_Classes=nclasses, classFile=classfile, $
		spatCfile=spatCfile, windFile=windFile, max_classes=mc
	endfor
end

pro batchModel, fast=fast, REALLYFAST=REALLYFAST, nobasics=nobasics
	scenes=findfile('*')
	cd, current=masterDIR
	print, scenes
	for i=0, n_elements(scenes)-1 do begin
		print, scenes(i)

		current=strtok(scenes(i),':', /extract)
		if (strtok(scenes(i), '.'))(0) eq 0 and file_test(current(0), /directory) then begin
			cd, current(0)
			print, ' ', string(byte(13))
			print, '**************************',current(0),'**************************'
			if keyword_set(REALLYFAST) then begin
				doitREALLYFAST
			endif else if keyword_set(fast) then begin
				doitquick
			endif else if keyword_set(nobasics) then begin
				doitmoderate
			endif else HighPlainsDuneModel, 'ModelInput', /nocent
			
			close, /all
			cd, masterDIR
		endif
	endfor
	
end