Configuration peru12_nemo

• Objet: Vents Quikscat au lieu de erai
• demarrage avec le restart 2000 de erai70c au lieu de demmarrage Levitus en 1990

================================================================================ diff erai70c / erai72c ================================================================================ -------------------------------------------------------------------------------- diff -b -w erai70c/param_inputs/namelist.base.oce erai72c/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 191,192c191,192 < sn_wndi = 'u10' , 6 , 'u10' , .true. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' < sn_wndj = 'v10' , 6 , 'v10' , .true. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' --- > sn_wndi = 'uv10' , 24 , 'zonal_wind_speed' , .true. , .false. , 'yearly' ,'reshape_QUIKSCAT_PERU12_bicubic.nc', '' > sn_wndj = 'uv10' , 24 , 'meridional_wind_speed' , .true. , .false. , 'yearly' ,'reshape_QUIKSCAT_PERU12_bicubic.nc', '' -------------------------------------------------------------------------------- diff -b -w erai70c/scripts/namelist_exp.sh erai72c/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 32c32 < export LEVITUS=1 # for init state ... si demarrage climato (levitus ou autre) --- > export LEVITUS=0 # for init state ... 1 si demarrage climato (levitus ou autre) / 0 si restart 41c41 < export YEAR_BEGIN_EXP=1990 --- > export YEAR_BEGIN_EXP=2000 45c45 < export EXP_DUR_MTH=$(( 12 * 18 )) --- > export EXP_DUR_MTH=$(( 12 * 8 )) -------------------------------------------------------------------------------- diff -b -w erai70c/scripts/getfile_oce.sh erai72c/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 35,36c38,40 < ${io_getfile} ${INPUTDIR}/${BLK}/${YEAR_BEGIN_JOB}/u10.${YEAR_BEGIN_JOB}.ashei.PERU_075.nc u10_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BLK}/${YEAR_BEGIN_JOB}/v10.${YEAR_BEGIN_JOB}.ashei.PERU_075.nc v10_y${YEAR_BEGIN_JOB}.nc --- > ${io_getfile} ${DATA}/quikscat/quikscat_1d_${YEAR_BEGIN_JOB}_filled.nc uv10_y${YEAR_BEGIN_JOB}.nc

• Objet: OBC mercator ORCA025_LIM-T323 (au lieu de drakkar ORCA025-B83) + bathy trop12 (au lieu de Stephane) + Etat initial ad'hoc
• SBC solaire: Diurnal cycle sur mean 1 day au lieu de 6h
• SBC lw: interpolation true au lieu de false
• SSS On enlève le rappel en salinité
• OBC Mercator ORCA025_LIM-T323 (3 jours au lieu de 5 jours pour ORCA025-B83
• TKE ln_mxl0 = .true. au lieu de false ! surface mixing length scale = F(wind stress) (T) or not (F)
• RST Demarrage Etat moyen Mercator ORCA025_LIM-T323 au lieu de restart erai70c 2000
• Année 2008 en plus!!
• Bathy trop12 au lieu de Stephane
• OBC year + 1 car problème avec le découpage en 3 jours => va chercher l'année suivante pour les derniers jours...
• Compilation: Déplacement d'un write + ajout du flag -fno-alias
• W_damping=1

================================================================================ ===== diff erai72c / erai73c ================================================================================ -------------------------------------------------------------------------------- diff -b -w erai72c/param_inputs/namelist.base.oce erai73c/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 139c139 < ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave --- > ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave 193,194c193,194 < sn_qsr = 'radsw' , 6 , 'ssrd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' < sn_qlw = 'radlw' , 6 , 'strd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' --- > sn_qsr = 'radsw' , 24 , 'ssrd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' > sn_qlw = 'radlw' , 6 , 'strd' , .true. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' 300c300 < nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) --- > nn_sssr = 0 ! add a damping term in the surface freshwater flux (=2) 368,379c368,379 < sn_obcn(1) = 'obc_north_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(2) = 'obc_north_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(3) = 'obc_north_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(4) = 'obc_north_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(1) = 'obc_west_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(2) = 'obc_west_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(3) = 'obc_west_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(4) = 'obc_west_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(1) = 'obc_south_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(2) = 'obc_south_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(3) = 'obc_south_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(4) = 'obc_south_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' --- > sn_obcn(1) = 'obc_north_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(2) = 'obc_north_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(3) = 'obc_north_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(4) = 'obc_north_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(1) = 'obc_west_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(2) = 'obc_west_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(3) = 'obc_west_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(4) = 'obc_west_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(1) = 'obc_south_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(2) = 'obc_south_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(3) = 'obc_south_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(4) = 'obc_south_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' 700,701c700,701 < ln_mxl0 = .false. ! surface mixing length scale = F(wind stress) (T) or not (F) < rn_mxl0 = 0.4 ! surface buoyancy lenght scale minimum value --- > ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) > rn_mxl0 = 0.01 ! surface buoyancy lenght scale minimum value -------------------------------------------------------------------------------- diff -b -w erai72c/scripts/namelist_exp.sh erai73c/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 32c32 < export LEVITUS=0 # for init state ... 1 si demarrage climato (levitus ou autre) / 0 si restart --- > export LEVITUS=1 # for init state ... 1 si demarrage climato (levitus ou autre) / 0 si restart 45,46c45,46 < export EXP_DUR_MTH=$(( 12 * 8 )) < #export EXP_DUR_MTH=12 --- > export EXP_DUR_MTH=$(( 12 * 9 )) > #export EXP_DUR_MTH=0 -------------------------------------------------------------------------------- diff -b -w erai72c/scripts/getfile_oce.sh erai73c/scripts/getfile_oce.sh -------------------------------------------------------------------------------- > # bathy & coordinates 5c5,6 < ${io_getfile} ${INPUTDIR}/Stephane/PERU12-I/bathymetry_PERU12_mod2.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean.nc bathy_meter.nc 104c61,62 < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_S_y${YEAR_BEGIN_JOB}.nc . 108c66,67 < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_S_y${YEAR_BEGIN_JOB}.nc . 112c71,72 < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_S_y${YEAR_BEGIN_JOB}.nc . 113a74,89 > if [ ${YEAR_BEGIN_JOBp1} -ne ${YEAR_BEGIN_JOB} ] ; then > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_S_y${YEAR_BEGIN_JOBp1}.nc . > > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_S_y${YEAR_BEGIN_JOBp1}.nc . > > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_S_y${YEAR_BEGIN_JOBp1}.nc . > fi -------------------------------------------------------------------------------- diff -b -w erai72c/compile/nemo/arch-X64_CURIE.fcm erai73c/compile/nemo/arch-X64_CURIE.fcm -------------------------------------------------------------------------------- 22c22 < %FCFLAGS -i4 -r8 -O3 -xAVX -fp-model precise --- > %FCFLAGS -i4 -r8 -O3 -xAVX -fp-model precise -fno-alias -------------------------------------------------------------------------------- diff -b -w -r -x .svn erai72c/compile/nemo/MY_SRC erai73c/compile/nemo/MY_SRC -------------------------------------------------------------------------------- diff -b -w -r -x .svn erai72c/compile/nemo/MY_SRC/dynldf_lap.F90 erai73c/compile/nemo/MY_SRC/dynldf_lap.F90 76a77 > IF(lwp) WRITE(numout,*) 'Chris sponge sur la dynamique pour OBC' 79d79 < WRITE(numout,*) 'Clitandre sponge sur la dynamique' diff -b -w -r -x .svn erai72c/compile/nemo/MY_SRC/traldf_lap.F90 erai73c/compile/nemo/MY_SRC/traldf_lap.F90 86a87 > IF(lwp) WRITE(numout,*) 'Chris sponge sur les traceurs pour OBC' 100d100 < WRITE(numout,*) 'Chris sponge sur les traceurs' -------------------------------------------------------------------------------- cmp erai72c/compile/nemo/nemo.exe erai73c/compile/nemo/nemo.exe -------------------------------------------------------------------------------- erai72c/compile/nemo/nemo.exe erai73c/compile/nemo/nemo.exe differ: char 42, line 1

• SBC solaire: Retour à 6h
• SBC lw: Retour à pas d'interpolation
• SSS Retour au rappel en salinité
• TKE ln_mxl0 = retour à false au lieu de true
• Ne change pas grand chose => on reste sur erai73c car plus cohérent

================================================================================ diff erai73c / erai73c2 ================================================================================ -------------------------------------------------------------------------------- diff -b -w erai73c/param_inputs/namelist.base.oce erai73c2/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 139c139 < ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave --- > ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 193,194c193,194 < sn_qsr = 'radsw' , 24 , 'ssrd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' < sn_qlw = 'radlw' , 6 , 'strd' , .true. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' --- > sn_qsr = 'radsw' , 6 , 'ssrd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' > sn_qlw = 'radlw' , 6 , 'strd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' 300c300 < nn_sssr = 0 ! add a damping term in the surface freshwater flux (=2) --- > nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 700c700 < ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) --- > ln_mxl0 = .false. ! surface mixing length scale = F(wind stress) (T) or not (F) -------------------------------------------------------------------------------- diff -b -w erai73c/scripts/getfile_oce.sh erai73c2/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 23c23 < ${io_getfile} ${DATA}/erai_peru_mask/freq_1_day/${YEAR_BEGIN_JOB}/ssrd.${YEAR_BEGIN_JOB}.fshei.PERU_075.nc radsw_y${YEAR_BEGIN_JOB}.nc --- > ${io_getfile} ${DATA}/erai_peru_mask/${YEAR_BEGIN_JOB}/ssrd.${YEAR_BEGIN_JOB}.fshei.PERU_075.nc radsw_y${YEAR_BEGIN_JOB}.nc 44a45,50 > # &namsbc_ssr ! surface boundary condition : sea surface restoring > ${io_getfile} ${INPUTDIR}/Stephane/PERU12-I/sss_1m.nc sss_1m.nc

• Objet: bathy Stephane au lieu de trop12 sur OBC mercator ORCA025_LIM-T323

================================================================================ ===== diff erai73c / erai74c ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w erai73c/scripts/module.sh erai74c/scripts/module.sh -------------------------------------------------------------------------------- 19a20 > module switch intel intel/12.1.7.256 -------------------------------------------------------------------------------- diff -b -w erai73c/scripts/getfile_oce.sh erai74c/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 5c5 < ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIR}/Stephane/PERU12-I/bathymetry_PERU12_mod2.nc bathy_meter.nc 50c50 < ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/nomask/ORCA025_LIM-T323_clim_2000-2008_gridT_votemper_Perou.nc dta_temp.nc --- > ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/bathy_Stephane/ORCA025_LIM-T323_clim_2000-2008_gridT_votemper_Perou.nc dta_temp.nc 53c53 < ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/nomask/ORCA025_LIM-T323_clim_2000-2008_gridS_vosaline_Perou.nc dta_sal.nc --- > ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/bathy_Stephane/ORCA025_LIM-T323_clim_2000-2008_gridS_vosaline_Perou.nc dta_sal.nc -------------------------------------------------------------------------------- diff -b -w erai73c/scripts/header.CURIE erai74c/scripts/header.CURIE -------------------------------------------------------------------------------- 78c78 < export OBC="obc/ORCA025_LIM-T323" --- > export OBC="obc/ORCA025_LIM-T323/bathy_Stephane" ================================================================================

