Changeset 2136

Timestamp:

Oct 27, 2014, 6:47:03 PM (10 years ago)

Author:

lguez

Message:

Add the module from ecmwf model for water production by methane oxydation
and photolysis. Activated if def file option: ok_qch4=y.
LG, D.Cugnet

Location:

LMDZ5/trunk/libf/phylmd

Files:

• YOEMETH.F90 (added)
• clesphys.h (modified) (2 diffs)
• conf_phys_m.F90 (modified) (4 diffs)
• methox.F90 (added)
• phys_local_var_mod.F90 (modified) (4 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq.F90 (modified) (1 diff)
• sumethox.F90 (added)
• suphel.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/clesphys.h

@@ -77 +77 @@
        LOGICAL :: ok_strato
        LOGICAL :: ok_hines, ok_gwd_rando
+       LOGICAL :: ok_qch4
        LOGICAL :: ok_conserv_q
 
@@ -114 +115 @@
      &     , ip_ebil_phy                                                &
      &     , ok_lic_melt,           aer_type                            &
-     &     , iflag_rrtm, ok_strato,ok_hines                             &
+     &     , iflag_rrtm, ok_strato,ok_hines, ok_qch4                    &
      &     , iflag_ice_thermo, ok_gwd_rando, NSW                        &
      &     , ok_conserv_q, ok_all_xml

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -189 +189 @@
     LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
     real, SAVE:: GWD_RANDO_RUWMAX_omp, gwd_rando_sat_omp
+    LOGICAL,SAVE :: ok_qch4_omp
     LOGICAL,SAVE      :: carbon_cycle_tr_omp
     LOGICAL,SAVE      :: carbon_cycle_cpl_omp
@@ -1725 +1726 @@
     gwd_rando_sat_omp = 0.25
     CALL getin('gwd_rando_sat', gwd_rando_sat_omp)
+
+    !Config  key = ok_qch4
+    !Config  Desc = activation de la parametrisation du methane
+    !Config  Def  = .FALSE.
+    !Config  Help = Clef controlant l'activation de la parametrisation
+    !               de l'humidite due a oxydation+photolyse du methane strato
+
+    ok_qch4_omp=.FALSE.
+    CALL getin('ok_qch4',ok_qch4_omp)
 
     !Config Key  = OK_LES                                               
@@ -1955 +1965 @@
     GWD_RANDO_RUWMAX = GWD_RANDO_RUWMAX_omp
     gwd_rando_sat = gwd_rando_sat_omp
+    ok_qch4 = ok_qch4_omp
     ok_LES = ok_LES_omp
     callstats = callstats_omp
@@ -2147 +2158 @@
     write(lunout,*) 'ok_hines = ',  ok_hines
     write(lunout,*) 'ok_gwd_rando = ',  ok_gwd_rando
+    write(lunout,*) 'ok_qch4 = ',  ok_qch4
     write(lunout,*) 'GWD_RANDO_RUWMAX = ', GWD_RANDO_RUWMAX
     write(lunout,*) 'gwd_rando_sat = ', gwd_rando_sat

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -73 +73 @@
       REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:)
       !$OMP THREADPRIVATE(d_t_hin)
+! tendance due a l'oxydation du methane
+      REAL, SAVE, ALLOCATABLE :: d_q_ch4(:,:)
+      !$OMP THREADPRIVATE(d_q_ch4)
 
 ! tendance du a la conersion Ec -> E thermique
@@ -338 +341 @@
       allocate(topswcf_aero(klon,3), solswcf_aero(klon,3))
       allocate(d_u_hin(klon,klev),d_v_hin(klon,klev),d_t_hin(klon,klev))
+      allocate(d_q_ch4(klon,klev))
 !      allocate(tausum_aero(klon,nwave,naero_spc))
 !      allocate(tau3d_aero(klon,klev,nwave,naero_spc)) 
@@ -535 +539 @@
       deallocate(load_tmp7)
       deallocate(d_u_hin,d_v_hin,d_t_hin)
+      deallocate(d_q_ch4)
 
 !IM ajout variables CFMIP2/CMIP5
@@ -603 +608 @@
 END SUBROUTINE phys_local_var_end
 
-END MODULE phys_local_var_mod
+END MODULE phys_local_var_mod 

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1100 +1100 @@
   TYPE(ctrl_out), SAVE :: o_dthin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     'dthin', 'Hines GWD dT', 'K/s', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_dqch4 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'dqch4', 'H2O due to CH4 oxidation & photolysis', '(kg/kg)/s', (/ ('', i=1, 9) /))
 
   type(ctrl_out), save:: o_du_gwd_rando &

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -130 +130 @@
          o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
          o_dtoro, o_dulif, o_dvlif, o_dtlif, &
-         o_duhin, o_dvhin, o_dthin, o_rsu, &
+         o_duhin, o_dvhin, o_dthin, o_dqch4, o_rsu, &
          o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, &
          o_rlucs, o_rldcs, o_tnt, o_tntr, &
@@ -217 +217 @@
          d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
          d_v_lif, d_t_lif, d_u_hin, d_v_hin, d_t_hin, &
-         pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD
+         d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD
 
     USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &
@@ -1153 +1153 @@
        end IF
 
+       IF (ok_qch4) then
+          CALL histwrite_phy(o_dqch4, d_q_ch4 / pdtphys)
+       ENDIF
+
        CALL histwrite_phy(o_rsu, swup)
        CALL histwrite_phy(o_rsd, swdn)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -3495 +3495 @@
           , fs_bound, fq_bound )
   END IF
+
+  !DC Calcul de la tendance due au methane
+  IF(ok_qch4) THEN
+     CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
+  ! ajout de la tendance d'humidite due au methane
+     CALL add_phys_tend(du0,dv0,dt0,d_q_ch4*dtime,dql0,'q_ch4')
+  END IF
   !
   !

LMDZ5/trunk/libf/phylmd/suphel.F90

@@ -192 +192 @@
   ! ------------------------------------------------------------------
 
-  ! calculer les constantes pour les fonctions thermodynamiques
+  ! *       10.   CONSTANTS FOR THERMODYNAMICAL FUNCTIONS.
+  ! ----------------------------------------
 
   rvtmp2 = rcpv/rcpd - 1.
@@ -204 +205 @@
   r5ies = r3ies*(rtt-r4ies)
 
+  ! ------------------------------------------------------------------
+
+  ! *       10.   CONSTANTS FOR METHANE OXIDATION AND PHOTOLYSIS.
+  ! -----------------------------------------------
+
+  CALL SUMETHOX()
+
   RETURN
 END SUBROUTINE suphel