# Roughness\_mom flag in ACCESS-OM2 doesn't do anything?

**URL:** <https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585>\
**Category:** Ocean\
**Tags:** help\
**Created:** [21 November 2025 04:30 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585 "2025-11-21T04:30:44Z")\
**Posts on this page:** 9\
**Page:** 1

<div class="post-metadata">

**Author:** ![ongqingyee](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ongqingyee/32/126_2.png) [@ongqingyee](https://forum.access-hive.org.au/u/ongqingyee)\
**Post date:** [21 November 2025 04:30 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/1 "2025-11-21T04:30:44Z")

</div>

I’ve been trying to tweak roughness coefficient in OM2 as part of a larger project. I have identified the flag I [need](https://forum.access-hive.org.au/t/forcing-access-om2-using-esm1-5-data/4349/70) and ran some experiments. Changing this flag does alter the surface stresses (tau\_x output from OM2) in the way I want, the screenshot below shows the difference between my updated roughness (roughness 4) and the old experiment (fric\_mo), apologies for the naming.

 ![Screenshot 2025-11-21 at 3.21.48 pm](https://us1.discourse-cdn.com/flex020/uploads/access1/original/2X/9/9909033272d637e0ddc5a53285d2f258b86f91c9.png)

However, this difference does _not_ translate to temperature, ocean currents etc. This is because the tau\_x (and tau\_y) diagnostics are being updated but the stresses being read into the model are not. Below are the ice-ocean stresses from the fields\_i2o\_in\_ice.nc file I outputted, the models show no differences.

 ![image](https://us1.discourse-cdn.com/flex020/uploads/access1/original/2X/f/f33625d3f120c413194d824396cef09870a0eb80.jpeg)

There are also no differences in the roughness (ice/restart000/fields\_roughness.nc)

 ![Screenshot 2025-11-21 at 3.17.31 pm](https://us1.discourse-cdn.com/flex020/uploads/access1/original/2X/0/06149b09b414e0f9b652e146fe4e2fc53db61710.png)

So apparently the roughness\_mom coefficient in ice/input\_ice\_gfdl.nml is not being used and fed into the ocean model as the surface stress field. But the coefficient _is_ being used when calculating the tau\_x and tau\_y diagnostics that are output from OM2? Would anyone know where I should be changing the roughness instead?

My config files for these two experiments, where only the roughness\_mom parameter is different, is here:  
/home/561/qo9901/access-om2/1deg\_esm1p5\_fric\_mo\_roughness\_4\_i2o/ (roughness 4)  
/home/561/qo9901/access-om2/1deg\_esm1p5\_vgrid\_fric\_mo\_i2o/ (fric\_mo)

The i2o and roughness fields are outputted too. There are now big differences between these experiments and the NRI release of OM2, mostly in the mixing schemes in ocean/input.nml, and ncar\_ocean\_fluxes = false in &surface\_flux\_nml. But it is odd to me that the surface stress being felt/read into the ocean is different from the diagnostic output.

Thank you!

---

<div class="post-metadata">

**Author:** ![lawrenceabird](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/lawrenceabird/32/3271_2.png) [@lawrenceabird](https://forum.access-hive.org.au/u/lawrenceabird)\
**Post date:** [23 November 2025 21:14 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/2 "2025-11-23T21:14:44Z")

</div>



---

<div class="post-metadata">

**Author:** ![lawrenceabird](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/lawrenceabird/32/3271_2.png) [@lawrenceabird](https://forum.access-hive.org.au/u/lawrenceabird)\
**Post date:** [23 November 2025 21:15 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/3 "2025-11-23T21:15:36Z")

</div>

Hey @ongqingyee!

Thanks for raising this query. @dougiesquire is going to dig into this and he’ll get back to you soon!

Cheers,

Lawrence

---

<div class="post-metadata">

**Author:** ![dougiesquire](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/dougiesquire/32/16_2.png) [@dougiesquire](https://forum.access-hive.org.au/u/dougiesquire)\
**Post date:** [26 November 2025 03:00 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/4 "2025-11-26T03:00:17Z")

</div>

Hi @ongqingyee. Thanks for your query and sorry it took me a few days to get to it.

I’ve just taken a quick look and I think the value of `ocean_rough_nml::roughness_mom` set in `ice/input_ice_gfdl.nml` is only used if `ocean_rough_nml::rough_scheme = 'fixed'`. You are using `ocean_rough_nml::rough_scheme='beljaars'`, which is the default in ACCESS-OM2 configurations.

If you’re happy with a fixed roughness you could switch to `ocean_rough_nml::rough_scheme='fixed'` and then modify the roughness using `ocean_rough_nml::roughness_mom`. There is also a `'charnock'` roughness scheme that may satisfy your needs. [This](https://github.com/COSIMA/cice5/blob/edcfa6f9c76ed05b63196ce4b5355fa5a8f4fe3a/drivers/auscom/ocean_rough_mod.F90#L80-L175) is where the roughness is calculated if you’re interested or want to modify the code. It’s hard to say more without being across what you are trying to achieve.

One thing I’m a bit confused about is your statement:

> [@ongqingyee](#):
>
> So apparently the roughness\_mom coefficient in ice/input\_ice\_gfdl.nml is not being used and fed into the ocean model as the surface stress field. But the coefficient _is_ being used when calculating the tau\_x and tau\_y diagnostics that are output from OM2?

I just did two runs, only changing `ocean_rough_nml::roughness_mom`, and got exactly the same `tau_x` and `tau_y` for both as expected. Am I missing something about what you have done? Have you changed anything in your MOM5 code? Could you please check this again as I cannot recreate this behaviour.

---

<div class="post-metadata">

**Author:** ![ongqingyee](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ongqingyee/32/126_2.png) [@ongqingyee](https://forum.access-hive.org.au/u/ongqingyee)\
**Post date:** [26 November 2025 23:06 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/5 "2025-11-26T23:06:43Z")

</div>

Thank you very much Dougie, so does this mean if I want to tune for roughness, the ‘charnock’ parameter is what I would want?

> [@dougiesquire](#):
>
> Have you changed anything in your MOM5 code?

I have changed many many things in the MOM5 configuration, the configs are [here](https://github.com/ongqingyee/1deg_esm1p5_vgrid_fric_mo). I wrote some changes down but I think the `ncar_ocean_fluxes = false` in ice/input\_ice\_gfdl.nml might be the one that makes the biggest difference. I did this because I was hoping to replicate the ESM1.5 ocean in OM2, and the ncar\_ocean\_fluxes is not used in ESM1.5.

> [@ongqingyee](#):
>
> There are now big differences between these experiments and the NRI release of OM2, mostly in the mixing schemes in ocean/input.nml, and ncar\_ocean\_fluxes = false in &surface\_flux\_nml.

---

<div class="post-metadata">

**Author:** ![dougiesquire](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/dougiesquire/32/16_2.png) [@dougiesquire](https://forum.access-hive.org.au/u/dougiesquire)\
**Post date:** [27 November 2025 06:41 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/6 "2025-11-27T06:41:12Z")

</div>

> [@ongqingyee](#):
>
> Thank you very much Dougie, so does this mean if I want to tune for roughness, the ‘charnock’ parameter is what I would want?

If you set `ocean_rough_nml::rough_scheme='charnock'` in `ice/input_ice_gfdl.nml` then you can [modify the roughness by changing the value of `ocean_rough_nml::rough_scheme='charnock'`](https://github.com/COSIMA/cice5/blob/edcfa6f9c76ed05b63196ce4b5355fa5a8f4fe3a/drivers/auscom/ocean_rough_mod.F90#L104-L105). The roughness is obviously calculated in a different way if you do this. I don’t know whether this suits what you are wanting to do.

> [@ongqingyee](#):
>
> I have changed many many things in the MOM5 configuration, the configs are [here](https://github.com/ongqingyee/1deg_esm1p5_vgrid_fric_mo). I wrote some changes down but I think the `ncar_ocean_fluxes = false` in ice/input\_ice\_gfdl.nml might be the one that makes the biggest difference. I did this because I was hoping to replicate the ESM1.5 ocean in OM2, and the ncar\_ocean\_fluxes is not used in ESM1.5.

Sorry, I think my previous message wasn’t very clear. From looking at the MOM5 code, I don’t see how it’s possible that the MOM5 diagnostics `tau_x` and `tau_y` could be different from the zonal wind stresses received from the coupler. Could you please check this again as I cannot recreate this behaviour?

---

<div class="post-metadata">

**Author:** ![ongqingyee](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ongqingyee/32/126_2.png) [@ongqingyee](https://forum.access-hive.org.au/u/ongqingyee)\
**Post date:** [28 November 2025 01:27 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/7 "2025-11-28T01:27:17Z")

</div>

Hi @dougiesquire, I’ve rechecked this and turns out my experiments were of different durations, which resulted in the different stresses. 🤦🏻‍♀️ glad this wasn’t a crazy issue in the end.

Thank you for checking re the roughness. I now have something I can tune.

My aim in general is to get my OM run to match ESM1.5. However, ESM1.5 does [not](https://github.com/ACCESS-NRI/access-esm1.5-configs/blob/release-preindustrial%252Bconcentrations/ice/input_ice.nml) have an ice/input\_ice\_gfdl.nml. If these parameters are not prescribed how does ESM deal with roughness?

Otherwise I would probably tune these parameters to get ESM surface stresses to match my OM run.

---

<div class="post-metadata">

**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\
**Post date:** [28 November 2025 05:09 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/8 "2025-11-28T05:09:00Z")

</div>

Hi @ongqingyee, No we wouldnt expect the ESM(1.5) to have the same input files as the OM2. But the roughness as you have seen looking at the UM7.3 is more of a “charnock” approach. But it is variable, wind dependent, not fixed like you had been trying to set. So I would try the “charnock” setting in OM2 next. It might not work but it might get closer.

---

<div class="post-metadata">

**Author:** ![dougiesquire](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/dougiesquire/32/16_2.png) [@dougiesquire](https://forum.access-hive.org.au/u/dougiesquire)\
**Post date:** [21 January 2026 22:46 UTC](https://forum.access-hive.org.au/t/roughness-mom-flag-in-access-om2-doesnt-do-anything/5585/9 "2026-01-21T22:46:18Z")

</div>