Objet: implementation BDY dans little_nemo r85

passage r73 => r85 (on vire dynadv_ubs.F90 (bug UBS 1/32 1/8) et sbcblk_core.F90 (pour l'understress) de MY_SRC)
+ ajout dans la namelist de "rn_usecrt" pour gérer avec ou sans understress
Run de contrôle sur 15 jours pour retrouver erai74c => (cdo -diffv 0 errors found OK)

Merge des codes: Julien (NEMO_r4018 + MY_SRC) avec Little_nemo_r85

Reste le problème de la modif du calcul de flux de chaleur lié à emp (sbcblk_core.F90 trasbc.F90)
qui crée un bug alors qu'on devrait retrouver le même resultat en remplaçant T2m par tsurf...

- On écrase le BDY de little_NEMO par celui de NEMO_r4018 JULIEN_MY_SRC="/ccc/work/cont005/ra0542/hourdinc/sources/JULIEN_BDY/MY_SRC" NEMO_r4018="/ccc/work/cont005/ra0542/hourdinc/sources/NEMO_r4018_base_JULIEN/NEMO/OPA_SRC" little_NEMO_r85_BDY="/ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC" COMPDIR="/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/bdy/compile/nemo/MY_SRC" ------------------------------- A prendre dans MY_SRC Julien ------------------------------- BDY: bdy_oce.F90 bdydta.F90 bdydyn.F90 bdydyn2d.F90 bdyice_lim.F90 bdyini.F90 bdytides.F90 DOM: domvvl.F90 DOM: dommsk.F90 (identique entre Little_nemo_r85 et NEMO_r4018) dom_msk_nsa dans dom_msk QUESTCE QUE C'est??? dommsk.F90: !! dom_msk : compute land/ocean mask dommsk.F90: !! dom_msk_nsa : update land/ocean mask when no-slip accurate option is used. DOM: dom_oce.F90 DOM: domain.F90 DYN: dynnxt.F90 (identique entre Little_nemo_r85 et NEMO_r4018) DYN: dynspg.F90 DYN: dynspg_flt.F90 DYN: dynspg_ts.F90 DYN: dynzdf.F90 (identique entre Little_nemo_r85 et NEMO_r4018) DYN: dynzdf_imp.F90 LBC: lib_mpp.F90 (... to optimize the BDY and OBC communications) SBC: sbcblk_core.F90 DYN: sshwzv.F90 TRA: traadv_ubs.F90 (identique entre Little_nemo_r85 et NEMO_r4018) TRA: tradmp.F90 (identique entre Little_nemo_r85 et NEMO_r4018) TRA: traqsr.F90 TRD: trdmod.F90 (identique entre Little_nemo_r85 et NEMO_r4018) TRD: trdmod_oce.F90 (identique entre Little_nemo_r85 et NEMO_r4018) SBC: updtide.F90 ZDF: zdftke.F90 nemogcm.F90 oce.F90 step.F90 stpctl.F90 (set the minimum value of tke) (identique entre Little_nemo_r85 et NEMO_r4018) cd ${JULIEN_MY_SRC} cp bdy_oce.F90 bdydta.F90 bdydyn.F90 bdydyn2d.F90 bdyice_lim.F90 bdyini.F90 bdytides.F90 domvvl.F90 dommsk.F90 dom_oce.F90 domain.F90 dynnxt.F90 dynspg.F90 dynspg_flt.F90 dynspg_ts.F90 dynzdf.F90 dynzdf_imp.F90 lib_mpp.F90 sbcblk_core.F90 sshwzv.F90 traadv_ubs.F90 tradmp.F90 traqsr.F90 trdmod.F90 trdmod_oce.F90 updtide.F90 zdftke.F90 nemogcm.F90 oce.F90 step.F90 stpctl.F90 ${COMPDIR} ------------------------------- Ne pas prendre dans MY_SRC Julien ------------------------------- divcur.F90 (commentaire JMM) dynadv_ubs.F90 (shema 4eme ordre - diagnostique couteux) dynnept.F90 (ne marche pas et doit etre enlevee) domhgr.F90 dtatsd.F90 (separe rappel et initialisation LEVITUS) obc_par.F90 obc_par_GOLFO36.h90 obc_par_NATL025.h90 obc_par_TROPICAL.h90 obcdta.F90 obcini.F90 obcvol.F90 par_GOLFO36.h90 par_NATL025.h90 par_TROPICAL.h90 par_oce.F90 ------------------------------- Recuperer dans NEMO_r4018 ------------------------------- SBC: sbcapr.F90 sbctide.F90 tideini.F90 sbcrnf.F90 SBC: sbc_oce.F90 (pour nn_ice_embd mais ATTENTION!! vire emps qui est utilisé dans plusieurs routines!!!!) SBC: dynadv.F90 (identique entre Little_nemo_r85 et NEMO_r4018) LBC: lbclnk.F90 (pour lbc_bdy_lnk) DOM: domzgr_substitute.h90 (identique entre Little_nemo_r85 et NEMO_r4018) cd ${NEMO_r4018} cd SBC cp sbcapr.F90 sbctide.F90 tideini.F90 sbcrnf.F90 sbc_oce.F90 ${COMPDIR} cd ../LBC cp lbclnk.F90 ${COMPDIR} cd ../DOM cp domzgr_substitute.h90 ${COMPDIR} ------------------------------- N'existe pas dans Little_nemo_r85 et NEMO_r4018 et pas forcement utile??? ------------------------------- flxish.F90 !!========================================================================== !! *** MODULE flxish *** !! add net heat and fresh water flux from ice shelf melting !! into the adjacent ocean using the parameterisation by !! Beckmann and Goosse (2003), "A parameterization of ice shelf-ocean !! interaction for climate models", Ocean Modelling 5(2003) 157-170. !! (hereafter BG) !!========================================================================== limdmp_2.F90 !!====================================================================== !! *** MODULE limdmp_2 *** !! LIM-2 ice model : restoring Ice thickness and Fraction leads !!====================================================================== limsbc_2.F90 !!====================================================================== !! *** MODULE limsbc_2 *** !! LIM-2 : updates the fluxes at the ocean surface with ice-ocean fluxes !!====================================================================== ------------------------------- Doublons ------------------------------- nemogcm.F90.jc step.F90.jc ------------------------------- gmake: *** No rule to make target `asmbkg.o', needed by `nemogcm.o'. Stop. cd ${NEMO_r4018}/ASM; cp asmbkg.F90 ${COMPDIR} ------------------------------- gmake: *** No rule to make target `icbini.o', needed by `nemogcm.o'. Stop. new features in nemo_v3_5_beta: ICB: iceberg floats (Bigg et al 1997, Martin and Adcroft 2010) cd ${NEMO_r4018}; cp -r ICB ${little_NEMO_r85_BDY} ------------------------------- gmake: *** No rule to make target `cyclone.o', needed by `sbcblk_core.o'. Stop. cd ${NEMO_r4018}/SBC; cp cyclone.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/solmat.f90(232): error #6285: There is no matching specific subroutine for this generic subroutine call. [LBC_LNK_E] CALL lbc_lnk_e( gcp (:,:,1), c_solver_pt, 1. ) ! lateral boundary conditions dans LBC/lbclnk.F90: INTERFACE lbc_lnk_e MODULE PROCEDURE mpp_lnk_2d_e END INTERFACE dans LBC/lib_mpp.F90: SUBROUTINE mpp_lnk_2d_e( pt2d, cd_type, psgn, jpri, jprj ) incohérent avec le nombre d'arguments dans l'appel dans solmat.F90. => On récupère tout le répertoire SOL de NEMO_r4018 qui n'a visiblement changé que ce passage d'arguments dans toutes ses routines. rm -rf ${little_NEMO_r85_BDY}/SOL ; cp -r ${NEMO_r4018}/SOL ${little_NEMO_r85_BDY} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/closea.f90(248): error #6911: The syntax of this substring is invalid. [EMPS] emps(:,:) = emps(:,:) + zfwf(jc) / surf(jpncs+1) !! *** MODULE closea *** !! Closed Seas : specific treatments associated with closed seas cp ${NEMO_r4018}/DOM/closea.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/closea.f90(24): error #6580: Name in only-list does not exist. [DDPDD] USE lib_fortran, ONLY: glob_sum, DDPDD cp ${NEMO_r4018}/lib_fortran.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcana.f90(123): error #6911: The syntax of this substring is invalid. [EMPS] emps(:,:) = rn_emp0 !! *** MODULE sbcana *** !! Ocean forcing: analytical momentum, heat and freshwater forcings cp ${NEMO_r4018}/SBC/sbcana.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/domzgr.f90(524): error #6404: This name does not have a type, and must have an explicit type. [NCLOSEA] IF( nclosea == 0 ) THEN ! NO closed seas or lakes ! cp ${NEMO_r4018}/DOM/domzgr.F90 ${COMPDIR} !! 3.4 ! 2012-08 (J. Siddorn) added Siddorn and Furner stretching function ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/traqsr.f90(182): error #6404: This name does not have a type, and must have an explicit type. [R1_RAU0_RCP] qsr_hc(:,:,jk) = r1_rau0_rcp * ( etot3(:,:,jk) - etot3(:,:,jk+1) ) !! *** MODULE traqsr *** !! Ocean physics: solar radiation penetration in the top ocean levels r1_rau0_rcp définit dans phycst.F90 cp ${NEMO_r4018}/DOM/phycst.F90 ${COMPDIR} !! *** MODULE phycst *** !! Definition of of both ocean and ice parameters used in the code /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/trabbc.f90(178): error #6404: This name does not have a type, and must have an explicit type. [RO0CPR] qgh_trd0(:,:) = ro0cpr * rn_geoflx_cst en fait dans NEMO_r4018, ro0cpr est renommé r1_rau0_rcp (= 1. / ( rau0 * rcp )) On détruit donc traqsr.f90 récupéré dans MY_SRC de JULLIEN qui permettait d'être cohérent avec le phycst.F90 de NEMO_r4018 (avec r1_rau0_rcp) et on garde le traqsr.f90 de little_NEMO_r85 qui est cohérent avec le phycst.F90 de little_NEMO_r85 (avec ro0cpr) rm -f ${COMPDIR}/traqsr.f90 ${COMPDIR}/phycst.F90 ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcflx.f90(185): error #6911: The syntax of this substring is invalid. [EMPS] emps(:,:) = emp (:,:) ! Initialization of emps (needed when no ice model) cp ${NEMO_r4018}/SBC/sbcflx.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcblk_clio.f90(379): error #6911: The syntax of this substring is invalid. [EMPS] emps(:,:) = emp(:,:) finalement on vire toutes les routines qui font appel à emps à cause du remplacement de la routine sbc_oce.F90 pour avoir le mot clé nn_ice_embd. On le rajoutera quand nécessaire rm -f closea.F90 lib_fortran.F90 sbcana.F90 sbcflx.F90 sbc_oce.F90 ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/zdftke.f90(89): error #6401: The attributes of this name conflict with those made accessible by a USE statement. [AVT_K] /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/zdftke.f90(89): error #6401: The attributes of this name conflict with those made accessible by a USE statement. [AVM_K] /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/zdftke.f90(90): error #6401: The attributes of this name conflict with those made accessible by a USE statement. [AVMU_K] /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/zdftke.f90(90): error #6401: The attributes of this name conflict with those made accessible by a USE statement. [AVMV_K] REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avt_k , avm_k ! not enhanced Kz REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avt_k , avm_k ! not enhanced Kz REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmu_k, avmv_k ! not enhanced Kz REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmu_k, avmv_k ! not enhanced Kz avt_k avm_k avmu_k avmv_k sont définit dans zdf_oce.F90 dans little_NEMO_r85 et pas dans celui de NEMO_r4018 C'est la seule différence dans zdf_oce.F90 entre les 2 versions NEMO. On récupère donc celui de NEMO_r4018 dans little_NEMO_r85 cp ${NEMO_r4018}/ZDF/zdf_oce.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcrnf.f90(153): error #6404: This name does not have a type, and must have an explicit type. [R1_RAU0] rnf_tsc(:,:,jp_tem) = ( sf_t_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rau0 La version récupérée dans NEMO_r4018 a besoin de la version ad'hoc de phycst.F90 pour les définitions de constantes. On conserve la version de phycst.F90 de little_NEMO_r85 et on ajoute les définitions nécessaires pour les routines récupérées pour BDY. cp ${little_NEMO_r85_BDY}/DOM/phycst.F90 ${COMPDIR} ajout de r1_rau0: diff phycst.F90 /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90 56d55 < REAL(wp), PUBLIC :: r1_rau0 !: = 1. / rau0 [m3/kg] 116d114 < r1_rau0 = 1._wp / rau0 158d155 < WRITE(numout,*) ' 1. / rau0 r1_rau0 = ', r1_rau0, ' m^3/kg' ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/icbini.f90(350): error #6404: This name does not have a type, and must have an explicit type. [LN_ICEBERGS] NAMELIST/namberg/ ln_icebergs , ln_bergdia , nn_sample_rate , rn_initial_mass , cp ${little_NEMO_r85_BDY}/SBC/sbc_oce.F90 ${COMPDIR} On ajoute la définition de ln_icebergs et nn_ice_embd [hourdinc@curie51:/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/bdy/compile/nemo/MY_SRC] diff /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90 . 42c42,47 < INTEGER , PUBLIC :: nn_ice = 0 !: flag on ice in the surface boundary condition (=0/1/2/3) --- > LOGICAL , PUBLIC :: ln_icebergs = .FALSE. !: Icebergs > INTEGER , PUBLIC :: nn_ice = 0 !: flag for ice in the surface boundary condition (=0/1/2/3) > INTEGER , PUBLIC :: nn_ice_embd = 0 !: flag for levitating/embedding sea-ice in the ocean > ! !: =0 levitating ice (no mass exchange, concentration/dilution effect) > ! !: =1 levitating ice with mass and salt exchange but no presure effect > ! !: =2 embedded sea-ice (full salt and mass exchanges and pressure) ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcblk_core.f90(211): error #6911: The syntax of this substring is invalid. [SFX] sfx(:,:) = 0._wp ! salt flux; zero unless ice is present (computed in limsbc(_2).F90) On ajoute sfx qui est déclaré dans sbc_oce.F90 [hourdinc@curie51:/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/bdy/compile/nemo/MY_SRC] diff /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90 . 64a70 > REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: sfx , sfx_b !: salt flux [PSU/m2/s] 110a117 > & sfx (jpi,jpj) , & ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcblk_core.f90(405): error #6404: This name does not have a type, and must have an explicit type. [CPIC] & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic On ajoute cpic dans le phycst.F90 de little_NEMO_r85 aui est dans MY_SRC pour le cas ou ni defined key_lim3 ni defined key_cice [hourdinc@curie51:/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/bdy/compile/nemo/MY_SRC] diff /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90 . 72a74 > REAL(wp), PUBLIC :: cpic !: = rcpic / rhoic (specific heat for ice) [J/Kg/K] 168a173 > cpic = rcpic / rhoic ! specific heat for ice [J/Kg/K] ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/sbcmod.f90(312): error #6404: This name does not have a type, and must have an explicit type. [NCLOSEA] IF( nclosea == 1 ) CALL sbc_clo( kt ) ! treatment of closed sea in the model domain On ajoute dans MY_SRC/domain.F90: nclosea = nn_closea et dans MY_SRC/dom_oce.F90: INTEGER , PUBLIC :: nn_closea = 0 !: =0 suppress closed sea/lake from the ORCA domain or not (=1) ------------------------------- mpif90 -o nemo.exe /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/obj/nemo.o -L/ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/lib -l__fcm__nemo -L/usr/local/netcdf-4.2_hdf5_parallel/lib -L/usr/local/hdf5-1.8.9_parallel/lib -lnetcdff -lnetcdf -lhdf5_hl -lhdf5 -lhdf5 -lz -lcurl -lstdc++ -L /ccc/work/cont005/ra0542/hourdinc/sources/XIOS_r416/lib -lxios /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/lib/lib__fcm__nemo.a(step.o): In function `step_mp_stp_': /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/step.f90:(.text+0xc32): undefined reference to `asm_bkg_wri_' step.F90 appelle asm_bkg_wri qui est définit dans ASM/asmtrj.F90 dans la version little_NEMO_r85 et dans ASM/asmbkg.F90 dans la version NEMO_r4018. !! *** MODULE asmtrj -> asmbkg *** !! Assimilation trajectory interface: Write to file the background state and the model state trajectory On recupere tous le repertoire ASM de NEMO_r4018 dans little_NEMO_r85 (inutile pour nous) rm -rf ${little_NEMO_r85_BDY}/ASM ; cp -r ${NEMO_r4018}/ASM ${little_NEMO_r85_BDY} du coup on peut détruire MY_SRC/asmbkg.F90 rm -f ${COMPDIR}/asmbkg.F90 ------------------------------- Toujours le même problème: /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/step.f90:(.text+0xc32): undefined reference to `asm_bkg_wri_' cp ${little_NEMO_r85_BDY}/step_oce.F90 ${COMPDIR} On ajoute la ligne USE asmbkg pour que step.F90 appelle step_oce.F90 qui appelle la nouvelle routine ASM asmbkg.F90 qui elle contient la définition de asm_bkg_wri ------------------------------- a la compilation avec -g et -traceback: ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/step.f90:209: undefined reference to `bdy_tra_dmp_' /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/step.f90:251: undefined reference to `bdy_dyn3d_dmp_' On ajoute les lignes USE bdytra ! bdy cond. for tracers (bdy_tra routine) USE bdydyn3d ! bdy cond. for baroclinic vel. (bdy_dyn3d routine) pour que step.F90 appelle step_oce.F90 qui appelle les routines BDY bdytra.F90 et bdydyn3d.F90 qui contiennent les définitions de bdy_tra_dmp et bdy_dyn3d_dmp ------------------------------- === : E R R O R =========== iom_open ~~~ File ./snow_y2000.nc* not found on remet dans la version sbcblk_core.F90 de Julien les lignes de la version de little_NEMO_r85 qui permettent de ne pas prendre en compte les précips solides (neige et grêle) dans les forçages. Seb avait du virer tout ce qui était de l'ordre des precips solides contre l'avis de maître Gurvan. Mais à priori sans trop d'effet en forcé sur les tropiques. ! do we use snow information? jfld = jfld - COUNT( (/ nn_ice == 0 /) ) ------------------------------- Mais plante plus loin: "segmentation fault" ligne 403 de sbcblk_core.f90 dans la routine blk_oce_core car utilise snow dans le calcul de qns forrtl: severe (174): SIGSEGV, segmentation fault occurred Image PC Routine Line Source nemo.exe 0000000000AAB75A sbcblk_core_mp_bl 403 sbcblk_core.f90 qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) & ! Downward Non Solar flux & - sf(jp_snow)%fnow(:,:,1) * rn_pfac * lfus & ! remove latent melting heat for solid precip & - zevap(:,:) * pst(:,:) * rcp & ! remove evap heat content at SST & + ( sf(jp_prec)%fnow(:,:,1) - sf(jp_snow)%fnow(:,:,1) ) * rn_pfac & ! add liquid precip heat content at Tair & * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp & & + sf(jp_snow)%fnow(:,:,1) * rn_pfac & ! add solid precip heat content at min(Tair,Tsnow) & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic ! Chris: On enleve tout ce qui fait appel a snow (precip solide neige ou grele) qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) & ! Downward Non Solar flux & - zevap(:,:) * pst(:,:) * rcp & ! remove evap heat content at SST & + sf(jp_prec)%fnow(:,:,1) * rn_pfac & ! add liquid precip heat content at Tair & * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp réécriture dans la 3.5 du transport du flux de chaleur lié à evaporation/condensation. dans la 3.4: qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) ! Downward Non Solar flux Il faut donc prendre le pendant de cette formulation dans trasbc.F90 cp ${NEMO_r4018}/TRA/trasbc.F90 ${COMPDIR} ------------------------------- /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/trasbc.f90(198): error #6404: This name does not have a type, and must have an explicit type. [R1_RAU0_RCP] sbc_tsc(ji,jj,jp_tem) = r1_rau0_rcp * qns(ji,jj) ! non solar heat flux On ajoute dans MY_SRC/phycst.F90: REAL(wp), PUBLIC :: r1_rau0_rcp !: = 1. / ( rau0 * rcp ) ------------------------------- > Error [CObjectFactory::GetObject(const StdString & id)] : In file '/ccc/work/cont005/ra0542/hourdinc/sources/XIOS/src/object_factory_impl.hpp', line 79 -> [ id = emp_x_sst, U = field ] object is not referenced ! on commente ! CALL iom_put( "emp_x_sst", emp (:,:) * tsn(:,:,1,jp_tem) ) ! c/d term on sst ! CALL iom_put( "emp_x_sss", emp (:,:) * tsn(:,:,1,jp_sal) ) ! c/d term on sss ------------------------------- > Error [CObjectFactory::GetObject(const StdString & id)] : In file '/ccc/work/cont005/ra0542/hourdinc/sources/XIOS/src/object_factory_impl.hpp', line 79 -> [ id = qhc_oce, U = field ] object is not referenced ! On ajoute dans field_def.xml le nouveau champ qhc_oce lié à la nouvelle formulation <field id="qhc_oce" long_name="Downward Heat Content of E-P over open ocean" unit="W/m2" /> ------------------------------- > Error [CObjectFactory::GetObject(const StdString & id)] : In file '/ccc/work/cont005/ra0542/hourdinc/sources/XIOS/src/object_factory_impl.hpp', line 79 -> [ id = tem_ldfubs, U = field ] object is not referenced ! definitions de ces champs introuvables dans tous les exemples de field_def.xml des différentes configs => on commente les CALL iom_put (sans doute des termes de la version de développement de J. Chanut) traadv_ubs.F90: zwx zwy traadv_ubs.F90: tem_ldfubs trdmod.F90: tem_ldf tem_zdf tem_qsr tem_qns tem_xad tem_yad tem_zad ------------------------------- FINALEMENT on teste en virant sbcblk_core.F90 et trasbc.F90 pour rester avec l'ancienne formulation sur le flux de chaleur lié à emp. ------------------------------- Compilation avec key_bdy: /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/bdydta.f90(241): error #6627: This is an actual argument keyword name, and not a dummy argument name. [KIT] & kit=jit, kt_offset=time_offset ) /ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/CONFIG/PERU12/BLD/ppsrc/nemo/bdydta.f90(241): error #6627: This is an actual argument keyword name, and not a dummy argument name. [KT_OFFSET] & kit=jit, kt_offset=time_offset ) cp ${NEMO_r4018}/SBC/fldread.F90 ${COMPDIR} ------------------------------- JULIEN_MY_SRC="/ccc/work/cont005/ra0542/hourdinc/sources/JULIEN_BDY/MY_SRC" NEMO_r4018="/ccc/work/cont005/ra0542/hourdinc/sources/NEMO_r4018_base_JULIEN/NEMO/OPA_SRC" little_NEMO_r85_BDY="/ccc/work/cont005/ra0542/hourdinc/sources/little_NEMO_r85_BDY/NEMOGCM/NEMO/OPA_SRC" COMPDIR="/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/bdy/compile/nemo/MY_SRC" ------------------------------- TODO LISTE ------------------------------- AJOUTER ln_icebergs et nn_ice_embd dans la namelist PENSER a virer le laplacien dynamique de la sponge!!! dans la nameliste: nous: ln_sssr_bnd = .false. ! flag to bound erp term (associated with nn_sssr=2) julien: ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) faut il modifier : >nam_tide ! tide parameters (#ifdef key_tide) comparer les namelistes julien et la notre pour réglage des différents paramètres. Que fait on avec tide? Modifs t2 ou t10 h2 ou h10?!? dans le changement de version de sbcblk_core.F90?

On vire les routines de la sponge pour OBC (dom_oce.F90 dommsk.F90 domwri.F90 dynldf_lap.F90 traldf_lap.F90)

Objet: Run test de 5 jours avec bdy, mais seulement le fichier de coordonnées des points frontières comme info
Pas de prescription de T,S,U,V,ssh aux points frontières. Bdy lit l'état de démarrage pour imposer T et S aux frontières et mets U, V à zero.

================================================================================ ===== diff bdy / bdy2 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy/scripts/getfile_oce.sh bdy2/scripts/getfile_oce.sh -------------------------------------------------------------------------------- < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_U_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_V_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_S_y${YEAR_BEGIN_JOB}.nc . ... > ${io_getfile} ${INPUTDIR}/${BDY}/bdy_coordinates_PERU12_from_Stephane.nc bdy_coordinates.nc -------------------------------------------------------------------------------- diff -b -w bdy/scripts/header.CURIE bdy2/scripts/header.CURIE -------------------------------------------------------------------------------- 78c78 < export OBC="obc/ORCA025_LIM-T323/bathy_Stephane" --- > export BDY="bdy/ORCA025_LIM-T323/bathy_Stephane" -------------------------------------------------------------------------------- diff -b -w bdy/compile/nemo/cpp_PERU12.fcm bdy2/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=75 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput key_diahth --- > bld::tool::fppkeys key_peru12 key_depth=75 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_bdy key_iomput key_diahth -------------------------------------------------------------------------------- diff -b -w -r -x .svn bdy/compile/nemo/MY_SRC bdy2/compile/nemo/MY_SRC -------------------------------------------------------------------------------- Only in bdy2/compile/nemo/MY_SRC: fldread.F90 (pour pouvoir faire des sortiers au pdt)

Objet: Prescription de T,S,U,V,ssh aux frontières => erai74c en bdy

================================================================================ ===== diff bdy2 / bdy3 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy2/param_inputs/namelist.base.oce bdy3/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 454c454 < nn_dyn2d_dta = 0 ! = 0, bdy data are equal to the initial state --- > nn_dyn2d_dta = 1 ! = 0, bdy data are equal to the initial state 459c459 < nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state --- > nn_dyn3d_dta = 1 ! = 0, bdy data are equal to the initial state 462c462 < nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state --- > nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state 473,480c473,480 < bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true. , .false. , 'daily' , '' , '' < bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true. , .false. , 'daily' , '' , '' < bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true. , .false. , 'daily' , '' , '' < bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true. , .false. , 'daily' , '' , '' < bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true. , .false. , 'daily' , '' , '' < bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' < bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' < cn_dir = 'bdydta/' --- > bn_ssh = 'bdyT_u2d_PERU12' , 72 , 'sossheig' , .true. , .false. , 'yearly' , '' , '' > bn_u2d = 'bdyU_u2d_PERU12' , 72 , 'vobtcrtx' , .true. , .false. , 'yearly' , '' , '' > bn_v2d = 'bdyV_u2d_PERU12' , 72 , 'vobtcrty' , .true. , .false. , 'yearly' , '' , '' > bn_u3d = 'bdyU_u3d_PERU12' , 72 , 'vozocrtx' , .true. , .false. , 'yearly' , '' , '' > bn_v3d = 'bdyV_u3d_PERU12' , 72 , 'vomecrty' , .true. , .false. , 'yearly' , '' , '' > bn_tem = 'bdyT_tra_PERU12' , 72 , 'votemper' , .true. , .false. , 'yearly' , '' , '' > bn_sal = 'bdyT_tra_PERU12' , 72 , 'vosaline' , .true. , .false. , 'yearly' , '' , '' > cn_dir = '' -------------------------------------------------------------------------------- diff -b -w bdy2/scripts/getfile_oce.sh bdy3/scripts/getfile_oce.sh -------------------------------------------------------------------------------- > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOB}/bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . > > if [ ${YEAR_BEGIN_JOBp1} -ne ${YEAR_BEGIN_JOB} ] ; then > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyT_u2d_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyT_tra_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyU_u2d_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyU_u3d_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyV_u2d_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIR}/${BDY}/${YEAR_BEGIN_JOBp1}/bdyV_u3d_PERU12_y${YEAR_BEGIN_JOBp1}.nc . > fi

