Changeset 10111 for NEMO/trunk
- Timestamp:
- 2018-09-11T11:55:10+02:00 (6 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/cfgs/SHARED/namelist_pisces_ref
r10075 r10111 370 370 ln_hydrofe = .false. ! boolean for from hydrothermal vents 371 371 sedfeinput = 2.e-9 ! Coastal release of Iron 372 distcoast = 5.e3 ! Distance off the coast for Iron from sediments 372 373 dustsolub = 0.02 ! Solubility of the dusta 373 374 mfrac = 0.035 ! Fe mineral fraction of dust -
NEMO/trunk/src/TOP/PISCES/P4Z/p4zsbc.F90
r10068 r10111 37 37 REAL(wp), PUBLIC :: concfediaz !: Fe half-saturation Cste for diazotrophs 38 38 REAL(wp) :: hratio !: Fe:3He ratio assumed for vent iron supply 39 REAL(wp) :: distcoast !: Distance off the coast for Iron from sediments 39 40 REAL(wp), PUBLIC :: fep_rats !: Fep/Fer ratio from sed sources 40 41 REAL(wp), PUBLIC :: fep_rath !: Fep/Fer ratio from hydro sources … … 205 206 INTEGER :: ik50 ! last level where depth less than 50 m 206 207 INTEGER :: isrow ! index for ORCA1 starting row 207 REAL(wp) :: zexpide, zdenitide, zmaskt 208 REAL(wp) :: zexpide, zdenitide, zmaskt, zsurfc, zsurfp,ze3t, ze3t2, zcslp 208 209 REAL(wp) :: ztimes_dust, ztimes_riv, ztimes_ndep 209 210 REAL(wp), DIMENSION(nbtimes) :: zsteps ! times records … … 220 221 & sn_riverdip, sn_riverdop, sn_riverdsi, sn_ndepo, sn_ironsed, sn_hydrofe, & 221 222 & ln_dust, ln_solub, ln_river, ln_ndepo, ln_ironsed, ln_ironice, ln_hydrofe, & 222 & sedfeinput, d ustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, &223 & sedfeinput, distcoast, dustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, & 223 224 & hratio, fep_rats, fep_rath, lgw_rath 224 225 !!---------------------------------------------------------------------- … … 248 249 WRITE(numout,*) ' fe input from hydrothermal vents ln_hydrofe = ', ln_hydrofe 249 250 WRITE(numout,*) ' coastal release of iron sedfeinput = ', sedfeinput 251 WRITE(numout,*) ' distance off the coast distcoast = ', distcoast 250 252 WRITE(numout,*) ' solubility of the dust dustsolub = ', dustsolub 251 253 WRITE(numout,*) ' Mineral Fe content of the dust mfrac = ', mfrac … … 459 461 DO jj = 2, jpjm1 460 462 DO ji = fs_2, fs_jpim1 461 IF( tmask(ji,jj,jk) /= 0. ) THEN 462 zmaskt = tmask(ji+1,jj,jk) * tmask(ji-1,jj,jk) * tmask(ji,jj+1,jk) & 463 & * tmask(ji,jj-1,jk) * tmask(ji,jj,jk+1) 464 IF( zmaskt == 0. ) zcmask(ji,jj,jk ) = MAX( 0.1, zcmask(ji,jj,jk) ) 465 END IF 463 ze3t = e3t_0(ji,jj,jk) 464 zsurfc = e1u(ji,jj) * ( 1. - umask(ji ,jj ,jk) ) & 465 + e1u(ji,jj) * ( 1. - umask(ji-1,jj ,jk) ) & 466 + e2v(ji,jj) * ( 1. - vmask(ji ,jj ,jk) ) & 467 + e2v(ji,jj) * ( 1. - vmask(ji ,jj-1,jk) ) 468 zsurfp = zsurfc * ze3t / e1e2t(ji,jj) 469 ! estimation of the coastal slope : 5 km off the coast 470 ze3t2 = ze3t * ze3t 471 zcslp = SQRT( ( distcoast*distcoast + ze3t2 ) / ze3t2 ) 472 ! 473 zcmask(ji,jj,jk) = zcmask(ji,jj,jk) + zcslp * zsurfp 466 474 END DO 467 475 END DO
Note: See TracChangeset
for help on using the changeset viewer.