-
Notifications
You must be signed in to change notification settings - Fork 24
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
Highres jra bgc #34
Highres jra bgc #34
Conversation
Pointing to temp versions of BGC initial condition files and some forcings, and temporarily turning off other forcings just to get model running
Updated GIAF_JRA_HR and GIAFECO_JRA_HR to use CICE%CICE4 (requires cice5_20200430 or later)
Incorporates a few namelist changes from Alper's run that aren't on CICE master, and also uses xmlchange to turn off ideal age tracers (adding -trage 0 to CICE_CONFIG_OPTS)
With partial bottom cells, the k index of DZT runs from 0 to km+1, and the driver was written with the assumption that it was indexed from 1 to km... so now that indexing is explicit.
GIAFECO_JRA_HR -> G1850ECOIAF_JRA_HR; updated task counts for high res runs with BGC (or at least with BGC and CICE)
Default for G1850ECOIAF_JRA_HR w/ t13 grid is now 14666 procs
Updated init_ecosys_init_file to version that has all cocco tracers
Keith helped me generate the file with restoring time scale, and we are restoring back to initial conditions (which came from obs)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The highres_JRA_cice
testmod dir also highlights modifications we need to make for our experimental run; we should roll anything we can into the compset definition (such as changing OCN_BGC_CONFIG
).
Also, this PR will need updated tavg defaults for the high res.
|
||
<fesedflux_input%filename ocn_grid="gx3v7">ocn/pop/gx3v7/forcing/fesedflux_gx3v7_cesm1_97_2017.nc</fesedflux_input%filename> | ||
<fesedflux_input%filename ocn_grid="gx1v6">ocn/pop/gx1v6/forcing/fesedfluxTot_gx1v6_cesm2_2018_c180618.nc</fesedflux_input%filename> | ||
<fesedflux_input%filename ocn_grid="gx1v7">ocn/pop/gx1v6/forcing/fesedfluxTot_gx1v6_cesm2_2018_c180618.nc</fesedflux_input%filename> | ||
<fesedflux_input%filename ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/fesedflux_total_reduce_oxic_POP_tx0.1v3.c200420.nc</fesedflux_input%filename> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
@@ -1392,6 +1394,7 @@ | |||
<feventflux_input%filename ocn_grid="gx3v7">ocn/pop/gx3v7/forcing/feventflux_gx3v7_5gmol_cesm1_97_2017.nc</feventflux_input%filename> | |||
<feventflux_input%filename ocn_grid="gx1v6">ocn/pop/gx1v6/forcing/feventflux_gx1v6_5gmol_cesm1_97_2017.nc</feventflux_input%filename> | |||
<feventflux_input%filename ocn_grid="gx1v7">ocn/pop/gx1v6/forcing/feventflux_gx1v6_5gmol_cesm1_97_2017.nc</feventflux_input%filename> | |||
<feventflux_input%filename ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/fesedflux_total_reduce_oxic_POP_tx0.1v3.c200420.nc</feventflux_input%filename> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
<nox_flux_monthly_input%filename ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/ndep_ocn_1850_w_nhx_emis_tx0.1v3_c191115.nc</nox_flux_monthly_input%filename> | ||
<nhy_flux_monthly_input%filename ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/ndep_ocn_1850_w_nhx_emis_tx0.1v3_c191115.nc</nhy_flux_monthly_input%filename> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
@@ -1597,19 +1603,22 @@ | |||
<restore_data_filenames_derived ocn_grid="gx3v7">ocn/pop/gx3v7/forcing/ecosys_restore_POP_gx3v7_20170113.nc</restore_data_filenames_derived> | |||
<restore_data_filenames_derived ocn_grid="gx1v6">ocn/pop/gx1v6/forcing/ecosys_restore_POP_gx1v6_20170113.nc</restore_data_filenames_derived> | |||
<restore_data_filenames_derived ocn_grid="gx1v7">ocn/pop/gx1v6/forcing/ecosys_restore_POP_gx1v6_20170113.nc</restore_data_filenames_derived> | |||
<restore_data_filenames_derived ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/ecosys_jan_IC_omip_POP_tx0.1v3_c200602.nc</restore_data_filenames_derived> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
|
||
<restore_inv_tau_const>0.0</restore_inv_tau_const> | ||
|
||
<restore_inv_tau_input%filename ocn_grid="gx3v7">ocn/pop/gx3v7/forcing/ecosys_restore_inv_tau_POP_gx3v7_20170125.nc</restore_inv_tau_input%filename> | ||
<restore_inv_tau_input%filename ocn_grid="gx1v6">ocn/pop/gx1v6/forcing/ecosys_restore_inv_tau_POP_gx1v6_20170125.nc</restore_inv_tau_input%filename> | ||
<restore_inv_tau_input%filename ocn_grid="gx1v7">ocn/pop/gx1v6/forcing/ecosys_restore_inv_tau_POP_gx1v6_20170125.nc</restore_inv_tau_input%filename> | ||
<restore_inv_tau_input%filename ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/ecosys_restore_inv_tau_POP_tx0.1v3_20200603.nc</restore_inv_tau_input%filename> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
@@ -1643,6 +1652,8 @@ | |||
<init_ecosys_init_file ocn_grid="gx1v7" ocn_bgc_config="cesm2.1+cocco">ocn/pop/gx1v6/ic/ecosys_jan_IC_gx1v6_20200325.nc</init_ecosys_init_file> | |||
<init_ecosys_init_file ocn_grid="gx1v7" ocn_bgc_config="SPECTRA1.0">ocn/pop/gx1v6/ic/ecosys_jan_IC_gx1v6_9p6z_20200325_3impCalc.nc</init_ecosys_init_file> | |||
|
|||
<init_ecosys_init_file ocn_grid="tx0.1v3">ocn/pop/tx0.1v3/tmp/ecosys_jan_IC_omip_POP_tx0.1v3_c200602.nc</init_ecosys_init_file> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Need to either move this to proper directory (and import it to inputdata repo) or generate updated version
<lname>2000_DATM%JRA-1p4-2018_SLND_CICE%CICE4_POP2_DROF%JRA-1p4-2018_SGLC_SWAV</lname> | ||
</compset> | ||
|
||
<compset> | ||
<!-- latest JRA forcing, ecosys, high-res --> | ||
<alias>G1850ECOIAF_JRA_HR</alias> | ||
<lname>1850_DATM%JRA-1p4-2018_SLND_CICE%CICE4_POP2%ECO_DROF%JRA-1p4-2018_SGLC_SWAV</lname> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
make sure other G_HR
compsets are also using CICE%CICE4
<ntasks_lnd>72</ntasks_lnd> | ||
<ntasks_rof>72</ntasks_rof> | ||
<ntasks_ice>828</ntasks_ice> | ||
<ntasks_ocn>14666</ntasks_ocn> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
still testing to determine best value to use here
./xmlchange -a CICE_CONFIG_OPTS="-trage 0" | ||
fi | ||
|
||
./xmlchange OCN_BGC_CONFIG=cesm2.1+cocco |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
hopefully we'll have a cesm2.2+cocco
option before the release
Externals_POP.cfg
Outdated
@@ -6,7 +6,7 @@ local_path = externals/CVMix | |||
required = True | |||
|
|||
[MARBL] | |||
tag = marbl0.35.0 | |||
tag = marbl0.37.0 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Will want to use tag from marbl-ecosys/MARBL#347
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
actually want tag from marbl-ecosys/MARBL#353
Note that memory limitations prevent us from actually using this file with the MCT cap; hopefully it'll work better with CMEPS
riv_nut is new, fesedflux has been updated, feventflux is a new file (previously was using fesedflux)
When I had first created the cocco variables, the long_names still mentioned "Small Phyto". Also, I changed how I used nco so the global attributes are clearer about the process.
By setting tavg_contents_override_file in user_nl_pop, the [grid]_tavg_contents file generated by all the ocn.*.tavg.csh scripts is renamed autogen_tavg_contents in Buildconf/popconf and then tavg_contents_override_file is copied into Buildconf/popconf/[grid]_tavg_contents (at which point build-namelist and buildnml take over to use the override file)
First implementation was applying tavg_contents_override_file before adding tidal mixing variables to [grid]_tavg_contents (missed this in testing because I was looking at a high-res case that didn't use the new tidal scheme, I think)
Adding some _2 variables to allow variables in multiple streams and _z_t_150m variables to accumulate only the top 150 m of output. Full list: TEMP_z_t_150m TEMP_2 SALT_z_t_150m SALT_2 UVEL_z_t_150m UVEL_2 VVEL_z_t_150m VVEL_2 SHF_QSW_2
Also increased max_avail_tavg_field
Note that this duplicate variable follows the tracer module's instructions regarding whether to truncate at 150 meters or provide full depth. Given the addition of _z_t_150m tavg output in the previous commit, it may make sense to force the _2 tavg output to use full depth (in which case it should be renamed _z_t)
This namelist variable sets tadvect_ctype for ecosystem tracers that have virtual fluxes.
6d72b75 allows us to select a different advection scheme for ecosystem tracers with virtual fluxes. Running with
Results in
The default for |
Files are no longer in ocn/pop/tx0.1v3/tmp/ -- they are either in forcing/ or ic/ Files moved: fesedflux_total_reduce_oxic_POP_tx0.1v3.c200616.nc feventflux_5gmol_POP_tx0.1v3.c200616.nc ndep_ocn_1850_w_nhx_emis_tx0.1v3_c191115.nc riv_nut.gnews_gnm.JRA025m_to_tx0.1v3_nnsm_e333r100_190226.20200616.nc ecosys_restore_POP_tx0.1v3_20200604.nc ecosys_restore_inv_tau_POP_tx0.1v3_20200603.nc ecosys_jan_IC_omip_POP_tx0.1v3_c200617.nc
I've finished testing --
|
Description of changes:
Create a new compset set up to run a
G
compset withJRA
forcing and the MARBL tracers enabled on thetx0.1v3
grid. Will also automatically use CICE4 for existing_HR
compsets.Testing:
Test case/suite: I'll run
aux_pop
andaux_pop_MARBL
when this is ready to be merged and update the results in a later commentTest status: I expect this to be bit-for-bit, since the high-res compsets are not tested.
Fixes #30
Fixes #29
User interface (namelist or namelist defaults) changes? None