SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
diag_inland_watern.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 ! #########
6 SUBROUTINE diag_inland_water_n (DGF, DGL, DGW, U, &
7  hprogram, &
8  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs, &
9  pz0,pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m, &
10  pswd, pswu, pswbd, pswbu, plwd, plwu, pfmu, pfmv, &
11  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
12  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
13  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
14  pevap, pevapc, psubl, psublc )
15 ! ###############################################################################
16 !
17 !!**** *DIAG_INLAND_WATER_n * - Chooses the surface schemes for lakes diagnostics
18 !!
19 !! PURPOSE
20 !! -------
21 !
22 !!** METHOD
23 !! ------
24 !!
25 !! REFERENCE
26 !! ---------
27 !!
28 !!
29 !! AUTHOR
30 !! ------
31 !! V. Masson
32 !!
33 !! MODIFICATIONS
34 !! -------------
35 !! Original 01/2004
36 !! Modified 08/2009 : cumulated diag & t2m min/max
37 !! V.Masson 10/2013 Adds min and max 2m parameters
38 ! B. decharme 04/2013 : Add EVAP and SUBL diag
39 !!------------------------------------------------------------------
40 !
41 
42 !
43 !
44 !
45 !
46 !
47 USE modd_diag_flake_n, ONLY : diag_flake_t
48 USE modd_diag_ideal_n, ONLY : diag_ideal_t
49 USE modd_diag_watflux_n, ONLY : diag_watflux_t
50 USE modd_surf_atm_n, ONLY : surf_atm_t
51 !
52 USE modd_surf_par, ONLY : xundef
53 
54 USE modi_diag_watflux_n
55 USE modi_diag_flake_n
56 USE modi_diag_ideal_n
57 !
58 !
59 USE yomhook ,ONLY : lhook, dr_hook
60 USE parkind1 ,ONLY : jprb
61 !
62 IMPLICIT NONE
63 !
64 !* 0.1 declarations of arguments
65 !
66 !
67 TYPE(diag_flake_t), INTENT(INOUT) :: dgf
68 TYPE(diag_ideal_t), INTENT(INOUT) :: dgl
69 TYPE(diag_watflux_t), INTENT(INOUT) :: dgw
70 TYPE(surf_atm_t), INTENT(INOUT) :: u
71 !
72  CHARACTER(LEN=6), INTENT(IN) :: hprogram ! program calling surf. schemes
73 !
74 REAL, DIMENSION(:), INTENT(OUT) :: prn ! Net radiation (W/m2)
75 REAL, DIMENSION(:), INTENT(OUT) :: ph ! Sensible heat flux (W/m2)
76 REAL, DIMENSION(:), INTENT(OUT) :: ple ! Total latent heat flux (W/m2)
77 REAL, DIMENSION(:), INTENT(OUT) :: plei ! Sublimation latent heat flux (W/m2)
78 REAL, DIMENSION(:), INTENT(OUT) :: pgflux ! Storage flux (W/m2)
79 REAL, DIMENSION(:), INTENT(OUT) :: pevap ! Total evapotranspiration (kg/m2/s)
80 REAL, DIMENSION(:), INTENT(OUT) :: psubl ! Sublimation (kg/m2/s)
81 REAL, DIMENSION(:), INTENT(OUT) :: pri ! Richardson number (-)
82 REAL, DIMENSION(:), INTENT(OUT) :: pcd ! drag coefficient (W/s2)
83 REAL, DIMENSION(:), INTENT(OUT) :: pch ! transf. coef heat (W/s)
84 REAL, DIMENSION(:), INTENT(OUT) :: pce ! transf. coef vapor (W/s/K)
85 REAL, DIMENSION(:), INTENT(OUT) :: pqs
86 REAL, DIMENSION(:), INTENT(OUT) :: pz0 ! rough. length wind (m)
87 REAL, DIMENSION(:), INTENT(OUT) :: pz0h ! rough. length heat (m)
88 REAL, DIMENSION(:), INTENT(OUT) :: pts ! surface temperature (K)
89 REAL, DIMENSION(:), INTENT(OUT) :: pt2m ! temperature at 2m (K)
90 REAL, DIMENSION(:), INTENT(OUT) :: pq2m ! humidity at 2m (kg/kg)
91 REAL, DIMENSION(:), INTENT(OUT) :: phu2m ! relative humidity at 2m (-)
92 REAL, DIMENSION(:), INTENT(OUT) :: pzon10m ! zonal wind at 10m (m/s)
93 REAL, DIMENSION(:), INTENT(OUT) :: pmer10m ! meridian wind at 10m (m/s)
94 REAL, DIMENSION(:), INTENT(OUT) :: pswd ! incoming short wave radiation (W/m2)
95 REAL, DIMENSION(:), INTENT(OUT) :: pswu ! outgoing short wave radiation (W/m2)
96 REAL, DIMENSION(:,:), INTENT(OUT) :: pswbd ! incoming short wave radiation by spectral band(W/m2)
97 REAL, DIMENSION(:,:), INTENT(OUT) :: pswbu ! outgoing short wave radiation by spectral band(W/m2)
98 REAL, DIMENSION(:), INTENT(OUT) :: plwd ! incoming long wave radiation (W/m2)
99 REAL, DIMENSION(:), INTENT(OUT) :: plwu ! outgoing long wave radiation (W/m2)
100 REAL, DIMENSION(:), INTENT(OUT) :: pfmu ! zonal friction
101 REAL, DIMENSION(:), INTENT(OUT) :: pfmv ! meridian friction
102 REAL, DIMENSION(:), INTENT(OUT) :: prnc ! Net radiation (J/m2)
103 REAL, DIMENSION(:), INTENT(OUT) :: phc ! Sensible heat flux (J/m2)
104 REAL, DIMENSION(:), INTENT(OUT) :: plec ! Total latent heat flux (J/m2)
105 REAL, DIMENSION(:), INTENT(OUT) :: pleic ! Sublimation latent heat flux (J/m2)
106 REAL, DIMENSION(:), INTENT(OUT) :: pgfluxc ! Storage flux (J/m2)
107 REAL, DIMENSION(:), INTENT(OUT) :: pevapc ! Total evapotranspiration (kg/m2/s)
108 REAL, DIMENSION(:), INTENT(OUT) :: psublc ! Sublimation (kg/m2/s)
109 REAL, DIMENSION(:), INTENT(OUT) :: pswdc ! incoming short wave radiation (J/m2)
110 REAL, DIMENSION(:), INTENT(OUT) :: pswuc ! outgoing short wave radiation (J/m2)
111 REAL, DIMENSION(:), INTENT(OUT) :: plwdc ! incoming long wave radiation (J/m2)
112 REAL, DIMENSION(:), INTENT(OUT) :: plwuc ! outgoing long wave radiation (J/m2)
113 REAL, DIMENSION(:), INTENT(OUT) :: pfmuc ! zonal friction
114 REAL, DIMENSION(:), INTENT(OUT) :: pfmvc ! meridian friction
115 REAL, DIMENSION(:), INTENT(OUT) :: pt2m_min ! Minimum temperature at 2m (K)
116 REAL, DIMENSION(:), INTENT(OUT) :: pt2m_max ! Maximum temperature at 2m (K)
117 REAL, DIMENSION(:), INTENT(OUT) :: phu2m_min! Minimum relative humidity at 2m (-)
118 REAL, DIMENSION(:), INTENT(OUT) :: phu2m_max! Maximum relative humidity at 2m (-)
119 REAL, DIMENSION(:), INTENT(OUT) :: pwind10m ! wind at 10m (m/s)
120 REAL, DIMENSION(:), INTENT(OUT) :: pwind10m_max! Maximum wind at 10m (m/s)
121 REAL(KIND=JPRB) :: zhook_handle
122 !
123 !
124 !* 0.2 declarations of local variables
125 !
126 !-------------------------------------------------------------------------------------
127 !
128 IF (lhook) CALL dr_hook('DIAG_INLAND_WATER_N',0,zhook_handle)
129 IF (u%CWATER=='WATFLX') THEN
130  CALL diag_watflux_n(dgw, &
131  hprogram, &
132  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs,&
133  pz0, pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m,&
134  pswd, pswu, plwd, plwu, pswbd, pswbu, pfmu, pfmv, &
135  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
136  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
137  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
138  pevap, pevapc, psubl, psublc )
139 ELSE IF (u%CWATER=='FLAKE ') THEN
140  CALL diag_flake_n(dgf, &
141  hprogram, &
142  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs,&
143  pz0, pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m,&
144  pswd, pswu, plwd, plwu, pswbd, pswbu, pfmu, pfmv, &
145  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
146  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
147  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
148  pevap, pevapc, psubl, psublc )
149 !
150 ELSE IF (u%CWATER=='FLUX ') THEN
151  CALL diag_ideal_n(dgl, hprogram, &
152  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs,&
153  pz0, pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m,&
154  pswd, pswu, plwd, plwu, pswbd, pswbu, pfmu, pfmv, &
155  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
156  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
157  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
158  pevap, pevapc, psubl, psublc )
159 ELSE IF (u%CWATER=='NONE ') THEN
160  prn = xundef
161  ph = xundef
162  ple = xundef
163  plei = xundef
164  pevap = xundef
165  psubl = xundef
166  pgflux = xundef
167  pri = xundef
168  pcd = xundef
169  pch = xundef
170  pce = xundef
171  pqs = xundef
172  pz0 = xundef
173  pz0h = xundef
174  pts = xundef
175  pt2m = xundef
176  pq2m = xundef
177  phu2m = xundef
178  pzon10m = xundef
179  pmer10m = xundef
180  pswd = xundef
181  pswu = xundef
182  pswbd = xundef
183  pswbu = xundef
184  plwd = xundef
185  plwu = xundef
186  pfmu = xundef
187  pfmv = xundef
188  prnc = xundef
189  phc = xundef
190  plec = xundef
191  pleic = xundef
192  pevapc = xundef
193  psublc = xundef
194  pgfluxc = xundef
195  pswdc = xundef
196  pswuc = xundef
197  plwdc = xundef
198  plwuc = xundef
199  pfmuc = xundef
200  pfmvc = xundef
201  pt2m_min = xundef
202  pt2m_max = xundef
203  phu2m_min= xundef
204  phu2m_max= xundef
205  pwind10m = xundef
206  pwind10m_max = xundef
207 END IF
208 IF (lhook) CALL dr_hook('DIAG_INLAND_WATER_N',1,zhook_handle)
209 !
210 !-------------------------------------------------------------------------------------
211 !
212 END SUBROUTINE diag_inland_water_n
modd_surf_atm_n::surf_atm_t
Definition: modd_surf_atmn.F90:42
diag_inland_water_n
subroutine diag_inland_water_n(DGF, DGL, DGW, U, HPROGRAM, PRN, PH, PLE, PLEI, PGFLUX, PRI, PCD, PCH, PCE, PQS, PZ0, PZ0H, PT2M, PTS, PQ2M, PHU2M, PZON10M, PMER10M, PSWD, PSWU, PSWBD, PSWBU, PLWD, PLWU, PFMU, PFMV, PRNC, PHC, PLEC, PGFLUXC, PSWDC, PSWUC, PLWDC, PLWUC, PFMUC, PFMVC, PT2M_MIN, PT2M_MAX, PLEIC, PHU2M_MIN, PHU2M_MAX, PWIND10M, PWIND10M_MAX, PEVAP, PEVAPC, PSUBL, PSUBLC)
Definition: diag_inland_watern.F90:6
modd_diag_watflux_n::diag_watflux_t
Definition: modd_diag_watfluxn.F90:41
diag_watflux_n
subroutine diag_watflux_n(DGW, HPROGRAM, PRN, PH, PLE, PLEI, PGFLUX, PRI, PCD, PCH, PCE, PQS, PZ0, PZ0H, PT2M, PTS, PQ2M, PHU2M, PZON10M, PMER10M, PSWD, PSWU, PLWD, PLWU, PSWBD, PSWBU, PFMU, PFMV, PRNC, PHC, PLEC, PGFLUXC, PSWDC, PSWUC, PLWDC, PLWUC, PFMUC, PFMVC, PT2M_MIN, PT2M_MAX, PLEIC, PHU2M_MIN, PHU2M_MAX, PWIND10M, PWIND10M_MAX, PEVAP, PEVAPC, PSUBL, PSUBLC)
Definition: diag_watfluxn.F90:6
modd_diag_ideal_n::diag_ideal_t
Definition: modd_diag_idealn.F90:41
modd_diag_flake_n::diag_flake_t
Definition: modd_diag_flaken.F90:42
diag_ideal_n
subroutine diag_ideal_n(DGL, HPROGRAM, PRN, PH, PLE, PLEI, PGFLUX, PRI, PCD, PCH, PCE, PQS, PZ0, PZ0H, PT2M, PTS, PQ2M, PHU2M, PZON10M, PMER10M, PSWD, PSWU, PLWD, PLWU, PSWBD, PSWBU, PFMU, PFMV, PRNC, PHC, PLEC, PGFLUXC, PSWDC, PSWUC, PLWDC, PLWUC, PFMUC, PFMVC, PT2M_MIN, PT2M_MAX, PLEIC, PHU2M_MIN, PHU2M_MAX, PWIND10M, PWIND10M_MAX, PEVAP, PEVAPC, PSUBL, PSUBLC)
Definition: diag_idealn.F90:6
diag_flake_n
subroutine diag_flake_n(DGF, HPROGRAM, PRN, PH, PLE, PLEI, PGFLUX, PRI, PCD, PCH, PCE, PQS, PZ0, PZ0H, PT2M, PTS, PQ2M, PHU2M, PZON10M, PMER10M, PSWD, PSWU, PLWD, PLWU, PSWBD, PSWBU, PFMU, PFMV, PRNC, PHC, PLEC, PGFLUXC, PSWDC, PSWUC, PLWDC, PLWUC, PFMUC, PFMVC, PT2M_MIN, PT2M_MAX, PLEIC, PHU2M_MIN, PHU2M_MAX, PWIND10M, PWIND10M_MAX, PEVAP, PEVAPC, PSUBL, PSUBLC)
Definition: diag_flaken.F90:6
modd_diag_watflux_n
Definition: modd_diag_watfluxn.F90:6
modd_diag_flake_n
Definition: modd_diag_flaken.F90:6
modd_surf_par
Definition: modd_surf_par.F90:6
modd_diag_ideal_n
Definition: modd_diag_idealn.F90:6
modd_surf_atm_n
Definition: modd_surf_atmn.F90:6
Generated on Fri May 13 2016 09:31:03 for SURFEX by   doxygen 1.8.5