GOSI10 porting to NEMO5

[isabellaascione]

1. Testing GOSI10 with NEMO5 beta

• NEMO5 beta released in on https://forge.nemo-ocean.eu/devs/ukmo/gosi/gosi/-/commits/5.0-beta?ref_type=tags

• Create issue and branch on GOSI repository https://forge.nemo-ocean.eu/devs/ukmo/gosi/gosi/-/issues/19

• Port GOSI10 package changes to NEMO5 beta. GOSI10 will have different settings compared to standard GOSI10 at 4.2.2

• set up a rose suite to run with the modified version of GOSI10 at NEMO5 beta .

$${\color{red} ✗ }$$ run 30 years : the run crashed after 25 years with errors in writing the variable sfxice

$${\color{red} ✗ }$$ validation notes : these cannot be run for the last 10 years of the experiment due to the failure above described

2. Testing GOSI10 with NEMO5 release candidate

• NEMO5 RC released on 5th August 2024
• Update issue and code on GOSI repository : the branch for testing is named 19-upgrade-to-nemo5
• Set up rose suite to run with NEMO5 RC
• run 30 y experiment
• validation notes

[isabellaascione]

Changes identified with NEMO5 beta release:

NEMO-ocean parameters

	parameter status at NEMO5 beta	GOSI10 at NEMO5 beta	GOSI10 Trial at NEMO 4.2.2 , for comparison with GOSI10 at NEMO5 beta	GOSI10 at NEMO 4.2.2 (standard)
namrun
ln_top	$${\color{green}new}$$	TRUE (but used only with key_top)	not available	not available
ln_rst_eos	$${\color{gray} N/A }$$ (GOSI10 only)	TRUE (ported from GOSI10 package)	TRUE	TRUE
namdom
ln_linssh	$${\color{red}removed}$$	N/A	FALSE	FALSE
ln_shuman	$${\color{green}new}$$	FALSE	not available	not available
namlbc
ln_shlat2d	$${\color{gray} N/A}$$ (GOSI10 4.2.2 only)	N/A (not ported)	FALSE	TRUE
rn_shlat		0	0	-9999
namsbc
ln_dm2dc	$${\color{blue}preexisting}$$	FALSE	FALSE	TRUE
namsbc_blk
ln_MFS	$${\color{green}new}$$	FALSE	not available	not available
ln_prec_met	$${\color{green}new}$$	FALSE	not available	not available
rn_Cd_ia	$${\color{orange}renamed}$$	rn_Cd_ia=1.4e-3	rn_Cd_i=1.4e-3	rn_Cd_i=1.4e-3
rn_Ce_ia	$${\color{orange}renamed}$$	rn_Ce_ia=1.4e-3	rn_Ce_i=1.4e-3	rn_Ce_i=1.4e-3
rn_Ch_ia	$${\color{orange}renamed}$$	rn_Ch_ia=1.4e-3	rn_Ch_i=1.4e-3	rn_Ch_i=1.4e-3
ln_Cx_ice_frm	$${\color{green}new}$$	false	not available	not available
nn_frm	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cs_io	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cs_ia	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cr_ia	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cr_io	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cf_ia	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
rn_Cf_io	$${\color{green}new}$$	set to default as per reference namelist, but not used because ln_Cx_ice_frm=.false.)	not available	not available
namsbc_abl
ln_pga_abl	$${\color{green}new}$$	false (only for ln_abl = T)	not available	not available
rn_vfac	$${\color{green}new}$$	(only for ln_abl = T)	not available	not available
namsbc_blk
ln_crt_fbk	$${\color{blue}preexisting}$$	TRUE	TRUE	FALSE
rn_stau_a	$${\color{blue}preexisting}$$	-2.90E-03	-2.90E-03	not used
rn_stau_b	$${\color{blue}preexisting}$$	8.00E-03	8.00E-03	not used
rn_vfac	$${\color{gray} N/A}$$ (GOSI10 4.2.2 only)	N/A (not ported)	0	1
namtra_qsr
nn_chldta	$${\color{blue}preexisting}$$	1	1	0
rn_chl_conc	$${\color{gray} N/A }$$ (GOSI10 4.2.2 only)	N/A (not ported)	not used	0.1
nameos
rn_T0	$${\color{green}new}$$	10	hardcoded to 10	hardcoded to 10
rn_S0	$${\color{green}new}$$	35	hardcoded to 35	hardcoded to 35
namtra_adv
nn_fct_imp	$${\color{green}new}$$	it should be 2 but that does not work , initially set to 1. Then it was set to 2 with this fix	not available	not available
namtra_eiv
ln_eke_equ
namldf_eke	$${\color{green}new}$$ namelist with new parameters (not listed in this table)	unused because nn_aei_ijk_t=21	not available	not available
namvvl
nn_vvl_interp	$${\color{blue}preexisting}$$	2	2	0
namdyn_adv
ln_dynadv_ubs	$${\color{red}removed}$$	N/A	FALSE	FALSE
namdyn_vor
ln_dynvor_eeT	$${\color{red}removed}$$	N/A	FALSE	FALSE
namdyn_hpg
ln_hpg_zps	$${\color{red}removed}$$	N/A	FALSE	FALSE
namdyn_spg
ln_bt_av	$${\color{red}removed}$$ now it is hardcoded to TRUE (ll_bt_av in dynspg_ts.F90)	N/A	TRUE	TRUE
rn_htau_scaling	$${\color{gray} N/A}$$ (GOSI10 4.2.2 only)	ported	1	1
nammpp
ln_mppdelay	$${\color{green}new}$$	FALSE	not available	not available
nn_hls	$${\color{blue}preexisting}$$	2	1 (should be 2, but does not have much impact )	1
namzdf
ln_zdfiwm	$${\color{blue}preexisting}$$	TRUE	TRUE	FALSE
ln_zdftmx	$${\color{gray}N/A }$$ (GOSI10 4.2.2 only)	N/A (not ported)	FALSE	TRUE
rn_avm0	$${\color{blue}preexisting}$$	1.40E-06	1.40E-06	1.20E-04
rn_avt0	$${\color{blue}preexisting}$$	1.00E-10	1.00E-10	1.20E-05

