source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/rrtm/rrtm_taumol10.F90 @ 2174

Last change on this file since 2174 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM

Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 3.0 KB
Line 
1!*******************************************************************************
2SUBROUTINE RRTM_TAUMOL10 (KLEV,P_TAU,&
3 & P_TAUAERL,P_FAC00,P_FAC01,P_FAC10,P_FAC11,K_JP,K_JT,K_JT1,&
4 & P_COLH2O,K_LAYTROP,PFRAC) 
5
6!     BAND 10:  1390-1480 cm-1 (low - H2O; high - H2O)
7
8! Modifications
9!        M.Hamrud      01-Oct-2003 CY28 Cleaning
10
11!     D Salmond   2000-05-15 speed-up
12!     JJMorcrette 2000-05-17 speed-up
13
14USE PARKIND1  ,ONLY : JPIM     ,JPRB
15USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
16
17USE PARRRTM  , ONLY : JPLAY  ,JPBAND ,JPGPT  ,NG10   ,NGS9
18USE YOERRTWN , ONLY :      NSPA   ,NSPB
19USE YOERRTA10, ONLY : ABSA   ,ABSB   ,FRACREFA, FRACREFB
20
21!  Input
22!#include "yoeratm.h"
23
24!      REAL TAUAER(JPLAY)
25
26IMPLICIT NONE
27
28!  Output
29INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
30REAL(KIND=JPRB)   ,INTENT(OUT)   :: P_TAU(JPGPT,JPLAY)
31REAL(KIND=JPRB)   ,INTENT(IN)    :: P_TAUAERL(JPLAY,JPBAND)
32REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC00(JPLAY)
33REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC01(JPLAY)
34REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC10(JPLAY)
35REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC11(JPLAY)
36INTEGER(KIND=JPIM),INTENT(IN)    :: K_JP(JPLAY)
37INTEGER(KIND=JPIM),INTENT(IN)    :: K_JT(JPLAY)
38INTEGER(KIND=JPIM),INTENT(IN)    :: K_JT1(JPLAY)
39REAL(KIND=JPRB)   ,INTENT(IN)    :: P_COLH2O(JPLAY)
40INTEGER(KIND=JPIM),INTENT(IN)    :: K_LAYTROP
41REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFRAC(JPGPT,JPLAY)
42!- from AER
43!- from INTFAC     
44!- from INTIND
45!- from PROFDATA             
46!- from SP             
47INTEGER(KIND=JPIM) :: IND0(JPLAY),IND1(JPLAY)
48
49INTEGER(KIND=JPIM) :: IG, I_LAY
50REAL(KIND=JPRB) :: ZHOOK_HANDLE
51
52!      EQUIVALENCE (TAUAERL(1,10),TAUAER)
53
54!     Compute the optical depth by interpolating in ln(pressure) and
55!     temperature. 
56
57IF (LHOOK) CALL DR_HOOK('RRTM_TAUMOL10',0,ZHOOK_HANDLE)
58DO I_LAY = 1, K_LAYTROP
59  IND0(I_LAY) = ((K_JP(I_LAY)-1)*5+(K_JT(I_LAY)-1))*NSPA(10) + 1
60  IND1(I_LAY) = (K_JP(I_LAY)*5+(K_JT1(I_LAY)-1))*NSPA(10) + 1
61ENDDO
62
63!-- DS_000515
64DO IG = 1, NG10
65  DO I_LAY = 1, K_LAYTROP
66!-- DS_000515
67    P_TAU (NGS9+IG,I_LAY) = P_COLH2O(I_LAY) *&
68     & (P_FAC00(I_LAY) * ABSA(IND0(I_LAY)  ,IG) +&
69     & P_FAC10(I_LAY) * ABSA(IND0(I_LAY)+1,IG) +&
70     & P_FAC01(I_LAY) * ABSA(IND1(I_LAY)  ,IG) +&
71     & P_FAC11(I_LAY) * ABSA(IND1(I_LAY)+1,IG)) &
72     & + P_TAUAERL(I_LAY,10) 
73    PFRAC(NGS9+IG,I_LAY) = FRACREFA(IG)
74  ENDDO
75ENDDO
76
77DO I_LAY = K_LAYTROP+1, KLEV
78  IND0(I_LAY) = ((K_JP(I_LAY)-13)*5+(K_JT(I_LAY)-1))*NSPB(10) + 1
79  IND1(I_LAY) = ((K_JP(I_LAY)-12)*5+(K_JT1(I_LAY)-1))*NSPB(10) + 1
80ENDDO
81
82!-- JJM_000517
83DO IG = 1, NG10
84  DO I_LAY = K_LAYTROP+1, KLEV
85!-- JJM_000517
86    P_TAU (NGS9+IG,I_LAY) = P_COLH2O(I_LAY) *&
87     & (P_FAC00(I_LAY) * ABSB(IND0(I_LAY)  ,IG) +&
88     & P_FAC10(I_LAY) * ABSB(IND0(I_LAY)+1,IG) +&
89     & P_FAC01(I_LAY) * ABSB(IND1(I_LAY)  ,IG) +&
90     & P_FAC11(I_LAY) * ABSB(IND1(I_LAY)+1,IG)) &
91     & + P_TAUAERL(I_LAY,10) 
92    PFRAC(NGS9+IG,I_LAY) = FRACREFB(IG)
93  ENDDO
94ENDDO
95
96IF (LHOOK) CALL DR_HOOK('RRTM_TAUMOL10',1,ZHOOK_HANDLE)
97END SUBROUTINE RRTM_TAUMOL10
Note: See TracBrowser for help on using the repository browser.