Changeset 3447 for LMDZ6/trunk

Timestamp:

Jan 25, 2019, 5:05:38 PM (6 years ago)

Author:

oboucher

Message:

A first package with molar masses for radiatively active gases, a new switch for interactive radiative CO2
and minor changes to the carbon_cycle module.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• YOMCST.h (modified) (2 diffs)
• carbon_cycle_mod.F90 (modified) (3 diffs)
• conf_phys_m.F90 (modified) (6 diffs)
• suphel.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/YOMCST.h

@@ -20 +20 @@
       REAL RSIGMA
 ! A1.4 Thermodynamic gas phase
-      REAL R,RMD,RMO3,RMCO2,RMC,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
+      REAL RMO3,RMCO2,RMC,RMCH4,RMN2O,RMCFC11,RMCFC12
+      REAL R,RMD,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
       REAL RKAPPA,RETV
 ! A1.5,6 Thermodynamic liquid,solid phases
@@ -35 +36 @@
      &      ,RA    ,RG    ,R1SA                                         &
      &      ,RSIGMA                                                     &
-     &      ,R     ,RMD   ,RMO3  ,RMCO2, RMC, RMV   ,RD    ,RV    ,RCPD &
+     &      ,R     ,RMD   ,RMV   ,RD    ,RV    ,RCPD                    &
+     &      ,RMO3  ,RMCO2 ,RMC   ,RMCH4 ,RMN2O ,RMCFC11 ,RMCFC12        &
      &      ,RCPV  ,RCVD  ,RCVV  ,RKAPPA,RETV                           &
      &      ,RCW   ,RCS                                                 &

LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

@@ -38 +38 @@
 
 ! Variables read from parmeter file physiq.def
+  LOGICAL, PUBLIC :: carbon_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
+!$OMP THREADPRIVATE(carbon_cycle_cpl)
   LOGICAL, PUBLIC :: carbon_cycle_tr        ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys
 !$OMP THREADPRIVATE(carbon_cycle_tr)
-  LOGICAL, PUBLIC :: carbon_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
-!$OMP THREADPRIVATE(carbon_cycle_cpl)
+  LOGICAL, PUBLIC :: carbon_cycle_rad       ! CO2 interactive radiatively 
+!$OMP THREADPRIVATE(carbon_cycle_rad)
   INTEGER, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3
 !$OMP THREADPRIVATE(level_coupling_esm)
+  REAL, PUBLIC :: RCO2_glo, RCO2_tot
+!$OMP THREADPRIVATE(RCO2_glo, RCO2_tot)
 
   LOGICAL :: carbon_cycle_emis_comp_omp=.FALSE.
@@ -91 +95 @@
   REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: co2_send ! Field allocated in phyetat0
 !$OMP THREADPRIVATE(co2_send)
+
+  INTEGER, PARAMETER, PUBLIC :: id_CO2=1              !--temporaire OB -- to be changed
 
 ! nbfields : total number of fields 
@@ -552 +558 @@
           
           ! Transformation of atmospheric CO2 concentration for the radiation code
-          RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 
+          RCO2 = co2_ppm * 1.0e-06  * RMCO2/28.97 
           
           WRITE(lunout,*) 'RCO2 is now updated! RCO2 = ', RCO2

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -26 +26 @@
     USE surface_data
     USE phys_cal_mod
-    USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl, level_coupling_esm
+    USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl, carbon_cycle_rad, level_coupling_esm
     USE mod_grid_phy_lmdz, ONLY: klon_glo
     USE print_control_mod, ONLY: lunout
@@ -231 +231 @@
     LOGICAL, SAVE :: carbon_cycle_tr_omp
     LOGICAL, SAVE :: carbon_cycle_cpl_omp
+    LOGICAL, SAVE :: carbon_cycle_rad_omp
     INTEGER, SAVE :: level_coupling_esm_omp
     LOGICAL, SAVE :: adjust_tropopause_omp
@@ -2156 +2157 @@
     carbon_cycle_cpl_omp=.FALSE.
     CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp)
+
+    carbon_cycle_rad_omp=.FALSE.
+    CALL getin('carbon_cycle_rad',carbon_cycle_rad_omp)
 
     ! >> PC
@@ -2428 +2432 @@
     carbon_cycle_tr = carbon_cycle_tr_omp
     carbon_cycle_cpl = carbon_cycle_cpl_omp
+    carbon_cycle_rad = carbon_cycle_rad_omp
     level_coupling_esm = level_coupling_esm_omp
 
@@ -2529 +2534 @@
     IF (flag_bc_internal_mixture .AND. flag_aerosol.NE.6) THEN 
        CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1)
+    ENDIF
+
+    ! Test on carbon cycle
+    IF (carbon_cycle_tr .AND. .NOT. carbon_cycle_cpl) THEN
+       CALL abort_physic('conf_phys', 'carbon_cycle_cpl has to be TRUE if carbon_cycle_tr is on',1)
+    ENDIF
+    IF (carbon_cycle_rad .AND. .NOT. carbon_cycle_cpl) THEN
+       CALL abort_physic('conf_phys', 'carbon_cycle_cpl has to be TRUE if carbon_cycle_rad is on',1)
     ENDIF
 
@@ -2742 +2755 @@
     WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr
     WRITE(lunout,*) ' carbon_cycle_cpl = ', carbon_cycle_cpl
+    WRITE(lunout,*) ' carbon_cycle_rad = ', carbon_cycle_rad
     WRITE(lunout,*) ' level_coupling_esm = ', level_coupling_esm
 

LMDZ6/trunk/libf/phylmd/suphel.F90

@@ -118 +118 @@
   rmo3 = 47.9942
   rmco2 = 44.011
+  rmch4 = 16.043
+  rmn2o = 44.013
+  rmcfc11 = 137.3686
+  rmcfc12 = 120.9140
   rmc   = 12.0107
   rmv = 18.0153
@@ -133 +137 @@
   WRITE (UNIT=6, FMT='('' Ozone   mass = '',e13.7)') rmo3
   WRITE (UNIT=6, FMT='('' CO2     mass = '',e13.7)') rmco2
-  WRITE (UNIT=6, FMT='('' CO2     mass = '',e13.7)') rmc
+  WRITE (UNIT=6, FMT='('' C       mass = '',e13.7)') rmc
+  WRITE (UNIT=6, FMT='('' CH4     mass = '',e13.7)') rmch4
+  WRITE (UNIT=6, FMT='('' N2O     mass = '',e13.7)') rmn2o
+  WRITE (UNIT=6, FMT='('' CFC11   mass = '',e13.7)') rmcfc11
+  WRITE (UNIT=6, FMT='('' CFC12   mass = '',e13.7)') rmcfc12
   WRITE (UNIT=6, FMT='('' Vapour  mass = '',e13.7)') rmv
   WRITE (UNIT=6, FMT='('' Dry air cst. = '',e13.7)') rd