def PalumboPlots():
DataDirectory = "T://analysis_for_papers//Cosmo_paper//Palumbo_shield_plots//"
Filename = "SpawnedBasin_07C13-(Q8).bil"
Drapename = "SpawnedBasin_07C13-(Q8)_BASINS.bil"
ThisFile = DataDirectory+Filename
DrapeFile = DataDirectory+Drapename
LSDP.BasicDrapedPlotGridPlot(ThisFile, DrapeFile, thiscmap='terrain',drape_cmap='gray',
colorbarlabel='Elevation in meters',clim_val = (2000,5000),
drape_alpha = 0.4,FigFileName = 'Basin.pdf',FigFormat = 'pdf')
def TestNewMappingTools():
#DataDirectory = "C://basin_data//Chile//lat26p0//"
#Filename = "Site_lat26p0_UTM19_DEM_FILL.bil"
#DataDirectory = "C://basin_data//Model_results//June2015_Results//HighK//"
#DataDirectory = "T://test_clone//topodata//"
#DataDirectory = "T://analysis_for_papers//Cosmo_paper//Tibet//for_plotting//"
DataDirectory = "C://basin_data//CosmoPaper//DEMs//"
#Filename = "SanBern.bil"
Filename = "SpawnedBasin_07C13-(Q8)_SH.bil"
#Filename = "CRNvariable_long_0_0_1var128.asc"
#DrapeFileName = "CRNvariable_long_0_0_1var128_erosion.asc"
DrapeFileName = "SpawnedBasin_07C13-(Q8)_BASINS.bil"
DEMName = "SpawnedBasin_07C13-(Q8).bil"
ThisFile = DataDirectory + Filename
DrapeFile = DataDirectory + DrapeFileName
DEMFile = DataDirectory + DEMName
Shield2 = DataDirectory + "SpawnedBasin_Q8_phi30.bil"
FigFormat = 'svg'
FigFileN = 'Elevation.svg'
FigFileName = DataDirectory + FigFileN
ShieldFigName = DataDirectory + 'Shielding.svg'
NDV, xsize, ysize, GeoT, Projection, DataType = LSDP.GetGeoInfo(ThisFile)
print "NDV: " + str(NDV)
print "xsize: " + str(xsize)
print "ysize: " + str(ysize)
print "GeoT: " + str(GeoT)
print "Projection: " + str(Projection)
print "DataType: " + str(DataType)
CellSize, XMin, XMax, YMin, YMax = LSDP.GetUTMMaxMin(ThisFile)
print "CellSize: " + str(CellSize)
print "XMin: " + str(XMin)
print "XMax: " + str(XMax)
print "YMin: " + str(YMin)
print "YMax: " + str(YMax)
tcmap = 'jet'
drape_cmap = 'gray'
shield_map = 'summer'
drape_alpha = 0.4
tcmapcolorbarlabel = 'Topographic shielding'
clim_val = (0.72, 1)
auto_clim = (1000, 4500)
#LSDP.BasicDensityPlot(ThisFile,tcmap,tcmapcolorbarlabel,clim_val)
#LSDP.DrapedPlot(ThisFile,DrapeFile)
#LSDP.BasicDensityPlotGridPlot(ThisFile,tcmap,tcmapcolorbarlabel)
#Plot the basin over the elevation
cmap_label = "Elevation (m)"
LSDP.BasicDrapedPlotGridPlot(DEMFile, DrapeFile, tcmap, drape_cmap,
cmap_label, auto_clim, drape_alpha,
FigFileName, FigFormat)
# Now plot the two shielding rasters
LSDP.BasicDensityPlotGridPlot(ThisFile, shield_map, tcmapcolorbarlabel,
clim_val, ShieldFigName, FigFormat)
FigFormat = 'svg'
SNFileN = 'Shield30.svg'
ShieldFigName2 = DataDirectory + SNFileN
# Now plot the two shielding rasters
LSDP.BasicDensityPlotGridPlot(Shield2, shield_map, tcmapcolorbarlabel,
clim_val, ShieldFigName2, FigFormat)
