Changeset 3481 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM
- Timestamp:
- 2012-10-02T14:38:42+02:00 (12 years ago)
- Location:
- branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90
r3443 r3481 57 57 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: satmco2 !: atmospheric pco2 58 58 59 REAL(wp) :: t_oce_co2_flx !: Total ocean carbon flux60 REAL(wp) :: t_atm_co2_flx !: global mean of atmospheric pco261 REAL(wp) :: area !: ocean surface62 59 REAL(wp) :: xconv = 0.01_wp / 3600._wp !: coefficients for conversion 63 60 … … 200 197 201 198 t_oce_co2_flx = t_oce_co2_flx + glob_sum( oce_co2(:,:) ) ! Cumulative Total Flux of Carbon 202 IF( kt == nitend ) THEN 203 t_atm_co2_flx = glob_sum( satmco2(:,:) * patm(:,:) * e1e2t(:,:) ) ! Total atmospheric pCO2 204 ! 205 t_oce_co2_flx = (-1.) * t_oce_co2_flx * 12. / 1.e15 ! Conversion in PgC ; negative for out of the ocean 206 t_atm_co2_flx = t_atm_co2_flx / area ! global mean of atmospheric pCO2 207 ! 208 IF( lwp) THEN 209 WRITE(numout,*) 210 WRITE(numout,*) ' Global mean of atmospheric pCO2 (ppm) at it= ', kt, ' date= ', ndastp 211 WRITE(numout,*) '------------------------------------------------------- : ',t_atm_co2_flx 212 WRITE(numout,*) 213 WRITE(numout,*) ' Cumulative total Flux of Carbon out of the ocean (PgC) :' 214 WRITE(numout,*) '------------------------------------------------------- ',t_oce_co2_flx 215 ENDIF 216 ! 217 ENDIF 199 t_atm_co2_flx = glob_sum( satmco2(:,:) * patm(:,:) * e1e2t(:,:) ) ! Total atmospheric pCO2 218 200 219 201 IF(ln_ctl) THEN ! print mean trends (used for debugging) … … 299 281 ENDIF 300 282 ! 301 area = glob_sum( e1e2t(:,:) ) ! interior global domain surface302 !303 283 oce_co2(:,:) = 0._wp ! Initialization of Flux of Carbon 304 284 t_atm_co2_flx = 0._wp -
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90
r3475 r3481 52 52 REAL(wp) :: texcret !: 1 - excret 53 53 REAL(wp) :: texcret2 !: 1 - excret2 54 REAL(wp) :: tpp !: Total primary production55 54 56 55 … … 409 408 ! Total primary production per year 410 409 tpp = tpp + glob_sum( ( zprorca(:,:,:) + zprorcad(:,:,:) ) * cvol(:,:,:) ) 411 412 IF( kt == nitend .AND. jnt == nrdttrc ) THEN413 IF (lwp) THEN414 WRITE(numout,*) 'Total PP (Gtc) :'415 WRITE(numout,*) '-------------------- : ',tpp * 12. / 1.E12416 WRITE(numout,*)417 ENDIF418 ENDIF419 410 420 411 IF( ln_diatrc ) THEN -
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90
r3471 r3481 238 238 END DO 239 239 240 ! 241 IF( ln_diatrc ) THEN 240 ! Total primary production per year 241 t_oce_co2_exp = t_oce_co2_exp + glob_sum( ( sinking(:,:,iksed+1) + sinking2(:,:,iksed+1) ) * e1e2t(:,:) * tmask(:,:,1) ) 242 ! 243 IF( ln_diatrc ) THEN 242 244 zrfact2 = 1.e3 * rfact2r 243 245 ik1 = iksed + 1 … … 273 275 !! *** ROUTINE p4z_sink_init *** 274 276 !!---------------------------------------------------------------------- 277 278 t_oce_co2_exp = 0._wp 279 ! 275 280 END SUBROUTINE p4z_sink_init 276 281 … … 474 479 END DO 475 480 481 ! Total primary production per year 482 t_oce_co2_exp = t_oce_co2_exp + glob_sum( ( sinking(:,:,:) ) * cvol(:,:,:) ) 483 ! 476 484 IF( ln_diatrc ) THEN 477 485 ! … … 655 663 END DO 656 664 ! 665 t_oce_co2_exp = 0._wp 666 ! 657 667 IF( nn_timing == 1 ) CALL timing_stop('p4z_sink_init') 658 668 ! -
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
r3471 r3481 42 42 INTEGER :: numsil !: logical unit for Si budget 43 43 INTEGER :: numfer !: logical unit for Fer budget 44 INTEGER :: numco2 !: logical unit for co2 budget 44 45 45 46 !!---------------------------------------------------------------------- … … 147 148 END IF 148 149 ! 150 ! 151 IF( kt == nitend ) THEN 152 ! 153 t_atm_co2_flx = t_atm_co2_flx / glob_sum( e1e2t(:,:) ) 154 t_oce_co2_flx = t_oce_co2_flx * 12. / 1.e15 * (-1 ) 155 tpp = tpp * 1000. * 12. / 1.E15 156 t_oce_co2_exp = t_oce_co2_exp * 1000. * 12. / 1.E15 157 ! 158 IF(lwp) WRITE(numco2,9000) ndastp, t_atm_co2_flx, t_oce_co2_flx, tpp, t_oce_co2_exp 159 ! 160 ENDIF 161 149 162 IF( nn_timing == 1 ) CALL timing_stop('p4z_sms') 163 ! 164 9000 FORMAT(i8,f10.5,e18.10,f10.5,f10.5) 150 165 ! 151 166 END SUBROUTINE p4z_sms … … 368 383 CALL ctl_opn( numfer, 'iron.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 369 384 ENDIF 385 IF(lwp) CALL ctl_opn( numco2, 'carbon.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 370 386 ENDIF 371 387 … … 401 417 ENDIF 402 418 ENDIF 419 ! 403 420 9500 FORMAT(i10,e18.10) 404 421 ! -
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/sms_pisces.F90
r3446 r3481 46 46 REAL(wp) :: xkmort !: ??? 47 47 REAL(wp) :: ferat3 !: ??? 48 49 !!* diagnostic parameters 50 REAL(wp) :: tpp !: total primary production 51 REAL(wp) :: t_oce_co2_exp !: total carbon export 52 REAL(wp) :: t_oce_co2_flx !: Total ocean carbon flux 53 REAL(wp) :: t_atm_co2_flx !: global mean of atmospheric pco2 48 54 49 55 !!* Damping
Note: See TracChangeset
for help on using the changeset viewer.