Changeset 1298
- Timestamp:
- 2009-02-09T09:48:20+01:00 (15 years ago)
- Location:
- trunk/NEMO/TOP_SRC/PISCES
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90
r1152 r1298 24 24 USE p4zche 25 25 26 USE lib_mpp 27 26 28 IMPLICIT NONE 27 29 PRIVATE … … 34 36 35 37 REAL(wp) :: & 36 tco2flx = 0.!: Total flux of carbon per year38 tco2flx !: Total flux of carbon per year 37 39 38 40 !!* Substitution … … 184 186 185 187 IF( MOD( kt, nspyr ) == 0 ) THEN 188 IF( lk_mpp ) CALL mpp_sum( tco2flx ) ! sum over the global domain 189 186 190 WRITE(numout,*) ' Atmospheric pCO2 :' 187 191 WRITE(numout,*) '-------------------- : ',kt,' ',atcco2 188 192 WRITE(numout,*) '(ppm)' 189 193 WRITE(numout,*) 'Total Flux of Carbon :' 190 WRITE(numout,*) '-------------------- : ',tco2flx * 12. /1E15194 WRITE(numout,*) '-------------------- : ',tco2flx * 12. / 1e15 191 195 WRITE(numout,*) '(GtC/an)' 192 196 tco2flx = 0. … … 227 231 ENDIF 228 232 233 ! Initialization of total Flux of Carbon 234 tco2flx = 0. 235 229 236 END SUBROUTINE p4z_flx_init 230 237 -
trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90
r1263 r1298 47 47 texcret2 , & !: 1 - excret2 48 48 rpis180 , & !: rpi / 180 49 tpp = 0.!: Total primary production49 tpp !: Total primary production 50 50 51 51 !!* Substitution … … 233 233 END DO 234 234 235 ! Computation of the fractionnal day length236 ! -----------------------------------------237 ! zstrn(:,:) = strn(:,:)238 235 239 236 WHERE( zstrn(:,:) < 1.e0 ) zstrn(:,:) = 24. 240 237 zstrn(:,:) = 24. / zstrn(:,:) 241 ! DO jj = 1, jpj242 ! DO ji = 1, jpi243 !244 ! Computation of the maximum light intensity245 ! ------------------------------------------246 ! IF( zstrn(ji,jj) < 1.e0 ) zstrn(ji,jj) = 24.247 ! zstrn(ji,jj) = 24. / zstrn(ji,jj)248 ! END DO249 ! END DO250 238 251 239 !CDIR NOVERRCHK … … 257 245 258 246 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 259 ! Computation of the various production terms for nanophyto. 260 ! ---------------------------------------------------------- 247 ! Computation of the various production terms for nanophyto. 261 248 zetot2 = enano(ji,jj,jk) * zstrn(ji,jj) 262 249 zmax = MAX( 0.1, xlimphy(ji,jj,jk) ) … … 290 277 DO ji = 1, jpi 291 278 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 292 ! Computation of the various production terms for diatoms 293 ! ------------------------------------------------------- 279 ! Computation of the various production terms for diatoms 294 280 zetot2 = ediat(ji,jj,jk) * zstrn(ji,jj) 295 281 zmax = MAX( 0.1, xlimdia(ji,jj,jk) ) … … 319 305 ! 320 306 321 ! 322 ! Update the arrays TRA which contain the biological sources and sinks 323 ! -------------------------------------------------------------------- 324 307 ! Update the arrays TRA which contain the biological sources and sinks 325 308 DO jk = 1, jpkm1 326 309 DO jj = 1, jpj … … 361 344 zvol = zvol * facvol(ji,jj,jk) 362 345 #endif 363 tpp = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) ) * zvol 346 tpp = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) ) & 347 * zvol * tmask(ji,jj,jk) * tmask_i(ji,jj) 364 348 END DO 365 349 END DO 366 350 END DO 367 351 368 IF( lk_mpp ) CALL mpp_sum( tpp )369 352 370 353 IF( MOD( kt, nspyr ) == 0 ) THEN 354 IF( lk_mpp ) CALL mpp_sum( tpp ) 371 355 WRITE(numout,*) 'Total PP :' 372 356 WRITE(numout,*) '-------------------- : ',tpp * 12. / 1.E12 … … 377 361 #if defined key_trc_dia3d 378 362 zrfact2 = 1.e3 * rfact2r 379 ! Supplementary diagnostics 380 ! ------------------------- 363 ! Supplementary diagnostics 381 364 trc3d(:,:,:,jp_pcs0_3d + 3) = etot(:,:,:) 382 365 trc3d(:,:,:,jp_pcs0_3d + 4) = zprorca(:,:,:) * zrfact2 … … 437 420 texcret = 1.0 - excret 438 421 texcret2 = 1.0 - excret2 422 tpp = 0. 439 423 440 424 END SUBROUTINE p4z_prod_init
Note: See TracChangeset
for help on using the changeset viewer.