Resolve topographic meltpond bug in iceth_pnd.F90 on branch_4.2
Context
-
Branches impacted: branch_4.2
Analysis
Bug in topographic meltpond code in iceth_pnd.F90 on branch_4.2.
Fixed by David Schroeder in the main
by !476 (merged), but does not appear to have been fixed in branch_4.2.
Fix
- zdHui = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_i(ji,jj,jl) )**2 ) &
- - v_il(ji,jj,jl) / a_i(ji,jj,jl)
+ zdHui = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_ip(ji,jj,jl) )**2 ) &
+ - v_il(ji,jj,jl) / a_ip(ji,jj,jl)
See also: devs/ukmo/gosi/gosi@93f6eda1
Edited by Daley Calvert