New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6225 for branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90 – NEMO

Ignore:
Timestamp:
2016-01-08T10:35:19+01:00 (8 years ago)
Author:
jamesharle
Message:

Update MPP_BDY_UPDATE branch to be consistent with head of trunk

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

    r4624 r6225  
    1111   !!   'key_pisces'                                       PISCES bio-model 
    1212   !!---------------------------------------------------------------------- 
    13    !!   p4zsms        :  Time loop of passive tracers sms 
     13   !!   p4zsms         :  Time loop of passive tracers sms 
    1414   !!---------------------------------------------------------------------- 
    1515   USE oce_trc         !  shared variables between ocean and passive tracers 
    … …  
    2424   USE p4zsed          !  Sedimentation 
    2525   USE p4zint          !  time interpolation 
     26   USE p4zrem          !  remineralisation 
    2627   USE iom             !  I/O manager 
    27    USE trdmod_oce      !  Ocean trends variables 
    28    USE trdmod_trc      !  TOP trends variables 
     28   USE trd_oce         !  Ocean trends variables 
     29   USE trdtrc          !  TOP trends variables 
    2930   USE sedmodel        !  Sediment model 
    3031   USE prtctl_trc      !  print control for debugging 
    … …  
    3334   PRIVATE 
    3435 
    35    PUBLIC   p4z_sms_init    ! called in p4zsms.F90 
    36    PUBLIC   p4z_sms    ! called in p4zsms.F90 
    37  
    38    REAL(wp) :: alkbudget, no3budget, silbudget, ferbudget 
    39    INTEGER ::  numco2, numnut  !: logical unit for co2 budget 
     36   PUBLIC   p4z_sms_init   ! called in p4zsms.F90 
     37   PUBLIC   p4z_sms        ! called in p4zsms.F90 
     38 
     39   REAL(wp) :: alkbudget, no3budget, silbudget, ferbudget, po4budget 
     40   REAL(wp) :: xfact1, xfact2, xfact3 
     41   INTEGER ::  numco2, numnut, numnit  !: logical unit for co2 budget 
     42 
     43   !!* Array used to indicate negative tracer values 
     44   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   xnegtr     !: ??? 
     45 
    4046 
    4147   !!---------------------------------------------------------------------- 
    … …  
    6167      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    6268      !! 
    63       INTEGER ::   jnt, jn, jl 
     69      INTEGER ::   ji, jj, jk, jnt, jn, jl 
     70      REAL(wp) ::  ztra 
     71#if defined key_kriest 
     72      REAL(wp) ::  zcoef1, zcoef2 
     73#endif 
    6474      CHARACTER (len=25) :: charout 
    65       REAL(wp), POINTER, DIMENSION(:,:,:,:)  :: ztrdpis 
    6675      !!--------------------------------------------------------------------- 
    6776      ! 
    6877      IF( nn_timing == 1 )  CALL timing_start('p4z_sms') 
    6978      ! 
    70       IF( l_trdtrc )  THEN 
    71          CALL wrk_alloc( jpi, jpj, jpk, jp_pisces, ztrdpis )  
    72          DO jn = 1, jp_pisces 
    73             jl = jn + jp_pcs0 - 1 
    74             ztrdpis(:,:,:,jn) = trn(:,:,:,jl) 
    75          ENDDO 
    76       ENDIF 
    77       ! 
    7879      IF( kt == nittrc000 ) THEN 
     80        ! 
     81        ALLOCATE( xnegtr(jpi,jpj,jpk) ) 
    7982        ! 
    8083        CALL p4z_che                              ! initialize the chemical constants 
    … …  
    8891      IF( ln_pisdmp .AND. MOD( kt - nn_dttrc, nn_pisdmp ) == 0 )   CALL p4z_dmp( kt )      ! Relaxation of some tracers 
    8992      ! 
     93      !                                                                    !   set time step size (Euler/Leapfrog) 
     94      IF( ( neuler == 0 .AND. kt == nittrc000 ) .OR. ln_top_euler ) THEN   ;    rfact = rdttrc     !  at nittrc000 
     95      ELSEIF( kt <= nittrc000 + nn_dttrc )                          THEN   ;    rfact = 2. * rdttrc   ! at nittrc000 or nittrc000+nn_dttrc (Leapfrog) 
     96      ENDIF 
     97      ! 
     98      IF( ( ln_top_euler .AND. kt == nittrc000 )  .OR. ( .NOT.ln_top_euler .AND. kt <= nittrc000 + nn_dttrc ) ) THEN 
     99         rfactr  = 1. / rfact 
     100         rfact2  = rfact / FLOAT( nrdttrc ) 
     101         rfact2r = 1. / rfact2 
     102         xstep = rfact2 / rday         ! Time step duration for biology 
     103         IF(lwp) WRITE(numout,*)  
     104         IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rdt = ', rdt 
     105         IF(lwp) write(numout,*) '    PISCES  Biology time step    rfact2 = ', rfact2 
     106         IF(lwp) WRITE(numout,*) 
     107      ENDIF 
     108 
     109      IF( ( neuler == 0 .AND. kt == nittrc000 ) .OR. ln_top_euler ) THEN 
     110         DO jn = jp_pcs0, jp_pcs1              !   SMS on tracer without Asselin time-filter 
     111            trb(:,:,:,jn) = trn(:,:,:,jn) 
     112         END DO 
     113      ENDIF 
     114      ! 
    90115      IF( ndayflxtr /= nday_year ) THEN      ! New days 
    91116         ! 
    … …  
    105130      DO jnt = 1, nrdttrc          ! Potential time splitting if requested 
    106131         ! 
    107          CALL p4z_bio (kt, jnt)    ! Biology 
    108          CALL p4z_sed (kt, jnt)    ! Sedimentation 
    109          ! 
     132         CALL p4z_bio( kt, jnt )   ! Biology 
     133         CALL p4z_sed( kt, jnt )   ! Sedimentation 
     134         CALL p4z_lys( kt, jnt )   ! Compute CaCO3 saturation 
     135         CALL p4z_flx( kt, jnt )   ! Compute surface fluxes 
     136         ! 
     137         xnegtr(:,:,:) = 1.e0 
    110138         DO jn = jp_pcs0, jp_pcs1 
    111             trb(:,:,:,jn) = trn(:,:,:,jn) 
    112          ENDDO 
    113          ! 
     139            DO jk = 1, jpk 
     140               DO jj = 1, jpj 
     141                  DO ji = 1, jpi 
     142                     IF( ( trb(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 ) THEN 
     143                        ztra             = ABS( trb(ji,jj,jk,jn) ) / ( ABS( tra(ji,jj,jk,jn) ) + rtrn ) 
     144                        xnegtr(ji,jj,jk) = MIN( xnegtr(ji,jj,jk),  ztra ) 
     145                     ENDIF 
     146                 END DO 
     147               END DO 
     148            END DO 
     149         END DO 
     150         !                                ! where at least 1 tracer concentration becomes negative 
     151         !                                !  
     152         DO jn = jp_pcs0, jp_pcs1 
     153           trb(:,:,:,jn) = trb(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn) 
     154         END DO 
     155        ! 
     156         DO jn = jp_pcs0, jp_pcs1 
     157            tra(:,:,:,jn) = 0._wp 
     158         END DO 
     159         ! 
     160         IF( ln_top_euler ) THEN 
     161            DO jn = jp_pcs0, jp_pcs1 
     162               trn(:,:,:,jn) = trb(:,:,:,jn) 
     163            END DO 
     164         ENDIF 
    114165      END DO 
    115166 
    116       IF( l_trdtrc )  THEN 
    117          DO jn = 1, jp_pisces 
    118             jl = jn + jp_pcs0 - 1 
    119             ztrdpis(:,:,:,jn) = ( ztrdpis(:,:,:,jn) - trn(:,:,:,jl) ) * rfact2r 
    120          ENDDO 
    121       ENDIF 
    122  
    123       CALL p4z_lys( kt )             ! Compute CaCO3 saturation 
    124       CALL p4z_flx( kt )             ! Compute surface fluxes 
    125  
    126       DO jn = jp_pcs0, jp_pcs1 
    127         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. ) 
    128         CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. ) 
    129         CALL lbc_lnk( tra(:,:,:,jn), 'T', 1. ) 
     167#if defined key_kriest 
     168      !  
     169      zcoef1 = 1.e0 / xkr_massp  
     170      zcoef2 = 1.e0 / xkr_massp / 1.1 
     171      DO jk = 1,jpkm1 
     172         trb(:,:,jk,jpnum) = MAX(  trb(:,:,jk,jpnum), trb(:,:,jk,jppoc) * zcoef1 / xnumm(jk)  ) 
     173         trb(:,:,jk,jpnum) = MIN(  trb(:,:,jk,jpnum), trb(:,:,jk,jppoc) * zcoef2              ) 
    130174      END DO 
    131175      ! 
     176#endif 
     177      ! 
     178      ! 
     179      IF( l_trdtrc ) THEN 
     180         DO jn = jp_pcs0, jp_pcs1 
     181           CALL trd_trc( tra(:,:,:,jn), jn, jptra_sms, kt )   ! save trends 
     182         END DO 
     183      END IF 
     184      ! 
    132185      IF( lk_sed ) THEN  
    133186         ! 
    … …  
    135188         ! 
    136189         DO jn = jp_pcs0, jp_pcs1 
    137            CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. ) 
     190           CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. ) 
    138191         END DO 
    139192         ! 
    … …  
    142195      IF( lrst_trc )  CALL p4z_rst( kt, 'WRITE' )  !* Write PISCES informations in restart file  
    143196      ! 
    144       IF( l_trdtrc ) THEN 
    145          DO jn = 1, jp_pisces 
    146             jl = jn + jp_pcs0 - 1 
    147              ztrdpis(:,:,:,jn) = ztrdpis(:,:,:,jn) + tra(:,:,:,jl) 
    148              CALL trd_mod_trc( ztrdpis(:,:,:,jn), jn, jptra_trd_sms, kt )   ! save trends 
    149           END DO 
    150           CALL wrk_dealloc( jpi, jpj, jpk, jp_pisces, ztrdpis )  
    151       END IF 
    152       ! 
    153       CALL p4z_chk_mass( kt ) ! Mass conservation checking 
     197 
     198      IF( lk_iomput .OR. ln_check_mass )  CALL p4z_chk_mass( kt ) ! Mass conservation checking 
    154199 
    155200      IF ( lwm .AND. kt == nittrc000 ) CALL FLUSH    ( numonp )     ! flush output namelist PISCES 
    … …  
    281326               ztmas   = tmask(ji,jj,jk) 
    282327               ztmas1  = 1. - tmask(ji,jj,jk) 
    283                zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  ) 
    284                zco3    = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1 
    285                zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk ) 
     328               zcaralk = trb(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  ) 
     329               zco3    = ( zcaralk - trb(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1 
     330               zbicarb = ( 2. * trb(ji,jj,jk,jpdic) - zcaralk ) 
    286331               hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1 
    287332            END DO 
    … …  
    328373         ENDIF 
    329374         ! 
     375         IF( iom_varid( numrtr, 'tcflxcum', ldstop = .FALSE. ) > 0 ) THEN  ! cumulative total flux of carbon 
     376            CALL iom_get( numrtr, 'tcflxcum' , t_oce_co2_flx_cum  ) 
     377         ELSE 
     378            t_oce_co2_flx_cum = 0._wp 
     379         ENDIF 
     380         ! 
    330381      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN 
    331382         IF( kt == nitrst ) THEN 
    … …  
    337388         CALL iom_rstput( kt, nitrst, numrtw, 'Silicalim', xksi(:,:) ) 
    338389         CALL iom_rstput( kt, nitrst, numrtw, 'Silicamax', xksimax(:,:) ) 
     390         CALL iom_rstput( kt, nitrst, numrtw, 'tcflxcum', t_oce_co2_flx_cum ) 
    339391      ENDIF 
    340392      ! 
    … …  
    355407      REAL(wp) ::  silmean = 91.51     ! mean value of silicate 
    356408      ! 
    357       REAL(wp) :: zarea, zalksum, zpo4sum, zno3sum, zsilsum 
     409      REAL(wp) :: zarea, zalksumn, zpo4sumn, zno3sumn, zsilsumn 
     410      REAL(wp) :: zalksumb, zpo4sumb, zno3sumb, zsilsumb 
    358411      !!--------------------------------------------------------------------- 
    359412 
    … …  
    368421         zarea          = 1._wp / glob_sum( cvol(:,:,:) ) * 1e6               
    369422 
    370          zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
    371          zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
    372          zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
    373          zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
     423         zalksumn = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
     424         zpo4sumn = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
     425         zno3sumn = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
     426         zsilsumn = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
    374427  
    375          IF(lwp) WRITE(numout,*) '       TALK mean : ', zalksum 
    376          trn(:,:,:,jptal) = trn(:,:,:,jptal) * alkmean / zalksum 
    377  
    378          IF(lwp) WRITE(numout,*) '       PO4  mean : ', zpo4sum 
    379          trn(:,:,:,jppo4) = trn(:,:,:,jppo4) * po4mean / zpo4sum 
    380  
    381          IF(lwp) WRITE(numout,*) '       NO3  mean : ', zno3sum 
    382          trn(:,:,:,jpno3) = trn(:,:,:,jpno3) * no3mean / zno3sum 
    383  
    384          IF(lwp) WRITE(numout,*) '       SiO3 mean : ', zsilsum 
    385          trn(:,:,:,jpsil) = MIN( 400.e-6,trn(:,:,:,jpsil) * silmean / zsilsum ) 
    386          ! 
    387       ENDIF 
    388  
     428         IF(lwp) WRITE(numout,*) '       TALKN mean : ', zalksumn 
     429         trn(:,:,:,jptal) = trn(:,:,:,jptal) * alkmean / zalksumn 
     430 
     431         IF(lwp) WRITE(numout,*) '       PO4N  mean : ', zpo4sumn 
     432         trn(:,:,:,jppo4) = trn(:,:,:,jppo4) * po4mean / zpo4sumn 
     433 
     434         IF(lwp) WRITE(numout,*) '       NO3N  mean : ', zno3sumn 
     435         trn(:,:,:,jpno3) = trn(:,:,:,jpno3) * no3mean / zno3sumn 
     436 
     437         IF(lwp) WRITE(numout,*) '       SiO3N mean : ', zsilsumn 
     438         trn(:,:,:,jpsil) = MIN( 400.e-6,trn(:,:,:,jpsil) * silmean / zsilsumn ) 
     439         ! 
     440         ! 
     441         IF( .NOT. ln_top_euler ) THEN 
     442            zalksumb = glob_sum( trb(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
     443            zpo4sumb = glob_sum( trb(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
     444            zno3sumb = glob_sum( trb(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
     445            zsilsumb = glob_sum( trb(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
     446  
     447            IF(lwp) WRITE(numout,*) ' ' 
     448            IF(lwp) WRITE(numout,*) '       TALKB mean : ', zalksumb 
     449            trb(:,:,:,jptal) = trb(:,:,:,jptal) * alkmean / zalksumb 
     450 
     451            IF(lwp) WRITE(numout,*) '       PO4B  mean : ', zpo4sumb 
     452            trb(:,:,:,jppo4) = trb(:,:,:,jppo4) * po4mean / zpo4sumb 
     453 
     454            IF(lwp) WRITE(numout,*) '       NO3B  mean : ', zno3sumb 
     455            trb(:,:,:,jpno3) = trb(:,:,:,jpno3) * no3mean / zno3sumb 
     456 
     457            IF(lwp) WRITE(numout,*) '       SiO3B mean : ', zsilsumb 
     458            trb(:,:,:,jpsil) = MIN( 400.e-6,trb(:,:,:,jpsil) * silmean / zsilsumb ) 
     459        ENDIF 
     460        ! 
     461      ENDIF 
     462        ! 
    389463   END SUBROUTINE p4z_dmp 
    390464 
    … …  
    399473      ! 
    400474      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    401       !! 
     475      REAL(wp)             ::  zrdenittot, zsdenittot, znitrpottot 
     476      CHARACTER(LEN=100)   ::   cltxt 
     477      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zvol 
     478      INTEGER :: jk 
     479      !!---------------------------------------------------------------------- 
     480 
     481      ! 
    402482      !!--------------------------------------------------------------------- 
    403483 
    … …  
    406486            CALL ctl_opn( numco2, 'carbon.budget'  , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
    407487            CALL ctl_opn( numnut, 'nutrient.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
     488            CALL ctl_opn( numnit, 'nitrogen.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
     489            xfact1 = rfact2r * 12. / 1.e15 * ryyss    ! conversion molC/kt --> PgC/yr 
     490            xfact2 = 1.e+3 * rno3 * 14. / 1.e12 * ryyss   ! conversion molC/l/s ----> TgN/m3/yr 
     491            xfact3 = 1.e+3 * rfact2r * rno3   ! conversion molC/l/kt ----> molN/m3/s 
     492            cltxt='time-step   Alkalinity        Nitrate        Phosphorus         Silicate           Iron' 
     493            IF( lwp ) WRITE(numnut,*)  TRIM(cltxt) 
     494            IF( lwp ) WRITE(numnut,*)  
    408495         ENDIF 
    409496      ENDIF 
    410497 
    411       IF( ln_check_mass .AND. kt == nitend ) THEN      !   Compute the budget of NO3, ALK, Si, Fer 
     498      ! 
     499      IF( iom_use( "pno3tot" ) .OR. ( ln_check_mass .AND. kt == nitend )  ) THEN 
     500         !   Compute the budget of NO3, ALK, Si, Fer 
    412501         no3budget = glob_sum( (   trn(:,:,:,jpno3) + trn(:,:,:,jpnh4)  & 
    413502            &                    + trn(:,:,:,jpphy) + trn(:,:,:,jpdia)  & 
    … …  
    417506            &                    + trn(:,:,:,jpgoc)                     & 
    418507#endif 
    419             &                    + trn(:,:,:,jpdoc)                     ) * cvol(:,:,:)  )  
    420          !  
     508            &                    + trn(:,:,:,jpdoc)                     ) * cvol(:,:,:)  ) 
     509         ! 
     510         no3budget = no3budget / areatot 
     511         CALL iom_put( "pno3tot", no3budget ) 
     512      ENDIF 
     513      ! 
     514      IF( iom_use( "ppo4tot" ) .OR. ( ln_check_mass .AND. kt == nitend )  ) THEN 
     515         po4budget = glob_sum( (   trn(:,:,:,jppo4)                     & 
     516            &                    + trn(:,:,:,jpphy) + trn(:,:,:,jpdia)  & 
     517            &                    + trn(:,:,:,jpzoo) + trn(:,:,:,jpmes)  & 
     518            &                    + trn(:,:,:,jppoc)                     & 
     519#if ! defined key_kriest 
     520            &                    + trn(:,:,:,jpgoc)                     & 
     521#endif 
     522            &                    + trn(:,:,:,jpdoc)                     ) * cvol(:,:,:)  ) 
     523         po4budget = po4budget / areatot 
     524         CALL iom_put( "ppo4tot", po4budget ) 
     525      ENDIF 
     526      ! 
     527      IF( iom_use( "psiltot" ) .OR. ( ln_check_mass .AND. kt == nitend )  ) THEN 
    421528         silbudget = glob_sum( (   trn(:,:,:,jpsil) + trn(:,:,:,jpgsi)  & 
    422529            &                    + trn(:,:,:,jpdsi)                     ) * cvol(:,:,:)  ) 
    423          !  
     530         ! 
     531         silbudget = silbudget / areatot 
     532         CALL iom_put( "psiltot", silbudget ) 
     533      ENDIF 
     534      ! 
     535      IF( iom_use( "palktot" ) .OR. ( ln_check_mass .AND. kt == nitend )  ) THEN 
    424536         alkbudget = glob_sum( (   trn(:,:,:,jpno3) * rno3              & 
    425537            &                    + trn(:,:,:,jptal)                     & 
    426538            &                    + trn(:,:,:,jpcal) * 2.                ) * cvol(:,:,:)  ) 
    427          !  
     539         ! 
     540         alkbudget = alkbudget / areatot 
     541         CALL iom_put( "palktot", alkbudget ) 
     542      ENDIF 
     543      ! 
     544      IF( iom_use( "pfertot" ) .OR. ( ln_check_mass .AND. kt == nitend )  ) THEN 
    428545         ferbudget = glob_sum( (   trn(:,:,:,jpfer) + trn(:,:,:,jpnfe)  & 
    429546            &                    + trn(:,:,:,jpdfe)                     & 
    … …  
    434551            &                    + trn(:,:,:,jpzoo) * ferat3            & 
    435552            &                    + trn(:,:,:,jpmes) * ferat3            ) * cvol(:,:,:)  ) 
    436  
    437          ! 
     553         ! 
     554         ferbudget = ferbudget / areatot 
     555         CALL iom_put( "pfertot", ferbudget ) 
     556      ENDIF 
     557      ! 
     558 
     559      ! Global budget of N SMS : denitrification in the water column and in the sediment 
     560      !                          nitrogen fixation by the diazotrophs 
     561      ! -------------------------------------------------------------------------------- 
     562      IF( iom_use( "tnfix" ) .OR.  ( ln_check_mass .AND. kt == nitend )  ) THEN 
     563         znitrpottot  = glob_sum ( nitrpot(:,:,:) * nitrfix * cvol(:,:,:) ) 
     564         CALL iom_put( "tnfix"  , znitrpottot * 1.e+3 * rno3 )  ! Global  nitrogen fixation molC/l  to molN/m3  
     565      ENDIF 
     566      ! 
     567      IF( iom_use( "tdenit" ) .OR.  ( ln_check_mass .AND. kt == nitend )  ) THEN 
     568         zrdenittot   = glob_sum ( denitr(:,:,:) * rdenit * xnegtr(:,:,:) * cvol(:,:,:) ) 
     569         CALL iom_put( "tdenit"  , zrdenittot * 1.e+3 * rno3 )  ! Total denitrification molC/l to molN/m3  
     570      ENDIF 
     571      ! 
     572      IF( iom_use( "Sdenit" ) .OR.  ( ln_check_mass .AND. kt == nitend )  ) THEN 
     573         zsdenittot   = glob_sum ( sdenit(:,:) * e1e2t(:,:) ) 
     574         CALL iom_put( "Sdenit", sdenit(:,:) * xfact3 * tmask(:,:,1) )  ! Nitrate reduction in the sediments 
     575      ENDIF 
     576 
     577      IF( ln_check_mass .AND. kt == nitend ) THEN   ! Compute the budget of NO3, ALK, Si, Fer 
    438578         t_atm_co2_flx  = t_atm_co2_flx / glob_sum( e1e2t(:,:) ) 
    439          t_oce_co2_flx  = t_oce_co2_flx         * 12. / 1.e15 * (-1 ) 
    440          tpp            = tpp           * 1000. * 12. / 1.E15 
    441          t_oce_co2_exp  = t_oce_co2_exp * 1000. * 12. / 1.E15 
    442          ! 
    443          no3budget = no3budget / areatot 
    444          silbudget = silbudget / areatot 
    445          alkbudget = alkbudget / areatot 
    446          ferbudget = ferbudget / areatot 
    447          ! 
    448          IF(lwp) THEN 
    449             WRITE(numco2,9000) ndastp, t_atm_co2_flx, t_oce_co2_flx, tpp, t_oce_co2_exp 
    450             WRITE(numnut,9500) ndastp, alkbudget, no3budget, silbudget, ferbudget 
    451          ENDIF 
    452          ! 
    453       ENDIF 
    454        ! 
     579         t_oce_co2_flx  = t_oce_co2_flx         * xfact1 * (-1 ) 
     580         tpp            = tpp           * 1000. * xfact1 
     581         t_oce_co2_exp  = t_oce_co2_exp * 1000. * xfact1 
     582         IF( lwp ) WRITE(numco2,9000) ndastp, t_atm_co2_flx, t_oce_co2_flx, tpp, t_oce_co2_exp 
     583         IF( lwp ) WRITE(numnut,9100) ndastp, alkbudget        * 1.e+06, & 
     584             &                                no3budget * rno3 * 1.e+06, & 
     585             &                                po4budget * po4r * 1.e+06, & 
     586             &                                silbudget        * 1.e+06, & 
     587             &                                ferbudget        * 1.e+09 
     588         ! 
     589         IF( lwp ) WRITE(numnit,9200) ndastp, znitrpottot * xfact2  , & 
     590         &                             zrdenittot  * xfact2  , & 
     591         &                             zsdenittot  * xfact2 
     592 
     593      ENDIF 
     594      ! 
    455595 9000  FORMAT(i8,f10.5,e18.10,f10.5,f10.5) 
    456  9500  FORMAT(i8,4e18.10)      
     596 9100  FORMAT(i8,5e18.10) 
     597 9200  FORMAT(i8,3f10.5) 
     598 
    457599       ! 
    458600   END SUBROUTINE p4z_chk_mass 
Note: See TracChangeset for help on using the changeset viewer.