Changeset 12109 for NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/TOP/PISCES/P4Z/p4zflx.F90
- Timestamp:
- 2019-12-07T12:40:06+01:00 (4 years ago)
- 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 160 160 zfld = zfco2 * chemc(ji,jj,1) * zkgco2(ji,jj) ! (mol/L) * (m/s) 161 161 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) 163 163 ! 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) 165 165 166 166 ! Compute O2 flux … … 174 174 IF( iom_use("tcflx") .OR. iom_use("tcflxcum") .OR. kt == nitrst & 175 175 & .OR. (ln_check_mass .AND. kt == nitend) ) & 176 t_oce_co2_flx = glob_sum( 'p4zflx', oce_co2(:,:) ) ! Total Flux of Carbon176 t_oce_co2_flx = glob_sum( 'p4zflx', oce_co2(:,:) * e1e2t(:,:) * 1000. ) ! Total Flux of Carbon 177 177 t_oce_co2_flx_cum = t_oce_co2_flx_cum + t_oce_co2_flx ! Cumulative Total Flux of Carbon 178 178 ! t_atm_co2_flx = glob_sum( 'p4zflx', satmco2(:,:) * e1e2t(:,:) ) ! Total atmospheric pCO2 … … 186 186 187 187 IF( lk_iomput .AND. knt == nrdttrc ) THEN 188 CALL iom_put( "AtmCo2" , satmco2(:,:) * tmask(:,:,1) ) ! Atmospheric CO2 concentration 189 ! 188 190 ALLOCATE( zw2d(jpi,jpj) ) 189 191 IF( iom_use( "Cflx" ) ) THEN 190 zw2d(:,:) = oce_co2(:,:) / e1e2t(:,:) * rfact2r192 zw2d(:,:) = oce_co2(:,:) * 1000. ! conversion in molC/m2/s 191 193 CALL iom_put( "Cflx" , zw2d ) 192 194 ENDIF 193 195 IF( iom_use( "Oflx" ) ) THEN 194 zw2d(:,:) = zoflx(:,:) * 1000 * tmask(:,:,1)196 zw2d(:,:) = zoflx(:,:) * 1000. 195 197 CALL iom_put( "Oflx" , zw2d ) 196 198 ENDIF … … 203 205 CALL iom_put( "Dpco2" , zw2d ) 204 206 ENDIF 207 IF( iom_use( "pCO2sea" ) ) THEN 208 zw2d(:,:) = ( zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) 209 CALL iom_put( "pCO2sea" , zw2d ) 210 ENDIF 211 205 212 IF( iom_use( "Dpo2" ) ) THEN 206 213 zw2d(:,:) = ( atcox * patm(:,:) - atcox * trb(:,:,1,jpoxy) / ( chemo2(:,:,1) + rtrn ) ) * tmask(:,:,1) 207 214 CALL iom_put( "Dpo2" , zw2d ) 208 215 ENDIF 209 CALL iom_put( "tcflx" , t_oce_co2_flx * rfact2r) ! molC/s210 CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum 216 CALL iom_put( "tcflx" , t_oce_co2_flx ) ! molC/s 217 CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum ) ! molC 211 218 ! 212 219 DEALLOCATE( zw2d )
Note: See TracChangeset
for help on using the changeset viewer.