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

add cpld_control_wave with 1deg wave grid #385

Closed
Closed
Show file tree
Hide file tree
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 10 additions & 0 deletions tests/fv3_conf/cpld_control_run.IN
Original file line number Diff line number Diff line change
Expand Up @@ -50,6 +50,16 @@ cp @[INPUTDATA_ROOT]/CICE_FIX/@[OCNRES]/grid_cice_NEMS_mx@[OCNRES].nc .
cp @[INPUTDATA_ROOT]/CICE_FIX/@[OCNRES]/kmtu_cice_NEMS_mx@[OCNRES].nc .
cp @[INPUTDATA_ROOT]/CICE_FIX/@[OCNRES]/mesh.mx@[OCNRES].nc .

# WW3 fix/input
if [[ $CPLWAV == .T. && $CPLWAV2ATM == .T. ]]; then
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We would need this only if CPLWAV is true, cplwave2atm=f would still need mod_defs and the input file

cp @[INPUTDATA_ROOT_WW3]/mod_def.* .
if [[ $RT1DEG == .T. ]]; then
cp @[INPUTDATA_ROOT_WW3]/ww3_multi_1deg.inp ww3_multi.inp
else
cp @[INPUTDATA_ROOT_WW3]/ww3_multi.inp .
fi
fi

# No restart
if [ $WARM_START = .F. ]; then
# ICs
Expand Down
93 changes: 93 additions & 0 deletions tests/tests/cpld_control_wave
Original file line number Diff line number Diff line change
@@ -0,0 +1,93 @@
#
# cpld_control_wave test
#

export TEST_DESCR="Fully coupled FV3-CCPP-MOM6-CICE-CMEPS-WW3 system - C96MX100 with waves"

export CNTL_DIR="cpld_control_wave"

export LIST_FILES="phyf024.tile1.nc \
phyf024.tile2.nc \
phyf024.tile3.nc \
phyf024.tile4.nc \
phyf024.tile5.nc \
phyf024.tile6.nc \
dynf024.tile1.nc \
dynf024.tile2.nc \
dynf024.tile3.nc \
dynf024.tile4.nc \
dynf024.tile5.nc \
dynf024.tile6.nc \
RESTART/coupler.res \
RESTART/fv_core.res.nc \
RESTART/fv_core.res.tile1.nc \
RESTART/fv_core.res.tile2.nc \
RESTART/fv_core.res.tile3.nc \
RESTART/fv_core.res.tile4.nc \
RESTART/fv_core.res.tile5.nc \
RESTART/fv_core.res.tile6.nc \
RESTART/fv_srf_wnd.res.tile1.nc \
RESTART/fv_srf_wnd.res.tile2.nc \
RESTART/fv_srf_wnd.res.tile3.nc \
RESTART/fv_srf_wnd.res.tile4.nc \
RESTART/fv_srf_wnd.res.tile5.nc \
RESTART/fv_srf_wnd.res.tile6.nc \
RESTART/fv_tracer.res.tile1.nc \
RESTART/fv_tracer.res.tile2.nc \
RESTART/fv_tracer.res.tile3.nc \
RESTART/fv_tracer.res.tile4.nc \
RESTART/fv_tracer.res.tile5.nc \
RESTART/fv_tracer.res.tile6.nc \
RESTART/phy_data.tile1.nc \
RESTART/phy_data.tile2.nc \
RESTART/phy_data.tile3.nc \
RESTART/phy_data.tile4.nc \
RESTART/phy_data.tile5.nc \
RESTART/phy_data.tile6.nc \
RESTART/sfc_data.tile1.nc \
RESTART/sfc_data.tile2.nc \
RESTART/sfc_data.tile3.nc \
RESTART/sfc_data.tile4.nc \
RESTART/sfc_data.tile5.nc \
RESTART/sfc_data.tile6.nc \
RESTART/MOM.res.nc \
RESTART/iced.2016-10-04-00000.nc \
RESTART/ufs.cpld.cpl.r.2016-10-04-00000.nc \
out_grd.glo_1deg"

export_fv3
export_cpl

# set component and coupling timesteps
export DT_ATMOS='450'
export DT_CICE=${DT_ATMOS}
export DT_DYNAM_MOM6='900'
export DT_THERM_MOM6='1800'
export CPL_SLOW=${DT_THERM_MOM6}
export CPL_FAST=${DT_ATMOS}

export RESTART_INTERVAL="12"
export RESTART_N=${RESTART_INTERVAL}

export TASKS=$TASKS_cpl_wwav
export TPN=$TPN_cpl_wwav
export INPES=$INPES_cpl_wwav
export JNPES=$JNPES_cpl_wwav
export THRD=$THRD_cpl_wwav
export WRTTASK_PER_GROUP=$WPG_cpl_wwav

export med_petlist_bounds=$MPB_cpl_wwav
export atm_petlist_bounds=$APB_cpl_wwav
export ocn_petlist_bounds=$OPB_cpl_wwav
export ice_petlist_bounds=$IPB_cpl_wwav
export wav_petlist_bounds=$WPB_cpl_wwav
Copy link
Collaborator

@DeniseWorthen DeniseWorthen Jan 19, 2021

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This will set the tasking to the tasking when using waves, which is appropriate for the 1/4 deg model but not for 1deg. You'll need to create a new set of PE tasking for 'cpl_dflt' which is for the 1deg case, but which includes waves, maybe call it 'cpl_dflt_wwav'.


# nems.configure
export NEMS_CONFIGURE="nems.configure.cpld_wave.IN"
export coupling_interval_slow_sec=${CPL_SLOW}
export coupling_interval_fast_sec=${CPL_FAST}

export CPLWAV='.T.'
export CPLWAV2ATM='.T.'
export RT1DEG='.T.'
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

In default_vars, we set the default resolution for each component using "OCNRES" etc. I would suggest to create a default WW3 res which is 1deg. You wouldn't need the RT1DEG variable because the default tests would be 1 deg ww3, 1deg ocean, 1deg ice and c96 fv3.

The current ww3 tests which are all for 1/4 deg would then need to set a non-default resolution. This would be consistent w/ how we vary the resolution across the different tests.

Copy link
Collaborator Author

@aliabdolali aliabdolali Jan 19, 2021

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@DeniseWorthen Thanks for the suggestions. I modified the test but still get a failure.
Here are my changes:

  1. I added cpl_dflt_wwav. Do you have any suggestions for the numbers? I assumed the 1deg is lighter than 1/4, and heavier than dflt, so I set the total to 288. What each of these variables stand for? I mean where can I see the definitions?
  2. I removed export RT1DEG='.T.'
  3. I added WAVRES='1.00'

export FV3_RUN=cpld_control_run.IN