search

ICTP / RegCM

ICTP Regional Climate Model
https://www.ictp.it/esp/about
Other
54 stars 42 forks source link

floating divide by zero error #18

Closed DesmetQuentin closed 2 years ago

DesmetQuentin commented 2 years ago

Hello,

I'm the only one in my team to have this problem, so I guess I did something wrong... I get a floating divide by zero error at line 848 of _modtendendy.F90 for a restart run with idynamic = 2. I have checked the array by printing them in debug files, and the zeros appear on external points. The error disappears when I move the exchange call of sfs%psdota at line 800 in _modinit.F90 to between line 805 and 806. Is there some parameterization I missed, which would prevent doing this?

Cheers, Quentin

graziano-giuliani commented 2 years ago

Thanks Quentin, You are correct.

DesmetQuentin commented 2 years ago

I'm glad to help !