[MITgcm-support] Friction velocity ustar in S/R exf_wind.F
Martin Losch
Martin.Losch at awi.de
Mon Mar 5 04:18:44 EST 2018
Hi Stanislav,
I am no expert on this, but did you see lines 156/157:
wStress(i,j,bi,bj) = SQRT(usSq)
c ustar = SQRT(usSq/atmrho)
so ustar should be sqrt(usSq/atmrho). for algorithmic reasons wstress=sqrt(usSq) is computed first and then divided by sqrt(atmrho) later in lin 179. To this makes sense, although it’s not very easy to read.
Martin
> On 4. Mar 2018, at 13:58, Stanislav Martyanov <martyanov.sd at gmail.com> wrote:
>
> Hello!
>
> I think it should be
>
> ustar = SQRT (wStress(i,j,bi,bj) / atmrho)
>
> rather than
>
> ustar = wStress(i,j,bi,bj)*recip_sqrtRhoA
>
> in exf package in S/R exf_winf.F (line 179)?
>
> In the current version wStress is not affected by the function SQRT,
> but should be, in my opinion.
>
> Best regards,
> Stanislav
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list