Objet: Test de la clé de correction de flux aux frontières

================================================================================ ===== diff bdy3 / bdy3b ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy3/param_inputs/namelist.base.oce bdy3b/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 465c465 < ln_vol = .false. ! total volume correction (see nn_volctl parameter) --- > ln_vol = .true. ! total volume correction (see nn_volctl parameter)

Objet: bdy avec frontières extraites de trop12

• rn_usecrt=1 !!!: valeur par défaut parceque parametre oublié dans la nameliste => prise en compte de l'understress
• Etat de demarrage : T et S lu dans le 5 janvier 2000 du run trop 12 (u et v a zero mais mezo déjà présente via T et S)
• Changements entre namelist peru12 et namelist trop12
• Short Wave toutes les 6 h au lieu d'être interpolé depuis une moyenne journalière des 6 heures
• 2 bands au lieu de RGB
• + autres changements dans des paramètres (solver...)
• bathy trop12 au lieu de Stephane

================================================================================ ===== diff bdy3 / bdy4 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy3/param_inputs/namelist.base.oce bdy4/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 41a42 > 100,101c101,103 < sn_tem = 'dta_temp' , -1 , 'votemper' , .true. , .true. , 'yearly' , ' ' , ' ' < sn_sal = 'dta_sal' , -1 , 'vosaline' , .true. , .true. , 'yearly' , '' , ' ' --- > sn_tem = 'dta_temp' , -1 , 'thetao' , .false. , .true. , 'yearly' , ' ' , ' ' > sn_sal = 'dta_sal' , -1 , 'so' , .false. , .true. , 'yearly' , '' , ' ' > ! 139c141 < ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave --- > ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 193c195 < sn_qsr = 'radsw' , 24 , 'ssrd' , .false. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' --- > sn_qsr = 'radsw' , 6 , 'ssrd' , .true. , .false. , 'yearly' ,'reshape_ERAI_PERU12_bilinear.nc', '' 202d203 < rn_usecrt = 0. ! use surface currents to compute relative wind speed 252,253c253,254 < ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration < ln_qsr_2bd = .false. ! 2 bands light penetration --- > ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration > ln_qsr_2bd = .true. ! 2 bands light penetration 266c267 < Sn_cnf = 'runoff_1m_nomask' , 0 , 'socoefr' , .false. , .true. , 'yearly' , '' , '' --- > sn_cnf = 'runoff_1m_nomask' , 0 , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 307c308 < ln_sssr_bnd = .false. ! flag to bound erp term (associated with nn_sssr=2) --- > ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 421c422 < nn_cln_update = 3 ! baroclinic update frequency --- > nn_cln_update = 1 ! baroclinic update frequency 465c466 < ln_vol = .false. ! total volume correction (see nn_volctl parameter) --- > ln_vol = .true. ! total volume correction (see nn_volctl parameter) 473,479c474,480 < bn_ssh = 'bdyT_u2d_PERU12' , 72 , 'sossheig' , .true. , .false. , 'yearly' , '' , '' < bn_u2d = 'bdyU_u2d_PERU12' , 72 , 'vobtcrtx' , .true. , .false. , 'yearly' , '' , '' < bn_v2d = 'bdyV_u2d_PERU12' , 72 , 'vobtcrty' , .true. , .false. , 'yearly' , '' , '' < bn_u3d = 'bdyU_u3d_PERU12' , 72 , 'vozocrtx' , .true. , .false. , 'yearly' , '' , '' < bn_v3d = 'bdyV_u3d_PERU12' , 72 , 'vomecrty' , .true. , .false. , 'yearly' , '' , '' < bn_tem = 'bdyT_tra_PERU12' , 72 , 'votemper' , .true. , .false. , 'yearly' , '' , '' < bn_sal = 'bdyT_tra_PERU12' , 72 , 'vosaline' , .true. , .false. , 'yearly' , '' , '' --- > bn_ssh = 'bdyT_u2d_PERU12', 120 , 'zos' , .true. , .false. , 'yearly' , '' , '' > bn_u2d = 'bdyU_u2d_PERU12', 120 , 'uobtrope' , .true. , .false. , 'yearly' , '' , '' > bn_v2d = 'bdyV_u2d_PERU12', 120 , 'vobtrope' , .true. , .false. , 'yearly' , '' , '' > bn_u3d = 'bdyU_u3d_PERU12', 120 , 'uo' , .true. , .false. , 'yearly' , '' , '' > bn_v3d = 'bdyV_u3d_PERU12', 120 , 'vo' , .true. , .false. , 'yearly' , '' , '' > bn_tem = 'bdyT_tra_PERU12', 120 , 'thetao' , .true. , .false. , 'yearly' , '' , '' > bn_sal = 'bdyT_tra_PERU12', 120 , 'so' , .true. , .false. , 'yearly' , '' , '' 708c709 < nn_npc = 0 ! frequency of application of npc --- > nn_npc = 1 ! frequency of application of npc 731c732 < rn_ebb = 60 ! coef. of the surface input of tke (67.83 now usual value) --- > rn_ebb = 67.83 ! coef. of the surface input of tke (67.83 now usual value) 816c817 < nn_solv = 2 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) --- > nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 820c821 < nn_nmin = 200 ! minimum of iterations for the SOR solver --- > nn_nmin = 300 ! minimum of iterations for the SOR solver 824c825 < rn_sor = 1.973 ! optimal coefficient for SOR solver (to be adjusted with the domain) --- > rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 907c908 < ln_ptrcomp = .false. ! Add decomposition : overturning --- > ln_ptrcomp = .true. ! Add decomposition : overturning -------------------------------------------------------------------------------- diff -b -w bdy3/scripts/getfile_oce.sh bdy4/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 5c5 < ${io_getfile} ${INPUTDIR}/Stephane/PERU12-I/bathymetry_PERU12_mod2.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean.nc bathy_meter.nc 22,23d21 < # ERAI (freq 1 day) : radsw < ${io_getfile} ${DATA}/erai_peru_mask/freq_1_day/${YEAR_BEGIN_JOB}/ssrd.${YEAR_BEGIN_JOB}.fshei.PERU_075.nc radsw_y${YEAR_BEGIN_JOB}.nc 25c23,24 < # ERAI (freq 6 hours) : radlw, q2, t2, precip (freq 6 hours) --- > # ERAI (freq 6 hours) : radsw, radlw, q2, t2, precip (freq 6 hours) > ${io_getfile} ${DATA}/erai_peru_mask/${YEAR_BEGIN_JOB}/ssrd.${YEAR_BEGIN_JOB}.fshei.PERU_075.nc radsw_y${YEAR_BEGIN_JOB}.nc 50c49 < ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/bathy_Stephane/ORCA025_LIM-T323_clim_2000-2008_gridT_votemper_Perou.nc dta_temp.nc --- > ${io_getfile} ${INPUTDIR}/ini/tr12_quik/tr12_quik_5d_20000105_grid_T_3D.nc dta_temp.nc 53c52 < ${io_getfile} ${INPUTDIR}/ini/mercator/ORCA025_LIM-T323/bathy_Stephane/ORCA025_LIM-T323_clim_2000-2008_gridS_vosaline_Perou.nc dta_sal.nc --- > ${io_getfile} ${INPUTDIR}/ini/tr12_quik/tr12_quik_5d_20000105_grid_T_3D.nc dta_sal.nc 59d57 < #${io_getfile} ${INPUTDIR}/${BDY}/bdy_coordinates_PERU12_from_Stephane.nc bdy_coordinates.nc

