1 |
|
|
!OPTIONS XOPT(NOEVAL) |
2 |
|
|
SUBROUTINE SURDI15 |
3 |
|
|
|
4 |
|
|
!**** *SURDI15* - INITIALIZE COMMON YOMRDI15 CONTROLLING RADINT |
5 |
|
|
!**** FROZEN VERSION (CYCLE 15) OF SURDI |
6 |
|
|
|
7 |
|
|
! PURPOSE. |
8 |
|
|
! -------- |
9 |
|
|
! INITIALIZE YOMRDI15, THE COMMON THAT CONTROLS THE |
10 |
|
|
! RADIATION INTERFACE |
11 |
|
|
|
12 |
|
|
!** INTERFACE. |
13 |
|
|
! ---------- |
14 |
|
|
! CALL *SURDI15* FROM *SUECRAD* |
15 |
|
|
! ------- ------- |
16 |
|
|
|
17 |
|
|
! EXPLICIT ARGUMENTS : |
18 |
|
|
! -------------------- |
19 |
|
|
! NONE |
20 |
|
|
|
21 |
|
|
! IMPLICIT ARGUMENTS : |
22 |
|
|
! -------------------- |
23 |
|
|
! COMMON YOMRDI15 |
24 |
|
|
|
25 |
|
|
! METHOD. |
26 |
|
|
! ------- |
27 |
|
|
! SEE DOCUMENTATION |
28 |
|
|
|
29 |
|
|
! EXTERNALS. |
30 |
|
|
! ---------- |
31 |
|
|
! NONE |
32 |
|
|
|
33 |
|
|
! REFERENCE. |
34 |
|
|
! ---------- |
35 |
|
|
! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE |
36 |
|
|
! "IN CORE MODEL" |
37 |
|
|
|
38 |
|
|
! AUTHOR. |
39 |
|
|
! ------- |
40 |
|
|
! 96-11: Ph. Dandin. Meteo-France |
41 |
|
|
! ORIGINAL : 88-12-15 BY JEAN-JACQUES MORCRETTE *ECMWF* |
42 |
|
|
|
43 |
|
|
! MODIFICATIONS. |
44 |
|
|
! -------------- |
45 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
46 |
|
|
! A.Alias 05-12-2005 greenhouse gases variables (M.Deque) |
47 |
|
|
! A.Alias 13-06-2006 RI0 value can be changed via namscen.h |
48 |
|
|
|
49 |
|
|
! ------------------------------------------------------------------ |
50 |
|
|
|
51 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
52 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
53 |
|
|
! Ce qui concerne NULNAM commente par MPL le 15.04.09 |
54 |
|
|
!USE YOMLUN , ONLY : NULNAM |
55 |
|
|
USE YOMLUN_IFSAUX , ONLY : NULOUT |
56 |
|
|
|
57 |
|
|
USE YOMCST , ONLY : RI0 |
58 |
|
|
USE YOMRDI15 , ONLY : RSDTSN15 ,RRAE15 ,RMU0015 ,RALBICE15,& |
59 |
|
|
& RALBSEA15,RALBSNM15,RALBSNO15,RCARDI15 ,REMISS15 ,& |
60 |
|
|
& RSNOWAL15,RVLBDC15 ,RCH415 ,RN2O15 ,RCFC1115 ,& |
61 |
|
|
& RCFC1215 ,REPALB15 ,REPCLC15 ,REPH2O15 |
62 |
|
|
USE YOMRDU15 , ONLY : REPSEC15 |
63 |
|
|
|
64 |
|
|
IMPLICIT NONE |
65 |
|
|
|
66 |
|
|
REAL(KIND=JPRB) :: XCARDI, XCFC11, XCFC12, XCH4, XN2O |
67 |
|
|
REAL(KIND=JPRB) :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZSUPSAT |
68 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
69 |
|
|
|
70 |
|
|
#include "namscen.h" |
71 |
|
|
#include "posnam.intfb.h" |
72 |
|
|
! ---------------------------------------------------------------- |
73 |
|
|
|
74 |
|
|
!* 1. SET DEFAULT VALUES. |
75 |
|
|
! ------------------- |
76 |
|
|
|
77 |
|
|
IF (LHOOK) CALL DR_HOOK('SURDI15',0,ZHOOK_HANDLE) |
78 |
|
|
RSDTSN15 = 5.0_JPRB |
79 |
|
|
|
80 |
|
|
RRAE15 = 0.1277E-02_JPRB |
81 |
|
|
RMU0015 = RRAE15/SQRT(RRAE15*(RRAE15+2.0_JPRB)) |
82 |
|
|
|
83 |
|
|
RALBICE15 = 0.55_JPRB |
84 |
|
|
RALBSEA15 = 0.07_JPRB |
85 |
|
|
RALBSNO15 = 0.80_JPRB |
86 |
|
|
RALBSNM15 = 0.40_JPRB |
87 |
|
|
RSNOWAL15 = 0.01_JPRB |
88 |
|
|
!* Concentration of the various trace gases (IPCC/SACC values for 1990) |
89 |
|
|
! CO2 CH4 N2O CFC11 CFC12 |
90 |
|
|
! 353ppmv 1.72ppmv 310ppbv 280pptv 484pptv |
91 |
|
|
|
92 |
|
|
XCARDI = 353.E-06_JPRB |
93 |
|
|
XCH4 = 1.72E-06_JPRB |
94 |
|
|
XN2O = 310.E-09_JPRB |
95 |
|
|
XCFC11 = 280.E-12_JPRB |
96 |
|
|
XCFC12 = 484.E-12_JPRB |
97 |
|
|
|
98 |
|
|
ZAIRMWG = 28.970_JPRB |
99 |
|
|
ZCO2MWG = 44.011_JPRB |
100 |
|
|
ZCH4MWG = 16.043_JPRB |
101 |
|
|
ZN2OMWG = 44.013_JPRB |
102 |
|
|
ZC11MWG = 137.3686_JPRB |
103 |
|
|
ZC12MWG = 120.9140_JPRB |
104 |
|
|
|
105 |
|
|
! Ce qui concerne NAMSCEN commente par MPL le 15.04.09 |
106 |
|
|
!CALL POSNAM(NULNAM,'NAMSCEN') |
107 |
|
|
!READ (NULNAM, NAMSCEN) |
108 |
|
|
|
109 |
|
|
WRITE(NULOUT,'( & |
110 |
|
|
& '' CO2 = '',E14.7,'' CH4 = '',E14.7 & |
111 |
|
|
& ,'' N2O = '',E14.7,'' CFC11 = '',E14.7 & |
112 |
|
|
& ,'' CFC12 = '',E14.7,'' RI0 = '',E14.7 & |
113 |
|
|
& )') XCARDI,XCH4,XN2O,XCFC11,XCFC12,RI0 |
114 |
|
|
|
115 |
|
|
RCARDI15 = XCARDI*ZCO2MWG/ZAIRMWG |
116 |
|
|
RCH415 = XCH4*ZCH4MWG/ZAIRMWG |
117 |
|
|
RN2O15 = XN2O*ZN2OMWG/ZAIRMWG |
118 |
|
|
RCFC1115 = XCFC11*ZC11MWG/ZAIRMWG |
119 |
|
|
RCFC1215 = XCFC12*ZC12MWG/ZAIRMWG |
120 |
|
|
REMISS15 = 0.996_JPRB |
121 |
|
|
!ZSUPSAT = 0.01_JPRB |
122 |
|
|
RVLBDC15 = 0.5_JPRB |
123 |
|
|
|
124 |
|
|
REPSEC15=1.E-12_JPRB |
125 |
|
|
REPCLC15=1.E-12_JPRB |
126 |
|
|
REPH2O15=1.E-12_JPRB |
127 |
|
|
REPALB15=1.E-12_JPRB |
128 |
|
|
|
129 |
|
|
! ----------------------------------------------------------------- |
130 |
|
|
|
131 |
|
|
IF (LHOOK) CALL DR_HOOK('SURDI15',1,ZHOOK_HANDLE) |
132 |
|
|
END SUBROUTINE SURDI15 |