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 12109 for NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/TOP/PISCES/P4Z/p4zflx.F90 – NEMO

Ignore:
Timestamp:
2019-12-07T12:40:06+01:00 (4 years ago)
Author:
cetlod
Message:

check out & merge dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags branch onto dev_r12072_TOP-01_ENHANCE-11_CEthe

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/TOP/PISCES/P4Z/p4zflx.F90

    r11993 r12109  
    160160            zfld = zfco2 * chemc(ji,jj,1) * zkgco2(ji,jj)  ! (mol/L) * (m/s) 
    161161            zflu = zh2co3(ji,jj) * zkgco2(ji,jj)                                   ! (mol/L) (m/s) ? 
    162             oce_co2(ji,jj) = ( zfld - zflu ) * rfact2 * e1e2t(ji,jj) * tmask(ji,jj,1) * 1000. 
     162            oce_co2(ji,jj) = ( zfld - zflu ) * tmask(ji,jj,1)  
    163163            ! compute the trend 
    164             tra(ji,jj,1,jpdic) = tra(ji,jj,1,jpdic) + ( zfld - zflu ) * rfact2 / e3t_n(ji,jj,1) * tmask(ji,jj,1) 
     164            tra(ji,jj,1,jpdic) = tra(ji,jj,1,jpdic) + oce_co2(ji,jj) * rfact2 / e3t_n(ji,jj,1) 
    165165 
    166166            ! Compute O2 flux  
     
    174174      IF( iom_use("tcflx") .OR. iom_use("tcflxcum") .OR. kt == nitrst   & 
    175175         &                 .OR. (ln_check_mass .AND. kt == nitend) )    & 
    176          t_oce_co2_flx  = glob_sum( 'p4zflx', oce_co2(:,:) )                    !  Total Flux of Carbon 
     176         t_oce_co2_flx  = glob_sum( 'p4zflx', oce_co2(:,:) * e1e2t(:,:) * 1000. )                    !  Total Flux of Carbon 
    177177      t_oce_co2_flx_cum = t_oce_co2_flx_cum + t_oce_co2_flx       !  Cumulative Total Flux of Carbon 
    178178!      t_atm_co2_flx     = glob_sum( 'p4zflx', satmco2(:,:) * e1e2t(:,:) )       ! Total atmospheric pCO2 
     
    186186 
    187187      IF( lk_iomput .AND. knt == nrdttrc ) THEN 
     188         CALL iom_put( "AtmCo2" , satmco2(:,:) * tmask(:,:,1) )   ! Atmospheric CO2 concentration 
     189         !  
    188190         ALLOCATE( zw2d(jpi,jpj) )   
    189191         IF( iom_use( "Cflx"  ) )  THEN 
    190             zw2d(:,:) = oce_co2(:,:) / e1e2t(:,:) * rfact2r 
     192            zw2d(:,:) = oce_co2(:,:) * 1000.  ! conversion in molC/m2/s 
    191193            CALL iom_put( "Cflx"     , zw2d )  
    192194         ENDIF 
    193195         IF( iom_use( "Oflx"  ) )  THEN 
    194             zw2d(:,:) =  zoflx(:,:) * 1000 * tmask(:,:,1) 
     196            zw2d(:,:) =  zoflx(:,:) * 1000. 
    195197            CALL iom_put( "Oflx" , zw2d ) 
    196198         ENDIF 
     
    203205           CALL iom_put( "Dpco2" ,  zw2d ) 
    204206         ENDIF 
     207         IF( iom_use( "pCO2sea" ) ) THEN 
     208           zw2d(:,:) = ( zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) 
     209           CALL iom_put( "pCO2sea" ,  zw2d ) 
     210         ENDIF 
     211 
    205212         IF( iom_use( "Dpo2" ) )  THEN 
    206213           zw2d(:,:) = ( atcox * patm(:,:) - atcox * trb(:,:,1,jpoxy) / ( chemo2(:,:,1) + rtrn ) ) * tmask(:,:,1) 
    207214           CALL iom_put( "Dpo2"  , zw2d ) 
    208215         ENDIF 
    209          CALL iom_put( "tcflx"    , t_oce_co2_flx * rfact2r )   ! molC/s 
    210          CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum       )   ! molC 
     216         CALL iom_put( "tcflx"    , t_oce_co2_flx     )   ! molC/s 
     217         CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum )   ! molC 
    211218         ! 
    212219         DEALLOCATE( zw2d ) 
Note: See TracChangeset for help on using the changeset viewer.