Objet: volume correction=FALSE

================================================================================ ===== diff bdy4 / bdy4b ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy4/param_inputs/namelist.base.oce bdy4b/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 466c466 < ln_vol = .true. ! total volume correction (see nn_volctl parameter) --- > ln_vol = .false. ! total volume correction (see nn_volctl parameter)

Objet: Scripts à zéro

Objet: Mise à plat des scripts

Objet: obc au lieu de bdy

Attention: little_nemo r73 au lieu de r85_BDY, mais à priori sans impact parce qu'il y a toujours dans MY_SRC
dynadv_ubs.F90 pour la correction gamma 1/32 pour ubs et
sbcblk_core.F90 pour la non prise en compte de l'understress
(cf "diff erai74c/bdy" pour le passage de little_nemo r73 a r85_BDY ne changeait rien sur un run OBC

================================================================================ ===== diff bdy4 / erai75c ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy4/scripts/module.sh erai75c/scripts/module.sh -------------------------------------------------------------------------------- 8c8 < export NEMOGCM="${WORKDIR}/sources/little_NEMO_r85_BDY/NEMOGCM" --- > export NEMOGCM="${WORKDIR}/sources/little_NEMO/NEMOGCM"

Objet: Correction bug free-slip (Nicolas Ducousso)

• on remet le module.sh qui ne fonctionnait plus depuis la maintenance du 25 Mars 2014. Pour que les commandes modules fonctionnent dans les scripts, on remplace sh par bash.
• Nouvelle maintenance : 17 Juin 2014 . Même en remettant l'ancien environnement avec "module purge; module load bullxscs/4r2u2" on ne parvient pas à recréer le même executable nemo de bdy4b au bit près.
  En revanche, je refais tourner le run bdy4b recompilé sur 1 jours, et les 216 pdt du fichiers solver.stat sont identiques.
• Après cette vérif, on recompile nemo avec la correction de bug de Nicolas Ducousso et on relance le même run que bdy4b.
  http://forge.ipsl.jussieu.fr/nemo/ticket/1339
  The bugfix allows to recover a free-slip character for UBS scheme close to the boundary as illustrated in the joint movie cmp_evol_C4_vs_modified_UBS.mpeg.
  n relance le même run que bdy4b.

================================================================================ ===== diff bdy4b / bdy4c ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w bdy4b/scripts/module.sh bdy4c/scripts/module.sh -------------------------------------------------------------------------------- 1c1 < #!/bin/sh --- > #!/bin/bash 9a10 > 11,19c12 > module purge; module load bullxscs/4r2u2 20a14,16 > module load netcdf/4.2_hdf5_parallel > module load hdf5/1.8.9_parallel -------------------------------------------------------------------------------- diff -b -w -r -x .svn bdy4b/compile/nemo/MY_SRC bdy4c/compile/nemo/MY_SRC -------------------------------------------------------------------------------- Only in bdy4c/compile/nemo/MY_SRC: dynadv_ubs.F90 -------------------------------------------------------------------------------- cmp bdy4b/compile/nemo/nemo.exe bdy4c/compile/nemo/nemo.exe -------------------------------------------------------------------------------- bdy4b/compile/nemo/nemo.exe bdy4c/compile/nemo/nemo.exe differ: char 8, line 1

Objet: Correction bug dynkeg et dynvor (Nicolas Ducousso)

================================================================================ ===== diff bdy4c / bdy4d ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w -r -x .svn bdy4c/compile/nemo/MY_SRC bdy4d/compile/nemo/MY_SRC -------------------------------------------------------------------------------- Only in bdy4d/compile/nemo/MY_SRC: dynkeg.F90 Only in bdy4d/compile/nemo/MY_SRC: dynvor.F90

Objet: Test efficacité bdy (=> tout comme tr12_quik)

Réponse:Sur les diags actuels, bdy semble bien faire son travail... diag

• Version little_NEMO : tr12_quik r85 (donc sans les corrections de Nicolas Ducousso sur dynadv_ubs, dynkeg et dynvor)
• Version little_NEMO : on ramène toutes les routines différentes de little_NEMO_r85_BDY dans MY_SRC => retour à little_NEMO_r85
• Version little_NEMO : test sur 2 x 1 jours des 2 versions. Même solver.stat
• Version little_NEMO : r86: commit Print de rn_usecrt dans ocean.output : SBC/sbcblk_core.F90
• cpp : on vire key_dtatem key_dtasal obsolete de cpp_PERU12
• cpp : on vire key_traldf_c2d key_dynldf_c2d car on est en UBS
• cpp : on ne prend pas dans cpp_TROP12 key_ldfslp (car UBS) ni key_zdftmx (tide avec mask pour indonésie uniquement) => pas besoin des fichiers K1 M2...
• Forçages: on enlève le cycle diurne pour coller à trop12
• Forçages: on prend les fichiers erai et quikscat de Seb
• Forçages: on prend les fichiers de poids de PERU12 (bicubic pour les vents et bilin pour le reste comme dans trop12)
• 2 bands au lieu de RGB
• runoffs: on prend les fichiers de Seb
• Chlorophyll: on prend les fichiers de Seb (mais pas utilisé avec le 2 bands)
• bathy: TROP12_NEST_PERU12 + extraction de x=1-361 sur 382 + on bouche la mer des caraibes (idem a l'extraction de PERU12 sur la bathy de trop12)
• coordinates: TROP12 extraction de PERU12 (idem au coordinates de Stephane)
• le fichier restart de tr12_quik (tr12_eranocrt 19991231 n'existe plus) => T et S de tr12_quik_5d_20000105_grid_T_3D.nc
• Ajout des sorties mensuelles

Objet: 3 bug fix Nicolas Ducousso + 1 bug fix GM (no fmask)

Réponse: Tous ces bugs fix n'affectent pas beaucoup notre solution...

• Version little_NEMO r86 +
• Bug fix free slip de Nicolas Ducousso sur dynadv_ubs.F90 http://forge.ipsl.jussieu.fr/nemo/ticket/1339
  The bugfix allows to recover a free-slip character for UBS scheme close to the boundary as illustrated in the joint movie cmp_evol_C4_vs_modified_UBS.mpeg.
• Bug fix de Nicolas Ducousso sur dynkeg.F90
• Bug fix de Nicolas Ducousso sur dynvor.F90
• Bug fix de Gurvan Madec sur dynvor.F90 (no fmask)

================================================================================ ===== diff p12n_bdy_tr12n_quik / p12n_bdy_tr12n_quik_ND ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w -r -x .svn p12n_bdy_tr12n_quik/compile/nemo/MY_SRC p12n_bdy_tr12n_quik_ND/compile/nemo/MY_SRC -------------------------------------------------------------------------------- Only in p12n_bdy_tr12n_quik_ND/compile/nemo/MY_SRC: dynadv_ubs.F90 Only in p12n_bdy_tr12n_quik_ND/compile/nemo/MY_SRC: dynkeg.F90 Only in p12n_bdy_tr12n_quik_ND/compile/nemo/MY_SRC: dynvor.F90

Objet: Rendre plus visqueux avec Vector Form + bilap sur la dyn au lieu de UBS pour diminuer la var_ssh

Réponse: En rendant le milieu plus visqueux (surtout là où les courants sont faibles), avec vector from + bilap au lieu de UBS, on réduit bien la var de ssh, et on refroidit l'upwelling le long des côtes Chiliennes:

• Version little_NEMO r86 +

-------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik_ND/param_inputs/namelist.base.oce p12n_bdy_tr12n_quik_blp1/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 637c637 < ln_dynadv_vec = .false. ! vector form (T) or flux form (F) --- > ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 639c639 < ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme --- > ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 672c672 < ln_dynldf_bilap = .false. ! bilaplacian operator --- > ln_dynldf_bilap = .true. ! bilaplacian operator 680c680 < rn_ahm_0_blp = -1.8e11 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -2.5e10 ! horizontal bilaplacian eddy viscosity [m4/s]

Objet: UBS traceurs => fct + msc + ln_traldf_blp + nn_aht_ijk_t=31

Réponse:

• run tiré de tr025n_fm02: tr025n_fm00 + fct + msc + ln_traldf_blp + nn_aht_ijk_t=31
• nouveau schéma d'advection sur les traceurs: fct au lieu de tvd

Objet: sensibilité aux fontières (diff entre bdy tr12_quik et bdy soda2)

Réponse:les meilleurs resultats du run bdy tr12_quik par rapport au run bdy soda2 peut s'expliquer par le fait que les bdy tr12_quik sont générées par le même modèle que le modèle utilisé. Il y a une cohérence entre les 2. Ca peut inciter à faire du AGRIF plutôt que du BDY.

================================================================================ ===== diff p12n_bdy_tr12n_quik/ / p12n_bdy_soda2/ ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik//param_inputs/namelist.base.oce p12n_bdy_soda2//param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 471,477c471,477 < bn_ssh = 'bdyT_u2d_PERU12', 120 , 'zos' , .true. , .false. , 'yearly' , '' , '' < bn_u2d = 'bdyU_u2d_PERU12', 120 , 'uobtrope' , .true. , .false. , 'yearly' , '' , '' < bn_v2d = 'bdyV_u2d_PERU12', 120 , 'vobtrope' , .true. , .false. , 'yearly' , '' , '' < bn_u3d = 'bdyU_u3d_PERU12', 120 , 'uo' , .true. , .false. , 'yearly' , '' , '' < bn_v3d = 'bdyV_u3d_PERU12', 120 , 'vo' , .true. , .false. , 'yearly' , '' , '' < bn_tem = 'bdyT_tra_PERU12', 120 , 'thetao' , .true. , .false. , 'yearly' , '' , '' < bn_sal = 'bdyT_tra_PERU12', 120 , 'so' , .true. , .false. , 'yearly' , '' , '' --- > bn_ssh = 'bdyT_u2d_PERU12', 120 , 'SSH' , .true. , .false. , 'yearly' , '' , '' > bn_u2d = 'bdyU_u2d_PERU12', 120 , 'BTU' , .true. , .false. , 'yearly' , '' , '' > bn_v2d = 'bdyV_u2d_PERU12', 120 , 'BTV' , .true. , .false. , 'yearly' , '' , '' > bn_u3d = 'bdyU_u3d_PERU12', 120 , 'U' , .true. , .false. , 'yearly' , '' , '' > bn_v3d = 'bdyV_u3d_PERU12', 120 , 'V' , .true. , .false. , 'yearly' , '' , '' > bn_tem = 'bdyT_tra_PERU12', 120 , 'TEMP' , .true. , .false. , 'yearly' , '' , '' > bn_sal = 'bdyT_tra_PERU12', 120 , 'SALT' , .true. , .false. , 'yearly' , '' , '' -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik//scripts/getfile_oce.sh p12n_bdy_soda2//scripts/getfile_oce.sh -------------------------------------------------------------------------------- 78,83c78,83 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc 87,92c87,92 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_y${year_begin_job_m1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_y${year_begin_job_m1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_y${year_begin_job_m1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_y${year_begin_job_m1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_y${year_begin_job_m1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_y${year_begin_job_m1}.nc . --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${year_begin_job_m1}.nc bdyT_u2d_PERU12_y${year_begin_job_m1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda2_y${year_begin_job_m1}.nc bdyT_tra_PERU12_y${year_begin_job_m1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda2_y${year_begin_job_m1}.nc bdyU_u2d_PERU12_y${year_begin_job_m1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda2_y${year_begin_job_m1}.nc bdyU_u3d_PERU12_y${year_begin_job_m1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda2_y${year_begin_job_m1}.nc bdyV_2d_PERU12_y${year_begin_job_m1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda2_y${year_begin_job_m1}.nc bdyV_3d_PERU12_y${year_begin_job_m1}.nc 97,102c97,102 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_y${year_begin_job_p1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_y${year_begin_job_p1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_y${year_begin_job_p1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_y${year_begin_job_p1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_y${year_begin_job_p1}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_y${year_begin_job_p1}.nc . --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${year_begin_job_p1}.nc bdyT_u2d_PERU12_y${year_begin_job_p1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda2_y${year_begin_job_p1}.nc bdyT_tra_PERU12_y${year_begin_job_p1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda2_y${year_begin_job_p1}.nc bdyU_u2d_PERU12_y${year_begin_job_p1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda2_y${year_begin_job_p1}.nc bdyU_u3d_PERU12_y${year_begin_job_p1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda2_y${year_begin_job_p1}.nc bdyV_u2d_PERU12_y${year_begin_job_p1}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda2_y${year_begin_job_p1}.nc bdyV_u3d_PERU12_y${year_begin_job_p1}.nc -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik//scripts/CURIE.env p12n_bdy_soda2//scripts/CURIE.env -------------------------------------------------------------------------------- 14c14 < export BDY="bdy/tr12_quik/from_obc" --- > export BDY="bdy/soda2"

Objet: sensibilité aux conditions intiales (diff entre demarrage T et S de tr12_quik_5d_20000105_grid_T_3D.nc et T et S de soda2_TS_20000102.nc)

Réponse:

• Diff Etat de demarrage Soda2 et tr12_quik en Janvier 2000

-------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2/scripts/getfile_oce.sh p12n_bdy_soda2_B/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 67c67 < ${io_getfile} ${INPUTDIR}/ini/tr12_quik/tr12_quik_5d_20000105_grid_T_3D.nc dta_temp.nc --- > ${io_getfile} ${INPUTDIR}/ini/soda2/soda2_TS_20000102.nc dta_temp.nc 70c70 < ${io_getfile} ${INPUTDIR}/ini/tr12_quik/tr12_quik_5d_20000105_grid_T_3D.nc dta_sal.nc --- > ${io_getfile} ${INPUTDIR}/ini/soda2/soda2_TS_20000102.nc dta_sal.nc

Objet: Responsabilité des partial steps sue les pb d'EKE trop fort à la côte

Réponse:

Historique:

Pour passer de partial steps en steps, il faut fournir un fichier bathy en level et non en meter.
On relit donc le mbathy du mesh_mask de p12n_bdy_soda2_B qu'on écrit dans un fichier bathy_level.nc

================================================================================ ===== diff p12n_bdy_soda2_B / p12n_bdy_soda2_B_steps ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/param_inputs/namelist.base.oce p12n_bdy_soda2_B_steps/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 61,62c61,62 < ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) < ln_zps = .true. ! z-coordinate - partial steps (T/F) --- > ln_zco = .true. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) > ln_zps = .false. ! z-coordinate - partial steps (T/F) 652,653c652,653 < ln_hpg_zco = .false. ! z-coordinate - full steps < ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) --- > ln_hpg_zco = .true. ! z-coordinate - full steps > ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getfile_oce.sh p12n_bdy_soda2_B_steps/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 5c5 < ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean_361x505.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean_361x505_level.nc bathy_level.nc

