peru12 - Tablean des simulations

moyenne < 6 heures / an sur curie Fat nodes 64 nemo.exe + 1 xios_server.exe
(2 fat nodes 64 coeurs nemo.exe + 1 fat node 1 coeur xios_server.exe)

19900101_19901231/erai21c_19900101_19901231.jobid_795489.txt:795489 erai21c_1990 96 3 22/10/2012 18:47:41 5:54:10 RUNNING None 19910101_19911231/erai21c_19910101_19911231.jobid_795957.txt:795957 erai21c_1991 96 3 23/10/2012 00:42:37 5:54:26 RUNNING None 19920101_19921231/erai21c_19920101_19921231.jobid_796333.txt:796333 erai21c_1992 96 3 23/10/2012 06:37:24 5:52:27 RUNNING None 19930101_19931231/erai21c_19930101_19931231.jobid_797257.txt:797257 erai21c_1993 96 3 23/10/2012 15:34:22 5:57:30 RUNNING None 19940101_19941231/erai21c_19940101_19941231.jobid_798724.txt:798724 erai21c_1994 96 3 24/10/2012 09:10:38 5:54:00 RUNNING None

moyenne < 6 heures / an sur curie Fat nodes 64 nemo.exe + 1 xios_server.exe
(2 fat nodes 64 coeurs nemo.exe + 1 fat node 1 coeur xios_server.exe)

19900101_19901231/erai31c_19900101_19901231.jobid_798167.txt:798167 erai31c_1990 96 3 23/10/2012 17:59:15 5:59:17 RUNNING None 19910101_19911231/erai31c_19910101_19911231.jobid_798887.txt:798887 erai31c_1991 96 3 24/10/2012 09:10:38 5:57:42 RUNNING None 19920101_19921231/erai31c_19920101_19921231.jobid_800364.txt:800364 erai31c_1992 96 3 24/10/2012 15:13:28 5:57:21 RUNNING None 19930101_19931231/erai31c_19930101_19931231.jobid_801282.txt:801282 erai31c_1993 96 3 24/10/2012 21:14:34 5:52:19 RUNNING None 19940101_19941231/erai31c_19940101_19941231.jobid_801761.txt:801761 erai31c_1994 96 3 25/10/2012 03:10:02 5:52:29 RUNNING None

diff erai4l / erai21c:

-------------------------------------------------------------------------------- diff -b -w erai21c/param_inputs/namelist.base.oce erai4l/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 552,553c552,553 < ln_traldf_lap = .true. ! laplacian operator < ln_traldf_bilap = .false. ! bilaplacian operator --- > ln_traldf_lap = .false. ! laplacian operator > ln_traldf_bilap = .true. ! bilaplacian operator 556,557c556,557 < ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) < ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") --- > ln_traldf_hor = .true. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) > ln_traldf_iso = .false. ! iso-neutral (needs "key_ldfslp") 568c568 < rn_aht_0 = 300. ! horizontal eddy diffusivity for tracers [m2/s] --- > rn_aht_0 = -1.e9 ! horizontal eddy diffusivity for tracers [m2/s] 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -2.5e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai21c/compile/nemo/cpp_PERU12.fcm erai4l/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_ldfslp key_iomput --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

diff erai5f / erai21c:

-------------------------------------------------------------------------------- diff -b -w erai21c/param_inputs/namelist.base.oce erai5f/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai21c/compile/nemo/cpp_PERU12.fcm erai5f/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_ldfslp key_iomput --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

diff erai50c / erai31c:

diff erai51c / erai50c:

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:

Frontières boostées de 1° 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:

Frontières boostées de 1°5 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

moyenne < 6 heures / an sur curie Thin nodes 64 nemo.exe + 1 xios_server.exe + pdt=600s
(4 thin nodes 64 coeurs nemo.exe + 1 thin node 1 coeur xios_server.exe)

diff erai6a / erai21c:

-------------------------------------------------------------------------------- diff -b -w erai21c/param_inputs/namelist.base.oce erai6a/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 601c601 < ln_dynadv_vec = .true. ! vector form (T) or flux form (F) --- > ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 603c603 < ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme --- > ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s]

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt=169220 max abs(U): 21.95 , i j k: 13 376 1

diff erai6a / erai21c:

diff 19930101_19931231_1_PLANTE/namelist 19930101_19931231_3_PLANTE_sans_bilap_dyn/namelist 636c636 < ln_dynldf_bilap = .true. ! bilaplacian operator --- > ln_dynldf_bilap = .false. ! bilaplacian operator

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt=168314 max abs(U): 197.7 , i j k: 9 376 1

diff erai60c / erai31c:

