Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Remove unused inputs #13

Closed
aidanheerdegen opened this issue May 17, 2024 · 15 comments
Closed

Remove unused inputs #13

aidanheerdegen opened this issue May 17, 2024 · 15 comments

Comments

@aidanheerdegen
Copy link
Member

It appears the geothermal heating isn't turned on

https://github.com/ACCESS-NRI/access-esm1.5-configs/blob/pre-industrial/ocean/input.nml#L53

Should we remove the geothermal_heating.nc input file?

@aidanheerdegen
Copy link
Member Author

aidanheerdegen commented May 17, 2024

Also basin_mask.nc isn't a model input and was deleted from the OM2 1 degree configurations

COSIMA/access-om2#205

Edit: to be explicit, this is an input file, but it isn't used

https://github.com/mom-ocean/MOM5/blob/d7ba13a3f364ce130b6ad0ba813f01832cada7a2/src/mom5/ocean_core/ocean_thickness.F90#L567-L571

https://github.com/ACCESS-NRI/access-esm1.5-configs/blob/pre-industrial/ocean/input.nml#L469

@aidanheerdegen aidanheerdegen changed the title Remove unused geothermal heating input Remove unused inputs May 17, 2024
@aidanheerdegen
Copy link
Member Author

CABLE-AUX-1.4/core/biogeochem/pftlookup_csiro_v16_17tiles.csv is commented out for the historical configuration:

https://github.com/ACCESS-NRI/access-esm1.5-configs/blob/pre-industrial/atmosphere/cable.nml#L56

@aidanheerdegen
Copy link
Member Author

pftlookup_csiro_v16_17tiles_spinup.csv isn't referenced at all. Neither is poolcnpInTumbarumba.csv

@aidanheerdegen
Copy link
Member Author

cf_name_table.txt seems superfluous. As does SOURCES which just contains a bunch of paths.

@aidanheerdegen
Copy link
Member Author

We need some input from @MartinDix or Tilo on the CABLE inputs (or @ccarouge?).

@dougiesquire can you confirm I'm correct about the geothermal heating and basin mask

@dougiesquire
Copy link
Contributor

@dougiesquire can you confirm I'm correct about the #13 (comment) and #13 (comment)

Yup, looks right to me (for the pre-industrial config at least)

@ccarouge
Copy link
Member

@aidanheerdegen What do you mean by removing the input files? I guess I don't know what is happening around the ESM1.5 preparations so it's not clear what you mean.

For the pftlookup file, the namelist file points to a different file, CABLE only uses 1 of these files per configuration so it's ok to remove other pftlookup files.

For the pool file, I'm surprised there isn't a pool file specified in the CABLE namelist but that's reaching parts of CABLE I don't know much about. That said, the file you list poolcnpInTumbarumba.csv seems to be for a single point simulation (I don't know where this file is so I can't check) that isn't suitable for a global simulation so should be fine to remove. We'll see if @MartinDix knows better than me.

@aidanheerdegen
Copy link
Member Author

aidanheerdegen commented Jun 18, 2024

What do you mean by removing the input files? I guess I don't know what is happening around the ESM1.5 preparations so it's not clear what you mean.

There appears to be a number of input files that just aren't used by the model configuration. So we would like to not include them in the input directory (remove them)

@dougiesquire
Copy link
Contributor