Objet: Responsabilité de la pente sur les pb d'EKE trop fort à la côte

Réponse:

Historique:

la bathy de François est seuillée à -25m. On multiplie donc cette bathy (ETOPO1_NESTPERU_R12_mask_Caribbean_361x505_rfact.nc)
par l'ancien mesh_mask (mesh_mask_peru12.L75_from_trop12.nc) pour recréer une bathy lissée avec exactement la même côte
que l'expérience précédente (ETOPO1_NESTPERU_R12_mask_Caribbean_361x505_rfact_mask.nc)

================================================================================ ===== diff p12n_bdy_soda2_B / p12n_bdy_soda2_B_rfact ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getfile_oce.sh p12n_bdy_soda2_B_rfact/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 5c5 < ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean_361x505.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIR}/grd/ETOPO1_NESTPERU_R12_mask_Caribbean_361x505_rfact_mask.nc bathy_meter.nc

• little_NEMO : Corrections de Nicolas Ducousso sur dynadv_ubs, dynkeg et dynvor???
• Forçages : Remettre le cycle diurne? (pour trop12 aussi!)
• Forçages : ln_dm2dc = .true. + mettre a 24 heure au lieu de 6 h dans la namelist + pas d'interpolation (juste une valeur journalière)
• Forçages : les erai de Stephane sont différents de Seb. diff Seb/Steph
• Retour à RGB (au lieu du 2 bands)? Quel fichier de chlorophyll?
• Runoff : François est pour les virer
• Quel etat de demarrage? T et S de tr12_quik_5d_20000105_grid_T_3D.nc?

