source: LMDZ6/branches/DYNAMICO-conv/libf/phymar/surdi.F90 @ 3024

Last change on this file since 3024 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 3.7 KB
Line 
1SUBROUTINE SURDI
2
3
4!**** *SURDI*   - INITIALIZE COMMON YOERDI CONTROLLING RADINT
5
6!     PURPOSE.
7!     --------
8!           INITIALIZE YOERDI, THE COMMON THAT CONTROLS THE
9!           RADIATION INTERFACE
10
11!**   INTERFACE.
12!     ----------
13!        CALL *SURDI* FROM *SURAD*
14!              ------        -----
15
16!        EXPLICIT ARGUMENTS :
17!        --------------------
18!        NONE
19
20!        IMPLICIT ARGUMENTS :
21!        --------------------
22!        COMMON YOERDI
23
24!     METHOD.
25!     -------
26!        SEE DOCUMENTATION
27
28!     EXTERNALS.
29!     ----------
30!        NONE
31
32!     REFERENCE.
33!     ----------
34!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS MODEL
35
36!     AUTHOR.
37!     -------
38!        Original  JEAN-JACQUES MORCRETTE  *ECMWF*
39!        Modified   P. Viterbo   99-03-26    Tiling of the land surface
40
41!     MODIFICATIONS.
42!     --------------
43!        ORIGINAL : 88-12-15
44!     ------------------------------------------------------------------
45
46#include "tsmbkind.h"
47
48USE YOERDI   , ONLY : RRAE     ,RALBSEAD ,&
49            &RALBICEVS_AR,RALBICENI_AR,RALBICEVS_AN,RALBICENI_AN,&
50            &RALBSFO  ,REMISD   ,REMISL   ,REMISN   ,REMISS   ,&
51            &RCARDI   ,RCH4     ,RN2O     ,RO3      ,RCFC11   ,&
52            &RCFC12   ,REPALB   ,REPCLC   ,REPH2O   ,RSUNDUR
53
54IMPLICIT NONE
55
56
57!     LOCAL REAL SCALARS
58REAL_B :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZO3MWG
59INTEGER_M :: JM, IM
60
61!      ----------------------------------------------------------------
62
63!*       1.    SET DEFAULT VALUES.
64!              -------------------
65
66RRAE = 0.1277E-02_JPRB
67
68!* Threshold for computing sunshine duration (W/m2)
69RSUNDUR=120._JPRB
70
71!*  Ocean surface albedo for diffuse radiation (Taylor et al., 1997)     
72RALBSEAD = 0.06_JPRB
73!*  For sea ice, monthly values are based on Ebert and Curry, 1993, Table 2.
74!   We take dry snow albedo as the representative value for non-summer
75!   months, and bare sea-ice as the representative value for summer
76!   months. The values for Antarctic are shifted six-months.
77!*  Sea ice surf. albedo for visible rad. (snow covered; Ebert and Curry, 1993)
78RALBICEVS_AR(1:12) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
79                      &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
80                      &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
81!*  Sea ice surf. albedo for near-IR rad. (snow covered; Ebert and Curry, 1993)
82RALBICENI_AR(1:12) = (/0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB,&
83                      &0.664_JPRB,0.476_JPRB,0.288_JPRB,0.288_JPRB,&
84                      &0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB/)
85DO JM=1,12
86  IM=MOD(JM+5,12)+1
87  RALBICEVS_AN(JM)=RALBICEVS_AR(IM)
88  RALBICENI_AN(JM)=RALBICENI_AR(IM)
89ENDDO
90!*  Snow albedo in the presence of high vegetation
91RALBSFO=0.15_JPRB
92
93!- sea surface emissivity and other surfaces outside the window region     
94REMISS  = 0.99_JPRB
95!- snow window emissivity     
96REMISN  = 0.98_JPRB
97!- land window emissivity     
98REMISL  = 0.96_JPRB
99!- desert window emissivity (lower bound when dry)     
100REMISD  = 0.93_JPRB
101
102!*  Concentration of the various trace gases (IPCC/SACC values for 1990)
103!        CO2         CH4        N2O        CFC11       CFC12
104!      353ppmv     1.72ppmv   310ppbv     280pptv     484pptv
105
106ZAIRMWG = 28.970_JPRB
107ZCO2MWG = 44.011_JPRB
108ZCH4MWG = 16.043_JPRB
109ZN2OMWG = 44.013_JPRB
110ZO3MWG  = 47.9982_JPRB
111ZC11MWG = 137.3686_JPRB
112ZC12MWG = 120.9140_JPRB
113
114RCARDI  = 353.E-06_JPRB*ZCO2MWG/ZAIRMWG
115RCH4    = 1.72E-06_JPRB*ZCH4MWG/ZAIRMWG
116RN2O    = 310.E-09_JPRB*ZN2OMWG/ZAIRMWG
117RO3     =   1.E-06_JPRB*ZO3MWG /ZAIRMWG
118RCFC11  = 280.E-12_JPRB*ZC11MWG/ZAIRMWG
119RCFC12  = 484.E-12_JPRB*ZC12MWG/ZAIRMWG
120
121REPCLC=1.E-12_JPRB
122REPH2O=1.E-12_JPRB
123REPALB=1.E-12_JPRB
124
125!     -----------------------------------------------------------------
126
127RETURN
128END SUBROUTINE SURDI
Note: See TracBrowser for help on using the repository browser.