Sea-Ice parameters

parameter	parameter status in NEMO5 beta	GOSI10 at NEMO5 beta	GOSI10 trial at NEMO 4.2.2 (for comparison with GOSI10 at NEMO5 beta)	GOSI10 at NEMO 4.2.2 (standard)
rn_cio= 1.0e-2	$${\color{orange}renamed}$$	rn_Cd_io= 1.0e-2	rn_cio= 1.0e-2	rn_cio= 1.0e-2
ln_icedS	$${\color{red}removed}$$	N/A	ln_icedS=.true.	ln_icedS=.true.
ln_flushing	$${\color{green}new}$$	TRUE	Not available as option in the namelist but included in the code	Not available as option in the namelist but included in the code
ln_drainage	$${\color{green}new}$$	TRUE	Not available as option in the namelist but included in the code	Not available as option in the namelist but included in the code
rn_sinew	$${\color{green}new}$$	when using nn_icesal=2, initially set to 0.3 as per recommendation for NEMO5beta. then set to 0.25 as per reccomendation for NEMO5rc. rn_sinew *sss_1d is used instead of rn_simax, so will give different result from 4.2.2	N/A	N/A
rn_simax	$${\color{red}removed}$$	N/A	20	20
rn_sal_himin	$${\color{green}new}$$	0	Hardcoded to 0.	Hardcoded to 0.
ln_pnd_rain	$${\color{green}new}$$	FALSE	N/A	N/A
rn_pnd_hl_min	$${\color{green}new}$$	0.005	Not available as parameter, but hardcoded to 0.005	Not available as parameter, but hardcoded to 0.005
rn_pnd_hl_max	$${\color{green}new}$$	0.015	Not available as parameter, but hardcoded to 0.015	Not available as parameter, but hardcoded to 0.015
rn_alb_lpnd	$${\color{green}new}$$	set to 0.36, it should do the same as 4.2.2	N/A	N/A
nn_liquidus	$${\color{green}new}$$	1	N/A	N/A
nn_pnd_brsal	$${\color{gray} N/A }$$ (GOSI10 4.2.2 only )	N/A, use nn_liquidus instead	0	0
nn_drainage	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
nn_flushing	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_flushrate	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_vbrc	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
nn_sal_scheme	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_alpha_RJW	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_Rc_RJW	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_alpha_GN	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_Rc_GN	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
rn_alpha_CW	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A
ln_sal_chk	$${\color{green}new}$$	not used (only used with nn_icesal=4)	N/A	N/A

[isabellaascione]

BUG

the parameter nn_fct_imp has been introduced with the 5-beta release. The default value is 1. To replicate the behaviour of 4.2.2, it should be set to 2. However this option seems to be broken. I let Daley know about this, and an issue should be opened on the main nemo repository.

update: Daley proposed a fix:
https://forge.nemo-ocean.eu/nemo/nemo/-/issues/499

[dcalve]

The r5282_nemo5rc_changes_and_fix_intel_xios MOCI branch should be used to run GOSI10 after it has been rebased onto the NEMO 5 release candidate (currently it is a branch of the beta).

This is required because the NEMO 5 release candidate has removed lib_cray.f90, necessitating an update to the fcm-make configurations in MOCI. The above branch contains these changes as well as those from r4684_fix_intel_xios, which we currently use.

Suites using GOSI10 based on the NEMO 5 release candidate will need the following update in app/fcm_make_ocean/rose-app.conf:

config_root_path=fcm:MOCI.xm/branches/dev/daleycalvert/r5282_nemo5rc_changes_and_fix_intel_xios
nemo_version_file=nemo5.0-version.cfg

[isabellaascione]

The NEMO 5 RC is released.

moving from NEMO5-beta to NEMO5-RC ,

the new reference namelist says to use 0.25 instead of 0.3 (as in NEMO5 beta) when using nn_icesal=2

Additionally, the following changes happened to the namelists:

namhsb	removed namelist
ln_diahsb	removed
namdom
ln_crs	removed
namcfg
ln_use_jattr	removed
namcrs	removed
nn_factx	removed
nn_facty	removed
nn_binref	removed
ln_msh_crs	removed
nn_crs_kz	removed
ln_crs_wn	removed
namc1d	it now includes the parameters from the namelists namc1d_dyndmp and namc1d_uvd, which were present in 5-beta. namc1d_dyndmp and namc1d_uvd no longer exist as namelists
nam_tide
ln_tide_dia	removed
namtra_eiv
ln_ldfeiv_dia	removed
namflo	removed
ln_floats	removed
jpnfl	removed
jpnnewflo	removed
ln_rstflo	removed
nn_writefl	removed
nn_stockfl	removed
ln_argo	removed
ln_flork4	removed
ln_ariane	removed
ln_flo_ascii	removed
nam_dia25h	removed
ln_dia25h	removed
namnc4
ln_nc4zip	default is .false.