Changeset 12551 for NEMO/branches
- Timestamp:
- 2020-03-13T15:51:01+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback
- Files:
-
- 1 deleted
- 12 edited
-
cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg (modified) (2 diffs)
-
cfgs/SHARED/namelist_ref (modified) (6 diffs)
-
src/ABL/sbcabl.F90 (modified) (1 diff)
-
src/ICE/icectl.F90 (modified) (1 diff)
-
src/OCE/BDY/bdydta.F90 (modified) (9 diffs)
-
src/OCE/LBC/lib_mpp.F90 (modified) (2 diffs)
-
src/OCE/SBC/cpl_oasis3.F90 (modified) (2 diffs)
-
src/OCE/ZDF/zdfiwm.F90 (modified) (3 diffs)
-
tests/BENCH/EXPREF/namelist_cfg_orca025_like (modified) (1 diff)
-
tests/BENCH/EXPREF/namelist_cfg_orca12_like (modified) (1 diff)
-
tests/BENCH/EXPREF/namelist_cfg_orca1_like (modified) (1 diff)
-
tests/BENCH/EXPREF/namelist_pisces_cfg (modified) (1 diff)
-
tests/BENCH/MY_SRC/zdfiwm.F90 (deleted)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r12495 r12551 28 28 &namdom ! time and space domain 29 29 !----------------------------------------------------------------------- 30 rn_Dt = 5400. ! time step for the dynamics and tracer30 rn_Dt = 5400. ! time step for the dynamics and tracer 31 31 / 32 32 !----------------------------------------------------------------------- … … 381 381 ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency 382 382 ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) 383 384 cn_dir = './' ! root directory for the iwm data location 385 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 386 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 387 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 388 sn_mpb = 'mixing_power_bot' , -12 , 'field' , .false. , .true. , 'yearly' , '' , '' , '' 389 sn_mpp = 'mixing_power_pyc' , -12 , 'field' , .false. , .true. , 'yearly' , '' , '' , '' 390 sn_mpc = 'mixing_power_cri' , -12 , 'field' , .false. , .true. , 'yearly' , '' , '' , '' 391 sn_dsb = 'decay_scale_bot' , -12 , 'field' , .false. , .true. , 'yearly' , '' , '' , '' 392 sn_dsc = 'decay_scale_cri' , -12 , 'field' , .false. , .true. , 'yearly' , '' , '' , '' 383 393 / 384 394 !!====================================================================== -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/cfgs/SHARED/namelist_ref
r12499 r12551 72 72 ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time 73 73 ! 74 rn_Dt = 5400. ! time step for the dynamics and tracer74 rn_Dt = 5400. ! time step for the dynamics and tracer 75 75 rn_atfp = 0.1 ! asselin time filter parameter 76 76 ! … … 217 217 nn_ice = 0 ! =0 no ice boundary condition 218 218 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 219 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 220 ! ! except in AGRIF zoom where it has to be specified 219 ! ! =2 or 3 for SI3 and CICE, respectively 221 220 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 222 221 ! ! =F levitating ice (no pressure, mass and salt exchanges) … … 272 271 rn_pfac = 1. ! multipl. factor for precipitation (total & snow) 273 272 rn_efac = 1. ! multipl. factor for evaporation (0. or 1.) 274 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 273 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 275 274 ! ! used to calculate the wind stress 276 275 ! ! (0. => absolute or 1. => relative winds) … … 292 291 sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 293 292 sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 294 sn_hpgi = 'NO NE', 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , ''295 sn_hpgj = 'NO NE', 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , ''293 sn_hpgi = 'NOT USED' , 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , '' 294 sn_hpgj = 'NOT USED' , 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , '' 296 295 sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 297 296 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' … … 999 998 ln_bt_auto = .true. ! Number of sub-step defined from: 1000 999 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 1001 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds1000 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds 1002 1001 rn_bt_alpha = 0. ! Temporal diffusion parameter (if ln_bt_av=F) 1003 1002 / … … 1193 1192 ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency 1194 1193 ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) 1195 / 1196 1194 1195 cn_dir = './' ! root directory for the iwm data location 1196 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 1197 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 1198 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 1199 sn_mpb = 'NOT USED' , -12 , 'mixing_power_bot' , .false. , .true. , 'yearly' , '' , '' , '' 1200 sn_mpp = 'NOT USED' , -12 , 'mixing_power_pyc' , .false. , .true. , 'yearly' , '' , '' , '' 1201 sn_mpc = 'NOT USED' , -12 , 'mixing_power_cri' , .false. , .true. , 'yearly' , '' , '' , '' 1202 sn_dsb = 'NOT USED' , -12 , 'decay_scale_bot' , .false. , .true. , 'yearly' , '' , '' , '' 1203 sn_dsc = 'NOT USED' , -12 , 'decay_scale_cri' , .false. , .true. , 'yearly' , '' , '' , '' 1204 / 1197 1205 !!====================================================================== 1198 1206 !! *** Diagnostics namelists *** !! -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/ABL/sbcabl.F90
r12495 r12551 75 75 !!--------------------------------------------------------------------- 76 76 77 REWIND( numnam_ref )! Namelist namsbc_abl in reference namelist : ABL parameters77 ! Namelist namsbc_abl in reference namelist : ABL parameters 78 78 READ ( numnam_ref, namsbc_abl, IOSTAT = ios, ERR = 901 ) 79 79 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in reference namelist' ) 80 ! 81 REWIND( numnam_cfg ) ! Namelist namsbc_abl in configuration namelist : ABL parameters 80 ! Namelist namsbc_abl in configuration namelist : ABL parameters 82 81 READ ( numnam_cfg, namsbc_abl, IOSTAT = ios, ERR = 902 ) 83 82 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in configuration namelist' ) -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/ICE/icectl.F90
r12495 r12551 725 725 726 726 CALL prt_ctl_info(' ') 727 CALL prt_ctl_info(' - Heat / FW fluxes : ')728 CALL prt_ctl_info(' ~~~~~~~~~~~~~~~~~~ ')729 CALL prt_ctl(tab2d_1=sst_m , clinfo1= ' sst : ', tab2d_2=sss_m , clinfo2= ' sss : ')730 CALL prt_ctl(tab2d_1=qsr , clinfo1= ' qsr : ', tab2d_2=qns , clinfo2= ' qns : ')731 CALL prt_ctl(tab2d_1=emp , clinfo1= ' emp : ', tab2d_2=sfx , clinfo2= ' sfx : ')732 733 CALL prt_ctl_info(' ')734 727 CALL prt_ctl_info(' - Stresses : ') 735 728 CALL prt_ctl_info(' ~~~~~~~~~~ ') -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/BDY/bdydta.F90
r12396 r12551 92 92 INTEGER :: ii, ij, ik, igrd, ipl ! local integers 93 93 INTEGER, DIMENSION(jpbgrd) :: ilen1 94 INTEGER, DIMENSION(:), POINTER :: nblen, nblenrim ! short cuts95 94 TYPE(OBC_DATA) , POINTER :: dta_alias ! short cut 96 95 TYPE(FLD), DIMENSION(:), POINTER :: bf_alias … … 108 107 DO jbdy = 1, nb_bdy 109 108 ! 110 nblen => idx_bdy(jbdy)%nblen111 nblenrim => idx_bdy(jbdy)%nblenrim112 !113 109 IF( nn_dyn2d_dta(jbdy) == 0 ) THEN 114 ilen1(:) = nblen(:)115 110 IF( dta_bdy(jbdy)%lneed_ssh ) THEN 116 111 igrd = 1 117 DO ib = 1, i len1(igrd)112 DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) ! ssh is allocated and used only on the rim 118 113 ii = idx_bdy(jbdy)%nbi(ib,igrd) 119 114 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 123 118 IF( dta_bdy(jbdy)%lneed_dyn2d) THEN 124 119 igrd = 2 125 DO ib = 1, ilen1(igrd)120 DO ib = 1, SIZE(dta_bdy(jbdy)%u2d) ! u2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 126 121 ii = idx_bdy(jbdy)%nbi(ib,igrd) 127 122 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 129 124 END DO 130 125 igrd = 3 131 DO ib = 1, ilen1(igrd)126 DO ib = 1, SIZE(dta_bdy(jbdy)%v2d) ! v2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 132 127 ii = idx_bdy(jbdy)%nbi(ib,igrd) 133 128 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 138 133 ! 139 134 IF( nn_dyn3d_dta(jbdy) == 0 ) THEN 140 ilen1(:) = nblen(:)141 135 IF( dta_bdy(jbdy)%lneed_dyn3d ) THEN 142 136 igrd = 2 143 DO ib = 1, i len1(igrd)137 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 144 138 DO ik = 1, jpkm1 145 139 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 149 143 END DO 150 144 igrd = 3 151 DO ib = 1, i len1(igrd)145 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 152 146 DO ik = 1, jpkm1 153 147 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 160 154 161 155 IF( nn_tra_dta(jbdy) == 0 ) THEN 162 ilen1(:) = nblen(:)163 156 IF( dta_bdy(jbdy)%lneed_tra ) THEN 164 157 igrd = 1 165 DO ib = 1, i len1(igrd)158 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 166 159 DO ik = 1, jpkm1 167 160 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 176 169 #if defined key_si3 177 170 IF( nn_ice_dta(jbdy) == 0 ) THEN ! set ice to initial values 178 ilen1(:) = nblen(:)179 171 IF( dta_bdy(jbdy)%lneed_ice ) THEN 180 172 igrd = 1 181 173 DO jl = 1, jpl 182 DO ib = 1, i len1(igrd)174 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 183 175 ii = idx_bdy(jbdy)%nbi(ib,igrd) 184 176 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 338 330 DO jbdy = 1, nb_bdy ! Tidal component added in ts loop 339 331 IF ( nn_dyn2d_dta(jbdy) .GE. 2 ) THEN 340 nblen => idx_bdy(jbdy)%nblen 341 nblenrim => idx_bdy(jbdy)%nblenrim 342 IF( cn_dyn2d(jbdy) == 'frs' ) THEN ; ilen1(:)=nblen(:) 343 ELSE ; ilen1(:)=nblenrim(:) 332 IF( cn_dyn2d(jbdy) == 'frs' ) THEN ; ilen1(:)=idx_bdy(jbdy)%nblen(:) 333 ELSE ; ilen1(:)=idx_bdy(jbdy)%nblenrim(:) 344 334 ENDIF 345 335 IF ( dta_bdy(jbdy)%lneed_ssh ) dta_bdy_s(jbdy)%ssh(1:ilen1(1)) = dta_bdy(jbdy)%ssh(1:ilen1(1)) -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/LBC/lib_mpp.F90
r12377 r12551 402 402 # if defined key_mpi2 403 403 IF( ln_timing ) CALL tic_tac( .TRUE., ld_global = .TRUE.) 404 CALL mpi_allreduce( y_in(:), todelay(idvar)%y1d(:), isz, MPI_DOUBLE_COMPLEX, mpi_sumdd, ilocalcomm, ndelayid(idvar), ierr ) 404 CALL mpi_allreduce( y_in(:), todelay(idvar)%y1d(:), isz, MPI_DOUBLE_COMPLEX, mpi_sumdd, ilocalcomm, ierr ) 405 ndelayid(idvar) = 1 405 406 IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.) 406 407 # else … … 469 470 # if defined key_mpi2 470 471 IF( ln_timing ) CALL tic_tac( .TRUE., ld_global = .TRUE.) 471 CALL mpi_allreduce( p_in(:), todelay(idvar)%z1d(:), isz, MPI_DOUBLE_PRECISION, mpi_max, ilocalcomm, ndelayid(idvar), ierr ) 472 CALL mpi_allreduce( p_in(:), todelay(idvar)%z1d(:), isz, MPI_DOUBLE_PRECISION, mpi_max, ilocalcomm, ierr ) 473 ndelayid(idvar) = 1 472 474 IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.) 473 475 # else -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/SBC/cpl_oasis3.F90
r12377 r12551 365 365 WRITE(numout,*) ' - Minimum value is ', MINVAL(pdata(nldi:nlei,nldj:nlej,jc)) 366 366 WRITE(numout,*) ' - Maximum value is ', MAXVAL(pdata(nldi:nlei,nldj:nlej,jc)) 367 WRITE(numout,*) ' - Sum value is ', SUM(pdata(nldi:nlei,nldj:nlej,jc))367 WRITE(numout,*) ' - Sum value is ', SUM(pdata(nldi:nlei,nldj:nlej,jc)) 368 368 WRITE(numout,*) '****************' 369 369 ENDIF … … 444 444 WRITE(numout,*) 'oasis_get: kstep', kstep 445 445 WRITE(numout,*) 'oasis_get: info ', kinfo 446 WRITE(numout,*) ' - Minimum value is ', MINVAL(pdata( :,:,jc))447 WRITE(numout,*) ' - Maximum value is ', MAXVAL(pdata( :,:,jc))448 WRITE(numout,*) ' - Sum value is ', SUM(pdata(:,:,jc))446 WRITE(numout,*) ' - Minimum value is ', MINVAL(pdata(nldi:nlei,nldj:nlej,jc)) 447 WRITE(numout,*) ' - Maximum value is ', MAXVAL(pdata(nldi:nlei,nldj:nlej,jc)) 448 WRITE(numout,*) ' - Sum value is ', SUM(pdata(nldi:nlei,nldj:nlej,jc)) 449 449 WRITE(numout,*) '****************' 450 450 ENDIF -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/ZDF/zdfiwm.F90
r12495 r12551 23 23 USE phycst ! physical constants 24 24 ! 25 USE fldread ! field read 25 26 USE prtctl ! Print control 26 27 USE in_out_manager ! I/O manager … … 383 384 !! de Lavergne et al. in prep., 2017 384 385 !!---------------------------------------------------------------------- 385 INTEGER :: inum ! local integer 386 INTEGER :: ifpr ! dummy loop indices 387 INTEGER :: inum ! local integer 386 388 INTEGER :: ios 387 389 REAL(wp) :: zbot, zpyc, zcri ! local scalars 388 !! 389 NAMELIST/namzdf_iwm/ nn_zpyc, ln_mevar, ln_tsdiff 390 ! 391 CHARACTER(len=256) :: cn_dir ! Root directory for location of ssr files 392 INTEGER, PARAMETER :: jpiwm = 5 ! maximum number of files to read 393 INTEGER, PARAMETER :: jp_mpb = 1 394 INTEGER, PARAMETER :: jp_mpp = 2 395 INTEGER, PARAMETER :: jp_mpc = 3 396 INTEGER, PARAMETER :: jp_dsb = 4 397 INTEGER, PARAMETER :: jp_dsc = 5 398 ! 399 TYPE(FLD_N), DIMENSION(jpiwm) :: slf_iwm ! array of namelist informations 400 TYPE(FLD_N) :: sn_mpb, sn_mpp, sn_mpc ! informations about Mixing Power field to be read 401 TYPE(FLD_N) :: sn_dsb, sn_dsc ! informations about Decay Scale field to be read 402 TYPE(FLD ), DIMENSION(jpiwm) :: sf_iwm ! structure of input fields (file informations, fields read) 403 ! 404 NAMELIST/namzdf_iwm/ nn_zpyc, ln_mevar, ln_tsdiff, & 405 & cn_dir, sn_mpb, sn_mpp, sn_mpc, sn_dsb, sn_dsc 390 406 !!---------------------------------------------------------------------- 391 407 ! … … 422 438 IF( zdf_iwm_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'zdf_iwm_init : unable to allocate iwm arrays' ) 423 439 ! 440 ! store namelist information in an array 441 slf_iwm(jp_mpb) = sn_mpb ; slf_iwm(jp_mpp) = sn_mpp ; slf_iwm(jp_mpc) = sn_mpc 442 slf_iwm(jp_dsb) = sn_dsb ; slf_iwm(jp_dsc) = sn_dsc 443 ! 444 DO ifpr= 1, jpiwm 445 ALLOCATE( sf_iwm(ifpr)%fnow(jpi,jpj,1) ) 446 IF( slf_iwm(ifpr)%ln_tint )ALLOCATE( sf_iwm(ifpr)%fdta(jpi,jpj,1,2) ) 447 END DO 448 449 ! fill sf_iwm with sf_iwm and control print 450 CALL fld_fill( sf_iwm, slf_iwm , cn_dir, 'zdfiwm_init', 'iwm input file', 'namiwm' ) 451 452 ! ! hard-coded default definition (to be defined in namelist ?) 453 sf_iwm(jp_mpb)%fnow(:,:,1) = 1.e-6 454 sf_iwm(jp_mpp)%fnow(:,:,1) = 1.e-6 455 sf_iwm(jp_mpc)%fnow(:,:,1) = 1.e-10 456 sf_iwm(jp_dsb)%fnow(:,:,1) = 100. 457 sf_iwm(jp_dsc)%fnow(:,:,1) = 100. 458 424 459 ! ! read necessary fields 425 CALL iom_open('mixing_power_bot',inum) ! energy flux for high-mode wave breaking [W/m2] 426 CALL iom_get (inum, jpdom_data, 'field', ebot_iwm, 1 ) 427 CALL iom_close(inum) 428 ! 429 CALL iom_open('mixing_power_pyc',inum) ! energy flux for pynocline-intensified wave breaking [W/m2] 430 CALL iom_get (inum, jpdom_data, 'field', epyc_iwm, 1 ) 431 CALL iom_close(inum) 432 ! 433 CALL iom_open('mixing_power_cri',inum) ! energy flux for critical slope wave breaking [W/m2] 434 CALL iom_get (inum, jpdom_data, 'field', ecri_iwm, 1 ) 435 CALL iom_close(inum) 436 ! 437 CALL iom_open('decay_scale_bot',inum) ! spatially variable decay scale for high-mode wave breaking [m] 438 CALL iom_get (inum, jpdom_data, 'field', hbot_iwm, 1 ) 439 CALL iom_close(inum) 440 ! 441 CALL iom_open('decay_scale_cri',inum) ! spatially variable decay scale for critical slope wave breaking [m] 442 CALL iom_get (inum, jpdom_data, 'field', hcri_iwm, 1 ) 443 CALL iom_close(inum) 444 445 ebot_iwm(:,:) = ebot_iwm(:,:) * ssmask(:,:) 446 epyc_iwm(:,:) = epyc_iwm(:,:) * ssmask(:,:) 447 ecri_iwm(:,:) = ecri_iwm(:,:) * ssmask(:,:) 460 CALL fld_read( nit000, 1, sf_iwm ) 461 462 ebot_iwm(:,:) = sf_iwm(1)%fnow(:,:,1) * ssmask(:,:) ! energy flux for high-mode wave breaking [W/m2] 463 epyc_iwm(:,:) = sf_iwm(2)%fnow(:,:,1) * ssmask(:,:) ! energy flux for pynocline-intensified wave breaking [W/m2] 464 ecri_iwm(:,:) = sf_iwm(3)%fnow(:,:,1) * ssmask(:,:) ! energy flux for critical slope wave breaking [W/m2] 465 hbot_iwm(:,:) = sf_iwm(4)%fnow(:,:,1) ! spatially variable decay scale for high-mode wave breaking [m] 466 hcri_iwm(:,:) = sf_iwm(5)%fnow(:,:,1) ! spatially variable decay scale for critical slope wave breaking [m] 448 467 449 468 zbot = glob_sum( 'zdfiwm', e1e2t(:,:) * ebot_iwm(:,:) ) 450 469 zpyc = glob_sum( 'zdfiwm', e1e2t(:,:) * epyc_iwm(:,:) ) 451 470 zcri = glob_sum( 'zdfiwm', e1e2t(:,:) * ecri_iwm(:,:) ) 471 452 472 IF(lwp) THEN 453 473 WRITE(numout,*) ' High-mode wave-breaking energy: ', zbot * 1.e-12_wp, 'TW' -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca025_like
r12495 r12551 50 50 !----------------------------------------------------------------------- 51 51 ln_usr = .true. ! user defined formulation (T => check usrdef_sbc) 52 nn_ice = 2 ! =0 no ice boundary condition 53 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 54 ! ! =2 or 3 for SI3 and CICE, respectively 52 55 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 53 56 / -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca12_like
r12495 r12551 50 50 !----------------------------------------------------------------------- 51 51 ln_usr = .true. ! user defined formulation (T => check usrdef_sbc) 52 nn_ice = 2 ! =0 no ice boundary condition 53 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 54 ! ! =2 or 3 for SI3 and CICE, respectively 52 55 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 53 56 / -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca1_like
r12495 r12551 50 50 !----------------------------------------------------------------------- 51 51 ln_usr = .true. ! user defined formulation (T => check usrdef_sbc) 52 nn_ice = 2 ! =0 no ice boundary condition 53 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 54 ! ! =2 or 3 for SI3 and CICE, respectively 52 55 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 53 56 / -
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_pisces_cfg
r10343 r12551 81 81 / 82 82 !----------------------------------------------------------------------- 83 &nampis sbc! parameters for inputs deposition83 &nampisbc ! parameters for inputs deposition 84 84 !----------------------------------------------------------------------- 85 ln_dust = .false. ! boolean for dust input from the atmosphere86 ln_solub = .false. ! boolean for variable solubility of atm. Iron87 ln_river = .false. ! boolean for river input of nutrients88 ln_ndepo = .false. ! boolean for atmospheric deposition of N89 85 ln_ironsed = .false. ! boolean for Fe input from sediments 90 86 ln_ironice = .false. ! boolean for Fe input from sea ice
Note: See TracChangeset
for help on using the changeset viewer.
