module_surface_driver.F
References to this file elsewhere.
1 !WRF:MEDIATION_LAYER:PHYSICS
2 !
3 MODULE module_surface_driver
4 CONTAINS
5
6 SUBROUTINE surface_driver( &
7 & acsnom,acsnow,akhs,akms,albedo,br,canwat &
8 & ,chklowq,dt,dx,dz8w,dzs,glw &
9 & ,grdflx,gsw,swdown,gz1oz0,hfx,ht,ifsnow,isfflx &
10 & ,isltyp,itimestep,ivgtyp,lowlyr,mavail,rmol &
11 & ,num_soil_layers,p8w,pblh,pi_phy,pshltr,psih &
12 & ,psim,p_phy,q10,q2,qfx,qsfc,qshltr,qz0 &
13 & ,raincv,rho,sfcevp,sfcexc,sfcrunoff &
14 & ,smois,smstav,smstot,snoalb,snow,snowc,snowh,stepbl &
15 & ,th10,th2,thz0,th_phy,tmn,tshltr,tsk,tslb &
16 & ,t_phy,u10,udrunoff,ust,uz0,u_frame,u_phy,v10,vegfra &
17 & ,vz0,v_frame,v_phy,warm_rain,wspd,xice,xland,z,znt,zs &
18 & ,ct,tke_myj &
19 & ,albbck,lh,sh2o,shdmax,shdmin,z0 &
20 & ,flqc,flhc,psfc,sst,sst_update,t2,emiss &
21 & ,sf_sfclay_physics,sf_surface_physics,ra_lw_physics &
22 ! Optional urban
23 & ,declin_urb,cosz_urb2d,omg_urb2d,xlat_urb2d & !I urban
24 & ,num_roof_layers, num_wall_layers & !I urban
25 & ,num_road_layers, dzr, dzb, dzg & !I urban
26 & ,tr_urb2d,tb_urb2d,tg_urb2d,tc_urb2d,qc_urb2d & !H urban
27 & ,uc_urb2d & !H urban
28 & ,xxxr_urb2d,xxxb_urb2d,xxxg_urb2d,xxxc_urb2d & !H urban
29 & ,trl_urb3d,tbl_urb3d,tgl_urb3d & !H urban
30 & ,sh_urb2d,lh_urb2d,g_urb2d,rn_urb2d,ts_urb2d & !H urban
31 & ,frc_urb2d, utype_urb2d & !H urban
32 & ,ucmcall & ! urban
33 & , ids,ide,jds,jde,kds,kde &
34 & , ims,ime,jms,jme,kms,kme &
35 & , i_start,i_end,j_start,j_end,kts,kte,num_tiles &
36 ! Optional moisture tracers
37 & ,qv_curr, qc_curr, qr_curr &
38 & ,qi_curr, qs_curr, qg_curr &
39 ! Optional moisture tracer flags
40 & ,f_qv,f_qc,f_qr &
41 & ,f_qi,f_qs,f_qg &
42 ! Other optionals (more or less em specific)
43 & ,capg,hol,mol &
44 & ,rainncv,rainbl,regime,thc &
45 & ,qsg,qvg,qcg,soilt1,tsnav &
46 & ,smfr3d,keepfr3dflag &
47 ! Other optionals (more or less nmm specific)
48 & ,potevp,snopcx,soiltb,sr &
49 ! Optional observation nudging
50 & ,uratx,vratx,tratx &
51 )
52
53 #if ( ! NMM_CORE == 1 )
54 USE module_state_description, ONLY : SFCLAYSCHEME &
55 ,MYJSFCSCHEME &
56 ,GFSSFCSCHEME &
57 ,SLABSCHEME &
58 ,LSMSCHEME &
59 ,RUCLSMSCHEME
60 #else
61 USE module_state_description, ONLY : SFCLAYSCHEME &
62 ,MYJSFCSCHEME &
63 ,GFSSFCSCHEME &
64 ,SLABSCHEME &
65 ,NMMLSMSCHEME &
66 ,LSMSCHEME &
67 ,RUCLSMSCHEME
68 #endif
69 USE module_model_constants
70 ! *** add new modules of schemes here
71
72 USE module_sf_sfclay
73 USE module_sf_myjsfc
74 USE module_sf_gfs
75 USE module_sf_noahlsm
76 USE module_sf_ruclsm
77 #if ( NMM_CORE == 1 )
78 USE module_sf_lsm_nmm
79 #endif
80
81 USE module_sf_slab
82 !
83 USE module_sf_sfcdiags
84 !
85
86 ! This driver calls subroutines for the surface parameterizations.
87 !
88 ! surface layer: (between surface and pbl)
89 ! 1. sfclay
90 ! 2. myjsfc
91 ! surface: ground temp/lsm scheme:
92 ! 1. slab
93 ! 2. Noah LSM
94 ! 99. NMM LSM (NMM core only)
95 !------------------------------------------------------------------
96 IMPLICIT NONE
97 !======================================================================
98 ! Grid structure in physics part of WRF
99 !----------------------------------------------------------------------
100 ! The horizontal velocities used in the physics are unstaggered
101 ! relative to temperature/moisture variables. All predicted
102 ! variables are carried at half levels except w, which is at full
103 ! levels. Some arrays with names (*8w) are at w (full) levels.
104 !
105 !----------------------------------------------------------------------
106 ! In WRF, kms (smallest number) is the bottom level and kme (largest
107 ! number) is the top level. In your scheme, if 1 is at the top level,
108 ! then you have to reverse the order in the k direction.
109 !
110 ! kme - half level (no data at this level)
111 ! kme ----- full level
112 ! kme-1 - half level
113 ! kme-1 ----- full level
114 ! .
115 ! kms+2 - half level
116 ! kms+2 ----- full level
117 ! kms+1 - half level
118 ! kms+1 ----- full level
119 ! kms - half level
120 ! kms ----- full level
121 !
122 !======================================================================
123 ! Definitions
124 !-----------
125 ! Theta potential temperature (K)
126 ! Qv water vapor mixing ratio (kg/kg)
127 ! Qc cloud water mixing ratio (kg/kg)
128 ! Qr rain water mixing ratio (kg/kg)
129 ! Qi cloud ice mixing ratio (kg/kg)
130 ! Qs snow mixing ratio (kg/kg)
131 !-----------------------------------------------------------------
132 !-- itimestep number of time steps
133 !-- GLW downward long wave flux at ground surface (W/m^2)
134 !-- GSW net short wave flux at ground surface (W/m^2)
135 !-- SWDOWN downward short wave flux at ground surface (W/m^2)
136 !-- EMISS surface emissivity (between 0 and 1)
137 !-- TSK surface temperature (K)
138 !-- TMN soil temperature at lower boundary (K)
139 !-- XLAND land mask (1 for land, 2 for water)
140 !-- ZNT time-varying roughness length (m)
141 !-- Z0 background roughness length (m)
142 !-- MAVAIL surface moisture availability (between 0 and 1)
143 !-- UST u* in similarity theory (m/s)
144 !-- MOL T* (similarity theory) (K)
145 !-- HOL PBL height over Monin-Obukhov length
146 !-- PBLH PBL height (m)
147 !-- CAPG heat capacity for soil (J/K/m^3)
148 !-- THC thermal inertia (Cal/cm/K/s^0.5)
149 !-- SNOWC flag indicating snow coverage (1 for snow cover)
150 !-- HFX net upward heat flux at the surface (W/m^2)
151 !-- QFX net upward moisture flux at the surface (kg/m^2/s)
152 !-- LH net upward latent heat flux at surface (W/m^2)
153 !-- REGIME flag indicating PBL regime (stable, unstable, etc.)
154 !-- tke_myj turbulence kinetic energy from Mellor-Yamada-Janjic (MYJ) (m^2/s^2)
155 !-- akhs sfc exchange coefficient of heat/moisture from MYJ
156 !-- akms sfc exchange coefficient of momentum from MYJ
157 !-- thz0 potential temperature at roughness length (K)
158 !-- uz0 u wind component at roughness length (m/s)
159 !-- vz0 v wind component at roughness length (m/s)
160 !-- qsfc specific humidity at lower boundary (kg/kg)
161 !-- uratx ratio of u over u10 (Added for obs-nudging)
162 !-- vratx ratio of v over v10 (Added for obs-nudging)
163 !-- tratx ratio of t over th2 (Added for obs-nudging)
164 !-- u10 diagnostic 10-m u component from surface layer
165 !-- v10 diagnostic 10-m v component from surface layer
166 !-- th2 diagnostic 2-m theta from surface layer and lsm
167 !-- t2 diagnostic 2-m temperature from surface layer and lsm
168 !-- q2 diagnostic 2-m mixing ratio from surface layer and lsm
169 !-- tshltr diagnostic 2-m theta from MYJ
170 !-- th10 diagnostic 10-m theta from MYJ
171 !-- qshltr diagnostic 2-m specific humidity from MYJ
172 !-- q10 diagnostic 10-m specific humidity from MYJ
173 !-- lowlyr index of lowest model layer above ground
174 !-- rr dry air density (kg/m^3)
175 !-- u_phy u-velocity interpolated to theta points (m/s)
176 !-- v_phy v-velocity interpolated to theta points (m/s)
177 !-- th_phy potential temperature (K)
178 !-- moist moisture array (4D - last index is species) (kg/kg)
179 !-- p_phy pressure (Pa)
180 !-- pi_phy exner function (dimensionless)
181 !-- pshltr diagnostic shelter (2m) pressure from MYJ (Pa)
182 !-- p8w pressure at full levels (Pa)
183 !-- t_phy temperature (K)
184 !-- dz8w dz between full levels (m)
185 !-- z height above sea level (m)
186 !-- DX horizontal space interval (m)
187 !-- DT time step (second)
188 !-- PSFC pressure at the surface (Pa)
189 !-- SST sea-surface temperature (K)
190 !-- TSLB
191 !-- ZS
192 !-- DZS
193 !-- num_soil_layers number of soil layer
194 !-- IFSNOW ifsnow=1 for snow-cover effects
195 !
196 !-- ids start index for i in domain
197 !-- ide end index for i in domain
198 !-- jds start index for j in domain
199 !-- jde end index for j in domain
200 !-- kds start index for k in domain
201 !-- kde end index for k in domain
202 !-- ims start index for i in memory
203 !-- ime end index for i in memory
204 !-- jms start index for j in memory
205 !-- jme end index for j in memory
206 !-- kms start index for k in memory
207 !-- kme end index for k in memory
208 !-- its start index for i in tile
209 !-- ite end index for i in tile
210 !-- jts start index for j in tile
211 !-- jte end index for j in tile
212 !-- kts start index for k in tile
213 !-- kte end index for k in tile
214 !
215 !******************************************************************
216 !------------------------------------------------------------------
217
218 INTEGER, INTENT(IN) :: &
219 & ids,ide,jds,jde,kds,kde &
220 & ,ims,ime,jms,jme,kms,kme &
221 & ,kts,kte,num_tiles
222
223 INTEGER, INTENT(IN) :: sf_sfclay_physics,sf_surface_physics,ra_lw_physics,sst_update
224
225 INTEGER, INTENT(IN) :: ucmcall !urban
226
227 INTEGER, DIMENSION(num_tiles), INTENT(IN) :: &
228 & i_start,i_end,j_start,j_end
229
230 INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: ISLTYP
231 INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: IVGTYP
232 INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: LOWLYR
233 INTEGER, INTENT(IN ):: IFSNOW
234 INTEGER, INTENT(IN ):: ISFFLX
235 INTEGER, INTENT(IN ):: ITIMESTEP
236 INTEGER, INTENT(IN ):: NUM_SOIL_LAYERS
237 INTEGER, INTENT(IN ):: STEPBL
238 LOGICAL, INTENT(IN ):: WARM_RAIN
239 REAL , INTENT(IN ):: U_FRAME
240 REAL , INTENT(IN ):: V_FRAME
241 REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT):: SMOIS
242 REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT):: TSLB
243 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: GLW
244 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: GSW,SWDOWN
245 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: HT
246 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: RAINCV
247 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: SST
248 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: TMN
249 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: VEGFRA
250 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: XICE
251 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: XLAND
252 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: MAVAIL
253 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: SNOALB
254 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: ACSNOW
255 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: AKHS
256 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: AKMS
257 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: ALBEDO
258 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: CANWAT
259
260
261 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: GRDFLX
262 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: HFX
263 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: RMOL
264 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: PBLH
265 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: Q2
266 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: QFX
267 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: QSFC
268 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: QZ0
269 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SFCRUNOFF
270 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SMSTAV
271 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SMSTOT
272 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SNOW
273 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SNOWC
274 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: SNOWH
275 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: TH2
276 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: THZ0
277 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: TSK
278 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: UDRUNOFF
279 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: UST
280 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: UZ0
281 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: VZ0
282 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: WSPD
283 REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT):: ZNT
284 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: BR
285 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: CHKLOWQ
286 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: GZ1OZ0
287 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: PSHLTR
288 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: PSIH
289 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: PSIM
290 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: Q10
291 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: QSHLTR
292 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: TH10
293 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: TSHLTR
294 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: U10
295 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: V10
296 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT):: PSFC
297 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: ACSNOM
298 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: SFCEVP
299 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: SFCEXC
300 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: FLHC
301 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: FLQC
302 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CT
303 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: DZ8W
304 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: P8W
305 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: PI_PHY
306 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: P_PHY
307 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: RHO
308 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: TH_PHY
309 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: T_PHY
310 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: U_PHY
311 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: V_PHY
312 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ):: Z
313 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: TKE_MYJ
314 REAL, DIMENSION(1:num_soil_layers), INTENT(IN):: DZS
315 REAL, DIMENSION(1:num_soil_layers), INTENT(IN):: ZS
316 REAL, INTENT(IN ):: DT
317 REAL, INTENT(IN ):: DX
318
319 ! arguments for NCAR surface physics
320
321 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ):: ALBBCK ! INOUT needed for NMM
322 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ):: LH
323 REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT):: SH2O
324 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: SHDMAX
325 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: SHDMIN
326 REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN ):: Z0
327
328 !
329 ! Optional
330 !
331
332 !
333 ! Observation nudging
334 !
335 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT):: uratx !Added for obs-nudging
336 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT):: vratx !Added for obs-nudging
337 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT):: tratx !Added for obs-nudging
338 !
339 ! Flags relating to the optional tendency arrays declared above
340 ! Models that carry the optional tendencies will provdide the
341 ! optional arguments at compile time; these flags all the model
342 ! to determine at run-time whether a particular tracer is in
343 ! use or not.
344 !
345 LOGICAL, INTENT(IN), OPTIONAL :: &
346 f_qv &
347 ,f_qc &
348 ,f_qr &
349 ,f_qi &
350 ,f_qs &
351 ,f_qg
352
353 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
354 OPTIONAL, INTENT(INOUT) :: &
355 ! optional moisture tracers
356 ! 2 time levels; if only one then use CURR
357 qv_curr, qc_curr, qr_curr &
358 ,qi_curr, qs_curr, qg_curr
359
360 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: capg
361 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: emiss
362 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: hol
363 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: mol
364 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: regime
365 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN ):: rainncv
366 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: RAINBL
367 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: t2
368 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN ):: thc
369 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: qsg
370 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: qvg
371 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: qcg
372 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: soilt1
373 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: tsnav
374 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: potevp ! NMM LSM
375 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: snopcx ! NMM LSM
376 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: soiltb ! NMM LSM
377 REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: sr ! NMM and RUC LSM
378 REAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), OPTIONAL, INTENT(INOUT):: smfr3d
379 REAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), OPTIONAL, INTENT(INOUT):: keepfr3dflag
380
381 ! LOCAL VAR
382
383 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ) ::v_phytmp
384 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ) ::u_phytmp
385
386 REAL, DIMENSION( ims:ime, jms:jme ) :: ZOL
387
388 REAL, DIMENSION( ims:ime, jms:jme ) :: &
389 QGH, &
390 CHS, &
391 CPM, &
392 CHS2, &
393 CQS2
394
395 REAL :: DTMIN,DTBL
396 !
397 INTEGER :: i,J,K,NK,jj,ij
398 LOGICAL :: radiation, myj, frpcpn
399 !-------------------------------------------------
400 ! urban related variables are added to declaration
401 !-------------------------------------------------
402 REAL, OPTIONAL, INTENT(IN) :: DECLIN_URB !urban
403 REAL, OPTIONAL , DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: COSZ_URB2D !urban
404 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: OMG_URB2D !urban
405 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: XLAT_URB2D !urban
406 INTEGER, OPTIONAL, INTENT(IN) :: num_roof_layers !urban
407 INTEGER, OPTIONAL, INTENT(IN) :: num_wall_layers !urban
408 INTEGER, OPTIONAL, INTENT(IN) :: num_road_layers !urban
409 REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(IN) :: DZR !urban
410 REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(IN) :: DZB !urban
411 REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(IN) :: DZG !urban
412
413 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TR_URB2D !urban
414 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TB_URB2D !urban
415 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TG_URB2D !urban
416 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: TC_URB2D !urban
417 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: QC_URB2D !urban
418 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: UC_URB2D !urban
419 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: XXXR_URB2D !urban
420 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: XXXB_URB2D !urban
421 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: XXXG_URB2D !urban
422 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: XXXC_URB2D !urban
423 REAL, OPTIONAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), & !urban
424 INTENT(INOUT) :: TRL_URB3D !urban
425 REAL, OPTIONAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), & !urban
426 INTENT(INOUT) :: TBL_URB3D !urban
427 REAL, OPTIONAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), & !urban
428 INTENT(INOUT) :: TGL_URB3D !urban
429 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SH_URB2D !urban
430 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: LH_URB2D !urban
431 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: G_URB2D !urban
432 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: RN_URB2D !urban
433 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT):: TS_URB2D !urban
434 !
435 REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: FRC_URB2D !urban
436 INTEGER, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: UTYPE_URB2D !urban
437
438 REAL, DIMENSION( ims:ime, jms:jme ) :: PSIM_URB2D !urban local var
439 REAL, DIMENSION( ims:ime, jms:jme ) :: PSIH_URB2D !urban local var
440 REAL, DIMENSION( ims:ime, jms:jme ) :: GZ1OZ0_URB2D !urban local var
441 !m REAL, DIMENSION( ims:ime, jms:jme ) :: AKHS_URB2D !urban local var
442 REAL, DIMENSION( ims:ime, jms:jme ) :: AKMS_URB2D !urban local var
443 REAL, DIMENSION( ims:ime, jms:jme ) :: U10_URB2D !urban local var
444 REAL, DIMENSION( ims:ime, jms:jme ) :: V10_URB2D !urban local var
445 REAL, DIMENSION( ims:ime, jms:jme ) :: TH2_URB2D !urban local var
446 REAL, DIMENSION( ims:ime, jms:jme ) :: Q2_URB2D !urban local var
447 REAL, DIMENSION( ims:ime, jms:jme ) :: UST_URB2D !urban local var
448
449 !------------------------------------------------------------------
450 CHARACTER*256 :: message
451 !------------------------------------------------------------------
452 !
453
454 if (sf_sfclay_physics .eq. 0) return
455
456 v_phytmp = 0.
457 u_phytmp = 0.
458 ZOL = 0.
459 QGH = 0.
460 CHS = 0.
461 CPM = 0.
462 CHS2 = 0.
463 DTMIN = 0.
464 DTBL = 0.
465
466 ! RAINBL in mm (Accumulation between PBL calls)
467
468 IF ( PRESENT( rainncv ) .AND. PRESENT( rainbl ) ) THEN
469 !$OMP PARALLEL DO &
470 !$OMP PRIVATE ( ij, i, j, k )
471 DO ij = 1 , num_tiles
472 DO j=j_start(ij),j_end(ij)
473 DO i=i_start(ij),i_end(ij)
474 RAINBL(i,j) = RAINBL(i,j) + RAINCV(i,j) + RAINNCV(i,j)
475 RAINBL(i,j) = MAX (RAINBL(i,j), 0.0)
476 ENDDO
477 ENDDO
478 ENDDO
479 !$OMP END PARALLEL DO
480 ELSE IF ( PRESENT( rainbl ) ) THEN
481 !$OMP PARALLEL DO &
482 !$OMP PRIVATE ( ij, i, j, k )
483 DO ij = 1 , num_tiles
484 DO j=j_start(ij),j_end(ij)
485 DO i=i_start(ij),i_end(ij)
486 RAINBL(i,j) = RAINBL(i,j) + RAINCV(i,j)
487 RAINBL(i,j) = MAX (RAINBL(i,j), 0.0)
488 ENDDO
489 ENDDO
490 ENDDO
491 !$OMP END PARALLEL DO
492 ENDIF
493 ! Update SST
494 IF (sst_update .EQ. 1) THEN
495 !$OMP PARALLEL DO &
496 !$OMP PRIVATE ( ij, i, j, k )
497 DO ij = 1 , num_tiles
498 DO j=j_start(ij),j_end(ij)
499 DO i=i_start(ij),i_end(ij)
500 IF(XLAND(i,j) .GT. 1.5)TSK(i,j)=SST(i,j)
501 ENDDO
502 ENDDO
503 ENDDO
504 !$OMP END PARALLEL DO
505 ENDIF
506
507 IF (itimestep .eq. 1 .or. mod(itimestep,STEPBL) .eq. 0) THEN
508
509 radiation = .false.
510 myj = .false.
511 frpcpn = .false.
512
513 IF (ra_lw_physics .gt. 0) radiation = .true.
514
515 !----
516 ! CALCULATE CONSTANT
517
518 DTMIN=DT/60.
519 ! Surface schemes need PBL time step for updates and accumulations
520 ! Assume these schemes provide no tendencies
521 DTBL=DT*STEPBL
522
523 ! SAVE OLD VALUES
524
525
526 !$OMP PARALLEL DO &
527 !$OMP PRIVATE ( ij, i, j, k )
528 DO ij = 1 , num_tiles
529 DO j=j_start(ij),j_end(ij)
530 DO i=i_start(ij),i_end(ij)
531 ! PSFC : in Pa
532 PSFC(I,J)=p8w(I,kts,J)
533 ! REVERSE ORDER IN THE VERTICAL DIRECTION
534 DO k=kts,kte
535 v_phytmp(i,k,j)=v_phy(i,k,j)+v_frame
536 u_phytmp(i,k,j)=u_phy(i,k,j)+u_frame
537 ENDDO
538 ENDDO
539 ENDDO
540 ENDDO
541 !$OMP END PARALLEL DO
542
543 !$OMP PARALLEL DO &
544 !$OMP PRIVATE ( ij, i, j, k )
545 DO ij = 1 , num_tiles
546 sfclay_select: SELECT CASE(sf_sfclay_physics)
547
548 CASE (SFCLAYSCHEME)
549 #if (NMM_CORE != 1)
550 ! DX varies spatially in NMM, therefore, SFCLAY cannot be called
551 ! because it takes a scalar DX. NMM passes in a dummy value for this
552 ! scalar. NEEDS FURTHER ATTENTION. JM 20050215
553 IF (PRESENT(qv_curr) .AND. &
554 PRESENT(mol) .AND. PRESENT(regime) .AND. &
555 .TRUE. ) THEN
556 CALL wrf_debug( 100, 'in SFCLAY' )
557 CALL SFCLAY(u_phytmp,v_phytmp,t_phy,qv_curr,&
558 p_phy,dz8w,cp,g,rcp,r_d,xlv,psfc,chs,chs2,cqs2,cpm, &
559 znt,ust,pblh,mavail,zol,mol,regime,psim,psih, &
560 xland,hfx,qfx,lh,tsk,flhc,flqc,qgh,qsfc,rmol, &
561 uratx,vratx,tratx, &
562 u10,v10,th2,t2,q2, &
563 gz1oz0,wspd,br,isfflx,dx, &
564 svp1,svp2,svp3,svpt0,ep_1,ep_2,karman,eomeg,stbolt, &
565 ids,ide, jds,jde, kds,kde, &
566 ims,ime, jms,jme, kms,kme, &
567 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
568 ELSE
569 CALL wrf_error_fatal('Lacking arguments for SFCLAY in surface driver')
570 ENDIF
571
572 #else
573 CALL wrf_error_fatal('SFCLAY cannot be used with NMM')
574 #endif
575 CASE (MYJSFCSCHEME)
576 IF (PRESENT(qv_curr) .AND. PRESENT(qc_curr) .AND. &
577 .TRUE. ) THEN
578
579 myj =.true.
580
581 CALL wrf_debug(100,'in MYJSFC')
582 CALL MYJSFC(itimestep,ht,dz8w, &
583 p_phy,p8w,th_phy,t_phy, &
584 qv_curr,qc_curr, &
585 u_phy,v_phy,tke_myj, &
586 tsk,qsfc,thz0,qz0,uz0,vz0, &
587 lowlyr, &
588 xland, &
589 ust,znt,z0,pblh,mavail,rmol, &
590 akhs,akms, &
591 chs,chs2,cqs2,hfx,qfx,lh,flhc,flqc,qgh,cpm,ct, &
592 u10,v10,t2,th2,tshltr,th10,q2,qshltr,q10,pshltr, &
593 ids,ide, jds,jde, kds,kde, &
594 ims,ime, jms,jme, kms,kme, &
595 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
596 ELSE
597 CALL wrf_error_fatal('Lacking arguments for MYJSFC in surface driver')
598 ENDIF
599
600 CASE (GFSSFCSCHEME)
601 IF (PRESENT(qv_curr) .AND. .TRUE. ) THEN
602 CALL wrf_debug( 100, 'in GFSSFC' )
603 CALL SF_GFS(u_phytmp,v_phytmp,t_phy,qv_curr, &
604 p_phy,CP,RCP,R_d,XLV,PSFC,CHS,CHS2,CQS2,CPM, &
605 ZNT,UST,PSIM,PSIH, &
606 XLAND,HFX,QFX,LH,TSK,FLHC,FLQC, &
607 QGH,QSFC,U10,V10, &
608 GZ1OZ0,WSPD,BR,ISFFLX, &
609 EP_1,EP_2,KARMAN,itimestep, &
610 ids,ide, jds,jde, kds,kde, &
611 ims,ime, jms,jme, kms,kme, &
612 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
613 CALL wrf_debug(100,'in SFCDIAGS')
614 ELSE
615 CALL wrf_error_fatal('Lacking arguments for SF_GFS in surface driver')
616 ENDIF
617
618 CASE DEFAULT
619
620 WRITE( message , * ) &
621 'The sfclay option does not exist: sf_sfclay_physics = ', sf_sfclay_physics
622 CALL wrf_error_fatal ( message )
623
624 END SELECT sfclay_select
625 ENDDO
626 !$OMP END PARALLEL DO
627
628 IF (ISFFLX.EQ.0 ) GOTO 430
629 !$OMP PARALLEL DO &
630 !$OMP PRIVATE ( ij, i, j, k )
631 DO ij = 1 , num_tiles
632
633 sfc_select: SELECT CASE(sf_surface_physics)
634
635 CASE (SLABSCHEME)
636
637 IF (PRESENT(qv_curr) .AND. &
638 PRESENT(capg) .AND. &
639 .TRUE. ) THEN
640 DO j=j_start(ij),j_end(ij)
641 DO i=i_start(ij),i_end(ij)
642 ! CQS2 ACCOUNTS FOR MAVAIL FOR SFCDIAGS 2M Q
643 CQS2(I,J)= CQS2(I,J)*MAVAIL(I,J)
644 ENDDO
645 ENDDO
646
647 CALL wrf_debug(100,'in SLAB')
648 CALL SLAB(t_phy,qv_curr,p_phy,flhc,flqc, &
649 psfc,xland,tmn,hfx,qfx,lh,tsk,qsfc,chklowq, &
650 gsw,glw,capg,thc,snowc,emiss,mavail, &
651 dtbl,rcp,xlv,dtmin,ifsnow, &
652 svp1,svp2,svp3,svpt0,ep_2,karman,eomeg,stbolt, &
653 tslb,zs,dzs,num_soil_layers,radiation, &
654 ids,ide, jds,jde, kds,kde, &
655 ims,ime, jms,jme, kms,kme, &
656 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
657
658 DO j=j_start(ij),j_end(ij)
659 DO i=i_start(ij),i_end(ij)
660 SFCEVP(I,J)= SFCEVP(I,J) + QFX(I,J)*DTBL
661 ENDDO
662 ENDDO
663
664 CALL wrf_debug(100,'in SFCDIAGS')
665 CALL SFCDIAGS(hfx,qfx,tsk,qsfc,chs2,cqs2,t2,th2,q2, &
666 psfc,cp,r_d,rcp, &
667 ids,ide, jds,jde, kds,kde, &
668 ims,ime, jms,jme, kms,kme, &
669 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
670
671 ELSE
672 CALL wrf_error_fatal('Lacking arguments for SLAB in surface driver')
673 ENDIF
674
675 #if ( NMM_CORE == 1 )
676 CASE (NMMLSMSCHEME)
677 IF (PRESENT(qv_curr) .AND. PRESENT(rainbl) .AND. &
678 PRESENT(potevp) .AND. PRESENT(snopcx) .AND. &
679 PRESENT(soiltb) .AND. PRESENT(sr) .AND. &
680 .TRUE. ) THEN
681 CALL wrf_debug(100,'in NMM LSM')
682 CALL nmmlsm(dz8w,qv_curr,p8w,rho, &
683 t_phy,th_phy,tsk,chs, &
684 hfx,qfx,qgh,swdown,glw,lh,rmol, &
685 smstav,smstot,sfcrunoff, &
686 udrunoff,ivgtyp,isltyp,vegfra,sfcevp,potevp, &
687 grdflx,sfcexc,acsnow,acsnom,snopcx, &
688 albbck,tmn,xland,xice,qz0, &
689 th2,q2,snowc,cqs2,qsfc,soiltb,chklowq,rainbl, &
690 num_soil_layers,dtbl,dzs,itimestep, &
691 smois,tslb,snow,canwat,cpm,rcp,sr, & !tslb
692 albedo,snoalb,sh2o,snowh, &
693 ids,ide, jds,jde, kds,kde, &
694 ims,ime, jms,jme, kms,kme, &
695 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
696 CALL wrf_debug(100,'back from NMM LSM')
697 ELSE
698 CALL wrf_error_fatal('Lacking arguments for NMMLSM in surface driver')
699 ENDIF
700 #endif
701
702 CASE (LSMSCHEME)
703
704 IF (PRESENT(qv_curr) .AND. PRESENT(rainbl) .AND. &
705 ! PRESENT(emiss) .AND. PRESENT(t2) .AND. &
706 ! PRESENT(declin_urb) .AND. PRESENT(cosz_urb2d) .AND. &
707 ! PRESENT(omg_urb2d) .AND. PRESENT( xlat_urb2d) .AND. &
708 ! PRESENT(dzr) .AND. &
709 ! PRESENT( dzb) .AND. PRESENT(dzg) .AND. &
710 ! PRESENT(tr_urb2d) .AND. PRESENT(tb_urb2d) .AND. &
711 ! PRESENT(tg_urb2d) .AND. PRESENT(tc_urb2d) .AND. &
712 ! PRESENT(qc_urb2d) .AND. PRESENT(uc_urb2d) .AND. &
713 ! PRESENT(xxxr_urb2d) .AND. PRESENT(xxxb_urb2d) .AND. &
714 ! PRESENT(xxxg_urb2d) .AND. &
715 ! PRESENT(xxxc_urb2d) .AND. PRESENT(trl_urb3d) .AND. &
716 ! PRESENT(tbl_urb3d) .AND. PRESENT(tgl_urb3d) .AND. &
717 ! PRESENT(sh_urb2d) .AND. PRESENT(lh_urb2d) .AND. &
718 ! PRESENT(g_urb2d) .AND. PRESENT(rn_urb2d) .AND. &
719 ! PRESENT(ts_urb2d) .AND. &
720 ! PRESENT(frc_urb2d) .AND. PRESENT(utype_urb2d) .AND. &
721 .TRUE. ) THEN
722 !------------------------------------------------------------------
723 CALL wrf_debug(100,'in NOAH LSM')
724 CALL lsm(dz8w,qv_curr,p8w,t_phy,tsk, &
725 hfx,qfx,lh,grdflx,qgh,gsw,glw,smstav,smstot, &
726 sfcrunoff,udrunoff,ivgtyp,isltyp,vegfra, &
727 albedo,albbck,znt,z0, tmn,xland,xice, emiss, &
728 snowc,qsfc,rainbl, &
729 num_soil_layers,dtbl,dzs,itimestep, &
730 smois,tslb,snow,canwat, &
731 chs, chs2, cqs2, cpm,rcp, &
732 sh2o,snowh, & !h
733 u_phy,v_phy, & !I
734 snoalb,shdmin,shdmax, & !i
735 acsnom,acsnow, & !o
736 ids,ide, jds,jde, kds,kde, &
737 ims,ime, jms,jme, kms,kme, &
738 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte, &
739 ucmcall &
740 !Optional urban
741 ,tr_urb2d,tb_urb2d,tg_urb2d,tc_urb2d,qc_urb2d, & !H urban
742 uc_urb2d, & !H urban
743 xxxr_urb2d,xxxb_urb2d,xxxg_urb2d,xxxc_urb2d, & !H urban
744 trl_urb3d,tbl_urb3d,tgl_urb3d, & !H urban
745 sh_urb2d,lh_urb2d,g_urb2d,rn_urb2d,ts_urb2d, & !H urban
746 psim_urb2d,psih_urb2d,u10_urb2d,v10_urb2d, & !O urban
747 GZ1OZ0_urb2d, AKMS_URB2D, & !O urban
748 th2_urb2d,q2_urb2d,ust_urb2d, & !O urban
749 declin_urb,cosz_urb2d,omg_urb2d, & !I urban
750 xlat_urb2d, & !I urban
751 num_roof_layers, num_wall_layers, & !I urban
752 num_road_layers, DZR, DZB, DZG, & !I urban
753 FRC_URB2D, UTYPE_URB2D & ! urban
754 )
755
756
757 DO j=j_start(ij),j_end(ij)
758 DO i=i_start(ij),i_end(ij)
759 CHKLOWQ(I,J)= 1.0
760 SFCEVP(I,J)= SFCEVP(I,J) + QFX(I,J)*DTBL
761 ENDDO
762 ENDDO
763
764 CALL SFCDIAGS(HFX,QFX,TSK,QSFC,CHS2,CQS2,T2,TH2,Q2, &
765 PSFC,CP,R_d,RCP, &
766 ids,ide, jds,jde, kds,kde, &
767 ims,ime, jms,jme, kms,kme, &
768 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
769
770 !urban
771 IF(UCMCALL.eq.1) THEN
772 DO j=j_start(ij),j_end(ij) !urban
773 DO i=i_start(ij),i_end(ij) !urban
774 IF( IVGTYP(I,J) == 1 .or. IVGTYP(I,J) == 31 .or. & !urban
775 IVGTYP(I,J) == 32 .or. IVGTYP(I,J) == 33 ) THEN !urban
776 ! TH2(I,J) = TH2_URB2D(I,J) !urban
777 ! T2(I,J) = TH2_URB2D(I,J)/(1.E5/PSFC(I,J))**RCP !urban
778 !m T2(I,J) = TH2_URB2D(I,J) !urban
779 T2(I,J) = FRC_URB2D(i,j)*TH2_URB2D(I,J) + (1-FRC_URB2D(i,j))*T2(I,J) !urban
780 TH2(I,J) = T2(I,J)*(1.E5/PSFC(I,J))**RCP !urban
781 !m Q2(I,J) = Q2_URB2D(I,J) !urban
782 Q2(I,J) = FRC_URB2D(i,j)*Q2_URB2D(I,J) +(1-FRC_URB2D(i,j))* Q2(I,J) !urban
783 U10(I,J) = U10_URB2D(I,J) !urban
784 V10(I,J) = V10_URB2D(I,J) !urban
785 PSIM(I,J) = PSIM_URB2D(I,J) !urban
786 PSIH(I,J) = PSIH_URB2D(I,J) !urban
787 GZ1OZ0(I,J) = GZ1OZ0_URB2D(I,J) !urban
788 !m AKHS(I,J) = AKHS_URB2D(I,J) !urban
789 AKHS(I,J) = CHS(I,J) !urban
790 AKMS(I,J) = AKMS_URB2D(I,J) !urban
791 END IF !urban
792 ENDDO !urban
793 ENDDO !urban
794 ENDIF
795 !------------------------------------------------------------------
796
797 ELSE
798 CALL wrf_error_fatal('Lacking arguments for LSM in surface driver')
799 ENDIF
800
801 CASE (RUCLSMSCHEME)
802 IF (PRESENT(qv_curr) .AND. PRESENT(qc_curr) .AND. &
803 ! PRESENT(emiss) .AND. PRESENT(t2) .AND. &
804 PRESENT(qsg) .AND. PRESENT(qvg) .AND. &
805 PRESENT(qcg) .AND. PRESENT(soilt1) .AND. &
806 PRESENT(tsnav) .AND. PRESENT(smfr3d) .AND. &
807 PRESENT(keepfr3dflag) .AND. PRESENT(rainbl) .AND. &
808 .TRUE. ) THEN
809
810 IF( PRESENT(sr) ) THEN
811 frpcpn=.true.
812 ELSE
813 SR = 1.
814 ENDIF
815
816 CALL wrf_debug(100,'in RUC LSM')
817 CALL LSMRUC(dtbl,itimestep,num_soil_layers, &
818 zs,rainbl,snow,snowh,snowc,sr,frpcpn, &
819 dz8w,p8w,t_phy,qv_curr,qc_curr,rho, & !p8w in [pa]
820 glw,gsw,emiss,chklowq, &
821 flqc,flhc,mavail,canwat,vegfra,albedo,znt, &
822 snoalb, albbck, & !new
823 qsfc,qsg,qvg,qcg,soilt1,tsnav, &
824 tmn,ivgtyp,isltyp,xland,xice, &
825 cp,g,xlv,stbolt, &
826 smois,smstav,smstot,tslb,tsk,hfx,qfx,lh, &
827 sfcrunoff,udrunoff,sfcexc, &
828 sfcevp,grdflx,acsnow, &
829 smfr3d,keepfr3dflag, &
830 myj, &
831 ids,ide, jds,jde, kds,kde, &
832 ims,ime, jms,jme, kms,kme, &
833 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
834
835 IF(.not. MYJ) then
836
837 CALL SFCDIAGS(HFX,QFX,TSK,QVG,CHS2,CQS2,T2,TH2,Q2, &
838 PSFC,CP,R_d,RCP, &
839 ids,ide, jds,jde, kds,kde, &
840 ims,ime, jms,jme, kms,kme, &
841 i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte )
842 ENDIF
843
844
845 ELSE
846 CALL wrf_error_fatal('Lacking arguments for RUCLSM in surface driver')
847 ENDIF
848
849 CASE DEFAULT
850
851 WRITE( message , * ) &
852 'The surface option does not exist: sf_surface_physics = ', sf_surface_physics
853 CALL wrf_error_fatal ( message )
854
855 END SELECT sfc_select
856 ENDDO
857 !$OMP END PARALLEL DO
858
859 430 CONTINUE
860
861
862 ! Reset RAINBL in mm (Accumulation between PBL calls)
863
864 IF ( PRESENT( rainbl ) ) THEN
865 !$OMP PARALLEL DO &
866 !$OMP PRIVATE ( ij, i, j, k )
867 DO ij = 1 , num_tiles
868 DO j=j_start(ij),j_end(ij)
869 DO i=i_start(ij),i_end(ij)
870 RAINBL(i,j) = 0.
871 ENDDO
872 ENDDO
873 ENDDO
874 !$OMP END PARALLEL DO
875 ENDIF
876
877 ENDIF
878
879 END SUBROUTINE surface_driver
880
881 END MODULE module_surface_driver
882