Changeset 3581 for LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90
- Timestamp:
- Oct 10, 2019, 2:35:59 PM (5 years ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90 (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90
r3459 r3581 7 7 ! ----------------------- 8 8 ! Control module for the carbon CO2 tracers : 9 ! - Identification 10 ! - Get concentrations comming from coupled model or read from file to tracers 11 ! - Calculate new RCO2 for radiation scheme 12 ! - Calculate new carbon flux for sending to coupled models (PISCES and ORCHIDEE) 13 ! 14 ! Module permettant de mettre a jour les champs (puits et sources) pour le 15 ! transport de CO2 en online (IPSL-CM et LMDZOR) et offline (lecture de carte) 9 ! - Initialisation of carbon cycle fields 10 ! - Definition of fluxes to be exchanged 11 ! 12 ! Rest of code is in tracco2i.F90 16 13 ! 17 14 ! Le cas online/offline est defini par le flag carbon_cycle_cpl (y/n) … … 35 32 SAVE 36 33 PRIVATE 37 PUBLIC :: carbon_cycle_init, carbon_cycle,infocfields_init34 PUBLIC :: carbon_cycle_init, infocfields_init 38 35 39 36 ! Variables read from parmeter file physiq.def … … 46 43 INTEGER, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3 47 44 !$OMP THREADPRIVATE(level_coupling_esm) 48 REAL, PUBLIC :: RCO2_glo, RCO2_tot 49 !$OMP THREADPRIVATE(RCO2_glo, RCO2_tot) 45 REAL, PUBLIC :: RCO2_glo 46 !$OMP THREADPRIVATE(RCO2_glo) 47 REAL, PUBLIC :: RCO2_tot 48 !$OMP THREADPRIVATE(RCO2_tot) 50 49 51 50 LOGICAL :: carbon_cycle_emis_comp_omp=.FALSE. … … 84 83 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 85 84 !$OMP THREADPRIVATE(fco2_land) 85 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_nbp ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 86 !$OMP THREADPRIVATE(fco2_land_nbp) 87 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_nep ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 88 !$OMP THREADPRIVATE(fco2_land_nep) 89 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_fLuc ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 90 !$OMP THREADPRIVATE(fco2_land_fLuc) 91 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_fwoodharvest ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 92 !$OMP THREADPRIVATE(fco2_land_fwoodharvest) 93 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_fHarvest ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 94 !$OMP THREADPRIVATE(fco2_land_fHarvest) 86 95 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocean ! Net flux from ocean [kgCO2/m2/s] 87 96 !$OMP THREADPRIVATE(fco2_ocean) … … 191 200 CONTAINS 192 201 193 SUBROUTINE carbon_cycle_init( tr_seri, pdtphys, aerosol, radio)202 SUBROUTINE carbon_cycle_init() 194 203 ! This subroutine is called from traclmdz_init, only at first timestep. 195 204 ! - Read controle parameters from .def input file … … 199 208 200 209 USE dimphy 201 USE geometry_mod, ONLY : cell_area202 USE mod_phys_lmdz_transfert_para203 USE infotrac_phy, ONLY: nbtr, nqo, niadv, tname204 210 USE IOIPSL 205 USE surface_data, ONLY : ok_veget, type_ocean206 USE phys_cal_mod, ONLY : mth_len207 211 USE print_control_mod, ONLY: lunout 208 212 … … 210 214 INCLUDE "clesphys.h" 211 215 212 ! Input argument213 REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: tr_seri ! Concentration Traceur [U/KgA]214 REAL,INTENT(IN) :: pdtphys ! length of time step in physiq (sec)215 216 ! InOutput arguments217 LOGICAL,DIMENSION(nbtr), INTENT(INOUT) :: aerosol218 LOGICAL,DIMENSION(nbtr), INTENT(INOUT) :: radio219 220 216 ! Local variables 221 INTEGER :: ierr, it, iiq, itc 222 INTEGER :: teststop 223 224 ! 1) Read controle parameters from .def input file 225 ! ------------------------------------------------ 226 ! Read fosil fuel value if no transport 227 IF (.NOT. carbon_cycle_tr) THEN 228 !$OMP MASTER 229 fos_fuel_s_omp = 0. 230 CALL getin ('carbon_cycle_fos_fuel',fos_fuel_s_omp) 231 !$OMP END MASTER 232 !$OMP BARRIER 233 fos_fuel_s=fos_fuel_s_omp 234 WRITE(lunout,*) 'carbon_cycle_fos_fuel = ', fos_fuel_s 235 END IF 236 237 ! Read parmeter for calculation compatible emission 238 IF (.NOT. carbon_cycle_tr) THEN 239 !$OMP MASTER 240 carbon_cycle_emis_comp_omp=.FALSE. 241 CALL getin('carbon_cycle_emis_comp',carbon_cycle_emis_comp_omp) 242 !$OMP END MASTER 243 !$OMP BARRIER 244 carbon_cycle_emis_comp=carbon_cycle_emis_comp_omp 245 WRITE(lunout,*) 'carbon_cycle_emis_comp = ',carbon_cycle_emis_comp 246 IF (carbon_cycle_emis_comp) THEN 247 CALL abort_physic('carbon_cycle_init', 'carbon_cycle_emis_comp option not yet implemented!!',1) 248 END IF 249 END IF 250 251 ! Read parameter for interactive calculation of the CO2 value for the radiation scheme 252 !$OMP MASTER 253 RCO2_inter_omp=.FALSE. 254 CALL getin('RCO2_inter',RCO2_inter_omp) 255 !$OMP END MASTER 256 !$OMP BARRIER 257 RCO2_inter=RCO2_inter_omp 258 WRITE(lunout,*) 'RCO2_inter = ', RCO2_inter 259 IF (RCO2_inter) THEN 260 WRITE(lunout,*) 'RCO2 will be recalculated once a day' 261 WRITE(lunout,*) 'RCO2 initial = ', RCO2 262 END IF 263 264 265 ! 2) Search for carbon tracers and set default values 266 ! --------------------------------------------------- 267 itc=0 268 DO it=1,nbtr 269 !! iiq=niadv(it+2) ! jyg 270 iiq=niadv(it+nqo) ! jyg 271 272 SELECT CASE(tname(iiq)) 273 CASE("fCO2_ocn") 274 itc = itc + 1 275 co2trac(itc)%name='fCO2_ocn' 276 co2trac(itc)%id=it 277 co2trac(itc)%file='fl_co2_ocean.nc' 278 IF (carbon_cycle_cpl .AND. type_ocean=='couple') THEN 279 co2trac(itc)%cpl=.TRUE. 280 co2trac(itc)%updatefreq = 86400 ! Once a day as the coupling with OASIS/PISCES 281 ELSE 282 co2trac(itc)%cpl=.FALSE. 283 co2trac(itc)%updatefreq = 86400*mth_len ! Once a month 284 END IF 285 CASE("fCO2_land") 286 itc = itc + 1 287 co2trac(itc)%name='fCO2_land' 288 co2trac(itc)%id=it 289 co2trac(itc)%file='fl_co2_land.nc' 290 IF (carbon_cycle_cpl .AND. ok_veget) THEN 291 co2trac(itc)%cpl=.TRUE. 292 co2trac(itc)%updatefreq = INT(pdtphys) ! Each timestep as the coupling with ORCHIDEE 293 ELSE 294 co2trac(itc)%cpl=.FALSE. 295 ! co2trac(itc)%updatefreq = 10800 ! 10800sec = 3H 296 co2trac(itc)%updatefreq = 86400*mth_len ! Once a month 297 END IF 298 CASE("fCO2_land_use") 299 itc = itc + 1 300 co2trac(itc)%name='fCO2_land_use' 301 co2trac(itc)%id=it 302 co2trac(itc)%file='fl_co2_land_use.nc' 303 IF (carbon_cycle_cpl .AND. ok_veget) THEN 304 co2trac(it)%cpl=.TRUE. 305 co2trac(itc)%updatefreq = INT(pdtphys) ! Each timestep as the coupling with ORCHIDEE 306 ELSE 307 co2trac(itc)%cpl=.FALSE. 308 co2trac(itc)%updatefreq = 10800 ! 10800sec = 3H 309 END IF 310 CASE("fCO2_fos_fuel") 311 itc = itc + 1 312 co2trac(itc)%name='fCO2_fos_fuel' 313 co2trac(itc)%id=it 314 co2trac(itc)%file='fossil_fuel.nc' 315 co2trac(itc)%cpl=.FALSE. ! This tracer always read from file 316 ! co2trac(itc)%updatefreq = 86400 ! 86400sec = 24H Cadule case 317 co2trac(itc)%updatefreq = 86400*mth_len ! Once a month 318 CASE("fCO2_bbg") 319 itc = itc + 1 320 co2trac(itc)%name='fCO2_bbg' 321 co2trac(itc)%id=it 322 co2trac(itc)%file='fl_co2_bbg.nc' 323 co2trac(itc)%cpl=.FALSE. ! This tracer always read from file 324 co2trac(itc)%updatefreq = 86400*mth_len ! Once a month 325 CASE("fCO2") 326 ! fCO2 : One tracer transporting the total CO2 flux 327 itc = itc + 1 328 co2trac(itc)%name='fCO2' 329 co2trac(itc)%id=it 330 co2trac(itc)%file='fl_co2.nc' 331 IF (carbon_cycle_cpl) THEN 332 co2trac(itc)%cpl=.TRUE. 333 ELSE 334 co2trac(itc)%cpl=.FALSE. 335 END IF 336 co2trac(itc)%updatefreq = 86400 337 ! DOES THIS WORK ???? Problematic due to implementation of the coupled fluxes... 338 CALL abort_physic('carbon_cycle_init','transport of total CO2 has to be implemented and tested',1) 339 END SELECT 340 END DO 341 342 ! Total number of carbon CO2 tracers 343 ntr_co2 = itc 344 345 ! Definition of control varaiables for the tracers 346 DO it=1,ntr_co2 347 aerosol(co2trac(it)%id) = .FALSE. 348 radio(co2trac(it)%id) = .FALSE. 349 END DO 350 351 ! Vector indicating which timestep to read for each tracer 352 ! Always start read in the beginning of the file 353 co2trac(:)%readstep = 0 354 355 356 ! 3) Allocate variables 357 ! --------------------- 358 ! Allocate vector for storing fluxes to inject 359 ALLOCATE(dtr_add(klon,maxco2trac), stat=ierr) 360 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 11',1) 361 362 ! Allocate variables for cumulating fluxes from ORCHIDEE 363 IF (RCO2_inter) THEN 364 IF (.NOT. carbon_cycle_tr .AND. carbon_cycle_cpl) THEN 365 ALLOCATE(fco2_land_day(klon), stat=ierr) 366 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 2',1) 367 fco2_land_day(1:klon) = 0. 368 369 ALLOCATE(fco2_lu_day(klon), stat=ierr) 370 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 3',1) 371 fco2_lu_day(1:klon) = 0. 372 END IF 373 END IF 374 375 376 ! 4) Test for compatibility 377 ! ------------------------- 378 ! IF (carbon_cycle_cpl .AND. type_ocean/='couple') THEN 379 ! WRITE(lunout,*) 'Coupling with ocean model is needed for carbon_cycle_cpl' 380 ! CALL abort_physic('carbon_cycle_init', 'coupled ocean is needed for carbon_cycle_cpl',1) 381 ! END IF 382 ! 383 ! IF (carbon_cycle_cpl .AND..NOT. ok_veget) THEN 384 ! WRITE(lunout,*) 'Coupling with surface land model ORCHDIEE is needed for carbon_cycle_cpl' 385 ! CALL abort_physic('carbon_cycle_init', 'ok_veget is needed for carbon_cycle_cpl',1) 386 ! END IF 387 388 ! Compiler test : following should never happen 389 teststop=0 390 DO it=1,teststop 391 CALL abort_physic('carbon_cycle_init', 'Entering loop from 1 to 0',1) 392 END DO 393 394 IF (ntr_co2==0) THEN 395 ! No carbon tracers found in tracer.def. It is not possible to do carbon cycle 396 WRITE(lunout,*) 'No carbon tracers found in tracer.def. Not ok with carbon_cycle_tr and/or carbon_cycle_cp' 397 CALL abort_physic('carbon_cycle_init', 'No carbon tracers found in tracer.def',1) 398 END IF 399 400 ! 5) Calculate total area of the earth surface 401 ! -------------------------------------------- 402 CALL reduce_sum(SUM(cell_area),airetot) 403 CALL bcast(airetot) 217 INTEGER :: ierr 218 219 IF (carbon_cycle_cpl) THEN 220 221 ierr=0 222 223 IF (.NOT.ALLOCATED(fco2_land)) ALLOCATE(fco2_land(klon), stat=ierr) 224 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land',1) 225 fco2_land(1:klon) = 0. 226 227 IF (.NOT.ALLOCATED(fco2_land_nbp)) ALLOCATE(fco2_land_nbp(klon), stat=ierr) 228 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_nbp',1) 229 fco2_land_nbp(1:klon) = 0. 230 231 IF (.NOT.ALLOCATED(fco2_land_nep)) ALLOCATE(fco2_land_nep(klon), stat=ierr) 232 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_nep',1) 233 fco2_land_nep(1:klon) = 0. 234 235 IF (.NOT.ALLOCATED(fco2_land_fLuc)) ALLOCATE(fco2_land_fLuc(klon), stat=ierr) 236 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_fLuc',1) 237 fco2_land_fLuc(1:klon) = 0. 238 239 IF (.NOT.ALLOCATED(fco2_land_fwoodharvest)) ALLOCATE(fco2_land_fwoodharvest(klon), stat=ierr) 240 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_fwoodharvest',1) 241 fco2_land_fwoodharvest(1:klon) = 0. 242 243 IF (.NOT.ALLOCATED(fco2_land_fHarvest)) ALLOCATE(fco2_land_fHarvest(klon), stat=ierr) 244 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_fHarvest',1) 245 fco2_land_fHarvest(1:klon) = 0. 246 247 IF (.NOT.ALLOCATED(fco2_ff)) ALLOCATE(fco2_ff(klon), stat=ierr) 248 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_ff',1) 249 fco2_ff(1:klon) = 0. 250 251 IF (.NOT.ALLOCATED(fco2_bb)) ALLOCATE(fco2_bb(klon), stat=ierr) 252 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_bb',1) 253 fco2_bb(1:klon) = 0. 254 255 IF (.NOT.ALLOCATED(fco2_ocean)) ALLOCATE(fco2_ocean(klon), stat=ierr) 256 IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_ocean',1) 257 fco2_bb(1:klon) = 0. 258 ENDIF 404 259 405 260 END SUBROUTINE carbon_cycle_init 406 261 407 SUBROUTINE carbon_cycle(nstep, pdtphys, pctsrf, tr_seri, source)408 ! Subroutine for injection of co2 in the tracers409 !410 ! - Find out if it is time to update411 ! - Get tracer from coupled model or from file412 ! - Calculate new RCO2 value for the radiation scheme413 ! - Calculate CO2 flux to send to ocean and land models (PISCES and ORCHIDEE)414 415 USE infotrac_phy, ONLY: nbtr416 USE dimphy417 USE mod_phys_lmdz_transfert_para418 USE phys_cal_mod, ONLY : mth_cur, mth_len419 USE phys_cal_mod, ONLY : day_cur420 USE indice_sol_mod421 USE print_control_mod, ONLY: lunout422 USE geometry_mod, ONLY : cell_area423 424 IMPLICIT NONE425 426 INCLUDE "clesphys.h"427 INCLUDE "YOMCST.h"428 429 ! In/Output arguments430 INTEGER,INTENT(IN) :: nstep ! time step in physiq431 REAL,INTENT(IN) :: pdtphys ! length of time step in physiq (sec)432 REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Surface fraction433 REAL, DIMENSION(klon,klev,nbtr), INTENT(INOUT) :: tr_seri ! All tracers434 REAL, DIMENSION(klon,nbtr), INTENT(INOUT) :: source ! Source for all tracers435 436 ! Local variables437 INTEGER :: it438 LOGICAL :: newmonth ! indicates if a new month just started439 LOGICAL :: newday ! indicates if a new day just started440 LOGICAL :: endday ! indicated if last time step in a day441 442 REAL, PARAMETER :: fact=1.E-15/2.12 ! transformation factor from gC/m2/day => ppm/m2/day443 REAL, DIMENSION(klon) :: fco2_tmp444 REAL :: sumtmp445 REAL :: delta_co2_ppm446 447 448 ! 1) Calculate logicals indicating if it is a new month, new day or the last time step in a day (end day)449 ! -------------------------------------------------------------------------------------------------------450 451 newday = .FALSE.; endday = .FALSE.; newmonth = .FALSE.452 453 IF (MOD(nstep,INT(86400./pdtphys))==1) newday=.TRUE.454 IF (MOD(nstep,INT(86400./pdtphys))==0) endday=.TRUE.455 IF (newday .AND. day_cur==1) newmonth=.TRUE.456 457 ! 2) For each carbon tracer find out if it is time to inject (update)458 ! --------------------------------------------------------------------459 DO it = 1, ntr_co2460 IF ( MOD(nstep,INT(co2trac(it)%updatefreq/pdtphys)) == 1 ) THEN461 co2trac(it)%updatenow = .TRUE.462 ELSE463 co2trac(it)%updatenow = .FALSE.464 END IF465 END DO466 467 ! 3) Get tracer update468 ! --------------------------------------469 DO it = 1, ntr_co2470 IF ( co2trac(it)%updatenow ) THEN471 IF ( co2trac(it)%cpl ) THEN472 ! Get tracer from coupled model473 SELECT CASE(co2trac(it)%name)474 CASE('fCO2_land') ! from ORCHIDEE475 dtr_add(:,it) = fco2_land_inst(:)*pctsrf(:,is_ter)*fact ! [ppm/m2/day]476 CASE('fCO2_land_use') ! from ORCHIDEE477 dtr_add(:,it) = fco2_lu_inst(:) *pctsrf(:,is_ter)*fact ! [ppm/m2/day]478 CASE('fCO2_ocn') ! from PISCES479 dtr_add(:,it) = fco2_ocn_day(:) *pctsrf(:,is_oce)*fact ! [ppm/m2/day]480 CASE DEFAULT481 WRITE(lunout,*) 'Error with tracer ',co2trac(it)%name482 CALL abort_physic('carbon_cycle', 'No coupling implemented for this tracer',1)483 END SELECT484 ELSE485 ! Read tracer from file486 co2trac(it)%readstep = co2trac(it)%readstep + 1 ! increment time step in file487 ! Patricia CALL read_map2D(co2trac(it)%file,'fco2',co2trac(it)%readstep,.FALSE.,dtr_add(:,it))488 CALL read_map2D(co2trac(it)%file,'fco2',co2trac(it)%readstep,.TRUE.,dtr_add(:,it))489 490 ! Converte from kgC/m2/h to kgC/m2/s491 dtr_add(:,it) = dtr_add(:,it)/3600492 ! Add individual treatment of values read from file493 SELECT CASE(co2trac(it)%name)494 CASE('fCO2_land')495 dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_ter)496 CASE('fCO2_land_use')497 dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_ter)498 CASE('fCO2_ocn')499 dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_oce)500 ! Patricia :501 ! CASE('fCO2_fos_fuel')502 ! dtr_add(:,it) = dtr_add(:,it)/mth_len503 ! co2trac(it)%readstep = 0 ! Always read same value for fossil fuel(Cadule case)504 END SELECT505 END IF506 END IF507 END DO508 509 ! 4) Update co2 tracers :510 ! Loop over all carbon tracers and add source511 ! ------------------------------------------------------------------512 IF (carbon_cycle_tr) THEN513 DO it = 1, ntr_co2514 IF (.FALSE.) THEN515 tr_seri(1:klon,1,co2trac(it)%id) = tr_seri(1:klon,1,co2trac(it)%id) + dtr_add(1:klon,it)516 source(1:klon,co2trac(it)%id) = 0.517 ELSE518 source(1:klon,co2trac(it)%id) = dtr_add(1:klon,it)519 END IF520 END DO521 END IF522 523 524 ! 5) Calculations for new CO2 value for the radiation scheme(instead of reading value from .def)525 ! ----------------------------------------------------------------------------------------------526 IF (RCO2_inter) THEN527 ! Cumulate fluxes from ORCHIDEE at each timestep528 IF (.NOT. carbon_cycle_tr .AND. carbon_cycle_cpl) THEN529 IF (newday) THEN ! Reset cumulative variables once a day530 fco2_land_day(1:klon) = 0.531 fco2_lu_day(1:klon) = 0.532 END IF533 fco2_land_day(1:klon) = fco2_land_day(1:klon) + fco2_land_inst(1:klon) ![gC/m2/day]534 fco2_lu_day(1:klon) = fco2_lu_day(1:klon) + fco2_lu_inst(1:klon) ![gC/m2/day]535 END IF536 537 ! At the end of a new day, calculate a mean scalare value of CO2538 ! JG : Ici on utilise uniquement le traceur du premier couche du modele. Est-ce que c'est correcte ?539 IF (endday) THEN540 541 IF (carbon_cycle_tr) THEN542 ! Sum all co2 tracers to get the total delta CO2 flux543 fco2_tmp(:) = 0.544 DO it = 1, ntr_co2545 fco2_tmp(1:klon) = fco2_tmp(1:klon) + tr_seri(1:klon,1,co2trac(it)%id)546 END DO547 548 ELSE IF (carbon_cycle_cpl) THEN ! no carbon_cycle_tr549 ! Sum co2 fluxes comming from coupled models and parameter for fossil fuel550 fco2_tmp(1:klon) = fos_fuel_s + ((fco2_lu_day(1:klon) + fco2_land_day(1:klon))*pctsrf(1:klon,is_ter) &551 + fco2_ocn_day(:)*pctsrf(:,is_oce)) * fact552 END IF553 554 ! Calculate a global mean value of delta CO2 flux555 fco2_tmp(1:klon) = fco2_tmp(1:klon) * cell_area(1:klon)556 CALL reduce_sum(SUM(fco2_tmp),sumtmp)557 CALL bcast(sumtmp)558 delta_co2_ppm = sumtmp/airetot559 560 ! Add initial value for co2_ppm and delta value561 co2_ppm = co2_ppm0 + delta_co2_ppm562 563 ! Transformation of atmospheric CO2 concentration for the radiation code564 RCO2 = co2_ppm * 1.0e-06 * RMCO2 / RMD565 566 WRITE(lunout,*) 'RCO2 is now updated! RCO2 = ', RCO2567 END IF ! endday568 569 END IF ! RCO2_inter570 571 572 ! 6) Calculate CO2 flux to send to ocean and land models : PISCES and ORCHIDEE573 ! ----------------------------------------------------------------------------574 IF (carbon_cycle_cpl) THEN575 576 IF (carbon_cycle_tr) THEN577 ! Sum all co2 tracers to get the total delta CO2 flux at first model layer578 fco2_tmp(:) = 0.579 DO it = 1, ntr_co2580 fco2_tmp(1:klon) = fco2_tmp(1:klon) + tr_seri(1:klon,1,co2trac(it)%id)581 END DO582 co2_send(1:klon) = fco2_tmp(1:klon) + co2_ppm0583 ELSE584 ! Send a scalare value in 2D variable to ocean and land model (PISCES and ORCHIDEE)585 co2_send(1:klon) = co2_ppm586 END IF587 588 END IF589 590 END SUBROUTINE carbon_cycle591 592 262 SUBROUTINE infocfields_init 593 263
Note: See TracChangeset
for help on using the changeset viewer.
