Changeset 4482 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/tracco2i_mod.F90
- Timestamp:
- Mar 29, 2023, 3:14:27 PM (15 months ago)
- Location:
- LMDZ6/branches/LMDZ_ECRad
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ_ECRad
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ_ECRad/libf/phylmd/tracco2i_mod.F90
r4179 r4482 103 103 IF (debutphy) THEN 104 104 105 ! Initiali sation de tr_seri(id_CO2) si pas initialise105 ! Initialization of tr_seri(id_CO2) If it is not initialized 106 106 IF (MAXVAL(tr_seri(:,:,id_CO2)).LT.1.e-15) THEN 107 107 tr_seri(:,:,id_CO2)=co2_ppm0*1.e-6/RMD*RMCO2 !--initialised from co2_ppm0 in rdem … … 226 226 227 227 DO i=1, klon 228 fco2_ocean_cor(i)=(var_fco2_ocean_cor*(RMCO2/RMC)*(pctsrf(i,is_oce)+pctsrf(i,is_sic))/ocean_area_tot/(secinday*year_len))*1.e12 228 fco2_ocean_cor(i)=(var_fco2_ocean_cor*(RMCO2/RMC) & 229 *(pctsrf(i,is_oce)+pctsrf(i,is_sic))/ocean_area_tot & 230 /(secinday*year_len))*1.e12 229 231 ENDDO 230 232 … … 305 307 306 308 !--compute a global mean CO2 value and print its value in ppm 307 IF (is_mpi_root ) THEN309 IF (is_mpi_root .AND. is_omp_root) THEN 308 310 RCO2_tot=SUM(co2_glo*m_air_glo) !--unit kg CO2 309 311 RCO2_glo=RCO2_tot/SUM(m_air_glo) !--unit kg CO2 / kg air 312 ! the following operation is only to maintain precision consistency 313 ! of RCO2_glo which differs whether it is directly computed or read from 314 ! a restart file (after having been computed) 315 RCO2_glo = FLOAT(INT(RCO2_glo * 1e8))/1e8 310 316 PRINT *,'tracco2i: global CO2 in ppm =', RCO2_glo*1.e6*RMD/RMCO2 311 317 PRINT *,'tracco2i: total CO2 in kg =', RCO2_tot … … 314 320 CALL bcast(RCO2_glo) 315 321 day_pre=day_cur 322 316 323 !--if not carbon_cycle_tr, then we reinitialize the CO2 each day to its global mean value 317 324 IF (.NOT.carbon_cycle_tr) THEN
Note: See TracChangeset
for help on using the changeset viewer.