New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
WorkingGroups/LIM/lim3_monocat – NEMO
wiki:WorkingGroups/LIM/lim3_monocat

First stesp on LIM3 monocategory

nemo v3_4, revision 4695: http://forge.ipsl.jussieu.fr/nemo/svn/branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM

  • which test and which errors ?

CHANGES NEEDED :

  1. copied cards 1 driver of ORCA1_LIM3 of claire, but changes needed for OPA_RDT, input files, etc...
  2. cp namelist_ref in namelist_gfg (if not error in supergrep opa9.driver in libigcm)
  3. keep modification of original namelist_cfg : jpidta, jpjdta, jperio, etc ... (if not error in tradmp)
  4. OPTION IOSERVER=TRUE it does not work (job still running but no output, solver.stat :0 it, in ocean.output: " * Euler time step switch is 0, * AAAAAAAA )
  5. in lim3.card restart: * [RestartFiles]
  • List= (${config_UserChoices_JobName}_${NEMO_END}_restart_ice.nc, restart_icemod.nc, restart_ice.nc)

' 1. * * '

TEST NAME CHANGES in MY_SRC ERRORS
ORCA2_LIM3_jpl1 par_ice.F90 (jpl = 1) NO
ORCA2_LIM3_jpl5_nlay2 par_ice.F90 (jpl = 1, nlay= 2) NaN in solver.start at kt=2 , exists files output.abort (to be rebuild)
ORCA2_LIM3_nolim_itme_zghe par_ice.F90 (jpl = 1)
limthd_dif.F90 (modify conductivity)
sbcice_lim.F90 (rm 'CALL lim_itd_me')
thd_ice.F90 (add thth variable: thick. for comp. of eq. thermal conduct)
after 8 years of simulation: zonal velocity max > 20 m/s

.... ongoing work

nemo v3_6, revision 4779: http://forge.ipsl.jussieu.fr/nemo/svn/trunk/NEMOGCM

  • which test and which errors ?
directory on ada TEST NAME CHANGES in /workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/modeles/NEMOGCM/CONFIG/ORCA2_LIM3_xxxxx/MY_SRC ERRORS COMMENT OUTPUT
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_jpl1 par_ice.F90 (jpl = 1) NO, 1Y ok
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_jpl5_nlay2 par_ice.F90 (jpl = 1, nlay= 2) NO, 1Y ok
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_nolim_itme_zghe/IGCM00 par_ice.F90 (jpl = 1)
limthd_dif.F90 (modify conductivity)
sbcice_lim.F90 (rm 'CALL lim_itd_me')
thd_ice.F90 (add thth variable: thick. for comp. of eq. thermal conduct)
'(1Y = 5475 iterations)
after 1 year of simulation:
KO at it: 5912 iter: 800 r: NaN b: NaN

see abort file here:
/workgpfs/rech/eee/!reee605/RUN_DIR/l3zdif.17201/REBUILD_20020101/
output.abort.nc'
test: move CALL lim_var_glo2eqv
CALL lim_itd_me_zapsmall in sbcice_lim.F90 and removed from limitd_me.F90,
.... compiling
.... running
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_jpl1 par_ice.F90 (jpl = 1) TEST OF 2 years after 1 year of simulation: KO at it: 5477 iter: 800 r: NaN b: NaN
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_nolim_itme/IGCM_limdyn_false/jpl1nody sbcice_lim.F90 (rm 'CALL lim_itd_me')
ln_limdyn = false (only thermo)
Error [CONetCDF4::CheckError(int status)] : In file '/workgpfs/rech/eee/!reee605/XIOS/src/output/onetcdf4.cpp', line 79 -> [ status = -60 ] NetCDF: Numeric conversion not representable put icemod enable=false in iodef.xml
OK: RUN 3 years without ERRORS
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_nolim_itme/IGCM_himax_init_5m sbcice_lim.F90 (rm 'CALL lim_itd_me')
namiceini = hti_ini_n = 5.0
after 1 year of simulation: KO at it: 5817 iter: 800 r: NaN b: NaN to see if hi_max is still ~3metres or not
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config ORCA2_LIM3_nolim_itme/IGCM_bulk_sal_4 sbcice_lim.F90 (rm 'CALL lim_itd_me')
num_sal = 1
bulk_sal = 4.0
after 1 year of simulation: KO at it: 5882 iter: 800 r: NaN b: NaN
/workgpfs/rech/eee/!reee605/lim3_monocat_nemo36/modipsl/config

NOTA: ORCA2_LIM3_nolim_itme_zghe version 4779, try to update to revision 4788, see ticket #1394: "Bug related to the porosity of ridges", re-compiled, and re-run

re-run after bug-fix et run for 5697 iterations.

run again 1 time: ...... ongoing

stop again at 5697 iterations

Simona Flavoni last update: Timestamp?

Last modified 9 years ago Last modified on 2014-10-28T10:53:28+01:00