#R-code
#Precipitation Fields creation

baseDirectory<-"/hidrosigDataBases/Whitewater_database/Rasters/Hydrology/storms/"
#baseDirectory<-"/hidrosigDataBases/Walnut_Gulch_AZ_database/Rasters/Hydrology/storms/"
#baseDirectory<-"/hidrosigDataBases/FractalTrees_database/Rasters/Hydrology/storms/"

stormDirectory<-"simulated_events/"
fileLabel<-"uniform_100_10"
nCols<-"20"
nRows<-"20"
dir.create(paste(baseDirectory,stormDirectory,sep=""))

thisFile<-paste(baseDirectory,stormDirectory,fileLabel,".metaVHC",sep="")

parameters=c(	"[Name]",
	    	"[Southernmost Latitude]",
		"[Westernmost Longitude]",
		"[Longitudinal Resolution (ArcSec)]",
		"[Latitudinal Resolution (ArcSec)]",
		"[# Columns]",
		"[# Rows]",
		"[Format]",
		"[Missing]",
		"[Temporal Resolution]",
		"[Units]",
		"[Information]")
	    
#Whitewater Responses
reponses=c(	"Simulated Precipitation Field",
	    	"37:39:28.50 N",
		"97:18:28.49 W",
		"120",
		"120",
		nCols,
		nRows,
		"float",
		"-9999",
		"10-minutes",
		"mm/h",
		"Simulated fields using R-program")

#Walnut Gulch Responses
#reponses=c(	"Simulated Precipitation Field",
#	    	"31:37:30.00 N",
#		"110:15:00.00 W",
#		"120",
#		"120",
#		nCols,
#		nRows,
#		"float",
#		"-9999",
#		"10-minutes",
#		"mm/h",
#		"Simulated fields using R-program")


#Test DEM Responses
#reponses=c(	"Simulated Precipitation Field",
#	    	"37:32:39.00 N",
#		"96:41:47.98 W",
#		"120",
#		"120",
#		nCols,
#		nRows,
#		"float",
#		"-9999",
#		"1-minutes",
#		"mm/h",
#		"Simulated fields using R-program")

#Fractal DEM Responses
#reponses=c(	"Simulated Precipitation Field",
#	    	"00:00:00.00 N",
#		"00:00:00.00 W",
#		"1200",
#		"1200",
#		nCols,
#		nRows,
#		"float",
#		"-9999",
#		"1-minutes",
#		"mm/h",
#		"Simulated fields using R-program")

filePipe<-file(thisFile,"w")

for (i in 1:length(parameters)){
	cat(parameters[i], file = filePipe,sep="\n")
	cat(reponses[i],"\n", file = filePipe,sep="\n")
}

close(filePipe)

#Event Parameters
Duration<-1.0		#minutes
Intensity<-100.0 	#mm/h


thisFile<-paste(baseDirectory,stormDirectory,fileLabel,".180000.01.July.1971.vhc",sep="")

filePipe<-file(thisFile,"wb")
rain<-Intensity*rep(1.0,as.integer(nCols)*as.integer(nRows))
#rain<-Intensity*runif(29*38)
writeBin(rain,filePipe,size=4,endian="big")
close(filePipe)