'
'----------------------------
'--- Creation information ---
'----------------------------
'
'Name: loadgrid.ave
'Version: 1.0
'Date: 02/28/97
'Author: Ann Quenzer
'        Center for Research in Water Resources
'        The University of Texas at Austin
'        quenzer@mail.utexas.edu
'
'---------------------------
'--- Purpose/Description ---
'---------------------------
'
'Makes a land surface loading grid (Kg/day) to be used in a weighted 
'flow accumulation which determines the loadings to the bay.
'
'----------------
'--- Get view ---
'----------------
'
theview = av.getactivedoc
'
'-------------------------------
'--- Get themes for loadgrid ---
'-------------------------------
'
thethemes = theview.getthemes
if (thethemes.count = 0) then
    msgbox.error("No themes found", "LOAD GRID")
    exit
end
thegthemes = list.make
for each thetheme in thethemes
    if (thetheme.getclass.getclassname = "gtheme") then
        thegthemes.add(thetheme)
    end
end
if (thegthemes.count = 0) then
    msgbox.error("No grid themes found", "LOAD GRID")
    exit
end
conctheme = msgbox.listasstring(thegthemes, "Concentration Grid", "LOAD GRID")
if (conctheme = nil) then
    exit
end
'
'
thethemes = theview.getthemes
if (thethemes.count = 0) then
    msgbox.error("No themes found", "LOAD GRID")
    exit
end
thegthemes = list.make
for each thetheme in thethemes
    if (thetheme.getclass.getclassname = "gtheme") then
        thegthemes.add(thetheme)
    end
end
if (thegthemes.count = 0) then
    msgbox.error("No grid themes found", "LOAD GRID")
    exit
end
rotheme = msgbox.listasstring(thegthemes, "Runoff Grid", "LOAD GRID")
if (rotheme = nil) then
    exit
end
'
'
'
'-----------------
'--- Calculate ---
'-----------------
'
concgrid = conctheme.getgrid
rogrid = rotheme.getgrid
'
cellsize = number.makenull
dummmy = grid.getanalysiscellsize(cellsize)
extent = rect.makenull
dummy = grid.getanalysisextent(extent)
'
loadgrid = ((concgrid * rogrid) / 36525.asgrid)
loadfilename = av.getproject.makefilename("loadgrid", "")
loadgrid.savedataset("/home/ferdi/quenzer/trialrun/loadgrid".asfilename)
loadgtheme = gtheme.make(loadgrid)
theview.addtheme(loadgtheme)
loadgtheme.setvisible(true)
'
'final message to user
'
message = "Load grid calculated."
msgbox.info(message,"LOAD GRID")
'
'
'-----------
'--- End ---
'-----------
'