Ignore:
Timestamp:
2017-08-08T14:11:22+02:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — 08-08-2017 — clean trcbio_medusa split only - running this time

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_split_trcbiomedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r8395 r8419  
    3232      !!---------------------------------------------------------------------- 
    3333      USE bio_medusa_mod 
    34       USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt, e3t_n, tmask 
     34      USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt 
    3535      USE in_out_manager,    ONLY: lwp, numout 
    3636# if defined key_iomput 
     
    6161      INTEGER :: ji, jj 
    6262      INTEGER :: jn 
    63  
    64       REAL(wp) :: fq0,fq1,fq2,fq3 
    65  
    6663      !!---------------------------------------------------------------------- 
    6764      !! Process benthic in/out fluxes 
     
    238235         ENDIF 
    239236      endif 
    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 
    323237       
    324238      IF( ln_diatrc ) THEN 
Note: See TracChangeset for help on using the changeset viewer.