To make the PCRaster calculation in NGMS the following xml files have to be created:
New xml files:
• Module configuration files:
WMW_interpolate_grid_to_branch_grid
WMW_PCR_generate_branch_file
WMW_QTRIB_ref(Historic)
• workflow file:
WMW_configuration_only
Modules directory:
• Fixed:
Ldd.map
Strdir.map
Scalar.map
• fromNGMS:
branch0000.000
• toNGMS:
brances_grid.map
cross_points.map
stream_directions.map
In order to create the above mentioned files it is necessary to follow the instructions below:
1. Import the branch locations into PCRaster. This is done with the PCRaster general adaptor export, which produces branch0000.000 file.
1.1. Copy the X and Y locations from the Branches.xml and make an new irregular grid (XXX_branch_grid) in Grids.xml
1.2. Make a new location ID for XXX_branch_grid.
1.3. Make XXX_interpolate_grid_to_branch_grid file. This file converts the locationid from XXX_merged to the XXX_branch_grid. And back to the XXX_merged grid, making use of inverse distance interpolation technique. Adapt the search radius to grid size. (= smaller than half a grid cell)
1.4. Make general adaptor run file: XXX_PCR_generate_branch_file. This file exports the branch0000.000 file.
2. Make stream flow direction map.
2.1. Export the stream flow direction map with the Modflow general adaptor
2.2. Rename stream flow direction asci file into "strdir.asc"
3. Make initial scalar and ldd map
3.1. use mapattr command to produce scalar.map mapattr scalar.map
number of columns =360
number of rows=240
data type = scalar
projection= y increases from bottom to top
x upperleft corner=400000
y upperleft corner=170000
celllength = 250
3.2. use mapattr command to produce ldd.map mapattr ldd.map
number of colums =360
number of rows=240
data type = ldd
projection= y increases from bottom to top
x upperleft corner=400000
y upperleft corner=170000
celllength = 250
4. Make a new workflow and adapt the WorkflowDescriptorDescriptors and the ModuleinstanceDescriptors file
5. Run workflow XXX_PCR_generate_branch_file. This run produces the following files: brances_grid.map, cross_points.map and stream_directions.map
6. make XXX_QTRIB_ref(Historic) transformation file
7. adapt interpolation and compare files
8. adapt Moduledescriptor, ModuleinstanceDescriptors and the Parameter file