Changeset 3453 for LMDZ6/trunk/libf
- Timestamp:
- Jan 28, 2019, 7:31:11 PM (6 years ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 4 edited
-
carbon_cycle_mod.F90 (modified) (1 diff)
-
phys_output_ctrlout_mod.F90 (modified) (3 diffs)
-
phys_output_write_mod.F90 (modified) (10 diffs)
-
tracco2i_mod.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90
r3447 r3453 82 82 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_bb ! Emission from biomass burning [kgCO2/m2/s] 83 83 !$OMP THREADPRIVATE(fco2_bb) 84 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land ! Net flux from terrestrial ecocsystems [kgCO2/m2/s] 85 !$OMP THREADPRIVATE(fco2_land) 86 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocean ! Net flux from ocean [kgCO2/m2/s] 87 !$OMP THREADPRIVATE(fco2_ocean) 84 88 85 89 REAL, DIMENSION(:,:), ALLOCATABLE :: dtr_add ! Tracer concentration to be injected -
LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
r3332 r3453 7 7 USE indice_sol_mod 8 8 USE aero_mod 9 10 11 9 12 10 IMPLICIT NONE … … 1014 1012 !FC 1015 1013 1016 1017 1014 TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_l_mixmin = (/ & 1018 1015 ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11, 11/),'l_mixmin_ter', & … … 1280 1277 !--end add ThL 1281 1278 1279 !---CO2 fluxes for interactive CO2 configuration 1280 TYPE(ctrl_out), SAVE :: o_flx_co2_ff = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), & 1281 'flx_co2_ff', 'CO2 flux from fossil fuel and cement', '1', (/ ('', i=1, 10) /)) 1282 TYPE(ctrl_out), SAVE :: o_flx_co2_bb = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), & 1283 'flx_co2_bb', 'CO2 flux from biomass burning', '1', (/ ('', i=1, 10) /)) 1284 TYPE(ctrl_out), SAVE :: o_flx_co2_ocean = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), & 1285 'flx_co2_ocean', 'CO2 flux from the ocean', '1', (/ ('', i=1, 10) /)) 1286 TYPE(ctrl_out), SAVE :: o_flx_co2_land = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), & 1287 'flx_co2_land', 'CO2 flux from the land', '1', (/ ('', i=1, 10) /)) 1288 1282 1289 #ifdef CPP_StratAer 1283 1290 !--extinction coefficient -
LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
r3435 r3453 196 196 ! Tropopause 197 197 o_p_tropopause, o_z_tropopause, o_t_tropopause, & 198 o_col_O3_strato, o_col_O3_tropo ! Added ThL 198 o_col_O3_strato, o_col_O3_tropo, & 199 !--interactive CO2 200 o_flx_co2_ocean, o_flx_co2_land, o_flx_co2_ff, o_flx_co2_bb 199 201 200 202 … … 327 329 USE CHEM_REP, ONLY : nas, nbnas, tnamenas, ttextnas 328 330 #endif 331 332 USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean 329 333 330 334 USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, & … … 1253 1257 IF (slab_gm) THEN 1254 1258 CALL histwrite_phy(o_slab_gm, dt_gm(:,1:nslay)) 1255 END 1259 ENDIF 1256 1260 IF (slab_hdiff) THEN 1257 1261 IF (nslay.EQ.1) THEN … … 1663 1667 ENDIF 1664 1668 CALL histwrite_phy(o_dtcon, zx_tmp_fi3d) 1665 if(iflag_thermals.eq.0)then1669 IF (iflag_thermals.EQ.0) THEN 1666 1670 IF (vars_defined) THEN 1667 1671 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + & … … 1669 1673 ENDIF 1670 1674 CALL histwrite_phy(o_tntc, zx_tmp_fi3d) 1671 else if(iflag_thermals.ge.1.and.iflag_wake.EQ.1)then1675 ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN 1672 1676 IF (vars_defined) THEN 1673 1677 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + & … … 1676 1680 ENDIF 1677 1681 CALL histwrite_phy(o_tntc, zx_tmp_fi3d) 1678 endif1682 ENDIF 1679 1683 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys 1680 1684 CALL histwrite_phy(o_ducon, zx_tmp_fi3d) … … 1727 1731 CALL histwrite_phy(o_plulst, plul_st) 1728 1732 IF (vars_defined) THEN 1729 doi=1,klon1733 DO i=1,klon 1730 1734 zx_tmp_fi2d(1:klon)=lmax_th(:) 1731 enddo1735 ENDDO 1732 1736 ENDIF 1733 1737 CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d) … … 2192 2196 IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN 2193 2197 !--3D fields 2194 !CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))2195 !CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))2196 !CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))2197 !CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))2198 !CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))2199 !CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))2200 !CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))2201 !CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))2202 !CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))2203 !CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))2204 !CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))2205 !CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))2206 !CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))2207 !CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))2208 2209 !CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))2198 CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo)) 2199 CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo)) 2200 CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo)) 2201 CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo)) 2202 CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo)) 2203 CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo)) 2204 CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo)) 2205 CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo)) 2206 CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo)) 2207 CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo)) 2208 CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo)) 2209 CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo)) 2210 CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo)) 2211 CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo)) 2212 !--2D fields 2213 CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo)) 2210 2214 zx_tmp_fi2d=0. 2211 2215 IF (vars_defined) THEN … … 2217 2221 CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d) 2218 2222 #endif 2219 ENDIF 2220 ENDDO 2221 ENDIF 2223 ENDIF !--type_trac 2224 ! 2225 IF (type_trac == 'co2i') THEN 2226 !--3D fields 2227 CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo)) 2228 CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo)) 2229 CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo)) 2230 CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo)) 2231 !--2D fields 2232 !--CO2 burden 2233 zx_tmp_fi2d=0. 2234 IF (vars_defined) THEN 2235 DO k=1,klev 2236 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo) 2237 ENDDO 2238 ENDIF 2239 CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d) 2240 !--CO2 net fluxes 2241 CALL histwrite_phy(o_flx_co2_land, fco2_land) 2242 CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean) 2243 CALL histwrite_phy(o_flx_co2_ff, fco2_ff) 2244 CALL histwrite_phy(o_flx_co2_bb, fco2_bb) 2245 ENDIF !--type_trac co2i 2246 2247 ENDDO !--iq 2248 ENDIF !--nqtot > nqo+1 2222 2249 2223 2250 IF (type_trac == 'repr') THEN … … 2230 2257 2231 2258 ENDIF !(iflag_phytrac==1) 2232 2233 2259 2234 2260 IF (.NOT.vars_defined) THEN -
LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90
r3450 r3453 12 12 USE infotrac_phy 13 13 USE geometry_mod, ONLY: cell_area 14 USE carbon_cycle_mod, ONLY: id_CO2, nbcf_in, fields_in, cfname_in, fco2_ocn_day, fco2_ff, fco2_bb 14 USE carbon_cycle_mod, ONLY: id_CO2, nbcf_in, fields_in, cfname_in 15 USE carbon_cycle_mod, ONLY: fco2_ocn_day, fco2_ff, fco2_bb, fco2_land, fco2_ocean 15 16 USE carbon_cycle_mod, ONLY: carbon_cycle_tr, carbon_cycle_rad, RCO2_glo, RCO2_tot 16 17 USE mod_grid_phy_lmdz … … 50 51 INTEGER :: it, k, i, nb 51 52 REAL, DIMENSION(klon,klev) :: m_air ! mass of air in every grid box [kg] 52 REAL, DIMENSION(klon) :: co2land ! surface land CO2 emissions [kg CO2/m2/s]53 REAL, DIMENSION(klon) :: co2ocean ! surface ocean CO2 emissions [kg CO2/m2/s]54 53 REAL, DIMENSION(klon_glo,klev) :: co2_glo ! variable temporaire sur la grille global 55 54 REAL, DIMENSION(klon_glo,klev) :: m_air_glo ! variable temporaire sur la grille global … … 85 84 86 85 !--retrieving land and ocean CO2 flux 87 co2land(:)=0.088 co2ocean(:)=0.086 fco2_land(:)=0.0 87 fco2_ocean(:)=0.0 89 88 DO nb=1, nbcf_in 90 89 !--fCO2_nep comes in unit of kg C m-2 s-1 91 90 !--converting to kg CO2 m-2 s-1 92 IF (cfname_in(nb) == "fCO2_nep" ) co2land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter)91 IF (cfname_in(nb) == "fCO2_nep" ) fco2_land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter) 93 92 !--fCO2_fgco2 comes in unit of mol C02 m-2 s-1 94 93 !--converting to kg CO2 m-2 s-1 + change sign 95 IF (cfname_in(nb) == "fCO2_fgco2" ) co2ocean(:)=-1.*fco2_ocn_day(:)*RMCO2/1.e3*(pctsrf(:,is_oce)+pctsrf(:,is_sic))94 IF (cfname_in(nb) == "fCO2_fgco2" ) fco2_ocean(:)=-1.*fco2_ocn_day(:)*RMCO2/1.e3*(pctsrf(:,is_oce)+pctsrf(:,is_sic)) 96 95 ENDDO 97 96 98 !--preparing the net anthropogenic flux at the surface for mixing layer 97 !--preparing the net anthropogenic flux at the surface for mixing layer 99 98 !--unit kg CO2 / m2 / s 100 source(:,id_CO2)=fco2_ff(:)+fco2_bb(:)+ co2land(:)+co2ocean(:)99 source(:,id_CO2)=fco2_ff(:)+fco2_bb(:)+fco2_land(:)+fco2_ocean(:) 101 100 102 101 !--computing global mean CO2 for radiation
Note: See TracChangeset
for help on using the changeset viewer.
