1 |
|
|
! |
2 |
|
|
! $Id: surdi.F90 2626 2016-09-15 14:20:56Z musat $ |
3 |
|
|
! |
4 |
|
1 |
SUBROUTINE SURDI |
5 |
|
|
|
6 |
|
|
!**** *SURDI* - INITIALIZE COMMON YOERDI CONTROLLING RADINT |
7 |
|
|
|
8 |
|
|
! PURPOSE. |
9 |
|
|
! -------- |
10 |
|
|
! INITIALIZE YOERDI, THE COMMON THAT CONTROLS THE |
11 |
|
|
! RADIATION INTERFACE |
12 |
|
|
|
13 |
|
|
!** INTERFACE. |
14 |
|
|
! ---------- |
15 |
|
|
! CALL *SURDI* FROM *SURAD* |
16 |
|
|
! ----- ----- |
17 |
|
|
|
18 |
|
|
! EXPLICIT ARGUMENTS : |
19 |
|
|
! -------------------- |
20 |
|
|
! NONE |
21 |
|
|
|
22 |
|
|
! IMPLICIT ARGUMENTS : |
23 |
|
|
! -------------------- |
24 |
|
|
! COMMON YOERDI |
25 |
|
|
|
26 |
|
|
! METHOD. |
27 |
|
|
! ------- |
28 |
|
|
! SEE DOCUMENTATION |
29 |
|
|
|
30 |
|
|
! EXTERNALS. |
31 |
|
|
! ---------- |
32 |
|
|
! NONE |
33 |
|
|
|
34 |
|
|
! REFERENCE. |
35 |
|
|
! ---------- |
36 |
|
|
! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS MODEL |
37 |
|
|
|
38 |
|
|
! AUTHOR. |
39 |
|
|
! ------- |
40 |
|
|
! Original JEAN-JACQUES MORCRETTE *ECMWF* |
41 |
|
|
! Modified P. Viterbo 99-03-26 Tiling of the land surface |
42 |
|
|
! Modified P. Viterbo 24-05-2004 surf library |
43 |
|
|
|
44 |
|
|
! MODIFICATIONS. |
45 |
|
|
! -------------- |
46 |
|
|
! ORIGINAL : 88-12-15 |
47 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
48 |
|
|
! JJMorcrette 2004-10-07 Gas concentrations |
49 |
|
|
! ------------------------------------------------------------------ |
50 |
|
|
|
51 |
|
|
USE PARKIND1 ,ONLY : JPRB |
52 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
53 |
|
|
|
54 |
|
|
USE YOERDI , ONLY : RRAE ,& |
55 |
|
|
! & RCARDI ,RCH4 ,RN2O ,RO3 ,RCFC11 ,& |
56 |
|
|
! & RCFC12, & |
57 |
|
|
& REPCLC ,REPH2O ,RSUNDUR ,& |
58 |
|
|
& RCCO2 ,RCCH4 ,RCN2O ,RCCFC11 ,RCCFC12 |
59 |
|
|
|
60 |
|
|
|
61 |
|
|
IMPLICIT NONE |
62 |
|
|
|
63 |
|
|
!MPL/IM 20160915 on prend GES de phylmd |
64 |
|
|
!#include "clesphys.h" |
65 |
|
|
REAL(KIND=JPRB) :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZO3MWG |
66 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
67 |
|
|
|
68 |
|
|
! ---------------------------------------------------------------- |
69 |
|
|
|
70 |
|
|
!* 1. SET DEFAULT VALUES. |
71 |
|
|
! ------------------- |
72 |
|
|
|
73 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SURDI',0,ZHOOK_HANDLE) |
74 |
|
1 |
RRAE = 0.1277E-02_JPRB |
75 |
|
|
|
76 |
|
|
!* Threshold for computing sunshine duration (W/m2) |
77 |
|
1 |
RSUNDUR=120._JPRB |
78 |
|
|
|
79 |
|
|
!* For sea ice, monthly values are based on Ebert and Curry, 1993, Table 2. |
80 |
|
|
! We take dry snow albedo as the representative value for non-summer |
81 |
|
|
! months, and bare sea-ice as the representative value for summer |
82 |
|
|
! months. The values for Antarctic are shifted six-months. |
83 |
|
|
! All computations brought back to *SUSWN* |
84 |
|
|
|
85 |
|
|
!* Concentration of the various trace gases (IPCC/SACC values for 1990) |
86 |
|
|
! CO2 CH4 N2O CFC11 CFC12 |
87 |
|
|
! 353ppmv 1.72ppmv 310ppbv 280pptv 484pptv |
88 |
|
|
|
89 |
|
|
ZAIRMWG = 28.970_JPRB |
90 |
|
|
ZCO2MWG = 44.011_JPRB |
91 |
|
|
ZCH4MWG = 16.043_JPRB |
92 |
|
|
ZN2OMWG = 44.013_JPRB |
93 |
|
|
ZO3MWG = 47.9982_JPRB |
94 |
|
|
ZC11MWG = 137.3686_JPRB |
95 |
|
|
ZC12MWG = 120.9140_JPRB |
96 |
|
|
|
97 |
|
|
!RCARDI = 353.E-06_JPRB*ZCO2MWG/ZAIRMWG |
98 |
|
|
!RCH4 = 1.72E-06_JPRB*ZCH4MWG/ZAIRMWG |
99 |
|
|
!RN2O = 310.E-09_JPRB*ZN2OMWG/ZAIRMWG |
100 |
|
|
!RO3 = 1.E-06_JPRB*ZO3MWG /ZAIRMWG |
101 |
|
|
!RCFC11 = 280.E-12_JPRB*ZC11MWG/ZAIRMWG |
102 |
|
|
!RCFC12 = 484.E-12_JPRB*ZC12MWG/ZAIRMWG |
103 |
|
|
|
104 |
|
|
!RCARDI = RCCO2 * ZCO2MWG/ZAIRMWG |
105 |
|
|
!RCH4 = RCCH4 * ZCH4MWG/ZAIRMWG |
106 |
|
|
!RN2O = RCN2O * ZN2OMWG/ZAIRMWG |
107 |
|
|
!RO3 = 1.E-06_JPRB*ZO3MWG /ZAIRMWG |
108 |
|
|
!RCFC11 = RCCFC11 * ZC11MWG/ZAIRMWG |
109 |
|
|
!RCFC12 = RCCFC12 * ZC12MWG/ZAIRMWG |
110 |
|
|
|
111 |
|
|
|
112 |
|
1 |
REPCLC=1.E-12_JPRB |
113 |
|
1 |
REPH2O=1.E-12_JPRB |
114 |
|
|
|
115 |
|
|
|
116 |
|
|
! ----------------------------------------------------------------- |
117 |
|
|
|
118 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SURDI',1,ZHOOK_HANDLE) |
119 |
|
1 |
END SUBROUTINE SURDI |