Changeset 8423 for branches/NERC/dev_r5518_GO6_split_trcbiomedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
- Timestamp:
- 2017-08-08T16:14:00+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_split_trcbiomedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
r8419 r8423 32 32 !!---------------------------------------------------------------------- 33 33 USE bio_medusa_mod 34 USE dom_oce, ONLY: atfp, atfp1, neuler, rdt 34 USE dom_oce, ONLY: atfp, atfp1, neuler, rdt, e3t_n, tmask 35 35 USE in_out_manager, ONLY: lwp, numout 36 36 # if defined key_iomput … … 61 61 INTEGER :: ji, jj 62 62 INTEGER :: jn 63 64 REAL(wp) :: fq0,fq1,fq2,fq3 65 63 66 !!---------------------------------------------------------------------- 64 67 !! Process benthic in/out fluxes … … 235 238 ENDIF 236 239 endif 237 240 241 # if defined key_debug_medusa 242 !! AXY (12/07/17) 243 !!----------------------------------------------------------------- 244 !! Check conservation of MEDUSA's sinks-minus-sources using fflx_X 245 !! diagnostics (i.e. biogeochemical processes only) 246 !! - fflx_X diagnostics *should* include all transfers between 247 !! modelled components 248 !! - they should also include gains / losses due to air-sea 249 !! fluxes of C and O2, aeolian and seafloor inputs of Fe, and 250 !! inputs from seafloor "benthic buckets" (N, Si, Fe, C and 251 !! alkalinity) 252 !! - however, they do not include the transfer of material to 253 !! "benthic buckets" by sedimenting slow- and fast-sinking 254 !! detritus since these are separate 2D reservoirs 255 !! - consequently, for a given water column, the integrated 256 !! fluxes should sum to the "loss" of material to the "benthic 257 !! buckets" 258 !! - if they do not, this suggests that MEDUSA contains errors 259 !! in its accounting (e.g. processes omitted from calculated 260 !! fluxes) 261 !! - here, the local integrated fluxes and benthic inputs (plus 262 !! air-sea fluxes in the case of C) are reported together with 263 !! the resulting error 264 !! - only N, Si, C and alkalinity inventories considered; Fe and 265 !! O2 overlooked because of wholesale loss (and addition, in 266 !! the case of O2) of these tracers within the water column 267 !!----------------------------------------------------------------- 268 !! 269 !! nitrogen 270 DO jj = 2,jpjm1 271 DO ji = 2,jpim1 272 if (tmask(ji,jj,1) == 1) then 273 fq0 = fflx_n(ji,jj) 274 fq1 = f_sbenin_n(ji,jj) + f_fbenin_n(ji,jj) 275 fq2 = fq0 + fq1 276 IF (lwp) write (numout,'(a,2i3,a,3f15.10)') & 277 'AXY N cons: (i,j)=',ji,jj,', (flx,ben,err)=', & 278 fq0,fq1,fq2 279 ENDIF 280 ENDDO 281 ENDDO 282 !! silicon 283 DO jj = 2,jpjm1 284 DO ji = 2,jpim1 285 if (tmask(ji,jj,1) == 1) then 286 fq0 = fflx_si(ji,jj) 287 fq1 = f_fbenin_si(ji,jj) 288 fq2 = fq0 + fq1 289 IF (lwp) write (numout,'(a,2i3,a,3f15.10)') & 290 'AXY Si cons: (i,j)=',ji,jj,', (flx,ben,err)=', & 291 fq0,fq1,fq2 292 ENDIF 293 ENDDO 294 ENDDO 295 !! carbon 296 DO jj = 2,jpjm1 297 DO ji = 2,jpim1 298 if (tmask(ji,jj,1) == 1) then 299 fq0 = fflx_c(ji,jj) 300 fq1 = f_sbenin_c(ji,jj) + f_fbenin_c(ji,jj) + f_fbenin_ca(ji,jj) 301 fq2 = f_co2flux(ji,jj) * e3t_n(ji,jj,1) 302 fq3 = fq0 + fq1 303 IF (lwp) write (numout,'(a,2i3,a,4f15.10)') & 304 'AXY C cons: (i,j)=',ji,jj,', (flx,ben,asf,err)=', & 305 fq0,fq1,fq2,fq3 306 ENDIF 307 ENDDO 308 ENDDO 309 !! alkalinity 310 DO jj = 2,jpjm1 311 DO ji = 2,jpim1 312 if (tmask(ji,jj,1) == 1) then 313 fq0 = fflx_a(ji,jj) 314 fq1 = 2.0 * f_fbenin_ca(ji,jj) 315 fq2 = fq0 + fq1 316 IF (lwp) write (numout,'(a,2i3,a,3f15.10)') & 317 'AXY alk cons: (i,j)=',ji,jj,', (flx,ben,err)=', & 318 fq0,fq1,fq2 319 ENDIF 320 ENDDO 321 ENDDO 322 # endif 323 238 324 IF( ln_diatrc ) THEN 239 325 !!-----------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.