source: LMDZ6/trunk/libf/phymar/rrtm_kgb15_00.F90 @ 3956

Last change on this file since 3956 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
Line 
1!OCL SCALAR
2SUBROUTINE RRTM_KGB15_00
3
4!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
5!     BAND 15:  2380-2600 cm-1 (low - N2O,CO2; 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 YOERRTO15, ONLY : KAO     ,SELFREFO   ,FRACREFAO
14USE YOERRTA15, ONLY : STRRAT
15
16!     ------------------------------------------------------------------
17
18
19IMPLICIT NONE
20STRRAT = 0.2883201_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.72146E-03_JPRB, 1.41213E-03_JPRB, 1.15839E-03_JPRB, 9.50236E-04_JPRB, 7.79489E-04_JPRB,&
32&6.39423E-04_JPRB, 5.24525E-04_JPRB, 4.30274E-04_JPRB, 3.52958E-04_JPRB, 2.89535E-04_JPRB/)
33SELFREFO( :, 2) = (/&
34&1.83092E-03_JPRB, 1.50082E-03_JPRB, 1.23023E-03_JPRB, 1.00843E-03_JPRB, 8.26611E-04_JPRB,&
35&6.77578E-04_JPRB, 5.55414E-04_JPRB, 4.55276E-04_JPRB, 3.73192E-04_JPRB, 3.05908E-04_JPRB/)
36SELFREFO( :, 3) = (/&
37&1.92620E-03_JPRB, 1.57871E-03_JPRB, 1.29391E-03_JPRB, 1.06048E-03_JPRB, 8.69169E-04_JPRB,&
38&7.12370E-04_JPRB, 5.83857E-04_JPRB, 4.78528E-04_JPRB, 3.92201E-04_JPRB, 3.21447E-04_JPRB/)
39SELFREFO( :, 4) = (/&
40&2.11224E-03_JPRB, 1.72966E-03_JPRB, 1.41638E-03_JPRB, 1.15984E-03_JPRB, 9.49760E-04_JPRB,&
41&7.77734E-04_JPRB, 6.36867E-04_JPRB, 5.21514E-04_JPRB, 4.27055E-04_JPRB, 3.49704E-04_JPRB/)
42SELFREFO( :, 5) = (/&
43&2.31367E-03_JPRB, 1.88988E-03_JPRB, 1.54371E-03_JPRB, 1.26095E-03_JPRB, 1.02999E-03_JPRB,&
44&8.41325E-04_JPRB, 6.87221E-04_JPRB, 5.61344E-04_JPRB, 4.58523E-04_JPRB, 3.74536E-04_JPRB/)
45SELFREFO( :, 6) = (/&
46&2.47614E-03_JPRB, 2.03107E-03_JPRB, 1.66599E-03_JPRB, 1.36654E-03_JPRB, 1.12091E-03_JPRB,&
47&9.19433E-04_JPRB, 7.54170E-04_JPRB, 6.18612E-04_JPRB, 5.07419E-04_JPRB, 4.16213E-04_JPRB/)
48SELFREFO( :, 7) = (/&
49&2.55137E-03_JPRB, 2.10353E-03_JPRB, 1.73430E-03_JPRB, 1.42987E-03_JPRB, 1.17889E-03_JPRB,&
50&9.71958E-04_JPRB, 8.01350E-04_JPRB, 6.60688E-04_JPRB, 5.44717E-04_JPRB, 4.49103E-04_JPRB/)
51SELFREFO( :, 8) = (/&
52&2.59317E-03_JPRB, 2.13957E-03_JPRB, 1.76532E-03_JPRB, 1.45653E-03_JPRB, 1.20175E-03_JPRB,&
53&9.91539E-04_JPRB, 8.18098E-04_JPRB, 6.74996E-04_JPRB, 5.56926E-04_JPRB, 4.59508E-04_JPRB/)
54SELFREFO( :, 9) = (/&
55&2.83495E-03_JPRB, 2.34704E-03_JPRB, 1.94309E-03_JPRB, 1.60867E-03_JPRB, 1.33181E-03_JPRB,&
56&1.10259E-03_JPRB, 9.12826E-04_JPRB, 7.55721E-04_JPRB, 6.25656E-04_JPRB, 5.17975E-04_JPRB/)
57SELFREFO( :,10) = (/&
58&2.98196E-03_JPRB, 2.46746E-03_JPRB, 2.04173E-03_JPRB, 1.68946E-03_JPRB, 1.39796E-03_JPRB,&
59&1.15676E-03_JPRB, 9.57177E-04_JPRB, 7.92028E-04_JPRB, 6.55373E-04_JPRB, 5.42297E-04_JPRB/)
60SELFREFO( :,11) = (/&
61&2.98241E-03_JPRB, 2.46776E-03_JPRB, 2.04191E-03_JPRB, 1.68955E-03_JPRB, 1.39800E-03_JPRB,&
62&1.15675E-03_JPRB, 9.57141E-04_JPRB, 7.91974E-04_JPRB, 6.55308E-04_JPRB, 5.42226E-04_JPRB/)
63SELFREFO( :,12) = (/&
64&2.97837E-03_JPRB, 2.46464E-03_JPRB, 2.03952E-03_JPRB, 1.68773E-03_JPRB, 1.39662E-03_JPRB,&
65&1.15572E-03_JPRB, 9.56375E-04_JPRB, 7.91413E-04_JPRB, 6.54904E-04_JPRB, 5.41942E-04_JPRB/)
66SELFREFO( :,13) = (/&
67&2.97777E-03_JPRB, 2.46461E-03_JPRB, 2.03989E-03_JPRB, 1.68835E-03_JPRB, 1.39740E-03_JPRB,&
68&1.15659E-03_JPRB, 9.57271E-04_JPRB, 7.92305E-04_JPRB, 6.55767E-04_JPRB, 5.42759E-04_JPRB/)
69SELFREFO( :,14) = (/&
70&2.98331E-03_JPRB, 2.46947E-03_JPRB, 2.04413E-03_JPRB, 1.69205E-03_JPRB, 1.40061E-03_JPRB,&
71&1.15937E-03_JPRB, 9.59683E-04_JPRB, 7.94388E-04_JPRB, 6.57564E-04_JPRB, 5.44306E-04_JPRB/)
72SELFREFO( :,15) = (/&
73&2.99414E-03_JPRB, 2.47830E-03_JPRB, 2.05133E-03_JPRB, 1.69792E-03_JPRB, 1.40540E-03_JPRB,&
74&1.16328E-03_JPRB, 9.62863E-04_JPRB, 7.96978E-04_JPRB, 6.59673E-04_JPRB, 5.46023E-04_JPRB/)
75SELFREFO( :,16) = (/&
76&3.00006E-03_JPRB, 2.48297E-03_JPRB, 2.05501E-03_JPRB, 1.70081E-03_JPRB, 1.40766E-03_JPRB,&
77&1.16503E-03_JPRB, 9.64227E-04_JPRB, 7.98033E-04_JPRB, 6.60484E-04_JPRB, 5.46643E-04_JPRB/)
78
79FRACREFAO( :, 1) = (/&
80!     From P = 1053.6 mb.
81    &0.11287100_JPRB,0.12070200_JPRB,0.12729000_JPRB,0.12858100_JPRB,&
82    &0.12743001_JPRB,0.11961800_JPRB,0.10290400_JPRB,0.07888980_JPRB,&
83    &0.05900120_JPRB,0.00667979_JPRB,0.00552926_JPRB,0.00436993_JPRB,&
84    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
85FRACREFAO( :, 2) = (/&
86    &0.13918801_JPRB,0.16353001_JPRB,0.16155800_JPRB,0.14090499_JPRB,&
87    &0.11322300_JPRB,0.08757720_JPRB,0.07225720_JPRB,0.05173390_JPRB,&
88    &0.04731360_JPRB,0.00667979_JPRB,0.00552926_JPRB,0.00436993_JPRB,&
89    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
90FRACREFAO( :, 3) = (/&
91    &0.14687300_JPRB,0.17853101_JPRB,0.15664500_JPRB,0.13351700_JPRB,&
92    &0.10791200_JPRB,0.08684320_JPRB,0.07158090_JPRB,0.05198410_JPRB,&
93    &0.04340110_JPRB,0.00667979_JPRB,0.00552926_JPRB,0.00436993_JPRB,&
94    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
95FRACREFAO( :, 4) = (/&
96    &0.15760700_JPRB,0.17759100_JPRB,0.15158001_JPRB,0.13193300_JPRB,&
97    &0.10742800_JPRB,0.08693760_JPRB,0.07159490_JPRB,0.05196250_JPRB,&
98    &0.04065270_JPRB,0.00667979_JPRB,0.00552926_JPRB,0.00436993_JPRB,&
99    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
100FRACREFAO( :, 5) = (/&
101    &0.16646700_JPRB,0.17299300_JPRB,0.15018500_JPRB,0.13138700_JPRB,&
102    &0.10735900_JPRB,0.08713110_JPRB,0.07130330_JPRB,0.05279420_JPRB,&
103    &0.03766730_JPRB,0.00667979_JPRB,0.00552926_JPRB,0.00436993_JPRB,&
104    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
105FRACREFAO( :, 6) = (/&
106    &0.17546000_JPRB,0.16666500_JPRB,0.14969499_JPRB,0.13105400_JPRB,&
107    &0.10782500_JPRB,0.08718610_JPRB,0.07156770_JPRB,0.05308320_JPRB,&
108    &0.03753960_JPRB,0.00432465_JPRB,0.00509623_JPRB,0.00436993_JPRB,&
109    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
110FRACREFAO( :, 7) = (/&
111    &0.18378501_JPRB,0.16064601_JPRB,0.14940400_JPRB,0.13146400_JPRB,&
112    &0.10810300_JPRB,0.08775740_JPRB,0.07115360_JPRB,0.05400040_JPRB,&
113    &0.03689970_JPRB,0.00388333_JPRB,0.00323610_JPRB,0.00353414_JPRB,&
114    &0.00320611_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
115FRACREFAO( :, 8) = (/&
116    &0.18966800_JPRB,0.15744300_JPRB,0.14993000_JPRB,0.13152599_JPRB,&
117    &0.10899200_JPRB,0.08858690_JPRB,0.07142920_JPRB,0.05399600_JPRB,&
118    &0.03433460_JPRB,0.00374886_JPRB,0.00302066_JPRB,0.00240653_JPRB,&
119    &0.00199205_JPRB,0.00204765_JPRB,0.00077371_JPRB,0.00010894_JPRB/)
120FRACREFAO( :, 9) = (/&
121    &0.11887100_JPRB,0.12479600_JPRB,0.12569501_JPRB,0.12839900_JPRB,&
122    &0.12473500_JPRB,0.12012800_JPRB,0.11086700_JPRB,0.08493590_JPRB,&
123    &0.05063770_JPRB,0.00328723_JPRB,0.00266849_JPRB,0.00210232_JPRB,&
124    &0.00152114_JPRB,0.00095635_JPRB,0.00035374_JPRB,0.00004980_JPRB/)
125
126!     -----------------------------------------------------------------
127RETURN
128END SUBROUTINE RRTM_KGB15_00
Note: See TracBrowser for help on using the repository browser.