When I use this module, the software crashed if np > 4 for mpi. Then I try to find the reason. Finally, I located the error in Line 5794 call aexchange(nc,myid,nprocs,taub). After some test, I find that there is something wrong of taub. Fortunately, only one place change it.
In function Advance_Sed(DTin,Tin,taub_in), taub = taub_in*rho_water is set in line 358 and line 361. But the size of input value taub_in is m, while the size of taub is 0:MT.
It is simple to fix it. modify it to taub(1:m) = taub_in*rho_water.
I hope this report is helpful.
When I use this module, the software crashed if np > 4 for mpi. Then I try to find the reason. Finally, I located the error in Line 5794
call aexchange(nc,myid,nprocs,taub). After some test, I find that there is something wrong oftaub. Fortunately, only one place change it.In function
Advance_Sed(DTin,Tin,taub_in),taub = taub_in*rho_wateris set in line 358 and line 361. But the size of input valuetaub_inis m, while the size oftaubis 0:MT.It is simple to fix it. modify it to
taub(1:m) = taub_in*rho_water.I hope this report is helpful.