Some bugs are corrected in routing.f90 and testrouting.f90. The calculation of the Pfafstetter codes in river_desc.nc. To be seen later if this is needed.