Changeset 1737

Timestamp:

Mar 29, 2013, 6:58:35 PM (12 years ago)

Author:

idelkadi

Message:

Ajout de la cles logique ok_4xCO2atm controlant le 2e appel au code de transfert radiatif (calcule du forcage radiatif quand la concentration de CO2 est quadriplee) !

Location:

LMDZ5/trunk/libf/phylmd

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys_m.F90 (modified) (4 diffs)
• physiq.F (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/clesphys.h

@@ -19 +19 @@
 !IM ajout CFMIP2/CMIP5
        REAL co2_ppm_per
+       LOGICAL ok_4xCO2atm
        REAL(kind=8) RCO2_per,RCH4_per,RN2O_per,RCFC11_per,RCFC12_per
        REAL(kind=8) CH4_ppb_per,N2O_ppb_per,CFC11_ppt_per,CFC12_ppt_per
@@ -81 +82 @@
      &     , RCO2, RCH4, RN2O, RCFC11, RCFC12                           &
      &     , RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act       &
+     &     , ok_4xCO2atm                                                & 
      &     , RCO2_per, RCH4_per, RN2O_per, RCFC11_per, RCFC12_per       &
      &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -130 +130 @@
   REAL,SAVE :: CFC11_ppt_omp,RCFC11_omp,CFC11_ppt_per_omp,RCFC11_per_omp
   REAL,SAVE :: CFC12_ppt_omp,RCFC12_omp,CFC12_ppt_per_omp,RCFC12_per_omp
+  LOGICAL,SAVE :: ok_4xCO2atm_omp
   REAL,SAVE :: epmax_omp
   LOGICAL,SAVE :: ok_adj_ema_omp
@@ -532 +533 @@
   RCO2_per_omp = co2_ppm_per_omp * 1.0e-06  * 44.011/28.97
 
+!Config Key  = ok_4xCO2atm
+!Config Desc = Calcul ou non effet radiatif 4xco2
+!Config Def  = .false.
+!Config Help = 
+
+  ok_4xCO2atm_omp = .false.
+  call getin('ok_4xCO2atm',ok_4xCO2atm_omp)
+
 !Config Key  = RCH4_per
 !Config Desc = Concentration du CH4_per
@@ -1613 +1622 @@
     RCFC11_act = RCFC11
     RCFC12_act = RCFC12
+    ok_4xCO2atm = ok_4xCO2atm_omp
     RCO2_per = RCO2_per_omp
     RCH4_per = RCH4_per_omp
@@ -1821 +1831 @@
   write(lunout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
   write(lunout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
+  write(lunout,*)' ok_4xCO2atm=',ok_4xCO2atm
   write(lunout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
   write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -3316 +3316 @@
 cIM Par defaut on a les taux perturbes egaux aux taux actuels
 c
+      if (ok_4xCO2atm) then
        if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR.
      $RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR.
@@ -3352 +3353 @@
      o        topswcf_aerop, solswcf_aerop)
        endif
+      endif
 c
       ENDIF ! aerosol_couple