SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
diag_sean.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_sea_n (DGL, DGS, 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_SEA_n * - Chooses the surface schemes for sea 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 01/2006 : sea flux parameterization.
37 !! Modified 08/2009 : new diag
38 ! B. decharme 04/2013 : Add EVAP and SUBL diag
39 !!------------------------------------------------------------------
40 !
41 
42 !
43 !
44 !
45 !
46 USE modd_diag_ideal_n, ONLY : diag_ideal_t
47 USE modd_diag_seaflux_n, ONLY : diag_seaflux_t
48 USE modd_surf_atm_n, ONLY : surf_atm_t
49 !
50 USE modd_surf_par, ONLY : xundef
51 !
52 USE modi_diag_seaflux_n
53 USE modi_diag_ideal_n
54 !
55 USE yomhook ,ONLY : lhook, dr_hook
56 USE parkind1 ,ONLY : jprb
57 !
58 IMPLICIT NONE
59 !
60 !* 0.1 declarations of arguments
61 !
62 !
63 TYPE(diag_ideal_t), INTENT(INOUT) :: dgl
64 TYPE(diag_seaflux_t), INTENT(INOUT) :: dgs
65 TYPE(surf_atm_t), INTENT(INOUT) :: u
66 !
67  CHARACTER(LEN=6), INTENT(IN) :: hprogram ! program calling surf. schemes
68 !
69 REAL, DIMENSION(:), INTENT(OUT) :: prn ! Net radiation (W/m2)
70 REAL, DIMENSION(:), INTENT(OUT) :: ph ! Sensible heat flux (W/m2)
71 REAL, DIMENSION(:), INTENT(OUT) :: ple ! Total latent heat flux (W/m2)
72 REAL, DIMENSION(:), INTENT(OUT) :: plei ! Sublimation latent heat flux (W/m2)
73 REAL, DIMENSION(:), INTENT(OUT) :: pgflux ! Storage flux (W/m2)
74 REAL, DIMENSION(:), INTENT(OUT) :: pevap ! Total evapotranspiration (kg/m2/s)
75 REAL, DIMENSION(:), INTENT(OUT) :: psubl ! Sublimation (kg/m2/s)
76 REAL, DIMENSION(:), INTENT(OUT) :: pri ! Richardson number (-)
77 REAL, DIMENSION(:), INTENT(OUT) :: pcd ! drag coefficient (W/s2)
78 REAL, DIMENSION(:), INTENT(OUT) :: pch ! transf. coef heat (W/s)
79 REAL, DIMENSION(:), INTENT(OUT) :: pce ! transf. coef vapor (W/s/K)
80 REAL, DIMENSION(:), INTENT(OUT) :: pqs
81 REAL, DIMENSION(:), INTENT(OUT) :: pz0 ! rough. length wind (m)
82 REAL, DIMENSION(:), INTENT(OUT) :: pz0h ! rough. length heat (m)
83 REAL, DIMENSION(:), INTENT(OUT) :: pts ! surface temperature (K)
84 REAL, DIMENSION(:), INTENT(OUT) :: pt2m ! temperature at 2m (K)
85 REAL, DIMENSION(:), INTENT(OUT) :: pq2m ! humidity at 2m (kg/kg)
86 REAL, DIMENSION(:), INTENT(OUT) :: phu2m ! relative humidity at 2m (-)
87 REAL, DIMENSION(:), INTENT(OUT) :: pzon10m ! zonal wind at 10m (m/s)
88 REAL, DIMENSION(:), INTENT(OUT) :: pmer10m ! meridian wind at 10m(m/s)
89 REAL, DIMENSION(:), INTENT(OUT) :: pswd ! incoming short wave radiation (W/m2)
90 REAL, DIMENSION(:), INTENT(OUT) :: pswu ! outgoing short wave radiation (W/m2)
91 REAL, DIMENSION(:,:), INTENT(OUT) :: pswbd ! incoming short wave radiation by spectral band(W/m2)
92 REAL, DIMENSION(:,:), INTENT(OUT) :: pswbu ! outgoing short wave radiation by spectral band(W/m2)
93 REAL, DIMENSION(:), INTENT(OUT) :: plwd ! incoming long wave radiation (W/m2)
94 REAL, DIMENSION(:), INTENT(OUT) :: plwu ! outgoing long wave radiation (W/m2)
95 REAL, DIMENSION(:), INTENT(OUT) :: pfmu ! zonal friction
96 REAL, DIMENSION(:), INTENT(OUT) :: pfmv ! meridian friction
97 REAL, DIMENSION(:), INTENT(OUT) :: prnc ! Net radiation (J/m2)
98 REAL, DIMENSION(:), INTENT(OUT) :: phc ! Sensible heat flux (J/m2)
99 REAL, DIMENSION(:), INTENT(OUT) :: plec ! Total latent heat flux (J/m2)
100 REAL, DIMENSION(:), INTENT(OUT) :: pleic ! Sublimation latent heat flux (J/m2)
101 REAL, DIMENSION(:), INTENT(OUT) :: pgfluxc ! Storage flux (J/m2)
102 REAL, DIMENSION(:), INTENT(OUT) :: pevapc ! Total evapotranspiration (kg/m2/s)
103 REAL, DIMENSION(:), INTENT(OUT) :: psublc ! Sublimation (kg/m2/s)
104 REAL, DIMENSION(:), INTENT(OUT) :: pswdc ! incoming short wave radiation (J/m2)
105 REAL, DIMENSION(:), INTENT(OUT) :: pswuc ! outgoing short wave radiation (J/m2)
106 REAL, DIMENSION(:), INTENT(OUT) :: plwdc ! incoming long wave radiation (J/m2)
107 REAL, DIMENSION(:), INTENT(OUT) :: plwuc ! outgoing long wave radiation (J/m2)
108 REAL, DIMENSION(:), INTENT(OUT) :: pfmuc ! zonal friction
109 REAL, DIMENSION(:), INTENT(OUT) :: pfmvc ! meridian friction
110 REAL, DIMENSION(:), INTENT(OUT) :: pt2m_min ! Minimum temperature at 2m (K)
111 REAL, DIMENSION(:), INTENT(OUT) :: pt2m_max ! Maximum temperature at 2m (K)
112 REAL, DIMENSION(:), INTENT(OUT) :: phu2m_min! Minimum relative humidity at 2m (-)
113 REAL, DIMENSION(:), INTENT(OUT) :: phu2m_max! Maximum relative humidity at 2m (-)
114 REAL, DIMENSION(:), INTENT(OUT) :: pwind10m ! wind at 10m (m/s)
115 REAL, DIMENSION(:), INTENT(OUT) :: pwind10m_max! Maximum wind at 10m (m/s)
116 REAL(KIND=JPRB) :: zhook_handle
117 !
118 !
119 !* 0.2 declarations of local variables
120 !
121 !-------------------------------------------------------------------------------------
122 !
123 IF (lhook) CALL dr_hook('DIAG_SEA_N',0,zhook_handle)
124 IF (u%CSEA=='SEAFLX') THEN
125  CALL diag_seaflux_n(dgs, &
126  hprogram, &
127  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs,&
128  pz0, pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m,&
129  pswd, pswu, plwd, plwu, pswbd, pswbu, pfmu, pfmv, &
130  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
131  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
132  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
133  pevap, pevapc, psubl, psublc )
134 ELSEIF (u%CSEA=='FLUX') THEN
135  CALL diag_ideal_n(dgl, hprogram, &
136  prn, ph, ple, plei, pgflux, pri, pcd, pch, pce, pqs,&
137  pz0, pz0h, pt2m, pts, pq2m, phu2m, pzon10m, pmer10m,&
138  pswd, pswu, plwd, plwu, pswbd, pswbu, pfmu, pfmv, &
139  prnc, phc, plec, pgfluxc, pswdc, pswuc, plwdc, &
140  plwuc, pfmuc, pfmvc, pt2m_min, pt2m_max, pleic, &
141  phu2m_min, phu2m_max, pwind10m, pwind10m_max, &
142  pevap, pevapc, psubl, psublc )
143 ELSE IF (u%CSEA=='NONE ') THEN
144  prn = xundef
145  ph = xundef
146  ple = xundef
147  plei = xundef
148  pevap = xundef
149  psubl = xundef
150  pgflux = xundef
151  pri = xundef
152  pcd = xundef
153  pch = xundef
154  pce = xundef
155  pqs = xundef
156  pz0 = xundef
157  pz0h = xundef
158  pts = xundef
159  pt2m = xundef
160  pq2m = xundef
161  phu2m = xundef
162  pzon10m = xundef
163  pmer10m = xundef
164  pswd = xundef
165  pswu = xundef
166  pswbd = xundef
167  pswbu = xundef
168  plwd = xundef
169  plwu = xundef
170  pfmu = xundef
171  pfmv = xundef
172  prnc = xundef
173  phc = xundef
174  plec = xundef
175  pleic = xundef
176  pevapc = xundef
177  psublc = xundef
178  pgfluxc = xundef
179  pswdc = xundef
180  pswuc = xundef
181  plwdc = xundef
182  plwuc = xundef
183  pfmuc = xundef
184  pfmvc = xundef
185  pt2m_min = xundef
186  pt2m_max = xundef
187  phu2m_min= xundef
188  phu2m_max= xundef
189  pwind10m = xundef
190  pwind10m_max = xundef
191 END IF
192 IF (lhook) CALL dr_hook('DIAG_SEA_N',1,zhook_handle)
193 !
194 !-------------------------------------------------------------------------------------
195 !
196 END SUBROUTINE diag_sea_n
modd_surf_atm_n::surf_atm_t
Definition: modd_surf_atmn.F90:42
diag_sea_n
subroutine diag_sea_n(DGL, DGS, 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_sean.F90:6
modd_diag_ideal_n::diag_ideal_t
Definition: modd_diag_idealn.F90:41
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
modd_diag_seaflux_n
Definition: modd_diag_seafluxn.F90:6
modd_diag_seaflux_n::diag_seaflux_t
Definition: modd_diag_seafluxn.F90:42
modd_surf_par
Definition: modd_surf_par.F90:6
modd_diag_ideal_n
Definition: modd_diag_idealn.F90:6
diag_seaflux_n
subroutine diag_seaflux_n(DGS, 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_seafluxn.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