[MITgcm-support] some questions about the terms in the sea ice momentum equations

[Dimitris Menemenlis]

Before you do that, enable useDiagnostics, run the model for one time step in your particular configuration, dump "available_diagnostics.log" and look carefully at the list of terms that are made available from pkg/diagnostics. In the case of pkg/seaice they would tend to be preceded by "SI". What you need may already be listed there and a lot more convenient and safer to use. As a starting point you can also look at available_diagnostics.log from on of the verification experiments that is in a configuration similar to yours, e.g., MITgcm/verification/lab_sea/run/available_diagnostics.log

On Jun 30, 2014, at 5:54 PM, yang.zhang <yang.zhang at mun.ca> wrote:

> Hi, Wei
> 
>      I guess, you can look into the code where the output is written, and find the variable in the code standing for ocean-ice stress (maybe just linearly dependent on the difference between ice and water velocity), Coriolis force (maybe u just need the ice velocity to calculated this term), internal stress tensor , to write them out using fortran write statement (if there is no default output option for these terms u are interested, It is quite common that not every variable could be default output field  ). Good luck..
> 
> Zhang Yang
> 
> 
> On Sun, Jun 29, 2014 at 9:03 AM, wisteria0310 at gmail.com <wisteria0310 at gmail.com> wrote:
> Hello everyone,
> 
>  Recently I am studying the impacts of cyclones on the Arctic sea ice motion using the MITgcm model. As we know, the momentum equation of the sea ice motion is as follows:
> 
> <clip_image002(06(06-29-19-31-15).png>
> 
> According to the momentum equation, there are many terms that contribute to the values of sea ice velocity. Now I want to include them as the output fields of the model to analyze the variability of the sea ice motion in depth. However, I have only found the wind-ice stress (SIfu, SIfv) and the gradient (or tilt) of the sea surface height (phiSurfX, phiSurfY) in the model. What should I do for the ocean-ice stress, the Coriolis force and the divergence of the internal ice stress tensor? I am not sure what kind of forms they take in the model and how could I output them all?
> 
>   
> Best wishes,
> 
> Jennifer
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140630/c49d86fe/attachment-0001.htm>