phylmd/Radlwsw/lw.f

module lw_m

  IMPLICIT none

contains

  SUBROUTINE LW(PPMB, PDP, PDT0, PEMIS, PTL, PTAVE, PWV, POZON, PAER, PCLDLD, &
       PCLDLU, PVIEW, PCOLR, PCOLR0, PTOPLW, PSOLLW, PTOPLW0, PSOLLW0, &
       psollwdown, plwup, plwdn, plwup0, plwdn0)

    use LWU_m, only: LWU
    USE suphec_m, ONLY: md, rcpd, rday, rg, rmo3
    USE raddim, ONLY: kdlon, kflev
    USE raddimlw, ONLY: nua

    ! Method.

    ! 1. Computes the pressure and temperature weighted amounts of
    ! absorbers.

    ! 2. Computes the planck functions on the interfaces and the
    ! gradient of planck functions in the layers.

    ! 3. Performs the vertical integration distinguishing the con-
    ! tributions of the adjacent and distant layers and those from the
    ! boundaries.

    ! 4. Computes the clear-sky downward and upward emissivities.

    ! 5. Introduces the effects of the clouds on the fluxes.

    ! Reference: see radiation's part of the model's documentation and
    ! ECMWF research department documentation of the IFS

    ! Author:
    ! Jean-Jacques Morcrette *ECMWF*

    ! Original : 89-07-14

    DOUBLE PRECISION PCLDLD(KDLON, KFLEV) ! DOWNWARD EFFECTIVE CLOUD COVER
    DOUBLE PRECISION PCLDLU(KDLON, KFLEV) ! UPWARD EFFECTIVE CLOUD COVER
    DOUBLE PRECISION PDP(KDLON, KFLEV) ! LAYER PRESSURE THICKNESS (Pa)
    DOUBLE PRECISION PDT0(KDLON) ! SURFACE TEMPERATURE DISCONTINUITY (K)
    DOUBLE PRECISION PEMIS(KDLON) ! SURFACE EMISSIVITY
    DOUBLE PRECISION PPMB(KDLON, KFLEV+1) ! HALF LEVEL PRESSURE (mb)
    DOUBLE PRECISION POZON(KDLON, KFLEV) ! O3 CONCENTRATION (kg/kg)
    DOUBLE PRECISION PTL(KDLON, KFLEV+1) ! HALF LEVEL TEMPERATURE (K)
    DOUBLE PRECISION PAER(KDLON, KFLEV, 5) ! OPTICAL THICKNESS OF THE AEROSOLS
    DOUBLE PRECISION PTAVE(KDLON, KFLEV) ! LAYER TEMPERATURE (K)
    DOUBLE PRECISION PVIEW(KDLON) ! COSECANT OF VIEWING ANGLE
    DOUBLE PRECISION PWV(KDLON, KFLEV) ! SPECIFIC HUMIDITY (kg/kg)

    DOUBLE PRECISION PCOLR(KDLON, KFLEV) ! LONG-WAVE TENDENCY (K/day)
    DOUBLE PRECISION PCOLR0(KDLON, KFLEV) ! LONG-WAVE TENDENCY (K/day) clear-sky
    DOUBLE PRECISION PTOPLW(KDLON) ! LONGWAVE FLUX AT T.O.A.
    DOUBLE PRECISION PSOLLW(KDLON) ! LONGWAVE FLUX AT SURFACE
    DOUBLE PRECISION PTOPLW0(KDLON) ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
    DOUBLE PRECISION PSOLLW0(KDLON) ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
    ! Rajout LF
    double precision psollwdown(kdlon) ! LONGWAVE downwards flux at surface
    !IM
    DOUBLE PRECISION plwup(KDLON, KFLEV+1) ! LW up total sky
    DOUBLE PRECISION plwup0(KDLON, KFLEV+1) ! LW up clear sky
    DOUBLE PRECISION plwdn(KDLON, KFLEV+1) ! LW down total sky
    DOUBLE PRECISION plwdn0(KDLON, KFLEV+1) ! LW down clear sky

    DOUBLE PRECISION ZABCU(KDLON, NUA, 3*KFLEV+1)
    DOUBLE PRECISION ZOZ(KDLON, KFLEV)

    DOUBLE PRECISION ZFLUX(KDLON, 2, KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
    DOUBLE PRECISION ZFLUC(KDLON, 2, KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
    DOUBLE PRECISION ZBINT(KDLON, KFLEV+1) ! Intermediate variable
    DOUBLE PRECISION ZBSUI(KDLON) ! Intermediate variable
    DOUBLE PRECISION ZCTS(KDLON, KFLEV) ! Intermediate variable
    DOUBLE PRECISION ZCNTRB(KDLON, KFLEV+1, KFLEV+1) ! Intermediate variable
    SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB

    INTEGER ilim, i, k, kpl1

    INTEGER, PARAMETER:: lw0pas = 1 ! Every lw0pas steps, clear-sky is done
    INTEGER, PARAMETER:: lwpas = 1 ! Every lwpas steps, cloudy-sky is done
    ! In general, lw0pas and lwpas should be 1

    INTEGER:: itaplw0 = 0, itaplw = 0

    ! ------------------------------------------------------------------

    IF (MOD(itaplw0, lw0pas) == 0) THEN
       DO k = 1, KFLEV
          DO i = 1, KDLON
             ! convertir ozone de kg/kg en pa (modif MPL 100505)
             ZOZ(i, k) = POZON(i, k)*PDP(i, k) * MD/RMO3
          ENDDO
       ENDDO
       CALL LWU(PAER, PDP, PPMB, ZOZ, PTAVE, PVIEW, PWV, ZABCU)
       CALL LWBV(ILIM, PDP, PDT0, PEMIS, PPMB, PTL, PTAVE, ZABCU, &
            ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB)
       itaplw0 = 0
    ENDIF
    itaplw0 = itaplw0 + 1

    IF (MOD(itaplw, lwpas) == 0) THEN
       CALL LWC(ILIM, PCLDLD, PCLDLU, PEMIS, &
            ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB, &
            ZFLUX)
       itaplw = 0
    ENDIF
    itaplw = itaplw + 1

    DO k = 1, KFLEV
       kpl1 = k+1
       DO i = 1, KDLON
          PCOLR(i, k) = ZFLUX(i, 1, kpl1)+ZFLUX(i, 2, kpl1) &
               - ZFLUX(i, 1, k)- ZFLUX(i, 2, k)
          PCOLR(i, k) = PCOLR(i, k) * RDAY*RG/RCPD / PDP(i, k)
          PCOLR0(i, k) = ZFLUC(i, 1, kpl1)+ZFLUC(i, 2, kpl1) &
               - ZFLUC(i, 1, k)- ZFLUC(i, 2, k)
          PCOLR0(i, k) = PCOLR0(i, k) * RDAY*RG/RCPD / PDP(i, k)
       ENDDO
    ENDDO
    DO i = 1, KDLON
       PSOLLW(i) = -ZFLUX(i, 1, 1)-ZFLUX(i, 2, 1)
       PTOPLW(i) = ZFLUX(i, 1, KFLEV+1) + ZFLUX(i, 2, KFLEV+1)

       PSOLLW0(i) = -ZFLUC(i, 1, 1)-ZFLUC(i, 2, 1)
       PTOPLW0(i) = ZFLUC(i, 1, KFLEV+1) + ZFLUC(i, 2, KFLEV+1)
       psollwdown(i) = -ZFLUX(i, 2, 1)

       !IM attention aux signes !; LWtop >0, LWdn < 0
       DO k = 1, KFLEV+1
          plwup(i, k) = ZFLUX(i, 1, k)
          plwup0(i, k) = ZFLUC(i, 1, k)
          plwdn(i, k) = ZFLUX(i, 2, k)
          plwdn0(i, k) = ZFLUC(i, 2, k)
       ENDDO
    ENDDO

  END SUBROUTINE LW

end module lw_m
1	guez	71	module lw_m
2
3			IMPLICIT none
4
5			contains
6
7	guez	145	SUBROUTINE LW(PPMB, PDP, PDT0, PEMIS, PTL, PTAVE, PWV, POZON, PAER, PCLDLD, &
8			PCLDLU, PVIEW, PCOLR, PCOLR0, PTOPLW, PSOLLW, PTOPLW0, PSOLLW0, &
9	guez	71	psollwdown, plwup, plwdn, plwup0, plwdn0)
10
11			use LWU_m, only: LWU
12			USE suphec_m, ONLY: md, rcpd, rday, rg, rmo3
13			USE raddim, ONLY: kdlon, kflev
14			USE raddimlw, ONLY: nua
15
16			! Method.
17
18			! 1. Computes the pressure and temperature weighted amounts of
19			! absorbers.
20
21			! 2. Computes the planck functions on the interfaces and the
22			! gradient of planck functions in the layers.
23
24			! 3. Performs the vertical integration distinguishing the con-
25			! tributions of the adjacent and distant layers and those from the
26			! boundaries.
27
28			! 4. Computes the clear-sky downward and upward emissivities.
29
30			! 5. Introduces the effects of the clouds on the fluxes.
31
32			! Reference: see radiation's part of the model's documentation and
33			! ECMWF research department documentation of the IFS
34
35			! Author:
36			! Jean-Jacques Morcrette ECMWF
37
38			! Original : 89-07-14
39
40			DOUBLE PRECISION PCLDLD(KDLON, KFLEV) ! DOWNWARD EFFECTIVE CLOUD COVER
41			DOUBLE PRECISION PCLDLU(KDLON, KFLEV) ! UPWARD EFFECTIVE CLOUD COVER
42			DOUBLE PRECISION PDP(KDLON, KFLEV) ! LAYER PRESSURE THICKNESS (Pa)
43			DOUBLE PRECISION PDT0(KDLON) ! SURFACE TEMPERATURE DISCONTINUITY (K)
44			DOUBLE PRECISION PEMIS(KDLON) ! SURFACE EMISSIVITY
45			DOUBLE PRECISION PPMB(KDLON, KFLEV+1) ! HALF LEVEL PRESSURE (mb)
46			DOUBLE PRECISION POZON(KDLON, KFLEV) ! O3 CONCENTRATION (kg/kg)
47			DOUBLE PRECISION PTL(KDLON, KFLEV+1) ! HALF LEVEL TEMPERATURE (K)
48			DOUBLE PRECISION PAER(KDLON, KFLEV, 5) ! OPTICAL THICKNESS OF THE AEROSOLS
49			DOUBLE PRECISION PTAVE(KDLON, KFLEV) ! LAYER TEMPERATURE (K)
50			DOUBLE PRECISION PVIEW(KDLON) ! COSECANT OF VIEWING ANGLE
51			DOUBLE PRECISION PWV(KDLON, KFLEV) ! SPECIFIC HUMIDITY (kg/kg)
52
53			DOUBLE PRECISION PCOLR(KDLON, KFLEV) ! LONG-WAVE TENDENCY (K/day)
54			DOUBLE PRECISION PCOLR0(KDLON, KFLEV) ! LONG-WAVE TENDENCY (K/day) clear-sky
55			DOUBLE PRECISION PTOPLW(KDLON) ! LONGWAVE FLUX AT T.O.A.
56			DOUBLE PRECISION PSOLLW(KDLON) ! LONGWAVE FLUX AT SURFACE
57			DOUBLE PRECISION PTOPLW0(KDLON) ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
58			DOUBLE PRECISION PSOLLW0(KDLON) ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
59			! Rajout LF
60			double precision psollwdown(kdlon) ! LONGWAVE downwards flux at surface
61			!IM
62			DOUBLE PRECISION plwup(KDLON, KFLEV+1) ! LW up total sky
63			DOUBLE PRECISION plwup0(KDLON, KFLEV+1) ! LW up clear sky
64			DOUBLE PRECISION plwdn(KDLON, KFLEV+1) ! LW down total sky
65			DOUBLE PRECISION plwdn0(KDLON, KFLEV+1) ! LW down clear sky
66
67			DOUBLE PRECISION ZABCU(KDLON, NUA, 3*KFLEV+1)
68			DOUBLE PRECISION ZOZ(KDLON, KFLEV)
69
70			DOUBLE PRECISION ZFLUX(KDLON, 2, KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
71			DOUBLE PRECISION ZFLUC(KDLON, 2, KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
72			DOUBLE PRECISION ZBINT(KDLON, KFLEV+1) ! Intermediate variable
73			DOUBLE PRECISION ZBSUI(KDLON) ! Intermediate variable
74			DOUBLE PRECISION ZCTS(KDLON, KFLEV) ! Intermediate variable
75			DOUBLE PRECISION ZCNTRB(KDLON, KFLEV+1, KFLEV+1) ! Intermediate variable
76			SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB
77
78			INTEGER ilim, i, k, kpl1
79
80			INTEGER, PARAMETER:: lw0pas = 1 ! Every lw0pas steps, clear-sky is done
81			INTEGER, PARAMETER:: lwpas = 1 ! Every lwpas steps, cloudy-sky is done
82			! In general, lw0pas and lwpas should be 1
83
84			INTEGER:: itaplw0 = 0, itaplw = 0
85
86			! ------------------------------------------------------------------
87
88			IF (MOD(itaplw0, lw0pas) == 0) THEN
89			DO k = 1, KFLEV
90			DO i = 1, KDLON
91			! convertir ozone de kg/kg en pa (modif MPL 100505)
92			ZOZ(i, k) = POZON(i, k)PDP(i, k) MD/RMO3
93			ENDDO
94			ENDDO
95	guez	139	CALL LWU(PAER, PDP, PPMB, ZOZ, PTAVE, PVIEW, PWV, ZABCU)
96	guez	71	CALL LWBV(ILIM, PDP, PDT0, PEMIS, PPMB, PTL, PTAVE, ZABCU, &
97			ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB)
98			itaplw0 = 0
99			ENDIF
100			itaplw0 = itaplw0 + 1
101
102			IF (MOD(itaplw, lwpas) == 0) THEN
103			CALL LWC(ILIM, PCLDLD, PCLDLU, PEMIS, &
104			ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB, &
105			ZFLUX)
106			itaplw = 0
107			ENDIF
108			itaplw = itaplw + 1
109
110			DO k = 1, KFLEV
111			kpl1 = k+1
112			DO i = 1, KDLON
113			PCOLR(i, k) = ZFLUX(i, 1, kpl1)+ZFLUX(i, 2, kpl1) &
114			- ZFLUX(i, 1, k)- ZFLUX(i, 2, k)
115			PCOLR(i, k) = PCOLR(i, k) * RDAY*RG/RCPD / PDP(i, k)
116			PCOLR0(i, k) = ZFLUC(i, 1, kpl1)+ZFLUC(i, 2, kpl1) &
117			- ZFLUC(i, 1, k)- ZFLUC(i, 2, k)
118			PCOLR0(i, k) = PCOLR0(i, k) * RDAY*RG/RCPD / PDP(i, k)
119			ENDDO
120			ENDDO
121			DO i = 1, KDLON
122			PSOLLW(i) = -ZFLUX(i, 1, 1)-ZFLUX(i, 2, 1)
123			PTOPLW(i) = ZFLUX(i, 1, KFLEV+1) + ZFLUX(i, 2, KFLEV+1)
124
125			PSOLLW0(i) = -ZFLUC(i, 1, 1)-ZFLUC(i, 2, 1)
126			PTOPLW0(i) = ZFLUC(i, 1, KFLEV+1) + ZFLUC(i, 2, KFLEV+1)
127			psollwdown(i) = -ZFLUX(i, 2, 1)
128
129			!IM attention aux signes !; LWtop >0, LWdn < 0
130			DO k = 1, KFLEV+1
131			plwup(i, k) = ZFLUX(i, 1, k)
132			plwup0(i, k) = ZFLUC(i, 1, k)
133			plwdn(i, k) = ZFLUX(i, 2, k)
134			plwdn0(i, k) = ZFLUC(i, 2, k)
135			ENDDO
136			ENDDO
137
138			END SUBROUTINE LW
139
140			end module lw_m