NCAR/wrf_hydro_gis_preprocessor

Error while using python Build_Routing_Stack.py

Opened this issue · 0 comments

ved-ux commented

import Build_Routing_Stack

Define script input parameters using python variables

in_geogrid = os.path.join(data_folder, 'geo_em.d01.nc')
lakes = os.path.join(data_folder, 'lake_shapes', 'lakes.shp')
#csv = os.path.join(data_folder, 'croton_frxst_pts_FOSS.csv')
in_dem = os.path.join(data_folder, 'geo_em.d01_boundary_dem.tif')
regrid_factor = 4
routing_cells = 20
out_zip = os.path.join(output_folder, 'croton_test.zip')

Print information to screen for reference

print('Command to run:\n')
print('python Build_Routing_Stack.py \\n\t -i {0} \\n\t -l {1} \\n\t --CSV {2} \\n\t -d {3} \\n\t -R {4} \\n\t -t {5} \\n\t -o {6}\n'.format(in_geogrid, lakes, csv, in_dem, regrid_factor, routing_cells, out_zip))

Run the script with required parameters

#! python Build_Routing_Stack.py -i "{in_geogrid}" -l "{lakes}" --CSV "{csv}" -d "{in_dem}" -R "{regrid_factor}" -t "{routing_cells}" -o "{out_zip}"

Run the script with required parameters

! python Build_Routing_Stack.py -i "{in_geogrid}" -l "{lakes}" -d "{in_dem}" -R "{regrid_factor}" -t "{routing_cells}" -o "{out_zip}"

Terrain processing step initiated...
Using ./whitebox_tools: /lib/x86_64-linux-gnu/libm.so.6: version `GLIBC_2.35' not found (required by ./whitebox_tools)
Depression Filling algorithm: Whitebox Fill Depressions.
Traceback (most recent call last):
File "/media/wrf_hydro_fnl/New Volume/wrf-hydro-preprocessing_training/data/example_case/Gridded/Build_Routing_Stack.py", line 544, in
GEOGRID_STANDALONE(args.in_Geogrid,
File "/media/wrf_hydro_fnl/New Volume/wrf-hydro-preprocessing_training/data/example_case/Gridded/Build_Routing_Stack.py", line 305, in GEOGRID_STANDALONE
rootgrp2, fdir, fac, channelgrid, fill, order = wrfh.WB_functions(rootgrp2, outDEM,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/media/wrf_hydro_fnl/New Volume/wrf-hydro-preprocessing_training/data/example_case/Gridded/wrfhydro_functions.py", line 2106, in WB_functions
fill_arr, ndv = return_raster_array(fill_pits_file)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/media/wrf_hydro_fnl/New Volume/wrf-hydro-preprocessing_training/data/example_case/Gridded/wrfhydro_functions.py", line 1091, in return_raster_array
ds = gdal.Open(in_file, gdalconst.GA_ReadOnly)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/wrf_hydro_fnl/WRF_Intel/miniconda3/envs/wrf-python/lib/python3.11/site-packages/osgeo/gdal.py", line 4716, in Open
return _gdal.Open(*args)
^^^^^^^^^^^^^^^^^
RuntimeError: /media/wrf_hydro_fnl/New Volume/wrf-hydro-preprocessing_training/data/example_case/Gridded/Outputs/scratchdir/fill_depressions.tif: No such file or directory