source: LMDZ5/trunk/libf/phymar/rrtm_kgb16_00.F90

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

Inclusion de la physique de MAR


Integration of MAR physics

File size: 7.0 KB
RevLine 
[2089]1!OCL SCALAR
2SUBROUTINE RRTM_KGB16_00
3
4!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
5!     BAND 16:  2600-3000 cm-1 (low - H2O,CH4; high - nothing)
6!     Reformatted for F90 by JJMorcrette, ECMWF
7!     Reformatted for NEC by H.Gallée   , LGGE  (splitting)
8
9!     ------------------------------------------------------------------
10
11#include "tsmbkind.h"
12
13USE YOERRTO16, ONLY : KAO     ,SELFREFO   ,FRACREFAO
14USE YOERRTA16, ONLY : STRRAT
15
16!     ------------------------------------------------------------------
17
18
19IMPLICIT NONE
20STRRAT = 830.411_JPRB
21
22!     ------------------------------------------------------------------
23
24!     The array SELFREFO contains the coefficient of the water vapor
25!     self-continuum (including the energy term).  The first index
26!     refers to temperature in 7.2 degree increments.  For instance,
27!     JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8,
28!     etc.  The second index runs over the g-channel (1 to 16).
29
30SELFREFO( :, 1) = (/&
31&1.26476E-03_JPRB, 1.04962E-03_JPRB, 8.71080E-04_JPRB, 7.22907E-04_JPRB, 5.99939E-04_JPRB,&
32&4.97888E-04_JPRB, 4.13197E-04_JPRB, 3.42911E-04_JPRB, 2.84581E-04_JPRB, 2.36173E-04_JPRB/)
33SELFREFO( :, 2) = (/&
34&1.27380E-03_JPRB, 1.05867E-03_JPRB, 8.79870E-04_JPRB, 7.31270E-04_JPRB, 6.07767E-04_JPRB,&
35&5.05121E-04_JPRB, 4.19812E-04_JPRB, 3.48911E-04_JPRB, 2.89983E-04_JPRB, 2.41008E-04_JPRB/)
36SELFREFO( :, 3) = (/&
37&1.26332E-03_JPRB, 1.05867E-03_JPRB, 8.87172E-04_JPRB, 7.43454E-04_JPRB, 6.23018E-04_JPRB,&
38&5.22092E-04_JPRB, 4.37515E-04_JPRB, 3.66640E-04_JPRB, 3.07246E-04_JPRB, 2.57474E-04_JPRB/)
39SELFREFO( :, 4) = (/&
40&1.28078E-03_JPRB, 1.07802E-03_JPRB, 9.07363E-04_JPRB, 7.63722E-04_JPRB, 6.42820E-04_JPRB,&
41&5.41057E-04_JPRB, 4.55404E-04_JPRB, 3.83311E-04_JPRB, 3.22630E-04_JPRB, 2.71556E-04_JPRB/)
42SELFREFO( :, 5) = (/&
43&1.61188E-03_JPRB, 1.32375E-03_JPRB, 1.08712E-03_JPRB, 8.92794E-04_JPRB, 7.33203E-04_JPRB,&
44&6.02140E-04_JPRB, 4.94505E-04_JPRB, 4.06110E-04_JPRB, 3.33516E-04_JPRB, 2.73898E-04_JPRB/)
45SELFREFO( :, 6) = (/&
46&1.42549E-03_JPRB, 1.25964E-03_JPRB, 1.11309E-03_JPRB, 9.83595E-04_JPRB, 8.69162E-04_JPRB,&
47&7.68043E-04_JPRB, 6.78688E-04_JPRB, 5.99728E-04_JPRB, 5.29955E-04_JPRB, 4.68299E-04_JPRB/)
48SELFREFO( :, 7) = (/&
49&1.43253E-03_JPRB, 1.28535E-03_JPRB, 1.15330E-03_JPRB, 1.03480E-03_JPRB, 9.28488E-04_JPRB,&
50&8.33094E-04_JPRB, 7.47501E-04_JPRB, 6.70702E-04_JPRB, 6.01794E-04_JPRB, 5.39965E-04_JPRB/)
51SELFREFO( :, 8) = (/&
52&1.48973E-03_JPRB, 1.33348E-03_JPRB, 1.19363E-03_JPRB, 1.06844E-03_JPRB, 9.56378E-04_JPRB,&
53&8.56072E-04_JPRB, 7.66286E-04_JPRB, 6.85918E-04_JPRB, 6.13978E-04_JPRB, 5.49583E-04_JPRB/)
54SELFREFO( :, 9) = (/&
55&1.53002E-03_JPRB, 1.44607E-03_JPRB, 1.36674E-03_JPRB, 1.29175E-03_JPRB, 1.22088E-03_JPRB,&
56&1.15390E-03_JPRB, 1.09059E-03_JPRB, 1.03076E-03_JPRB, 9.74209E-04_JPRB, 9.20760E-04_JPRB/)
57SELFREFO( :,10) = (/&
58&1.78480E-03_JPRB, 1.68865E-03_JPRB, 1.59768E-03_JPRB, 1.51161E-03_JPRB, 1.43018E-03_JPRB,&
59&1.35313E-03_JPRB, 1.28023E-03_JPRB, 1.21126E-03_JPRB, 1.14601E-03_JPRB, 1.08427E-03_JPRB/)
60SELFREFO( :,11) = (/&
61&1.85121E-03_JPRB, 1.73903E-03_JPRB, 1.63364E-03_JPRB, 1.53464E-03_JPRB, 1.44164E-03_JPRB,&
62&1.35427E-03_JPRB, 1.27220E-03_JPRB, 1.19510E-03_JPRB, 1.12268E-03_JPRB, 1.05464E-03_JPRB/)
63SELFREFO( :,12) = (/&
64&2.10772E-03_JPRB, 1.96994E-03_JPRB, 1.84118E-03_JPRB, 1.72082E-03_JPRB, 1.60834E-03_JPRB,&
65&1.50321E-03_JPRB, 1.40495E-03_JPRB, 1.31311E-03_JPRB, 1.22728E-03_JPRB, 1.14705E-03_JPRB/)
66SELFREFO( :,13) = (/&
67&1.99323E-03_JPRB, 1.89279E-03_JPRB, 1.79741E-03_JPRB, 1.70684E-03_JPRB, 1.62083E-03_JPRB,&
68&1.53915E-03_JPRB, 1.46159E-03_JPRB, 1.38794E-03_JPRB, 1.31800E-03_JPRB, 1.25159E-03_JPRB/)
69SELFREFO( :,14) = (/&
70&2.05835E-03_JPRB, 1.93301E-03_JPRB, 1.81529E-03_JPRB, 1.70475E-03_JPRB, 1.60093E-03_JPRB,&
71&1.50344E-03_JPRB, 1.41189E-03_JPRB, 1.32591E-03_JPRB, 1.24516E-03_JPRB, 1.16934E-03_JPRB/)
72SELFREFO( :,15) = (/&
73&2.34120E-03_JPRB, 2.20012E-03_JPRB, 2.06754E-03_JPRB, 1.94295E-03_JPRB, 1.82587E-03_JPRB,&
74&1.71585E-03_JPRB, 1.61245E-03_JPRB, 1.51529E-03_JPRB, 1.42398E-03_JPRB, 1.33817E-03_JPRB/)
75SELFREFO( :,16) = (/&
76&3.03559E-03_JPRB, 2.80023E-03_JPRB, 2.58312E-03_JPRB, 2.38284E-03_JPRB, 2.19809E-03_JPRB,&
77&2.02767E-03_JPRB, 1.87045E-03_JPRB, 1.72543E-03_JPRB, 1.59165E-03_JPRB, 1.46825E-03_JPRB/)
78
79FRACREFAO( :, 1) = (/&
80!     From P = 862.6 mb.
81    &0.17356300_JPRB,0.18880001_JPRB,0.17704099_JPRB,0.13661300_JPRB,&
82    &0.10691600_JPRB,0.08222480_JPRB,0.05939860_JPRB,0.04230810_JPRB,&
83    &0.02526330_JPRB,0.00244532_JPRB,0.00193541_JPRB,0.00150415_JPRB,&
84    &0.00103528_JPRB,0.00067068_JPRB,0.00024951_JPRB,0.00003348_JPRB/)
85FRACREFAO( :, 2) = (/&
86    &0.17779499_JPRB,0.19837400_JPRB,0.16557600_JPRB,0.13470000_JPRB,&
87    &0.11013600_JPRB,0.08342720_JPRB,0.05987030_JPRB,0.03938700_JPRB,&
88    &0.02293650_JPRB,0.00238849_JPRB,0.00192400_JPRB,0.00149921_JPRB,&
89    &0.00103539_JPRB,0.00067150_JPRB,0.00024822_JPRB,0.00003348_JPRB/)
90FRACREFAO( :, 3) = (/&
91    &0.18535601_JPRB,0.19407199_JPRB,0.16053200_JPRB,0.13300700_JPRB,&
92    &0.10779000_JPRB,0.08408500_JPRB,0.06480450_JPRB,0.04070160_JPRB,&
93    &0.02203590_JPRB,0.00227779_JPRB,0.00189074_JPRB,0.00146888_JPRB,&
94    &0.00103147_JPRB,0.00066770_JPRB,0.00024751_JPRB,0.00003348_JPRB/)
95FRACREFAO( :, 4) = (/&
96    &0.19139200_JPRB,0.18917400_JPRB,0.15748601_JPRB,0.13240699_JPRB,&
97    &0.10557300_JPRB,0.08383260_JPRB,0.06724060_JPRB,0.04364450_JPRB,&
98    &0.02175820_JPRB,0.00225436_JPRB,0.00184421_JPRB,0.00143153_JPRB,&
99    &0.00103027_JPRB,0.00066066_JPRB,0.00024222_JPRB,0.00003148_JPRB/)
100FRACREFAO( :, 5) = (/&
101    &0.19547801_JPRB,0.18539500_JPRB,0.15442000_JPRB,0.13114899_JPRB,&
102    &0.10515600_JPRB,0.08350350_JPRB,0.06909780_JPRB,0.04671630_JPRB,&
103    &0.02168820_JPRB,0.00224400_JPRB,0.00182009_JPRB,0.00139098_JPRB,&
104    &0.00102582_JPRB,0.00065367_JPRB,0.00023202_JPRB,0.00003148_JPRB/)
105FRACREFAO( :, 6) = (/&
106    &0.19757500_JPRB,0.18266800_JPRB,0.15208900_JPRB,0.12897800_JPRB,&
107    &0.10637200_JPRB,0.08391220_JPRB,0.06989830_JPRB,0.04964120_JPRB,&
108    &0.02155800_JPRB,0.00224310_JPRB,0.00177358_JPRB,0.00138184_JPRB,&
109    &0.00101538_JPRB,0.00063370_JPRB,0.00023227_JPRB,0.00003148_JPRB/)
110FRACREFAO( :, 7) = (/&
111    &0.20145500_JPRB,0.17692900_JPRB,0.14940600_JPRB,0.12690400_JPRB,&
112    &0.10828800_JPRB,0.08553720_JPRB,0.07004940_JPRB,0.05153430_JPRB,&
113    &0.02268740_JPRB,0.00216943_JPRB,0.00178603_JPRB,0.00137754_JPRB,&
114    &0.00098344_JPRB,0.00063165_JPRB,0.00023218_JPRB,0.00003148_JPRB/)
115FRACREFAO( :, 8) = (/&
116    &0.20383500_JPRB,0.17047501_JPRB,0.14570600_JPRB,0.12679300_JPRB,&
117    &0.11043100_JPRB,0.08719150_JPRB,0.07045440_JPRB,0.05345420_JPRB,&
118    &0.02448340_JPRB,0.00215839_JPRB,0.00175893_JPRB,0.00138296_JPRB,&
119    &0.00098318_JPRB,0.00063188_JPRB,0.00023199_JPRB,0.00003148_JPRB/)
120FRACREFAO( :, 9) = (/&
121    &0.18680701_JPRB,0.15961801_JPRB,0.15092900_JPRB,0.13049100_JPRB,&
122    &0.11418400_JPRB,0.09380540_JPRB,0.07093450_JPRB,0.05664280_JPRB,&
123    &0.02938410_JPRB,0.00217751_JPRB,0.00176766_JPRB,0.00138275_JPRB,&
124    &0.00098377_JPRB,0.00063181_JPRB,0.00023193_JPRB,0.00003148_JPRB/)
125
126!     -----------------------------------------------------------------
127RETURN
128END SUBROUTINE RRTM_KGB16_00
Note: See TracBrowser for help on using the repository browser.