I was surprised/curious about all the ocean/biogeochemistry inputs so I had a deeper look:

  • cfc_auscom.nc: I can’t find any evidence that this is used by the ocean (I haven't checked other components)
  • co2_obs.nc: I presume this is anthropogenic atmospheric co2 concentration. But use_access_co2 = .true. in the field_table meaning the co2 concentration from the coupler is used. Also, aco2_file = INPUT/co2_b35_2D.nc in the field_table and this file is not provided as an input. To use the co2_obs.nc input for anthropogenic CO2 concentration, this should be aco2_file = INPUT/co2_obs.nc. Currently users would get an error if they tried to set use_access_co2 = .false..
  • ocmip2_abiotic_c14_atm_hist_om3_bc-1-9999.nc: I can’t find any evidence that this is used by the ocean (I haven't checked other components)
  • ocmip2_fice_monthly_om1p5_bc.nc: this is not used since ice_file4gasx = .false. in the field_table
  • ocmip2_siple_co2_atm_am2_bc-1-9999.nc: I can’t find any evidence that this is used by the ocean (I haven't checked other components)
  • ocmip2_xkw_monthly_om1p5_bc.nc: this is not used since gasx_from_file = .false. in the field_table

So...

If we only want to keep inputs that are used by the configuration without modification:

  • Remove cfc_auscom.nc
  • Remove co2_obs.nc
  • Remove ocmip2_abiotic_c14_atm_hist_om3_bc-1-9999.nc
  • Remove ocmip2_fice_monthly_om1p5_bc.nc
  • Remove ocmip2_siple_co2_atm_am2_bc-1-9999.nc
  • Remove ocmip2_xkw_monthly_om1p5_bc.nc

or, if we want to support users changing the WOMBAT field_table flags use_access_co2, gasx_from_file and ice_file4gasx:

(Note, someone who knows about these inputs should check all of this)

@blimlim
Copy link

blimlim commented Jul 1, 2024

I've had a go at removing the extra inputs from an esm1.5 configuration – most of the files could be removed without any issues, apart from:

  • ocmip2_xkw_monthly_om1p5_bc.nc
  • ocmip2_fice_monthly_om1p5_bc.nc

The model crashes when either of these are removed even though the related namelist settings are .false.. Changing both files' values to zero doesn't seem impact the simulation at all, so it looks like the model is just trying to read in their values without using them.

I'm wondering what a good way to categorise these files in the restructured input directories be would be. E.g. do you think it would it be more confusing to classify them as unused but still require them in the config.yaml file, or to instead include them with the other required files despite them being ignored in the config's namelist settings? Or if there are any other ideas that could be good too!

@aidanheerdegen
Copy link
Member Author

The model crashes when either of these are removed even though the related namelist settings are .false..

The version of csiro_bgc.F90 in the ACCESS-ESM1.5 fork doesn't do the check for gas_from_file before trying to open the file.

https://github.com/ACCESS-NRI/MOM5/blob/access-esm1.5/src/mom5/ocean_csiro_bgc/csiro_bgc.F90#L2027-L2029

Arguably then this should be back-ported to the ESM1.5 branch.

If it is just this commit then it may be straightforward to cherry-pick it

ACCESS-NRI/MOM5@3f7b6de

But I'll leave that up to @dougiesquire to make that call.

@dougiesquire
Copy link
Contributor

If it is just this commit then it may be straightforward to cherry-pick it

ACCESS-NRI/MOM5@3f7b6de

But I'll leave that up to @dougiesquire to make that call.

I agree. PR here

@dougiesquire
Copy link
Contributor

I agree. PR here

@blimlim, this PR has now been merged into the access-esm1.5 branch. When you get a chance, could you confirm that with this change you can now run without ocmip2_xkw_monthly_om1p5_bc.nc and ocmip2_fice_monthly_om1p5_bc.nc?

@blimlim
Copy link

blimlim commented Jul 4, 2024

When you get a chance, could you confirm that with this change you can now run without ocmip2_xkw_monthly_om1p5_bc.nc and ocmip2_fice_monthly_om1p5_bc.nc?

Can do! I've added an issue to the ESM1.5 repository to incorporate these changes into a new MOM5 executable, and will test this out once its is ready. @penguian mentioned that it would be useful to set up automated issues in https://github.com/ACCESS-NRI/ACCESS-ESM1.5/ that are automatically created whenever there are new commits to the esm1.5 components.

@aidanheerdegen
Copy link
Member Author

Closing as I think this is completed, but re-open if there is anything outstanding.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants