You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 5 Next »

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 1_interpolate_grid_to_branch_grid
    • XXX 1 _PCR_generate_branch_file
    • XXX 1 _QTRIB_ref(Historic)
  • workflow file:
    • XXX 1 _configuration_only
      1 = model abbreviation

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

  • No labels

Disclaimer