mRM initilization hangs for many hours at "call L11_stream_features" of mo_mrm_init.f90 file (L11: 1/64) whereas 1/32 L11 is quick

Hello,

mRM initilization of Elbe L1 (scale of 1/64) and L11 (1/32) takes several minutes with cygwin in my desktop but mRM hangs when both L1 and L11 are 1/64.

Considering the increase in cells size from 1/32 to 1/64 I would expect an increase in initialization time but here we have hours (more than 5 hours) of waiting.

I did a small search in the code and added some "write" to see where hanging happens.

Hanging is happening at the line of "call L11_stream_features" of mo_mrm_init.f90 file. There is a counter in “iD0” scanning at L0.

For your testing I prepared two mhm.nml files for 1/32 L11 and 1/64 L11.

1. mhm_6basin_ElbeCoarseROUTING_L11_1by32_shortRUNtime.nml

2. mhm_6EU_1by64_degrees_HANGS.nml

Download link: web.itu.edu.tr/demirelmc/mHM_EUscale_NEW_ONLYelbe.rar

I would be very grateful if @ottor and @thober could have time to see the issue.

The issue might be a similar case to #10 (closed)

do iBasin = 1, nBasins
  if (.not. read_restart) then
    call L11_flow_direction(iBasin)
    write(*,*) 'Position mrm init after call L11_flow_direction(iBasin)'

    call L11_flow_accumulation(iBasin)
    write(*,*) 'Position mrm init after call L11_flow_accumulation(iBasin)'

    call L11_set_network_topology(iBasin)
    write(*,*) 'Position mrm init after call L11_set_network_topology(iBasin)'

    call L11_routing_order(iBasin)
    write(*,*) 'Position mrm init after call L11_routing_order(iBasin)'
    call L11_link_location(iBasin)
    write(*,*) 'Position mrm init after call L11_link_location(iBasin)'
    call L11_set_drain_outlet_gauges(iBasin)
    write(*,*) 'Position mrm init after call L11_set_drain_outlet_gauges(iBasin)'
    ! stream characteristics
    call L11_stream_features(iBasin)
    write(*,*) 'Position mrm init after call L11_stream_features(iBasin)'
  end if
end do