Objet: Impact du refroidissement de 1° dans le coeur de l'EUC

Différences:

• Modif des frontières pour réduire la température de 1 degré dans le coeur de l'EUC
• Différence de démarrage

================================================================================ ===== diff p12n_bdy_soda2_B / p12n_bdy_soda2_B_T-1degre ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getfile_oce.sh p12n_bdy_soda2_B_T-1degre/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 78,83c78,83 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc ... --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_T-1degre_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getrst_oce.sh p12n_bdy_soda2_B_T-1degre/scripts/getrst_oce.sh -------------------------------------------------------------------------------- 13c13 < # for file in `${MACHINE_STOCKAGE} ls /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai70c/19991231/restart_oce_*_????.nc` --- > for file in `${MACHINE_STOCKAGE} ls /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/p12n_bdy_soda2_B/20021231/restart_oce_*_????.nc`

Objet: Impact du refroidissement de 1.5° dans le coeur de l'EUC

Différences:

• Modif des frontières pour réduire la température de 1.5 degré dans le coeur de l'EUC

================================================================================ ===== diff p12n_bdy_soda2_B / p12n_bdy_soda2_B_T-1.5degre ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getfile_oce.sh p12n_bdy_soda2_B_T-1.5degre/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 78,83c78,83 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc ... --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc

