source: LMDZ6/trunk/libf/phylmd/sumethox.F90 @ 3654

Last change on this file since 3654 was 2136, checked in by lguez, 10 years ago

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
File size: 1.5 KB
Line 
1SUBROUTINE SUMETHOX
2
3!**** *SUCLDP*   - INITIALIZE MODULE YOEMETH CONTROLLING *METHOX*
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOEMETH
8
9!**   INTERFACE.
10!     ----------
11!        CALL *SUMETHOX* FROM *SUPHEC*
12!              --------        ------
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        MODULE YOEMETH
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28!        NONE
29
30!     REFERENCE.
31!     ----------
32!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
33!     "INTEGRATED FORECASTING SYSTEM"
34
35!     AUTHOR.
36!     -------
37!        C.JAKOB   *ECMWF*
38
39!     MODIFICATIONS.
40!     --------------
41!        ORIGINAL : 98-04-07
42!        Modified : 02-01-29  A.Simmons: increase RQLIM from 3.75e-6
43!        M.Hamrud      01-Oct-2003 CY28 Cleaning
44!        D. Cugnet     24-Feb-2012 Adapted for LMDZ
45!     ------------------------------------------------------------------
46
47USE YOEMETH   , ONLY : RALPHA1 ,RALPHA2  ,RQLIM   ,&
48 & RPBOTOX,  RPBOTPH ,RPTOPOX  ,RPTOPPH ,&
49 & RALPHA3,  RLOGPPH 
50
51!*       1.    SET VALUES
52!              ----------
53
54IMPLICIT NONE
55RALPHA1=(19.*LOG(10.))/(LOG(20.)**4)
56RALPHA2=LOG(1.0/3.+0.01)
57RQLIM=4.25E-6
58RPBOTOX=10000.
59RPBOTPH=20.
60RPTOPOX=50.
61RPTOPPH=0.1
62RALPHA3=0.5*(LOG(100.)+RALPHA2)
63RLOGPPH=LOG(RPTOPPH/RPBOTPH)
64
65!     -----------------------------------------------------------------
66
67END SUBROUTINE SUMETHOX
68
69
Note: See TracBrowser for help on using the repository browser.