bug in Rothrock ice strength
When using Rothrock scheme for calculating the ice strength, there is a bug due to the call to rdgrft_prep. This subroutine calculates things necessary for the strength but also the opening and closing rates. These rates need rheology outputs and then they must not be calculated here for the strength.