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.
-
mhm_6basin_ElbeCoarseROUTING_L11_1by32_shortRUNtime.nml
-
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