<p>
        Hi everyone,
</p>
<p>
        <br>
</p>
<p>
        I am now doing ocean simulations on MITgcm with shelfice package, non-hydrostatic settings and flux-form momentum equations. Now I have some problems with the diagnostic output:
</p>
<p>
        <br>
</p>
<p>
        (1) For the advection term Um/Vm_Advec, the description is "<span style="font-family:Menlo;font-size:11px;">U/V momentum tendency from Advection terms", <span style="font-family:Arial;"><span style="font-size:16px;">but I checked the pkg/mom_fluxform/mom_fluxform.F and find that it may be wrong:</span></span></span>
</p>
<p>
        Line 950, the diagnostic output is filled by gV: 
        </p><p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                        CALL DIAGNOSTICS_FILL(gV(1-OLx,1-OLy,k,bi,bj),
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                     &                               'Vm_Advec',k,1,2,bi,bj,myThid)
        </p>

<p>
        <span style="white-space:normal;">Before that, gV is added by advection terms (line 689), metric terms (line 867 and 876) and Coriolis terms (line 907), so it seems that the final output</span> should be advection terms + metric terms + Coriolis terms.
</p>
<p>
        <br>
</p>
<p>
        (2) I also want to know how I can get the output of the non-hydrostatic pressure field (phi_NH) in the model.
</p>
<p>
        <br>
</p>
<p>
        (3) For the Ice shelf freshwater flux (<span style="font-family:Menlo;font-size:11px;">SHIfwFlx<span style="font-family:Arial;"><span style="font-size:16px;">), I get the output of all zeros, which is quite weird because there is salinity change so the ice should be melting/freezing. I'm not sure if this is caused by some settings for the shelfice package. My data.shelfice settings are as follows:</span></span></span>
</p>
<p>
        <span style="font-family:Menlo;font-size:11px;"><span style="font-family:Arial;"><span style="font-size:16px;">
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 &SHELFICE_PARM01
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEboundaryLayer = .FALSE.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEtopoFile='ice.topo',
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEloadAnomalyFile = 'p.anomaly',
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 useISOMIPTD = .FALSE.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEconserve = .TRUE.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEuseGammaFrict = .FALSE.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 no_slip_shelfice = .false.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEwriteState = .TRUE.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEthetaSurface = -81.2,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 rhoShelfIce = 926.,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEsaltTransCoeff = 4.E-8,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 SHELFICEheatTransCoeff = 7.E-7,
        </p>
        <p style="margin-bottom:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;">
                 &
        </p>
        <p>
                <span style="font-family:Menlo;font-size:11px;"><span style="font-family:Arial;"><span style="font-size:16px;"><br>
</span></span></span>
        </p>
Thanks,<br>
</span></span></span>
</p>
<p>
        Yaoxuan
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>