| [1227] | 1 | MODULE carbon_cycle_mod |
|---|
| [1454] | 2 | ! Controle module for the carbon CO2 tracers : |
|---|
| 3 | ! - Identification |
|---|
| 4 | ! - Get concentrations comming from coupled model or read from file to tracers |
|---|
| 5 | ! - Calculate new RCO2 for radiation scheme |
|---|
| 6 | ! - Calculate new carbon flux for sending to coupled models (PISCES and ORCHIDEE) |
|---|
| 7 | ! |
|---|
| [1227] | 8 | ! Author : Josefine GHATTAS, Patricia CADULE |
|---|
| 9 | |
|---|
| 10 | IMPLICIT NONE |
|---|
| 11 | SAVE |
|---|
| 12 | PRIVATE |
|---|
| 13 | PUBLIC :: carbon_cycle_init, carbon_cycle |
|---|
| 14 | |
|---|
| 15 | ! Variables read from parmeter file physiq.def |
|---|
| [1250] | 16 | LOGICAL, PUBLIC :: carbon_cycle_tr ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys |
|---|
| [1249] | 17 | !$OMP THREADPRIVATE(carbon_cycle_tr) |
|---|
| [1250] | 18 | LOGICAL, PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) |
|---|
| [1249] | 19 | !$OMP THREADPRIVATE(carbon_cycle_cpl) |
|---|
| [1454] | 20 | |
|---|
| [1795] | 21 | LOGICAL :: carbon_cycle_emis_comp_omp=.FALSE. |
|---|
| [1227] | 22 | LOGICAL :: carbon_cycle_emis_comp=.FALSE. ! Calculation of emission compatible |
|---|
| [1454] | 23 | !$OMP THREADPRIVATE(carbon_cycle_emis_comp) |
|---|
| [1227] | 24 | |
|---|
| [1795] | 25 | LOGICAL :: RCO2_inter_omp |
|---|
| [1454] | 26 | LOGICAL :: RCO2_inter ! RCO2 interactive : if true calculate new value RCO2 for the radiation scheme |
|---|
| 27 | !$OMP THREADPRIVATE(RCO2_inter) |
|---|
| 28 | |
|---|
| [1227] | 29 | ! Scalare values when no transport, from physiq.def |
|---|
| [1795] | 30 | REAL :: fos_fuel_s_omp |
|---|
| [1227] | 31 | REAL :: fos_fuel_s ! carbon_cycle_fos_fuel dans physiq.def |
|---|
| [1250] | 32 | !$OMP THREADPRIVATE(fos_fuel_s) |
|---|
| [1227] | 33 | REAL :: emis_land_s ! not yet implemented |
|---|
| [1250] | 34 | !$OMP THREADPRIVATE(emis_land_s) |
|---|
| [1227] | 35 | |
|---|
| [1454] | 36 | REAL :: airetot ! Total area of the earth surface |
|---|
| 37 | !$OMP THREADPRIVATE(airetot) |
|---|
| [1250] | 38 | |
|---|
| [1454] | 39 | INTEGER :: ntr_co2 ! Number of tracers concerning the carbon cycle |
|---|
| 40 | !$OMP THREADPRIVATE(ntr_co2) |
|---|
| 41 | |
|---|
| 42 | ! fco2_ocn_day : flux CO2 from ocean for 1 day (cumulated) [gC/m2/d]. Allocation and initalization done in cpl_mod |
|---|
| 43 | REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocn_day |
|---|
| [1250] | 44 | !$OMP THREADPRIVATE(fco2_ocn_day) |
|---|
| [1454] | 45 | |
|---|
| [1227] | 46 | REAL, DIMENSION(:), ALLOCATABLE :: fco2_land_day ! flux CO2 from land for 1 day (cumulated) [gC/m2/d] |
|---|
| [1250] | 47 | !$OMP THREADPRIVATE(fco2_land_day) |
|---|
| [1227] | 48 | REAL, DIMENSION(:), ALLOCATABLE :: fco2_lu_day ! Emission from land use change for 1 day (cumulated) [gC/m2/d] |
|---|
| [1250] | 49 | !$OMP THREADPRIVATE(fco2_lu_day) |
|---|
| [1227] | 50 | |
|---|
| [1454] | 51 | REAL, DIMENSION(:,:), ALLOCATABLE :: dtr_add ! Tracer concentration to be injected |
|---|
| 52 | !$OMP THREADPRIVATE(dtr_add) |
|---|
| 53 | |
|---|
| 54 | ! Following 2 fields will be allocated and initialized in surf_land_orchidee |
|---|
| [1227] | 55 | REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_inst ! flux CO2 from land at one time step |
|---|
| [1250] | 56 | !$OMP THREADPRIVATE(fco2_land_inst) |
|---|
| [1227] | 57 | REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_lu_inst ! Emission from land use change at one time step |
|---|
| [1250] | 58 | !$OMP THREADPRIVATE(fco2_lu_inst) |
|---|
| [1227] | 59 | |
|---|
| 60 | ! Calculated co2 field to be send to the ocean via the coupler and to ORCHIDEE |
|---|
| [1454] | 61 | REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: co2_send ! Field allocated in phyetat0 |
|---|
| [1250] | 62 | !$OMP THREADPRIVATE(co2_send) |
|---|
| [1227] | 63 | |
|---|
| [1454] | 64 | |
|---|
| 65 | TYPE, PUBLIC :: co2_trac_type |
|---|
| 66 | CHARACTER(len = 8) :: name ! Tracer name in tracer.def |
|---|
| 67 | INTEGER :: id ! Index in total tracer list, tr_seri |
|---|
| 68 | CHARACTER(len=30) :: file ! File name |
|---|
| 69 | LOGICAL :: cpl ! True if this tracers is coupled from ORCHIDEE or PISCES. |
|---|
| 70 | ! False if read from file. |
|---|
| 71 | INTEGER :: updatefreq ! Frequence to inject in second |
|---|
| 72 | INTEGER :: readstep ! Actual time step to read in file |
|---|
| 73 | LOGICAL :: updatenow ! True if this tracer should be updated this time step |
|---|
| 74 | END TYPE co2_trac_type |
|---|
| 75 | INTEGER,PARAMETER :: maxco2trac=5 ! Maximum number of different CO2 fluxes |
|---|
| 76 | TYPE(co2_trac_type), DIMENSION(maxco2trac) :: co2trac |
|---|
| 77 | |
|---|
| [1227] | 78 | CONTAINS |
|---|
| 79 | |
|---|
| [1454] | 80 | SUBROUTINE carbon_cycle_init(tr_seri, pdtphys, aerosol, radio) |
|---|
| 81 | ! This subroutine is called from traclmdz_init, only at first timestep. |
|---|
| 82 | ! - Read controle parameters from .def input file |
|---|
| 83 | ! - Search for carbon tracers and set default values |
|---|
| 84 | ! - Allocate variables |
|---|
| 85 | ! - Test for compatibility |
|---|
| 86 | |
|---|
| [1227] | 87 | USE dimphy |
|---|
| [2408] | 88 | USE geometry_mod, ONLY : cell_area |
|---|
| [1454] | 89 | USE mod_phys_lmdz_transfert_para |
|---|
| [2408] | 90 | USE infotrac_phy, ONLY: nbtr, nqo, niadv, tname |
|---|
| [1227] | 91 | USE IOIPSL |
|---|
| 92 | USE surface_data, ONLY : ok_veget, type_ocean |
|---|
| [1454] | 93 | USE phys_cal_mod, ONLY : mth_len |
|---|
| [2408] | 94 | USE print_control_mod, ONLY: lunout |
|---|
| [1227] | 95 | |
|---|
| 96 | IMPLICIT NONE |
|---|
| 97 | INCLUDE "clesphys.h" |
|---|
| 98 | |
|---|
| 99 | ! Input argument |
|---|
| 100 | REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: tr_seri ! Concentration Traceur [U/KgA] |
|---|
| [1454] | 101 | REAL,INTENT(IN) :: pdtphys ! length of time step in physiq (sec) |
|---|
| [1227] | 102 | |
|---|
| 103 | ! InOutput arguments |
|---|
| 104 | LOGICAL,DIMENSION(nbtr), INTENT(INOUT) :: aerosol |
|---|
| 105 | LOGICAL,DIMENSION(nbtr), INTENT(INOUT) :: radio |
|---|
| 106 | |
|---|
| 107 | ! Local variables |
|---|
| [1454] | 108 | INTEGER :: ierr, it, iiq, itc |
|---|
| 109 | INTEGER :: teststop |
|---|
| [1227] | 110 | |
|---|
| 111 | |
|---|
| 112 | |
|---|
| [1454] | 113 | ! 1) Read controle parameters from .def input file |
|---|
| 114 | ! ------------------------------------------------ |
|---|
| 115 | ! Read fosil fuel value if no transport |
|---|
| 116 | IF (.NOT. carbon_cycle_tr) THEN |
|---|
| [1795] | 117 | !$OMP MASTER |
|---|
| 118 | fos_fuel_s_omp = 0. |
|---|
| 119 | CALL getin ('carbon_cycle_fos_fuel',fos_fuel_s_omp) |
|---|
| 120 | !$OMP END MASTER |
|---|
| 121 | !$OMP BARRIER |
|---|
| 122 | fos_fuel_s=fos_fuel_s_omp |
|---|
| [1227] | 123 | WRITE(lunout,*) 'carbon_cycle_fos_fuel = ', fos_fuel_s |
|---|
| 124 | END IF |
|---|
| 125 | |
|---|
| 126 | |
|---|
| [1454] | 127 | ! Read parmeter for calculation compatible emission |
|---|
| 128 | IF (.NOT. carbon_cycle_tr) THEN |
|---|
| [1795] | 129 | !$OMP MASTER |
|---|
| 130 | carbon_cycle_emis_comp_omp=.FALSE. |
|---|
| 131 | CALL getin('carbon_cycle_emis_comp',carbon_cycle_emis_comp_omp) |
|---|
| 132 | !$OMP END MASTER |
|---|
| 133 | !$OMP BARRIER |
|---|
| 134 | carbon_cycle_emis_comp=carbon_cycle_emis_comp_omp |
|---|
| [1454] | 135 | WRITE(lunout,*) 'carbon_cycle_emis_comp = ',carbon_cycle_emis_comp |
|---|
| 136 | IF (carbon_cycle_emis_comp) THEN |
|---|
| [2408] | 137 | CALL abort_physic('carbon_cycle_init', 'carbon_cycle_emis_comp option not yet implemented!!',1) |
|---|
| [1454] | 138 | END IF |
|---|
| 139 | END IF |
|---|
| [1227] | 140 | |
|---|
| [1454] | 141 | ! Read parameter for interactive calculation of the CO2 value for the radiation scheme |
|---|
| [1795] | 142 | !$OMP MASTER |
|---|
| 143 | RCO2_inter_omp=.FALSE. |
|---|
| 144 | CALL getin('RCO2_inter',RCO2_inter_omp) |
|---|
| 145 | !$OMP END MASTER |
|---|
| 146 | !$OMP BARRIER |
|---|
| 147 | RCO2_inter=RCO2_inter_omp |
|---|
| [1454] | 148 | WRITE(lunout,*) 'RCO2_inter = ', RCO2_inter |
|---|
| 149 | IF (RCO2_inter) THEN |
|---|
| 150 | WRITE(lunout,*) 'RCO2 will be recalculated once a day' |
|---|
| 151 | WRITE(lunout,*) 'RCO2 initial = ', RCO2 |
|---|
| 152 | END IF |
|---|
| [1227] | 153 | |
|---|
| [1454] | 154 | |
|---|
| 155 | ! 2) Search for carbon tracers and set default values |
|---|
| 156 | ! --------------------------------------------------- |
|---|
| 157 | itc=0 |
|---|
| 158 | DO it=1,nbtr |
|---|
| [2298] | 159 | !! iiq=niadv(it+2) ! jyg |
|---|
| 160 | iiq=niadv(it+nqo) ! jyg |
|---|
| [1227] | 161 | |
|---|
| [1454] | 162 | SELECT CASE(tname(iiq)) |
|---|
| 163 | CASE("fCO2_ocn") |
|---|
| 164 | itc = itc + 1 |
|---|
| 165 | co2trac(itc)%name='fCO2_ocn' |
|---|
| 166 | co2trac(itc)%id=it |
|---|
| 167 | co2trac(itc)%file='fl_co2_ocean.nc' |
|---|
| 168 | IF (carbon_cycle_cpl .AND. type_ocean=='couple') THEN |
|---|
| 169 | co2trac(itc)%cpl=.TRUE. |
|---|
| 170 | co2trac(itc)%updatefreq = 86400 ! Once a day as the coupling with OASIS/PISCES |
|---|
| 171 | ELSE |
|---|
| 172 | co2trac(itc)%cpl=.FALSE. |
|---|
| 173 | co2trac(itc)%updatefreq = 86400*mth_len ! Once a month |
|---|
| [1227] | 174 | END IF |
|---|
| [1454] | 175 | CASE("fCO2_land") |
|---|
| 176 | itc = itc + 1 |
|---|
| 177 | co2trac(itc)%name='fCO2_land' |
|---|
| 178 | co2trac(itc)%id=it |
|---|
| 179 | co2trac(itc)%file='fl_co2_land.nc' |
|---|
| 180 | IF (carbon_cycle_cpl .AND. ok_veget) THEN |
|---|
| 181 | co2trac(itc)%cpl=.TRUE. |
|---|
| 182 | co2trac(itc)%updatefreq = INT(pdtphys) ! Each timestep as the coupling with ORCHIDEE |
|---|
| 183 | ELSE |
|---|
| 184 | co2trac(itc)%cpl=.FALSE. |
|---|
| 185 | ! co2trac(itc)%updatefreq = 10800 ! 10800sec = 3H |
|---|
| 186 | co2trac(itc)%updatefreq = 86400*mth_len ! Once a month |
|---|
| 187 | END IF |
|---|
| 188 | CASE("fCO2_land_use") |
|---|
| 189 | itc = itc + 1 |
|---|
| 190 | co2trac(itc)%name='fCO2_land_use' |
|---|
| 191 | co2trac(itc)%id=it |
|---|
| 192 | co2trac(itc)%file='fl_co2_land_use.nc' |
|---|
| 193 | IF (carbon_cycle_cpl .AND. ok_veget) THEN |
|---|
| 194 | co2trac(it)%cpl=.TRUE. |
|---|
| 195 | co2trac(itc)%updatefreq = INT(pdtphys) ! Each timestep as the coupling with ORCHIDEE |
|---|
| 196 | ELSE |
|---|
| 197 | co2trac(itc)%cpl=.FALSE. |
|---|
| 198 | co2trac(itc)%updatefreq = 10800 ! 10800sec = 3H |
|---|
| 199 | END IF |
|---|
| 200 | CASE("fCO2_fos_fuel") |
|---|
| 201 | itc = itc + 1 |
|---|
| 202 | co2trac(itc)%name='fCO2_fos_fuel' |
|---|
| 203 | co2trac(itc)%id=it |
|---|
| 204 | co2trac(itc)%file='fossil_fuel.nc' |
|---|
| 205 | co2trac(itc)%cpl=.FALSE. ! This tracer always read from file |
|---|
| 206 | ! co2trac(itc)%updatefreq = 86400 ! 86400sec = 24H Cadule case |
|---|
| 207 | co2trac(itc)%updatefreq = 86400*mth_len ! Once a month |
|---|
| 208 | CASE("fCO2_bbg") |
|---|
| 209 | itc = itc + 1 |
|---|
| 210 | co2trac(itc)%name='fCO2_bbg' |
|---|
| 211 | co2trac(itc)%id=it |
|---|
| 212 | co2trac(itc)%file='fl_co2_bbg.nc' |
|---|
| 213 | co2trac(itc)%cpl=.FALSE. ! This tracer always read from file |
|---|
| 214 | co2trac(itc)%updatefreq = 86400*mth_len ! Once a month |
|---|
| 215 | CASE("fCO2") |
|---|
| 216 | ! fCO2 : One tracer transporting the total CO2 flux |
|---|
| 217 | itc = itc + 1 |
|---|
| 218 | co2trac(itc)%name='fCO2' |
|---|
| 219 | co2trac(itc)%id=it |
|---|
| 220 | co2trac(itc)%file='fl_co2.nc' |
|---|
| 221 | IF (carbon_cycle_cpl) THEN |
|---|
| 222 | co2trac(itc)%cpl=.TRUE. |
|---|
| 223 | ELSE |
|---|
| 224 | co2trac(itc)%cpl=.FALSE. |
|---|
| 225 | END IF |
|---|
| 226 | co2trac(itc)%updatefreq = 86400 |
|---|
| 227 | ! DOES THIS WORK ???? Problematic due to implementation of the coupled fluxes... |
|---|
| [2408] | 228 | CALL abort_physic('carbon_cycle_init','transport of total CO2 has to be implemented and tested',1) |
|---|
| [1454] | 229 | END SELECT |
|---|
| 230 | END DO |
|---|
| [1227] | 231 | |
|---|
| [1454] | 232 | ! Total number of carbon CO2 tracers |
|---|
| 233 | ntr_co2 = itc |
|---|
| 234 | |
|---|
| 235 | ! Definition of control varaiables for the tracers |
|---|
| 236 | DO it=1,ntr_co2 |
|---|
| 237 | aerosol(co2trac(it)%id) = .FALSE. |
|---|
| 238 | radio(co2trac(it)%id) = .FALSE. |
|---|
| 239 | END DO |
|---|
| 240 | |
|---|
| 241 | ! Vector indicating which timestep to read for each tracer |
|---|
| 242 | ! Always start read in the beginning of the file |
|---|
| 243 | co2trac(:)%readstep = 0 |
|---|
| 244 | |
|---|
| [1227] | 245 | |
|---|
| [1454] | 246 | ! 3) Allocate variables |
|---|
| 247 | ! --------------------- |
|---|
| 248 | ! Allocate vector for storing fluxes to inject |
|---|
| 249 | ALLOCATE(dtr_add(klon,maxco2trac), stat=ierr) |
|---|
| [2408] | 250 | IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 11',1) |
|---|
| [1454] | 251 | |
|---|
| 252 | ! Allocate variables for cumulating fluxes from ORCHIDEE |
|---|
| 253 | IF (RCO2_inter) THEN |
|---|
| 254 | IF (.NOT. carbon_cycle_tr .AND. carbon_cycle_cpl) THEN |
|---|
| [1227] | 255 | ALLOCATE(fco2_land_day(klon), stat=ierr) |
|---|
| [2408] | 256 | IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 2',1) |
|---|
| [1454] | 257 | fco2_land_day(1:klon) = 0. |
|---|
| 258 | |
|---|
| [1227] | 259 | ALLOCATE(fco2_lu_day(klon), stat=ierr) |
|---|
| [2408] | 260 | IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 3',1) |
|---|
| [1454] | 261 | fco2_lu_day(1:klon) = 0. |
|---|
| [1227] | 262 | END IF |
|---|
| 263 | END IF |
|---|
| 264 | |
|---|
| [1454] | 265 | |
|---|
| 266 | ! 4) Test for compatibility |
|---|
| 267 | ! ------------------------- |
|---|
| 268 | ! IF (carbon_cycle_cpl .AND. type_ocean/='couple') THEN |
|---|
| 269 | ! WRITE(lunout,*) 'Coupling with ocean model is needed for carbon_cycle_cpl' |
|---|
| [2408] | 270 | ! CALL abort_physic('carbon_cycle_init', 'coupled ocean is needed for carbon_cycle_cpl',1) |
|---|
| [1454] | 271 | ! END IF |
|---|
| 272 | ! |
|---|
| 273 | ! IF (carbon_cycle_cpl .AND..NOT. ok_veget) THEN |
|---|
| 274 | ! WRITE(lunout,*) 'Coupling with surface land model ORCHDIEE is needed for carbon_cycle_cpl' |
|---|
| [2408] | 275 | ! CALL abort_physic('carbon_cycle_init', 'ok_veget is needed for carbon_cycle_cpl',1) |
|---|
| [1454] | 276 | ! END IF |
|---|
| 277 | |
|---|
| 278 | ! Compiler test : following should never happen |
|---|
| 279 | teststop=0 |
|---|
| 280 | DO it=1,teststop |
|---|
| [2408] | 281 | CALL abort_physic('carbon_cycle_init', 'Entering loop from 1 to 0',1) |
|---|
| [1454] | 282 | END DO |
|---|
| 283 | |
|---|
| 284 | IF (ntr_co2==0) THEN |
|---|
| 285 | ! No carbon tracers found in tracer.def. It is not possible to do carbon cycle |
|---|
| 286 | WRITE(lunout,*) 'No carbon tracers found in tracer.def. Not ok with carbon_cycle_tr and/or carbon_cycle_cp' |
|---|
| [2408] | 287 | CALL abort_physic('carbon_cycle_init', 'No carbon tracers found in tracer.def',1) |
|---|
| [1227] | 288 | END IF |
|---|
| [1454] | 289 | |
|---|
| 290 | ! 5) Calculate total area of the earth surface |
|---|
| 291 | ! -------------------------------------------- |
|---|
| [2408] | 292 | CALL reduce_sum(SUM(cell_area),airetot) |
|---|
| [1454] | 293 | CALL bcast(airetot) |
|---|
| [1227] | 294 | |
|---|
| 295 | END SUBROUTINE carbon_cycle_init |
|---|
| 296 | |
|---|
| [1454] | 297 | |
|---|
| 298 | SUBROUTINE carbon_cycle(nstep, pdtphys, pctsrf, tr_seri, source) |
|---|
| 299 | ! Subroutine for injection of co2 in the tracers |
|---|
| [1227] | 300 | ! |
|---|
| [1454] | 301 | ! - Find out if it is time to update |
|---|
| 302 | ! - Get tracer from coupled model or from file |
|---|
| 303 | ! - Calculate new RCO2 value for the radiation scheme |
|---|
| 304 | ! - Calculate CO2 flux to send to ocean and land models (PISCES and ORCHIDEE) |
|---|
| [1227] | 305 | |
|---|
| [2408] | 306 | USE infotrac_phy, ONLY: nbtr |
|---|
| [1227] | 307 | USE dimphy |
|---|
| [1454] | 308 | USE mod_phys_lmdz_transfert_para |
|---|
| [1227] | 309 | USE phys_cal_mod, ONLY : mth_cur, mth_len |
|---|
| 310 | USE phys_cal_mod, ONLY : day_cur |
|---|
| [1795] | 311 | USE indice_sol_mod |
|---|
| [2408] | 312 | USE print_control_mod, ONLY: lunout |
|---|
| 313 | USE geometry_mod, ONLY : cell_area |
|---|
| [1227] | 314 | |
|---|
| 315 | IMPLICIT NONE |
|---|
| 316 | |
|---|
| 317 | INCLUDE "clesphys.h" |
|---|
| [1454] | 318 | INCLUDE "YOMCST.h" |
|---|
| [1227] | 319 | |
|---|
| 320 | ! In/Output arguments |
|---|
| 321 | INTEGER,INTENT(IN) :: nstep ! time step in physiq |
|---|
| 322 | REAL,INTENT(IN) :: pdtphys ! length of time step in physiq (sec) |
|---|
| [1454] | 323 | REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Surface fraction |
|---|
| 324 | REAL, DIMENSION(klon,klev,nbtr), INTENT(INOUT) :: tr_seri ! All tracers |
|---|
| 325 | REAL, DIMENSION(klon,nbtr), INTENT(INOUT) :: source ! Source for all tracers |
|---|
| [1227] | 326 | |
|---|
| 327 | ! Local variables |
|---|
| [1454] | 328 | INTEGER :: it |
|---|
| [1227] | 329 | LOGICAL :: newmonth ! indicates if a new month just started |
|---|
| 330 | LOGICAL :: newday ! indicates if a new day just started |
|---|
| 331 | LOGICAL :: endday ! indicated if last time step in a day |
|---|
| 332 | |
|---|
| 333 | REAL, PARAMETER :: fact=1.E-15/2.12 ! transformation factor from gC/m2/day => ppm/m2/day |
|---|
| [1454] | 334 | REAL, DIMENSION(klon) :: fco2_tmp |
|---|
| [1227] | 335 | REAL :: sumtmp |
|---|
| 336 | REAL :: delta_co2_ppm |
|---|
| 337 | |
|---|
| 338 | |
|---|
| [1454] | 339 | ! 1) Calculate logicals indicating if it is a new month, new day or the last time step in a day (end day) |
|---|
| 340 | ! ------------------------------------------------------------------------------------------------------- |
|---|
| 341 | |
|---|
| [1227] | 342 | newday = .FALSE.; endday = .FALSE.; newmonth = .FALSE. |
|---|
| 343 | |
|---|
| 344 | IF (MOD(nstep,INT(86400./pdtphys))==1) newday=.TRUE. |
|---|
| 345 | IF (MOD(nstep,INT(86400./pdtphys))==0) endday=.TRUE. |
|---|
| 346 | IF (newday .AND. day_cur==1) newmonth=.TRUE. |
|---|
| [1454] | 347 | |
|---|
| 348 | ! 2) For each carbon tracer find out if it is time to inject (update) |
|---|
| 349 | ! -------------------------------------------------------------------- |
|---|
| 350 | DO it = 1, ntr_co2 |
|---|
| 351 | IF ( MOD(nstep,INT(co2trac(it)%updatefreq/pdtphys)) == 1 ) THEN |
|---|
| 352 | co2trac(it)%updatenow = .TRUE. |
|---|
| 353 | ELSE |
|---|
| 354 | co2trac(it)%updatenow = .FALSE. |
|---|
| [1227] | 355 | END IF |
|---|
| [1454] | 356 | END DO |
|---|
| [1227] | 357 | |
|---|
| [1454] | 358 | ! 3) Get tracer update |
|---|
| 359 | ! -------------------------------------- |
|---|
| 360 | DO it = 1, ntr_co2 |
|---|
| 361 | IF ( co2trac(it)%updatenow ) THEN |
|---|
| 362 | IF ( co2trac(it)%cpl ) THEN |
|---|
| 363 | ! Get tracer from coupled model |
|---|
| 364 | SELECT CASE(co2trac(it)%name) |
|---|
| 365 | CASE('fCO2_land') ! from ORCHIDEE |
|---|
| 366 | dtr_add(:,it) = fco2_land_inst(:)*pctsrf(:,is_ter)*fact ! [ppm/m2/day] |
|---|
| 367 | CASE('fCO2_land_use') ! from ORCHIDEE |
|---|
| 368 | dtr_add(:,it) = fco2_lu_inst(:) *pctsrf(:,is_ter)*fact ! [ppm/m2/day] |
|---|
| 369 | CASE('fCO2_ocn') ! from PISCES |
|---|
| 370 | dtr_add(:,it) = fco2_ocn_day(:) *pctsrf(:,is_oce)*fact ! [ppm/m2/day] |
|---|
| 371 | CASE DEFAULT |
|---|
| 372 | WRITE(lunout,*) 'Error with tracer ',co2trac(it)%name |
|---|
| [2408] | 373 | CALL abort_physic('carbon_cycle', 'No coupling implemented for this tracer',1) |
|---|
| [1454] | 374 | END SELECT |
|---|
| 375 | ELSE |
|---|
| 376 | ! Read tracer from file |
|---|
| 377 | co2trac(it)%readstep = co2trac(it)%readstep + 1 ! increment time step in file |
|---|
| 378 | ! Patricia CALL read_map2D(co2trac(it)%file,'fco2',co2trac(it)%readstep,.FALSE.,dtr_add(:,it)) |
|---|
| 379 | CALL read_map2D(co2trac(it)%file,'fco2',co2trac(it)%readstep,.TRUE.,dtr_add(:,it)) |
|---|
| [1227] | 380 | |
|---|
| [1454] | 381 | ! Converte from kgC/m2/h to kgC/m2/s |
|---|
| 382 | dtr_add(:,it) = dtr_add(:,it)/3600 |
|---|
| 383 | ! Add individual treatment of values read from file |
|---|
| 384 | SELECT CASE(co2trac(it)%name) |
|---|
| 385 | CASE('fCO2_land') |
|---|
| 386 | dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_ter) |
|---|
| 387 | CASE('fCO2_land_use') |
|---|
| 388 | dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_ter) |
|---|
| 389 | CASE('fCO2_ocn') |
|---|
| 390 | dtr_add(:,it) = dtr_add(:,it) *pctsrf(:,is_oce) |
|---|
| 391 | ! Patricia : |
|---|
| 392 | ! CASE('fCO2_fos_fuel') |
|---|
| 393 | ! dtr_add(:,it) = dtr_add(:,it)/mth_len |
|---|
| 394 | ! co2trac(it)%readstep = 0 ! Always read same value for fossil fuel(Cadule case) |
|---|
| 395 | END SELECT |
|---|
| [1227] | 396 | END IF |
|---|
| [1454] | 397 | END IF |
|---|
| 398 | END DO |
|---|
| [1227] | 399 | |
|---|
| [1454] | 400 | ! 4) Update co2 tracers : |
|---|
| 401 | ! Loop over all carbon tracers and add source |
|---|
| 402 | ! ------------------------------------------------------------------ |
|---|
| 403 | IF (carbon_cycle_tr) THEN |
|---|
| 404 | DO it = 1, ntr_co2 |
|---|
| 405 | IF (.FALSE.) THEN |
|---|
| 406 | tr_seri(1:klon,1,co2trac(it)%id) = tr_seri(1:klon,1,co2trac(it)%id) + dtr_add(1:klon,it) |
|---|
| 407 | source(1:klon,co2trac(it)%id) = 0. |
|---|
| 408 | ELSE |
|---|
| 409 | source(1:klon,co2trac(it)%id) = dtr_add(1:klon,it) |
|---|
| [1227] | 410 | END IF |
|---|
| [1454] | 411 | END DO |
|---|
| [1227] | 412 | END IF |
|---|
| 413 | |
|---|
| 414 | |
|---|
| [1454] | 415 | ! 5) Calculations for new CO2 value for the radiation scheme(instead of reading value from .def) |
|---|
| 416 | ! ---------------------------------------------------------------------------------------------- |
|---|
| 417 | IF (RCO2_inter) THEN |
|---|
| 418 | ! Cumulate fluxes from ORCHIDEE at each timestep |
|---|
| 419 | IF (.NOT. carbon_cycle_tr .AND. carbon_cycle_cpl) THEN |
|---|
| 420 | IF (newday) THEN ! Reset cumulative variables once a day |
|---|
| 421 | fco2_land_day(1:klon) = 0. |
|---|
| 422 | fco2_lu_day(1:klon) = 0. |
|---|
| 423 | END IF |
|---|
| 424 | fco2_land_day(1:klon) = fco2_land_day(1:klon) + fco2_land_inst(1:klon) ![gC/m2/day] |
|---|
| 425 | fco2_lu_day(1:klon) = fco2_lu_day(1:klon) + fco2_lu_inst(1:klon) ![gC/m2/day] |
|---|
| 426 | END IF |
|---|
| [1227] | 427 | |
|---|
| [1454] | 428 | ! At the end of a new day, calculate a mean scalare value of CO2 |
|---|
| 429 | ! JG : Ici on utilise uniquement le traceur du premier couche du modele. Est-ce que c'est correcte ? |
|---|
| 430 | IF (endday) THEN |
|---|
| [1227] | 431 | |
|---|
| [1454] | 432 | IF (carbon_cycle_tr) THEN |
|---|
| 433 | ! Sum all co2 tracers to get the total delta CO2 flux |
|---|
| 434 | fco2_tmp(:) = 0. |
|---|
| 435 | DO it = 1, ntr_co2 |
|---|
| 436 | fco2_tmp(1:klon) = fco2_tmp(1:klon) + tr_seri(1:klon,1,co2trac(it)%id) |
|---|
| 437 | END DO |
|---|
| [1227] | 438 | |
|---|
| [1454] | 439 | ELSE IF (carbon_cycle_cpl) THEN ! no carbon_cycle_tr |
|---|
| 440 | ! Sum co2 fluxes comming from coupled models and parameter for fossil fuel |
|---|
| 441 | fco2_tmp(1:klon) = fos_fuel_s + ((fco2_lu_day(1:klon) + fco2_land_day(1:klon))*pctsrf(1:klon,is_ter) & |
|---|
| 442 | + fco2_ocn_day(:)*pctsrf(:,is_oce)) * fact |
|---|
| [1227] | 443 | END IF |
|---|
| 444 | |
|---|
| [1454] | 445 | ! Calculate a global mean value of delta CO2 flux |
|---|
| [2408] | 446 | fco2_tmp(1:klon) = fco2_tmp(1:klon) * cell_area(1:klon) |
|---|
| [1227] | 447 | CALL reduce_sum(SUM(fco2_tmp),sumtmp) |
|---|
| [1454] | 448 | CALL bcast(sumtmp) |
|---|
| [1227] | 449 | delta_co2_ppm = sumtmp/airetot |
|---|
| 450 | |
|---|
| [1454] | 451 | ! Add initial value for co2_ppm and delta value |
|---|
| 452 | co2_ppm = co2_ppm0 + delta_co2_ppm |
|---|
| 453 | |
|---|
| 454 | ! Transformation of atmospheric CO2 concentration for the radiation code |
|---|
| 455 | RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97 |
|---|
| 456 | |
|---|
| 457 | WRITE(lunout,*) 'RCO2 is now updated! RCO2 = ', RCO2 |
|---|
| 458 | END IF ! endday |
|---|
| [1227] | 459 | |
|---|
| [1454] | 460 | END IF ! RCO2_inter |
|---|
| [1227] | 461 | |
|---|
| 462 | |
|---|
| [1454] | 463 | ! 6) Calculate CO2 flux to send to ocean and land models : PISCES and ORCHIDEE |
|---|
| 464 | ! ---------------------------------------------------------------------------- |
|---|
| 465 | IF (carbon_cycle_cpl) THEN |
|---|
| 466 | |
|---|
| [1227] | 467 | IF (carbon_cycle_tr) THEN |
|---|
| [1454] | 468 | ! Sum all co2 tracers to get the total delta CO2 flux at first model layer |
|---|
| 469 | fco2_tmp(:) = 0. |
|---|
| 470 | DO it = 1, ntr_co2 |
|---|
| 471 | fco2_tmp(1:klon) = fco2_tmp(1:klon) + tr_seri(1:klon,1,co2trac(it)%id) |
|---|
| 472 | END DO |
|---|
| 473 | co2_send(1:klon) = fco2_tmp(1:klon) + co2_ppm0 |
|---|
| [1227] | 474 | ELSE |
|---|
| 475 | ! Send a scalare value in 2D variable to ocean and land model (PISCES and ORCHIDEE) |
|---|
| [1454] | 476 | co2_send(1:klon) = co2_ppm |
|---|
| [1227] | 477 | END IF |
|---|
| 478 | |
|---|
| 479 | END IF |
|---|
| 480 | |
|---|
| 481 | END SUBROUTINE carbon_cycle |
|---|
| 482 | |
|---|
| 483 | END MODULE carbon_cycle_mod |
|---|