Objet: Impact des vents (erai au lieu de Quik)

Différences:

• Uniquement les vents

================================================================================ ===== diff p12n_bdy_tr12n_quik / p12n_bdy_tr12n_quik_wnd_erai ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/param_inputs/namelist.base.oce p12n_bdy_tr12n_quik_wnd_erai/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 193,194c193,194 < sn_wndi = 'quikscat' , 24, 'zonal_wind_speed' , .true. , .false., 'yearly' ,'wbicub.nc', '' < sn_wndj = 'quikscat' , 24, 'meridional_wind_speed', .true. , .false., 'yearly' ,'wbicub.nc', '' --- > sn_wndi = 'u10' , 6 , 'u10' , .true. , .false., 'yearly' ,'wbilin.nc', '' > sn_wndj = 'v10' , 6 , 'v10' , .true. , .false., 'yearly' ,'wbilin.nc', '' -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/scripts/getfile_oce.sh p12n_bdy_tr12n_quik_wnd_erai/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 13,14c13,15 < # QUIKSCAT (frep 1 day): u10, v10 < ${io_getfile} ${QUIKDIR}/quikscat_1d_${YEAR_BEGIN_JOB}_filled.nc quikscat_y${YEAR_BEGIN_JOB}.nc --- > # u10 v10 > ${io_getfile} ${BLKDIR}/u10.${YEAR_BEGIN_JOB}.ashei.GLOBAL_075.nc u10_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${BLKDIR}/v10.${YEAR_BEGIN_JOB}.ashei.GLOBAL_075.nc v10_y${YEAR_BEGIN_JOB}.nc 24c25,26 < ${io_getfile} ${QUIKDIR}/quikscat_1d_${year_begin_job_m1}_filled.nc quikscat_y${year_begin_job_m1}.nc --- > ${io_getfile} ${BLKDIR}/u10.${year_begin_job_m1}.ashei.GLOBAL_075.nc u10_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${BLKDIR}/v10.${year_begin_job_m1}.ashei.GLOBAL_075.nc v10_y${YEAR_BEGIN_JOB}.nc

Objet: Impact de soda3 plutôt que soda2 aux frontières sur nos solutions

Différences:

• Uniquement les frontières et l'état initial

================================================================================ ===== diff p12n_bdy_soda2_B / p12n_bdy_soda3 ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda2_B/scripts/getfile_oce.sh p12n_bdy_soda3/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 67c67 < ${io_getfile} ${INPUTDIR}/ini/soda2/soda2_TS_20000102.nc dta_temp.nc --- > ${io_getfile} ${INPUTDIR}/ini/soda3.3.1_5dy_ocean_or/soda3.3.1_5dy_ocean_or_TS_20000103_bathy_trop12.nc dta_temp.nc 70c70 < ${io_getfile} ${INPUTDIR}/ini/soda2/soda2_TS_20000102.nc dta_sal.nc --- > ${io_getfile} ${INPUTDIR}/ini/soda3.3.1_5dy_ocean_or/soda3.3.1_5dy_ocean_or_TS_20000103_bathy_trop12.nc dta_sal.nc 78,83c78,83 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda2_y${YEAR_BEGIN_JOB}.nc bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc ...

Objet: Impact du refroidissement de 1.5° dans le coeur de l'EUC

Différences:

• Modif des frontières pour réduire la température de 1.5 degré dans le coeur de l'EUC

================================================================================ ===== diff p12n_bdy_soda3/ / p12n_bdy_soda3_T-1.5degre/ ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda3//scripts/CURIE.env p12n_bdy_soda3_T-1.5degre//scripts/CURIE.env -------------------------------------------------------------------------------- 14c14 < export BDY="bdy/soda3.3.1_5dy_ocean_or" --- > export BDY="bdy/soda3.3.1_5dy_ocean_or_T-1.5degre" -------------------------------------------------------------------------------- diff -b -w p12n_bdy_soda3//scripts/getfile_oce.sh p12n_bdy_soda3_T-1.5degre//scripts/getfile_oce.sh -------------------------------------------------------------------------------- 78,83c78,83 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_soda3.3.1_5dy_ocean_or_y${YEAR_BEGIN_JOB}.nc bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc --- > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_soda3.3.1_5dy_ocean_or_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_soda3.3.1_5dy_ocean_or_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_soda3.3.1_5dy_ocean_or_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_soda3.3.1_5dy_ocean_or_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_soda3.3.1_5dy_ocean_or_T-1.5degre_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc ...

Objet: Impact du rechauffement de 1.5° dans le coeur de l'EUC

Différences:

• Modif des frontières pour augmenter la température de 1.5 degré dans le coeur de l'EUC
• Sorties T et S en 3D dans un seul fichier "T_3D" au lieu de 2 fichiers séparés
• modif de l'environnement machine par rapport à l'exp tr12_quik faite il y a 3 ans

================================================================================ ===== diff p12n_bdy_tr12n_quik / p12n_bdy_tr12n_quik_Tp1.5degre ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/param_inputs/iodef.xml p12n_bdy_tr12n_quik_Tp1.5degre/param_inputs/iodef.xml -------------------------------------------------------------------------------- 214c224 < <file id="5d_grid_T_T3D" name_suffix="_T3D" > --- > <file id="5d_grid_T_3D" name_suffix="_3D" > ... -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/scripts/getfile_oce.sh p12n_bdy_tr12n_quik_Tp1.5degre/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 78,83c78,84 < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc . < ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc . --- > > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_u2d_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyT_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyT_tra_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyT_tra_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u2d_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyU_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyU_u3d_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyU_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u2d_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyV_u2d_PERU12_y${YEAR_BEGIN_JOB}.nc > ${io_getfile} ${INPUTDIR}/${BDY}/bdyV_u3d_PERU12_tr12_quik_T+1.5degre_y${YEAR_BEGIN_JOB}.nc bdyV_u3d_PERU12_y${YEAR_BEGIN_JOB}.nc ... -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/scripts/CURIE.env p12n_bdy_tr12n_quik_Tp1.5degre/scripts/CURIE.env -------------------------------------------------------------------------------- 14c14 < export BDY="bdy/tr12_quik/from_obc" --- > export BDY="bdy/tr12_quik_T+1.5degre" ... -------------------------------------------------------------------------------- diff -b -w p12n_bdy_tr12n_quik/scripts/CURIE.header p12n_bdy_tr12n_quik_Tp1.5degre/scripts/CURIE.header -------------------------------------------------------------------------------- < module purge < module load bullxscs/4r2u2 --- > module unload nco cdo ncview netcdf hdf5 ...

================================================================================ ===== diff cpl00 / cpl01 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl00/param_inputs/namelist.base.oce cpl01/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 251,252c251,252 < ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration < ln_qsr_2bd = .true. ! 2 bands light penetration --- > ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration > ln_qsr_2bd = .false. ! 2 bands light penetration -------------------------------------------------------------------------------- diff -b -w cpl00/scripts/getrst_oce.sh cpl01/scripts/getrst_oce.sh -------------------------------------------------------------------------------- 7c7 < ${io_getfile} /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai70a/20061231/restart_oce_erai70a_20061231.nc restart_oce.nc --- > ${io_getfile} /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai70c/20061231/restart_oce_erai70c_20061231.nc restart_oce.nc ================================================================================ ===== ATMOSPHERE ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl00/param_inputs/namelist.base.atm cpl01/param_inputs/namelist.base.atm -------------------------------------------------------------------------------- 82c82 < ra_sw_physics = 2, --- > ra_sw_physics = 1, 84c84 < sf_sfclay_physics = 1, --- > sf_sfclay_physics = 5, 86c86 < bl_pbl_physics = 1, --- > bl_pbl_physics = 5, 91c91 < ifsnow = 1, --- > ifsnow = 0, ================================================================================ ===== COUPLAGE ===== ================================================================================