-------------------------------------------------------------------------------- diff -b -w erai31c/param_inputs/namelist.base.oce erai60c/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 601c601 < ln_dynadv_vec = .true. ! vector form (T) or flux form (F) --- > ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 603c603 < ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme --- > ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme 636c636 < ln_dynldf_bilap = .true. ! bilaplacian operator --- > ln_dynldf_bilap = .false. ! bilaplacian operator 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai31c/compile/nemo/cpp_PERU12.fcm erai60c/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput key_ldfslp --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt= 8532 max abs(U): 22.51 , i j k: 18 377 9

diff erai61c / erai31c:

-------------------------------------------------------------------------------- diff -b -w erai31c/param_inputs/namelist.base.oce erai61c/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 601c601 < ln_dynadv_vec = .true. ! vector form (T) or flux form (F) --- > ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 603c603 < ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme --- > ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme 636c636 < ln_dynldf_bilap = .true. ! bilaplacian operator --- > ln_dynldf_bilap = .false. ! bilaplacian operator 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai31c/scripts/namelist_exp.sh erai61c/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 23c23 < export TSP_OCE=800 --- > export TSP_OCE=400 -------------------------------------------------------------------------------- diff -b -w erai31c/compile/nemo/cpp_PERU12.fcm erai61c/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput key_ldfslp --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

diff erai62c / erai61c:

-------------------------------------------------------------------------------- diff -b -w erai61c/scripts/namelist_exp.sh erai62c/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 23c23 < export TSP_OCE=400 --- > export TSP_OCE=720

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt= 9972 max abs(U): 42.72 , i j k: 13 379 13

diff erai63c / erai62c:

-------------------------------------------------------------------------------- diff -b -w erai62c/scripts/namelist_exp.sh erai63c/scripts/namelist_exp.sh -------------------------------------------------------------------------------- 23c23 < export TSP_OCE=720 --- > export TSP_OCE=600

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt= 68824 max abs(U): 31.52 , i j k: 11 378 1

================================================================================ ===== 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

run	base	periode	restart	little_nemo	pdt	cpu		machine
run	base	periode	restart	little_nemo	pdt	proc	time	machine
erai60c	erai31c	1990-11 Mars 1990	levitus	r48	800			curie fat nodes
erai61c	erai60c	1990-1994	levitus	r48	400	64 nemo +1 xios	~6h/an	curie thin nodes
erai62c	erai61c	1990-24 Mars 1990	levitus	r48	720	64 nemo +1 xios		curie thin nodes
erai63c	erai62c	1990-22 Avril 1991	levitus	r48	600	64 nemo +1 xios		curie thin nodes
erai65c	erai61c + 46 niveaux => 75 niveaux (mêmes obc ORCA025.L46-B83-DFS4/ 46 => matlab => 75) (bug sur UBS traceur 2 x trop diffusif)
	erai61c	1990-1995	levitus	r62	600 400 400 400	64 nemo+1 xios 64 nemo+1 xios 63 nemo +1 xios 126 nemo + 2 xios	5,8 heures/an 9 heures/an 8,8 heures/an 6,6 heures/an	curie thin nodes
erai66c	erai65c + sans bug x2 sur UBS traceur => 2 x moins diffusif => plante sur obc west à l'équateur
	erai65c	1990-20 Avril 1993 ...10 Avril 1994	levitus	r62	600 400	126 nemo + 2 xios	4,5 heures/an 7.0 heures/an	curie thin nodes
erai67c	erai66c + Interpolation sur la frontière west au lieu de la dupliquer façon Stephane + Charles.
	erai66c	1990-30 Avril 1991 ...-26 Mai 1997	levitus	r62	600 400	126 nemo +2 xios	4,3 heures/an 6.6 heures/an	curie thin nodes
erai68c	erai67c + smooth la frontière west
	erai67c	1997-18 Fév 2000	restart 1997 erai67c	r62	400	126 nemo +2 xios	6,6 heures/an	curie thin nodes
erai69c	erai67c + sponge (laplacien + masque) sur la dynamique
	erai67c	1990-2007	levitus	r62	600	126 nemo +2 xios		curie thin nodes
erai70a	erai69c + Phasage avec la namelist trop12 de Seb pour Pulsation + update little_nemo r73 (div=0 pour OBC)
	erai69c	1990-2007	levitus	r73	600	126 nemo +2 xios		curie thin nodes
erai70b	erai70a + light penetration: retour a RGB (Red-Green-Blue) au lieu du 2 bands pour
	erai70a	1990-2007	levitus	r73	600 400	126 nemo +2 xios		curie thin nodes
erai70c	erai70b + Bug UBS dyn: gamma2 = 1/32 au lieu de 1/8
	erai70b	1990-2007	levitus	r73	400	126 nemo +2 xios		curie thin nodes
