Step 8. Derivation tributary flows (PCRtransformation)

In Step 7 a ModuleConfigFile is produced which uses a PCRaster transformation. This PCRaster transformation is used for the following 3 steps:

1. PCRaster determines the discharge of the tributaries; the flow of the branch rivers is substracted from the total surface water flow.
   only_Q_tributaries = Qflow * branches_grid;
2. Sums the discharges of all upstream river cells
   sum_Q_tributaries = upstream(stream_directions, only_Q_tributaries);
3. Make the map with discharges of all upstream river cells at the crossing points of rivers.
   Output = sum_Q_tributaries * cross_points;

This PCRaster tranformation needs input to be able to run and therefore the following xml files have to be created:
New xml files :

• Module configuration files:
  • XXX ¹_interpolate_grid_to_branch_grid
  • XXX ¹_PCR_generate_branch_file
  • XXX ¹_QTRIB_ref(Historic)
• workflow file:
  • XXX ¹_configuration_only
    ¹ = model abbreviation

Besides the xml files the following files in the Modules directory should be generated:

• 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.

• Copy the X and Y locations from the Branches.xml and make an new irregular grid (XXX_branch_grid) in Grids.xml
• Make a new location ID for XXX_branch_grid.
• 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)
• Make general adaptor run file: XXX_PCR_generate_branch_file. This file exports the branch0000.000 file.

1. Make stream flow direction map.

• Export the stream flow direction map with the Modflow general adaptor
• Rename stream flow direction asci file into "strdir.asc"

1. 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
2. Make a new workflow and adapt the WorkflowDescriptorDescriptors and the ModuleinstanceDescriptors file
3. Run workflow XXX_PCR_generate_branch_file. This run produces the following files: brances_grid.map, cross_points.map and stream_directions.map
4. make XXX_QTRIB_ref(Historic) transformation file
5. adapt interpolation and compare files
6. adapt Moduledescriptor and ModuleinstanceDescriptors if still applicable