[MITgcm-support] Problems about the diagnostic output of MITgcm

[曾耀萱]

Hi everyone,




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:




(1) For the advection term Um/Vm_Advec, the description is "U/V momentum tendency from Advection terms", but I checked the pkg/mom_fluxform/mom_fluxform.F and find that it may be wrong:

Line 950, the diagnostic output is filled by gV: 

        CALL DIAGNOSTICS_FILL(gV(1-OLx,1-OLy,k,bi,bj),

     &                               'Vm_Advec',k,1,2,bi,bj,myThid)

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 should be advection terms + metric terms + Coriolis terms.




(2) I also want to know how I can get the output of the non-hydrostatic pressure field (phi_NH) in the model.




(3) For the Ice shelf freshwater flux (SHIfwFlx), 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:

 &SHELFICE_PARM01

 SHELFICEboundaryLayer = .FALSE.,

 SHELFICEtopoFile='ice.topo',

 SHELFICEloadAnomalyFile = 'p.anomaly',

 useISOMIPTD = .FALSE.,

 SHELFICEconserve = .TRUE.,

 SHELFICEuseGammaFrict = .FALSE.,

 no_slip_shelfice = .false.,

 SHELFICEwriteState = .TRUE.,

 SHELFICEthetaSurface = -81.2,

 rhoShelfIce = 926.,

 SHELFICEsaltTransCoeff = 4.E-8,

 SHELFICEheatTransCoeff = 7.E-7,

 &




Thanks,


Yaoxuan














-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190906/d7835217/attachment.html>