3a_forcing/3_create_shapefile "create_ERA5_shapefile"

[bartvanosnabrugge]

• [ ] run test_bow_at_banff.py file to make sure current version is working
• [ ] copy/ change and adjust code to work as python function
• [ ] make sure to add docstring in the sphinx format (see for example util/util.py - read_summa_workflow_control_file or data_specific_processing/merit.py)
• [ ] Add to the description parts of the readme.rst file that come with the summa-cwarhm script that are relevant to understand what the function does precisely.
• [ ] point test_bow_at_banff.py to use the function instead of the wrapper
• [ ] test run test_bow_at_banff.py again

[ShervanGharari]

This can be done using easymore package in various forms:

1) if the user call the easymore nc_remapper the source shapefile is created in the temporary files.

2) calling the easymore functions based on the first example of the easymore: https://github.com/ShervanGharari/EASYMORE/blob/main/examples/Chapter1_E1.ipynb

instead of calling the esmr.nc_remapper() user can instead do the following:

esmr.NetCDF_SHP_lat_lon()
if hasattr(esmr, 'lat_expanded') and hasattr(esmr, 'lon_expanded'):
    esmr.lat_lon_SHP(esmr.lat_expanded, esmr.lon_expanded,\
        esmr.temp_dir+esmr.case_name+'_source_shapefile1.shp')
else:
    esmr.lat_lon_SHP(esmr.lat, self.lon,\
        esmr.temp_dir+esmr.case_name+'_source_shapefile1.shp')

note that the test.shp does not have a projection and that should be added

Note, the nc shapefile depending on the whether or not the nc file is regular or rotated the first and last row and first and last column might be missing in the shapefile. Therefore it is suggested to have a buffer for the domain of modelling.