When I run the model with the GRIDC=RECT the model works fine but the it crashes when I use GRIDC=TRAP. Is there a problem with this algorithm? I've checked the code the problems comes from a NaN value for the parameter "Z". In file "gate-spill-pipe.f90" / subroutine "GRID_AREA2" line 630.
One fix