• Compilation nemo: Déplacement d'un write + ajout du flag -fno-alias
• Only in cpl02/compile/nemo/MY_SRC: dynadv_ubs.F90 (Bug UBS dyn: gamma2 = 1/32 au lieu de 1/8)
• SSS On enlève le rappel en salinité
• OBC Mercator ORCA025_LIM-T323 (3 jours au lieu de 5 jours pour ORCA025-B83
• TKE ln_mxl0 = .true. au lieu de false ! surface mixing length scale = F(wind stress) (T) or not (F)
• RST Demarrage Etat moyen Mercator ORCA025_LIM-T323 au lieu de restart erai70c 2000
• Année 2008 en plus!!
• Bathy trop12 au lieu de Stephane
• OBC year + 1 car problème avec le découpage en 3 jours => va chercher l'année suivante pour les derniers jours...

================================================================================ ===== diff cpl01 / cpl02 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl01/param_inputs/namelist.base.oce cpl02/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 100,101c100,101 < sn_tem = 'dta_temp' , -1 , 'thetao' , .false. , .true. , 'yearly' , ' ' , ' ' < sn_sal = 'dta_sal' , -1 , 'so' , .false. , .true. , 'yearly' , '' , ' ' --- > sn_tem = 'dta_temp' , -1 , 'votemper' , .true. , .true. , 'yearly' , ' ' , ' ' > sn_sal = 'dta_sal' , -1 , 'vosaline' , .true. , .true. , 'yearly' , '' , ' ' 300c300 < nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) --- > nn_sssr = 0 ! add a damping term in the surface freshwater flux (=2) 368,379c368,379 < sn_obcn(1) = 'obc_north_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(2) = 'obc_north_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(3) = 'obc_north_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcn(4) = 'obc_north_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(1) = 'obc_west_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(2) = 'obc_west_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(3) = 'obc_west_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcw(4) = 'obc_west_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(1) = 'obc_south_U', 120 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(2) = 'obc_south_V', 120 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(3) = 'obc_south_TS', 120 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' < sn_obcs(4) = 'obc_south_TS', 120 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' --- > sn_obcn(1) = 'obc_north_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(2) = 'obc_north_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(3) = 'obc_north_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcn(4) = 'obc_north_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(1) = 'obc_west_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(2) = 'obc_west_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(3) = 'obc_west_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcw(4) = 'obc_west_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(1) = 'obc_south_U', 72 , 'vozocrtx' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(2) = 'obc_south_V', 72 , 'vomecrty' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(3) = 'obc_south_T', 72 , 'votemper' , .TRUE. , .FALSE. , 'yearly' , '' , '' > sn_obcs(4) = 'obc_south_S', 72 , 'vosaline' , .TRUE. , .FALSE. , 'yearly' , '' , '' 700,701c700,701 < ln_mxl0 = .false. ! surface mixing length scale = F(wind stress) (T) or not (F) < rn_mxl0 = 0.4 ! surface buoyancy lenght scale minimum value --- > ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) > rn_mxl0 = 0.01 ! surface buoyancy lenght scale minimum value -------------------------------------------------------------------------------- diff -b -w cpl01/scripts/module.sh cpl02/scripts/module.sh -------------------------------------------------------------------------------- 8,9c8,9 < export NEMOGCM="${WORKDIR}/sources/now/little_nemo" < export XIOS="${WORKDIR}/sources/now/xios_cpl" --- > export NEMOGCM="${WORKDIR}/sources/little_NEMO_r73/NEMOGCM" > export XIOS="${WORKDIR}/sources/now/xios_cpl_r430" 11c11 < export OASIS="${WORKDIR}/sources/now/oa3mct/CURIE" --- > export OASIS="${WORKDIR}/sources/now/oa3mct_r682/CURIE" -------------------------------------------------------------------------------- diff -b -w cpl01/scripts/getfile_oce.sh cpl02/scripts/getfile_oce.sh -------------------------------------------------------------------------------- 5,6c5,7 < ${io_getfile} ${INPUTDIRO}/${DATA}/bathymetry_PERU12_mod2.nc bathy_meter.nc --- > ${io_getfile} ${INPUTDIRO}/grd/ETOPO1_NESTPERU_R12_org_mask_Caribbean_Sean.nc bathy_meter.nc 22c25,26 < ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_north_S_y${YEAR_BEGIN_JOB}.nc . 26c30,31 < ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_west_S_y${YEAR_BEGIN_JOB}.nc . 30c35,53 < ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_TS_y${YEAR_BEGIN_JOB}.nc . --- > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_T_y${YEAR_BEGIN_JOB}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOB}/obc_south_S_y${YEAR_BEGIN_JOB}.nc . > > if [ ${YEAR_BEGIN_JOBp1} -ne ${YEAR_BEGIN_JOB} ] ; then > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_north_S_y${YEAR_BEGIN_JOBp1}.nc . > > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_west_S_y${YEAR_BEGIN_JOBp1}.nc . > > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_U_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_V_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_T_y${YEAR_BEGIN_JOBp1}.nc . > ${io_getfile} ${INPUTDIRO}/${OBC}/${YEAR_BEGIN_JOBp1}/obc_south_S_y${YEAR_BEGIN_JOBp1}.nc . > fi -------------------------------------------------------------------------------- diff -b -w cpl01/scripts/getrst_oce.sh cpl02/scripts/getrst_oce.sh -------------------------------------------------------------------------------- 7c7 < ${io_getfile} /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai70c/20061231/restart_oce_erai70c_20061231.nc restart_oce.nc --- > ${io_getfile} /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai73c/20041231/restart_oce_erai73c_20041231.nc restart_oce.nc -------------------------------------------------------------------------------- diff -b -w cpl01/scripts/header.CURIE cpl02/scripts/header.CURIE -------------------------------------------------------------------------------- 36,40c43,46 < export INPUTDIRO="/ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/inputs" < export DATA="Stephane/PERU12-I" < export OBC="obc/ORCA025-B83/D_PERU12.L75_west_interp" < export INPUTDIRA="/ccc/scratch/cont005/ra0542/hourdinc/peru12_wrf/inputs/l60" < export INPUTDIRC="/ccc/scratch/cont005/ra0542/hourdinc/peru12_now/inputs" --- > export INPUTDIRO="/ccc/work/cont005/ra0542/hourdinc/peru12_nemo/inputs" > export OBC="obc/ORCA025_LIM-T323/bathy_trop12" > export INPUTDIRA="/ccc/work/cont005/ra0542/hourdinc/peru12_wrf/inputs/l60" > export INPUTDIRC="/ccc/work/cont005/ra0542/hourdinc/peru12_now/inputs" -------------------------------------------------------------------------------- diff -b -w cpl01/compile/nemo/arch-X64_CURIE.fcm cpl02/compile/nemo/arch-X64_CURIE.fcm -------------------------------------------------------------------------------- 21c21,24 < %FCFLAGS -i4 -r8 -O3 -xAVX -fp-model precise --- > # for Curie Fat and Thin Nodes > #%FCFLAGS -i4 -r8 -O3 -fp-model precise > # for Curie Thin Node > %FCFLAGS -i4 -r8 -O3 -xAVX -fp-model precise -fno-alias -------------------------------------------------------------------------------- diff -b -w -r -x .svn cpl01/compile/nemo/MY_SRC cpl02/compile/nemo/MY_SRC -------------------------------------------------------------------------------- Only in cpl02/compile/nemo/MY_SRC: dynadv_ubs.F90 diff -b -w -r -x .svn cpl01/compile/nemo/MY_SRC/dynldf_lap.F90 cpl02/compile/nemo/MY_SRC/dynldf_lap.F90 76a77 > IF(lwp) WRITE(numout,*) 'Chris sponge sur la dynamique pour OBC' 78d78 < ! WRITE(numout,*) 'Clitandre sponge sur la dynamique: spomask' diff -b -w -r -x .svn cpl01/compile/nemo/MY_SRC/traldf_lap.F90 cpl02/compile/nemo/MY_SRC/traldf_lap.F90 86a87 > IF(lwp) WRITE(numout,*) 'Chris sponge sur les traceurs pour OBC' 100d100 < WRITE(numout,*) 'Chris sponge sur les traceurs' -------------------------------------------------------------------------------- cmp cpl01/compile/nemo/nemo.exe cpl02/compile/nemo/nemo.exe -------------------------------------------------------------------------------- cpl01/compile/nemo/nemo.exe cpl02/compile/nemo/nemo.exe differ: char 41, line 1 ================================================================================ ===== ATMOSPHERE ===== ================================================================================ ================================================================================ ===== COUPLAGE ===== ================================================================================

• modif de la namelist nemo et de la namelist couplage pour ne plus prendre en compte le courrant de surface dans le calcul du stress
• Modif du pas de temps atmosphérique??

================================================================================ ===== diff cpl01 / cpl04 ===== ================================================================================ ===== OCEAN ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl01/param_inputs/namelist.base.oce cpl04/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 228c228 < sn_snd_crt = 'oce only' , 'no' , 'spherical' ,'eastward-northward', 'T' --- > sn_snd_crt = 'none' , 'no' , 'spherical' ,'eastward-northward', 'T' -------------------------------------------------------------------------------- diff -b -w cpl01/scripts/namelist_exp.sh cpl04/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 84c84 < export TIMEsnd=50000 --- > export TIMEsnd=35000 87c87 < export TIMEsnd=36000 --- > export TIMEsnd=25000 ================================================================================ ===== ATMOSPHERE ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl01/param_inputs/namelist.base.atm /ccc/work/cont005/ra0542/oerderv/peru12_now/cpl04/param_inputs/namelist.base.atm -------------------------------------------------------------------------------- 35c35 < auxhist3_interval = 360, --- > auxhist3_interval = 1440, ================================================================================ ===== COUPLAGE ===== ================================================================================ -------------------------------------------------------------------------------- diff -b -w cpl01/param_inputs/namelist.base.cpl /ccc/work/cont005/ra0542/oerderv/peru12_now/cpl04/param_inputs/namelist.base.cpl -------------------------------------------------------------------------------- 6c6 < 9 --- > 7 52,55c52,55 < O_OCurx1 WRF_UOCE_d01 1 <cplfreq> 2 sstoc.nc EXPORTED < 361 505 486 510 n12t w12m LAG=<tsp_oce> < R 0 R 0 < LOCTRANS MAPPING --- > #### DEBUT COM VERA (POUR DESACTIVER COURANTS SURFACE (modifier aussi Nfields) ##### O_OCurx1 WRF_UOCE_d01 1 <cplfreq> 2 sstoc.nc EXPORTED > #361 505 486 510 n12t w12m LAG=<tsp_oce> > #R 0 R 0 > #LOCTRANS MAPPING 57,58c57,58 < AVERAGE < mozaic_o2a.nc --- > #AVERAGE > ###### ######mozaic_o2a.nc 64,67c64,67 < O_OCury1 WRF_VOCE_d01 1 <cplfreq> 2 sstoc.nc EXPORTED < 361 505 486 510 n12t w12m LAG=<tsp_oce> < R 0 R 0 < LOCTRANS MAPPING --- > #O_OCury1 WRF_VOCE_d01 1 <cplfreq> 2 sstoc.nc EXPORTED > #361 505 486 510 n12t w12m LAG=<tsp_oce> > #R 0 R 0 > #LOCTRANS MAPPING 69,70c69,70 < AVERAGE < mozaic_o2a.nc --- > #AVERAGE > ########## FIN COM VERA #####mozaic_o2a.nc