Running MOM6 for an idealized 2-layer solution

[xliu-ocean]

Hi everyone,

I am trying to use MOM6 to do some 2-layer idealized-basin tests. I revised the original MOM6-example "ocean_only" case to reproduce the classical Stommel solution. The parameters I used were given by R. Huang in his note (attached in the PDF file), with a sine-shape zonal wind, a reduced gravity of 0.015 and an initially 300m thermocline.

However, the results of my MOM6 test show a strong outcrop of the interface of the two layers near the northern boundary (shown in the figure). I checked the parameters I used in the model and cannot figure out the problem. The MOM_parameter_doc.all is also attached.

The wind stress was given by revising the code in user_surface_forcing.F90

forces%taux(I,j) = -0.15*G%mask2dCu(I,j) * & cos(PI/30.0*(G%geoLatCu(I,j)-20.0))* & US%kg_m3_to_R*US%m_s_to_L_T**2*US%L_to_Z

Does anyone have any idea about my setting of the model?

Thanks,
Xiaohui Liu

07ReduG_1.pdf

MOM_parameter_doc.all.txt

[xliu-ocean]

The figure of the layer interface at the 10 years might not so clear. I converted it to PDF type.
layer_interface_at10years.pdf

[marshallward]

Were you able to sort this out? @Hallberg-NOAA was hoping to respond to this when he had a chance.

[xliu-ocean]

I have done a series of tests and found that the bottom drag should be a critical parameter for the model results. However, the value of the bottom drag (kappa) was not given in Huang's note. I am now trying to figure out how to deal with this parameter. I closed this issue considering the complexity of this problem.
Another problem is that the bottom drag in a 1.5-layer model is actually at the bottom of the upper layer, while in MOM6, the bottom drag is at the lower layer.

Thanks so much for your reply @marshallward and @Hallberg-NOAA . I reopened this issue in case you are willing to respond.

Xiaohui