erai70d	erai70c + On retire la contrainte div=0 pour OBC (pas de div dans UBS => sans impact => erai70d=erai70c)
	erai70c	erai70d = erai70c	levitus	r73	400	126 nemo +2 xios		curie thin nodes
erai70e	erai70c + dynvor (celui utilisé par Gurvan et Andrew en z-tilde coord) - pour le pb damier à l'equateur (25/10/2013)
	erai70c	2000-2007	restart 2000 erai70c	r73	400	126 nemo +2 xios		curie thin nodes
erai71a	erai70a + Vertical diffusion: GLS (k-eps et CanutoA) au lieu de TKE ln_zdfevd=.true. au lieu de .false. NaN à partir de pdt=15
	erai70a	1990... Plante!	levitus	r73	600	126 nemo +2 xios		curie thin nodes
erai71b	erai71a + ln_zdfevd (enhanced vertical diffusion) = .false. ocean.abort a pdt=11
	erai71a	1990... Plante!	levitus	r73	600	126 nemo +2 xios		curie thin nodes
erai72c	erai70c + Vents quikscat au lieu de erai
		2000-2007	restart 2000 erai70c	r73	400	126 nemo +2 xios		curie thin nodes
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
erai73c	erai72c + OBC mercator ORCA025_LIM-T323 + bathy trop12 + Etat initial ad'hoc
		2000-2008	clim-ORCA025_LIM-T323	r73	400	126 nemo +2 xios		curie thin nodes
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
erai73c2	erai73c + retour solaire 6h + retour TKE ln_mxl0=.false.
		2000-2008	clim-ORCA025_LIM-T323	r73	400	126 nemo +2 xios		curie thin nodes
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
erai74c	erai73c + bathy Stephane au lieu de trop12 (OBC mercator ORCA025_LIM-T323)
		2000-2008	clim-ORCA025_LIM-T323	r73	400	126 nemo +2 xios		curie thin nodes
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" ================================================================================
bdy	erai74c + implementation BDY dans little_nemo r85 ... pour vérif du merge des codes (key_obc)
		2000-2003	clim-ORCA025_LIM-T323	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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)
bdy2	bdy + key_bdy (bdy coordinates, mais pas de data)
		2000 (test 5days)	clim-ORCA025_LIM-T323	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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)
bdy3	bdy2 + Prescription de T,S,U,V,ssh aux frontières => erai74c en bdy (frontières mercator ORCA025_LIM-T323)(Attention: sorties 5d incomplètes!!)
		2000-2008	clim-ORCA025_LIM-T323	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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
bdy3b	bdy3 + volume correction (perd 10cm par an)
		2000 (- day 365)	clim-ORCA025_LIM-T323	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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)
bdy4 (rn_usecrt=1 !!)	bdy3 + frontières extraites de trop12 - namelist de trop12 exactement (attention RGV/2band) + volume correction (perd 10cm par an)
		2000 - 2008	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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
bdy4b (rn_usecrt=1 !!)	bdy4 + volume correction=FALSE (frontières TROP12)
		2000 - 2008	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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)
newscript	Scripts à zéro

Objet: Scripts à zéro
bdy4b_sh1	bdy4b + Mise à plat des scripts
		2000 en cours...	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r85_BDY	400	126 nemo +2 xios		curie thin nodes
Objet: Mise à plat des scripts
erai75c	bdy4 + obc au lieu de bdy (little_nemo r73 au lieu de r85_BDY sans impact) sans understress comme tous les runs OBC
		2000 - 2008	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r73	400	126 nemo +2 xios		curie thin nodes
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"
bdy4c (rn_usecrt=1 !!)	bdy4b + Correction bug free-slip dans dynadv_ubs.F90 (Nicolas Ducousso)
		2000-2008	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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
bdy4d (rn_usecrt=1 !!)	bdy4c + Correction bug dynkeg et dynvor (Nicolas Ducousso) - dynvor déjà testée dans erai70e
		2000-2008	T et S de tr12_quik_5d_20000105_grid_T_3D.nc	r85_BDY	400	126 nemo +2 xios		curie thin nodes
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

Configuration peru12_nemo

Tablean des simulations

Où trouver les runs

Type de runs

"ERAI2..." (bilap/lap) avec understress

"ERAI3..." (bilap/lap) sans understress

"erai4..." (bilap/bilap)

"erai5..." (bilap/UBS)

"erai6... part I" (UBS/UBS)

"erai6... part II" (UBS/UBS)

Mise à plat + Convergence complète trop12 + bdy OK

"cpl"

Momentum				Tracer				ppacr	understress
advection	diffusion			advection	diffusion
UBS	none	none	none	UBS	none	none	none	9.000_wp	sans