Changeset 10368 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zsbc.F90
- Timestamp:
- 2018-12-03T12:45:01+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zsbc.F90
r10314 r10368 31 31 REAL(wp), PUBLIC :: dustsolub !: Solubility of the dust 32 32 REAL(wp), PUBLIC :: mfrac !: Mineral Content of the dust 33 REAL(wp), PUBLIC :: rdustfep !: Fraction of dust that is dissolvable 33 34 REAL(wp), PUBLIC :: icefeinput !: Iron concentration in sea ice 34 35 REAL(wp), PUBLIC :: wdust !: Sinking speed of the dust … … 134 135 DO ji = 1, jpi 135 136 zcoef = ryyss * e1e2t(ji,jj) * h_rnf(ji,jj) 136 rivalk(ji,jj) = sf_river(jr_dic)%fnow(ji,jj,1) 137 rivalk(ji,jj) = sf_river(jr_dic)%fnow(ji,jj,1) & 137 138 & * 1.E3 / ( 12. * zcoef + rtrn ) 138 rivdic(ji,jj) = ( sf_river(jr_dic)%fnow(ji,jj,1) + sf_river(jr_doc)%fnow(ji,jj,1) )&139 rivdic(ji,jj) = sf_river(jr_dic)%fnow(ji,jj,1) & 139 140 & * 1.E3 / ( 12. * zcoef + rtrn ) 140 rivdin(ji,jj) = ( sf_river(jr_din)%fnow(ji,jj,1) + sf_river(jr_don)%fnow(ji,jj,1) )&141 rivdin(ji,jj) = sf_river(jr_din)%fnow(ji,jj,1) & 141 142 & * 1.E3 / rno3 / ( 14. * zcoef + rtrn ) 142 rivdip(ji,jj) = ( sf_river(jr_dip)%fnow(ji,jj,1) + sf_river(jr_dop)%fnow(ji,jj,1) )&143 rivdip(ji,jj) = sf_river(jr_dip)%fnow(ji,jj,1) & 143 144 & * 1.E3 / po4r / ( 31. * zcoef + rtrn ) 144 rivdsi(ji,jj) = sf_river(jr_dsi)%fnow(ji,jj,1) 145 rivdsi(ji,jj) = sf_river(jr_dsi)%fnow(ji,jj,1) & 145 146 & * 1.E3 / ( 28.1 * zcoef + rtrn ) 147 rivdoc(ji,jj) = sf_river(jr_doc)%fnow(ji,jj,1) & 148 & * 1.E3 / ( 12. * zcoef + rtrn ) 146 149 END DO 147 150 END DO … … 158 161 rivdip(ji,jj) = ( sf_river(jr_dip)%fnow(ji,jj,1) ) & 159 162 & * 1.E3 / po4r / ( 31. * zcoef + rtrn ) * tmask(ji,jj,1) 160 rivdoc(ji,jj) = ( sf_river(jr_doc)%fnow(ji,jj,1) ) &161 & * 1.E3 / ( 12. * zcoef + rtrn ) * tmask(ji,jj,1)162 163 rivdon(ji,jj) = ( sf_river(jr_don)%fnow(ji,jj,1) ) & 163 164 & * 1.E3 / rno3 / ( 14. * zcoef + rtrn ) * tmask(ji,jj,1) 164 165 rivdop(ji,jj) = ( sf_river(jr_dop)%fnow(ji,jj,1) ) & 165 166 & * 1.E3 / po4r / ( 31. * zcoef + rtrn ) * tmask(ji,jj,1) 167 rivdsi(ji,jj) = sf_river(jr_dsi)%fnow(ji,jj,1) & 168 & * 1.E3 / ( 28.1 * zcoef + rtrn ) 169 rivdoc(ji,jj) = sf_river(jr_doc)%fnow(ji,jj,1) & 170 & * 1.E3 / ( 12. * zcoef + rtrn ) 166 171 END DO 167 172 END DO … … 223 228 & ln_dust, ln_solub, ln_river, ln_ndepo, ln_ironsed, ln_ironice, ln_hydrofe, & 224 229 & sedfeinput, distcoast, dustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, & 225 & hratio, fep_rats, fep_rath, lgw_rath230 & hratio, fep_rats, fep_rath, rdustfep, lgw_rath 226 231 !!---------------------------------------------------------------------- 227 232 ! … … 262 267 WRITE(numout,*) ' Fep/Fer ratio from sed sources fep_rats = ', fep_rats 263 268 WRITE(numout,*) ' Fep/Fer ratio from sed hydro sources fep_rath = ', fep_rath 269 WRITE(numout,*) ' Fraction of dust that is dissolvable rdustfep = ', rdustfep 264 270 WRITE(numout,*) ' Weak ligand ratio from sed hydro sources lgw_rath = ', lgw_rath 265 271 ENDIF … … 343 349 slf_river(jr_dsi) = sn_riverdsi 344 350 ! 345 ALLOCATE( rivdic(jpi,jpj), rivalk(jpi,jpj), rivdin(jpi,jpj), rivdip(jpi,jpj), rivdsi(jpi,jpj) )346 IF( ln_p5z ) ALLOCATE( rivdon(jpi,jpj), rivdop(jpi,jpj) , rivdoc(jpi,jpj))351 ALLOCATE( rivdic(jpi,jpj), rivalk(jpi,jpj), rivdin(jpi,jpj), rivdip(jpi,jpj), rivdsi(jpi,jpj), rivdoc(jpi,jpj) ) 352 IF( ln_p5z ) ALLOCATE( rivdon(jpi,jpj), rivdop(jpi,jpj) ) 347 353 ! 348 354 ALLOCATE( sf_river(jpriv), rivinput(jpriv), STAT=ierr1 ) !* allocate and fill sf_river (forcing structure) with sn_river_
Note: See TracChangeset
for help on using the changeset viewer.