[MITgcm-support] funky ice dynamics in doubly periodic domain
Jean Mensa
jean.mensa at yale.edu
Sun Mar 6 13:23:24 EST 2016
OK, I found the problem. In EXF, useRelativeWind has to be set to TRUE.
This keeps the ice from flying away...
thank you for the help!
j
On Tue, Feb 23, 2016 at 1:07 PM, Dimitris Menemenlis <dmenemenlis at gmail.com>
wrote:
> looking at available_diagnostics.log, the following surface stress
> variables are available:
>
> 79 |oceTAUX | 1 | 80 |UU U1|N/m^2 |zonal surface
> wind stress, >0 increases uVel
> 80 |oceTAUY | 1 | 79 |VV U1|N/m^2 |meridional surf.
> wind stress, >0 increases vVel
> 179 |EXFtaux | 1 | |UM U1|N/m^2 |zonal surface
> wind stress, >0 increases uVel
> 180 |EXFtauy | 1 | |VM U1|N/m^2 |meridional
> surface wind stress, >0 increases vVel
> 227 |SIfu | 1 | 228 |UU U1|N/m^2 |SEAICE zonal
> surface wind stress, >0 increases uVel
> 228 |SIfv | 1 | 227 |VV U1|N/m^2 |SEAICE merid.
> surface wind stress, >0 increases vVel
>
> and from seaice_diagnostics_init.F
>
> C pkg/diagnostics SIfu and oceTAUX, dumpfreq FU, and tavefreq FUtave
> C are identical but they differ from pkg/diagnostics EXFtaux, which
> C is stress before impact of ice. Also when using exf bulk
> C formulae, EXFtaux is defined on tracer rather than uvel points.
>
> C pkg/diagnostics SIfv and oceTAUY, dumpfreq FV, and tavefreq FVtave
> C are identical but they differ from pkg/diagnostics EXFtauy, which
> C is stress before impact of ice. Also when using exf bulk
> C formulae, EXFtauy is defined on tracer rather than vvel points.
>
> meaning that if siAREA=1, all of the above are identical and denote stress
> from ice to ocean,
> except that EXFtaux and EXFtauy are defined at the tracer point instead of
> the velocity
> points on the SouthWest C-grid.
>
> we don’t appear to have a diagnostic for air-ice stress, but that can be
> inferred
> from 10-m wind SEAICE_drag, and SEAICE_rhoAir
>
> since you have no shear in the surface forcing, I am assuming the internal
> stress will be negligible,
> so if you add the components of stress at top and bottom of the ice, you
> should be able to estimate
> expected acceleration of the sea ice, and that number should be much
> smaller than my 3 m/s/day
> back-of-the-envelope estimate from earlier today that neglected bottom
> friction.
>
> > On Feb 23, 2016, at 9:02 AM, Jean Mensa <jean.mensa at yale.edu> wrote:
> >
> > ok, for oceTAUX/oceTAUY I will have to rerun but I have KE series at
> surface. Velocities decay with depth. Ice velocity is not anomalously large
> and features seem reasonable.
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
>
> https://urldefense.proofpoint.com/v2/url?u=http-3A__mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=AwIGaQ&c=-dg2m7zWuuDZ0MUcV7Sdqw&r=5oFkuOCa8D2MEaB0wBLw8HmmtUZ0uoNntNFgQ9b_kVE&m=9kFa5U0NT4xR7lY6dvKOuQUO3n0I5ButKe8wIYC4D5U&s=QCopWjF-yiA7cZwNWDln90DaD19bTt9sb3vv3yFf0Lo&e=
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160306/85b0559b/attachment.htm>
More information about the MITgcm-support
mailing list