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	module lw_m
2
3	IMPLICIT none
4
5	contains
6
7	SUBROUTINE LW(PPMB, PDP, PDT0, PEMIS, PTL, PTAVE, PWV, POZON, PAER, PCLDLD, &
8	PCLDLU, PVIEW, PCOLR, PCOLR0, PTOPLW, PSOLLW, PTOPLW0, PSOLLW0, &
9	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	CALL LWU(PAER, PDP, PPMB, ZOZ, PTAVE, PVIEW, PWV, ZABCU)
96	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