ACCESS-coupled N48 for deep paleo

well, I have your old at Uni Sydney, but not your new BoM!

I am now starting to change the ancillary files for the coupled ACCESS-N48 run (u-cy630). I need to know:

(1) how do change the ancillary file names/path for the UM and MOM?

Currently I can’t even figure out what the path for the MOM ancillary files are. Some help here would be great.

@Aidan @atteggiani

Hi Dietmar,

The namelists are generated in roses/u-cy630/app/coupled/rose-app.conf.
Here there is a list with all the input files in the format:

[file:<reference file for the model>]
mode=<type of file> (usually symlink or omitted)
source=<the actual filepath>

For example, in your suite, the current grid_spec.nc is listed as:

[file:$OCN_RUNDIR/INPUT/grid_spec.nc]
source=~access/access-cm2/input_b/mom4/grid_spec.auscom.20110618.nc

If you want to change that, you have to change the source field.

Repeat the step for any input/ancillary file you want to change.

Cheers
Davide

1 Like

Dear ACCESS-hive community,
here is a short update on what I have been working on the last two month:

(*) Change land-sea distribution in coupled GFDL: following the previous advice I started with running the coupled GFDL model from @dkhutch. This works now with the help from @dkhutch. I can introduce new land and remove land. I am currently running a configuration that has only one tropical ocean basin. It has been running for 80yrs so far.

(*) Change land-sea distribution in ocean only GFDL model (mom): I followed the instructions from @rmholmes and it is fairly straight forward.

(*) ACCESS-N48 coupled model: I have done about 100 experiments with changing land and oceans, and worked on this for a few 100hrs. To make changes in the land-sea distributions, topographies and bathymetries I change the following ancillary files:

grid_spec.nc
masks.nc
qrparm.landfrac
qrparm.mask
qrparm.orog
rmp_momt_to_um_t_CONSERV_FRACAREA.nc
rmp_momt_to_um_t_PATCH.nc
rmp_momt_to_um_u_PATCH.nc
rmp_momt_to_um_v_PATCH.nc
rmp_um_t_to_momt_CONSERV_FRACAREA.nc
rmp_um_t_to_momt_PATCH.nc
rmp_um_u_to_momt_PATCH.nc
rmp_um_v_to_momt_PATCH.nc

This works sometimes: I can do the following:

() Close Drake passsage → OK (runs for 5yrs)
(
) introduce big island in tropical Pacific → OK (runs for 5yrs)
() remove australia → OK (runs for 5yrs)
(
) remove s-america → OK (runs for 2yrs)

Now the runs that do not work:

(*) remove most of the land in tropics:
data[ 90:190, 1:170] = 4000.0
data[ 90:190,220:360] = 4000.0

→ crash at start:
→ all kind of different warnings /last one:

? Warning code: -1
? Warning from routine: OASIS_INITA2O
? Warning message: Coupling field not enabled - 3D CO2 STASH code:0,252
? Warning from processor: 0
? Warning number: 14

→ Floating point exception

I did a few more tests:

() data[ 90:190, 5:330] = 4000.0 → crash at start
(
) data[ 90:190, 30:300] = 4000.0 → crash at start

I noticed that when I remove most of the land in tropics the land fraction is not zero at longitude at which the ocean map starts, see figures below. So there is clearly something wrong with the land fraction calculation. However, this I do not see in the other tests:

() data[ 90:190, 5:330] = 4000.0 → crash at start
(
) data[ 90:190, 30:300] = 4000.0 → crash at start

So, I am clueless on how to continue. Any ideas?

For now I have moved this project into my “lost cause” folder, as I have no more time to spend on this. May start to think about it again in October, after teaching.

best regards
Dietmar

@atteggiani @Scott @abhik @holger @Aidan

3 Likes

Thanks for the update Dietmar. I can appreciate that you don’t have time right now to keep working on it, but if you’d like to revisit it in October, please send us the work directories of the problematic runs and one of us (if we’re feeling brave) might poke around at it.