[MITgcm-support] FW: Diagnose momentum equation
Jean-Michel Campin
jmc at ocean.mit.edu
Sat Sep 26 11:19:28 EDT 2015
Hi Lianzhan,
It's likely that if Temp & Salt remains uniform, U/Vm_dPHdx will be zero
(no hydrostatic pressure anomaly horizontal gradient).
Not very clear why U/Vm_Cori would be zero everywhere except for surface level,
unless the velocity (uVel,vVel) is also zero everywhere below k=1.
What type of forcing do you have in this set-up ?
If you could send (in attachement) the 2 parameter files "data" and
"data.pkg", this would help to figure out what is happening.
Cheers,
Jean-Michel
On Sat, Sep 26, 2015 at 01:27:27PM +0000, lianzhan wrote:
> Hello everyone
>
> I want to diagnose the effect of all items in momentum equation (as below for zonal).
> dU/dt+UdU/dx+VdU/dy+WdU/dz=fv-1/rho*dP/dx+Fx+taux/rho
> Fx is dissipative stress, taux is wind stress. Is there any corresponding variable in pkg/diagnostics ? I have tried "U/Vm_Advec" for nonlinear term, but this variable is zero expect for the surface layer. Its value is as same as "U/Vm_Cori", the value of "U/Vm_Cori" is also zero below surface. And "U/Vm_dPHdx"(pressure gradient force ? ) is all zero.
> The settings of my exp is that there is an island located in an ideal rectangular bay, near which there is a sill. The model has 32 vertical layers and with free surface. The Temp/Salt is uniform.
> Do you have any recommendations or suggestions?
> Thanks in advance for your helpCheers,Lianzhan
>
>
> xiaober at outlook.com
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list