source: palm/trunk/SOURCE/ls_forcing_mod.f90 @ 1850

Last change on this file since 1850 was 1850, checked in by maronga, 5 years ago

added _mod string to several filenames to meet the naming convection for modules

  • Property svn:keywords set to Id
File size: 21.8 KB
Line 
1!> @file ls_forcing_mod.f90
2!--------------------------------------------------------------------------------!
3! This file is part of PALM.
4!
5! PALM is free software: you can redistribute it and/or modify it under the terms
6! of the GNU General Public License as published by the Free Software Foundation,
7! either version 3 of the License, or (at your option) any later version.
8!
9! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
10! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
11! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
12!
13! You should have received a copy of the GNU General Public License along with
14! PALM. If not, see <http://www.gnu.org/licenses/>.
15!
16! Copyright 1997-2016 Leibniz Universitaet Hannover
17!--------------------------------------------------------------------------------!
18!
19! Current revisions:
20! ------------------
21! Module renamed
22!
23!
24! Former revisions:
25! -----------------
26! $Id: ls_forcing_mod.f90 1850 2016-04-08 13:29:27Z maronga $
27!
28! 1682 2015-10-07 23:56:08Z knoop
29! Code annotations made doxygen readable
30!
31! 1602 2015-06-22 07:50:56Z heinze
32! PA0370 changed to PA0363
33!
34! 1382 2014-04-30 12:15:41Z boeske
35! Renamed variables which store large scale forcing tendencies
36! pt_lsa -> td_lsa_lpt, pt_subs -> td_sub_lpt,
37! q_lsa  -> td_lsa_q,   q_subs  -> td_sub_q,
38! high|lowpt_lsa -> high|low_td_lsa_lpt, ...
39!
40! 1365 2014-04-22 15:03:56Z boeske
41! Usage of large scale forcing for pt and q enabled:
42! Added new subroutine ls_advec for horizontal large scale advection and large
43! scale subsidence,
44! error message in init_ls_forcing specified,
45! variable t renamed nt
46!
47! 1353 2014-04-08 15:21:23Z heinze
48! REAL constants provided with KIND-attribute
49!
50! 1320 2014-03-20 08:40:49Z raasch
51! ONLY-attribute added to USE-statements,
52! kind-parameters added to all INTEGER and REAL declaration statements,
53! kinds are defined in new module kinds,
54! comment fields (!:) to be used for variable explanations added to
55! all variable declaration statements
56!
57! 1318 2014-03-17 13:35:16Z raasch
58! module interfaces removed
59!
60! 1299 2014-03-06 13:15:21Z heinze
61! Ensure a zero large scale vertical velocity at the surface
62! Bugfix: typo in case of boundary condition in if-clause
63!
64! 1276 2014-01-15 13:40:41Z heinze
65! Use LSF_DATA also in case of Dirichlet bottom boundary condition for scalars
66!
67! 1249 2013-11-06 10:45:47Z heinze
68! remove call of user module
69! reformatting
70!
71! 1241 2013-10-30 11:36:58Z heinze
72! Initial revision
73!
74! Description:
75! ------------
76!> Calculates large scale forcings (geostrophic wind and subsidence velocity) as
77!> well as surfaces fluxes dependent on time given in an external file (LSF_DATA).
78!> Code is based in parts on DALES and UCLA-LES.
79!--------------------------------------------------------------------------------!
80 MODULE ls_forcing_mod
81 
82
83    PRIVATE
84    PUBLIC init_ls_forcing, ls_forcing_surf, ls_forcing_vert, ls_advec
85    SAVE
86
87    INTERFACE ls_advec
88       MODULE PROCEDURE ls_advec
89       MODULE PROCEDURE ls_advec_ij
90    END INTERFACE ls_advec
91
92 CONTAINS
93
94!------------------------------------------------------------------------------!
95! Description:
96! ------------
97!> @todo Missing subroutine description.
98!------------------------------------------------------------------------------!
99    SUBROUTINE init_ls_forcing
100
101       USE arrays_3d,                                                          &
102           ONLY:  p_surf, pt_surf, q_surf, qsws_surf, shf_surf, td_lsa_lpt,    &
103                  td_lsa_q, td_sub_lpt, td_sub_q, time_surf, time_vert,        &
104                  ug_vert, vg_vert, wsubs_vert, zu
105
106       USE control_parameters,                                                 &
107           ONLY:  end_time, lsf_surf, lsf_vert, message_string, nlsf
108
109       USE indices,                                                            &
110           ONLY:  ngp_sums_ls, nzb, nz, nzt
111
112       USE kinds
113
114       USE statistics,                                                         &
115           ONLY:  sums_ls_l
116
117
118       IMPLICIT NONE
119
120       CHARACTER(100) ::  chmess      !<
121       CHARACTER(1)   ::  hash        !<
122
123       INTEGER(iwp) ::  ierrn         !<
124       INTEGER(iwp) ::  finput = 90   !<
125       INTEGER(iwp) ::  k             !<
126       INTEGER(iwp) ::  nt             !<
127
128       REAL(wp) ::  fac               !<
129       REAL(wp) ::  highheight        !<
130       REAL(wp) ::  highug_vert       !<
131       REAL(wp) ::  highvg_vert       !<
132       REAL(wp) ::  highwsubs_vert    !<
133       REAL(wp) ::  lowheight         !<
134       REAL(wp) ::  lowug_vert        !<
135       REAL(wp) ::  lowvg_vert        !<
136       REAL(wp) ::  lowwsubs_vert     !<
137       REAL(wp) ::  high_td_lsa_lpt   !<
138       REAL(wp) ::  low_td_lsa_lpt    !<
139       REAL(wp) ::  high_td_lsa_q     !<
140       REAL(wp) ::  low_td_lsa_q      !<
141       REAL(wp) ::  high_td_sub_lpt   !<
142       REAL(wp) ::  low_td_sub_lpt    !<
143       REAL(wp) ::  high_td_sub_q     !<
144       REAL(wp) ::  low_td_sub_q      !<
145       REAL(wp) ::  r_dummy           !<
146
147       ALLOCATE( p_surf(0:nlsf), pt_surf(0:nlsf), q_surf(0:nlsf),              &
148                 qsws_surf(0:nlsf), shf_surf(0:nlsf),                          &
149                 td_lsa_lpt(nzb:nzt+1,0:nlsf), td_lsa_q(nzb:nzt+1,0:nlsf),     &
150                 td_sub_lpt(nzb:nzt+1,0:nlsf), td_sub_q(nzb:nzt+1,0:nlsf),     &
151                 time_vert(0:nlsf), time_surf(0:nlsf),                         &
152                 ug_vert(nzb:nzt+1,0:nlsf), vg_vert(nzb:nzt+1,0:nlsf),         &
153                 wsubs_vert(nzb:nzt+1,0:nlsf) )
154
155       p_surf = 0.0_wp; pt_surf = 0.0_wp; q_surf = 0.0_wp; qsws_surf = 0.0_wp
156       shf_surf = 0.0_wp; time_vert = 0.0_wp; td_lsa_lpt = 0.0_wp
157       td_lsa_q = 0.0_wp; td_sub_lpt = 0.0_wp; td_sub_q = 0.0_wp
158       time_surf = 0.0_wp; ug_vert = 0.0_wp; vg_vert = 0.0_wp
159       wsubs_vert = 0.0_wp
160
161!
162!--    Array for storing large scale forcing and nudging tendencies at each
163!--    timestep for data output
164       ALLOCATE( sums_ls_l(nzb:nzt+1,0:7) )
165       sums_ls_l = 0.0_wp
166
167       ngp_sums_ls = (nz+2)*6
168
169       OPEN ( finput, FILE='LSF_DATA', STATUS='OLD', &
170              FORM='FORMATTED', IOSTAT=ierrn )
171
172       IF ( ierrn /= 0 )  THEN
173          message_string = 'file LSF_DATA does not exist'
174          CALL message( 'ls_forcing', 'PA0368', 1, 2, 0, 6, 0 )
175       ENDIF
176
177       ierrn = 0
178!
179!--    First three lines of LSF_DATA contain header
180       READ ( finput, FMT='(a100)', IOSTAT=ierrn ) chmess
181       READ ( finput, FMT='(a100)', IOSTAT=ierrn ) chmess
182       READ ( finput, FMT='(a100)', IOSTAT=ierrn ) chmess
183
184       IF ( ierrn /= 0 )  THEN
185          message_string = 'errors in file LSF_DATA'
186          CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
187       ENDIF
188
189!
190!--    Surface values are read in
191       nt     = 0
192       ierrn = 0
193
194       DO WHILE ( time_surf(nt) < end_time )
195          nt = nt + 1
196          READ ( finput, *, IOSTAT = ierrn ) time_surf(nt), shf_surf(nt),      &
197                                             qsws_surf(nt), pt_surf(nt),       &
198                                             q_surf(nt), p_surf(nt)
199
200          IF ( ierrn < 0 )  THEN
201            WRITE ( message_string, * ) 'No time dependent surface variables ',&
202                              'in&LSF_DATA for end of run found'
203
204             CALL message( 'ls_forcing', 'PA0363', 1, 2, 0, 6, 0 )
205          ENDIF
206       ENDDO
207
208
209       IF ( time_surf(1) > end_time )  THEN
210          WRITE ( message_string, * ) 'No time dependent surface variables in ',&
211                                     '&LSF_DATA for end of run found - ',      &
212                                     'lsf_surf is set to FALSE'
213          CALL message( 'ls_forcing', 'PA0371', 0, 0, 0, 6, 0 )
214          lsf_surf = .FALSE.
215       ENDIF
216
217!
218!--    Go to the end of the list with surface variables
219       DO WHILE ( ierrn == 0 )
220          READ ( finput, *, IOSTAT = ierrn ) r_dummy
221       ENDDO
222
223!
224!--    Profiles of ug, vg and w_subs are read in (large scale forcing)
225
226       nt = 0
227       DO WHILE ( time_vert(nt) < end_time )
228          nt = nt + 1
229          hash = "#"
230          ierrn = 1 ! not zero
231!
232!--       Search for the next line consisting of "# time",
233!--       from there onwards the profiles will be read
234          DO WHILE ( .NOT. ( hash == "#" .AND. ierrn == 0 ) ) 
235             READ ( finput, *, IOSTAT=ierrn ) hash, time_vert(nt)
236             IF ( ierrn < 0 )  THEN
237                WRITE( message_string, * ) 'No time dependent vertical profiles',&
238                                 ' in&LSF_DATA for end of run found'
239                CALL message( 'ls_forcing', 'PA0372', 1, 2, 0, 6, 0 )
240             ENDIF
241          ENDDO
242
243          IF ( nt == 1 .AND. time_vert(nt) > end_time ) EXIT
244
245          READ ( finput, *, IOSTAT=ierrn ) lowheight, lowug_vert, lowvg_vert,  &
246                                           lowwsubs_vert, low_td_lsa_lpt,      &
247                                           low_td_lsa_q, low_td_sub_lpt,       &
248                                           low_td_sub_q
249          IF ( ierrn /= 0 )  THEN
250             message_string = 'errors in file LSF_DATA'
251             CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
252          ENDIF
253
254          READ ( finput, *, IOSTAT=ierrn ) highheight, highug_vert,            &
255                                           highvg_vert, highwsubs_vert,        &
256                                           high_td_lsa_lpt, high_td_lsa_q,     &
257                                           high_td_sub_lpt, high_td_sub_q
258     
259          IF ( ierrn /= 0 )  THEN
260             message_string = 'errors in file LSF_DATA'
261             CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
262          ENDIF
263
264
265          DO  k = nzb, nzt+1
266             IF ( highheight < zu(k) )  THEN
267                lowheight      = highheight
268                lowug_vert     = highug_vert
269                lowvg_vert     = highvg_vert
270                lowwsubs_vert  = highwsubs_vert
271                low_td_lsa_lpt = high_td_lsa_lpt
272                low_td_lsa_q   = high_td_lsa_q
273                low_td_sub_lpt = high_td_sub_lpt
274                low_td_sub_q   = high_td_sub_q
275
276                ierrn = 0
277                READ ( finput, *, IOSTAT=ierrn ) highheight, highug_vert,      &
278                                                 highvg_vert, highwsubs_vert,  &
279                                                 high_td_lsa_lpt,              &
280                                                 high_td_lsa_q,                &
281                                                 high_td_sub_lpt, high_td_sub_q
282
283                IF ( ierrn /= 0 )  THEN
284                   WRITE( message_string, * ) 'zu(nzt+1) = ', zu(nzt+1), 'm ', &
285                        'is higher than the maximum height in LSF_DATA which ',&
286                        'is ', lowheight, 'm. Interpolation on PALM ',         &
287                        'grid is not possible.'
288                   CALL message( 'ls_forcing', 'PA0395', 1, 2, 0, 6, 0 )
289                ENDIF
290
291             ENDIF
292
293!
294!--          Interpolation of prescribed profiles in space
295             fac = (highheight-zu(k))/(highheight - lowheight)
296
297             ug_vert(k,nt)    = fac * lowug_vert                               &
298                                + ( 1.0_wp - fac ) * highug_vert
299             vg_vert(k,nt)    = fac * lowvg_vert                               &
300                                + ( 1.0_wp - fac ) * highvg_vert
301             wsubs_vert(k,nt) = fac * lowwsubs_vert                            &
302                                + ( 1.0_wp - fac ) * highwsubs_vert
303
304             td_lsa_lpt(k,nt) = fac * low_td_lsa_lpt                           &
305                                + ( 1.0_wp - fac ) * high_td_lsa_lpt
306             td_lsa_q(k,nt)   = fac * low_td_lsa_q                             &
307                                + ( 1.0_wp - fac ) * high_td_lsa_q
308             td_sub_lpt(k,nt) = fac * low_td_sub_lpt                           &
309                                + ( 1.0_wp - fac ) * high_td_sub_lpt
310             td_sub_q(k,nt)   = fac * low_td_sub_q                             &
311                                + ( 1.0_wp - fac ) * high_td_sub_q
312
313          ENDDO
314
315       ENDDO
316
317!
318!--    Large scale vertical velocity has to be zero at the surface
319       wsubs_vert(nzb,:) = 0.0_wp
320 
321       IF ( time_vert(1) > end_time )  THEN
322          WRITE ( message_string, * ) 'Time dependent large scale profile ',   &
323                             'forcing from&LSF_DATA sets in after end of ' ,   &
324                             'simulation - lsf_vert is set to FALSE'
325          CALL message( 'ls_forcing', 'PA0373', 0, 0, 0, 6, 0 )
326          lsf_vert = .FALSE.
327       ENDIF
328
329       CLOSE( finput )
330
331
332    END SUBROUTINE init_ls_forcing
333
334
335!------------------------------------------------------------------------------!
336! Description:
337! ------------
338!> @todo Missing subroutine description.
339!------------------------------------------------------------------------------!
340    SUBROUTINE ls_forcing_surf ( time )
341
342       USE arrays_3d,                                                          &
343           ONLY:  p_surf, pt_surf, q_surf, qsws, qsws_surf, shf, shf_surf,     &
344                  time_surf, time_vert, ug, ug_vert, vg, vg_vert
345
346       USE control_parameters,                                                 &
347           ONLY:  bc_q_b, ibc_pt_b, ibc_q_b, pt_surface, q_surface,            &
348                  surface_pressure
349
350       USE kinds
351
352       IMPLICIT NONE
353
354       INTEGER(iwp) ::  nt                     !<
355
356       REAL(wp)             :: fac            !<
357       REAL(wp), INTENT(in) :: time           !<
358
359!
360!--    Interpolation in time of LSF_DATA at the surface
361       nt = 1
362       DO WHILE ( time > time_surf(nt) )
363          nt = nt + 1
364       ENDDO
365       IF ( time /= time_surf(nt) )  THEN
366         nt = nt - 1
367       ENDIF
368
369       fac = ( time -time_surf(nt) ) / ( time_surf(nt+1) - time_surf(nt) )
370
371       IF ( ibc_pt_b == 0 )  THEN
372!
373!--       In case of Dirichlet boundary condition shf must not
374!--       be set - it is calculated via MOST in prandtl_fluxes
375          pt_surface = pt_surf(nt) + fac * ( pt_surf(nt+1) - pt_surf(nt) )
376
377       ELSEIF ( ibc_pt_b == 1 )  THEN
378!
379!--       In case of Neumann boundary condition pt_surface is needed for
380!--       calculation of reference density
381          shf        = shf_surf(nt) + fac * ( shf_surf(nt+1) - shf_surf(nt) )
382          pt_surface = pt_surf(nt) + fac * ( pt_surf(nt+1) - pt_surf(nt) )
383
384       ENDIF
385
386       IF ( ibc_q_b == 0 )  THEN
387!
388!--       In case of Dirichlet boundary condition qsws must not
389!--       be set - it is calculated via MOST in prandtl_fluxes
390          q_surface = q_surf(nt) + fac * ( q_surf(nt+1) - q_surf(nt) )
391
392       ELSEIF ( ibc_q_b == 1 )  THEN
393
394          qsws = qsws_surf(nt) + fac * ( qsws_surf(nt+1) - qsws_surf(nt) )
395
396       ENDIF
397
398       surface_pressure = p_surf(nt) + fac * ( p_surf(nt+1) - p_surf(nt) )
399
400    END SUBROUTINE ls_forcing_surf
401
402
403!------------------------------------------------------------------------------!
404! Description:
405! ------------
406!> @todo Missing subroutine description.
407!------------------------------------------------------------------------------!
408    SUBROUTINE ls_forcing_vert ( time )
409
410       USE arrays_3d,                                                          &
411           ONLY:  time_vert, ug, ug_vert, vg, vg_vert, w_subs, wsubs_vert
412
413       USE control_parameters,                                                 &
414           ONLY:  large_scale_subsidence
415
416       USE kinds
417
418       IMPLICIT NONE
419
420       INTEGER(iwp) ::  nt                     !<
421
422       REAL(wp)             ::  fac           !<
423       REAL(wp), INTENT(in) ::  time          !<
424
425!
426!--    Interpolation in time of LSF_DATA for ug, vg and w_subs
427       nt = 1
428       DO WHILE ( time > time_vert(nt) )
429          nt = nt + 1
430       ENDDO
431       IF ( time /= time_vert(nt) )  THEN
432         nt = nt - 1
433       ENDIF
434
435       fac = ( time-time_vert(nt) ) / ( time_vert(nt+1)-time_vert(nt) )
436
437       ug     = ug_vert(:,nt) + fac * ( ug_vert(:,nt+1) - ug_vert(:,nt) )
438       vg     = vg_vert(:,nt) + fac * ( vg_vert(:,nt+1) - vg_vert(:,nt) )
439
440       IF ( large_scale_subsidence )  THEN
441          w_subs = wsubs_vert(:,nt)                                            &
442                   + fac * ( wsubs_vert(:,nt+1) - wsubs_vert(:,nt) )
443       ENDIF
444
445    END SUBROUTINE ls_forcing_vert
446
447
448!------------------------------------------------------------------------------!
449! Description:
450! ------------
451!> Call for all grid points
452!------------------------------------------------------------------------------!
453    SUBROUTINE ls_advec ( time, prog_var )
454
455       USE arrays_3d,                                                          &
456           ONLY:  td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q, tend, time_vert       
457       
458       USE control_parameters,                                                 &
459           ONLY:  large_scale_subsidence, use_subsidence_tendencies
460       
461       USE indices
462       
463       USE kinds
464
465       IMPLICIT NONE
466
467       CHARACTER (LEN=*) ::  prog_var   !<
468
469       REAL(wp), INTENT(in)  :: time    !<
470       REAL(wp) :: fac                  !< 
471
472       INTEGER(iwp) ::  i               !<
473       INTEGER(iwp) ::  j               !<
474       INTEGER(iwp) ::  k               !<
475       INTEGER(iwp) ::  nt               !<
476
477!
478!--    Interpolation in time of LSF_DATA
479       nt = 1
480       DO WHILE ( time > time_vert(nt) )
481          nt = nt + 1
482       ENDDO
483       IF ( time /= time_vert(nt) )  THEN
484         nt = nt - 1
485       ENDIF
486
487       fac = ( time-time_vert(nt) ) / ( time_vert(nt+1)-time_vert(nt) )
488
489!
490!--    Add horizontal large scale advection tendencies of pt and q
491       SELECT CASE ( prog_var )
492
493          CASE ( 'pt' )
494
495             DO  i = nxl, nxr
496                DO  j = nys, nyn
497                   DO  k = nzb_u_inner(j,i)+1, nzt
498                      tend(k,j,i) = tend(k,j,i) + td_lsa_lpt(k,nt) + fac *     &
499                                    ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) )
500                   ENDDO
501                ENDDO
502             ENDDO
503
504          CASE ( 'q' )
505
506             DO  i = nxl, nxr
507                DO  j = nys, nyn
508                   DO  k = nzb_u_inner(j,i)+1, nzt
509                      tend(k,j,i) = tend(k,j,i) + td_lsa_q(k,nt) + fac *       &
510                                    ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) )
511                   ENDDO
512                ENDDO
513             ENDDO
514
515       END SELECT
516
517!
518!--    Subsidence of pt and q with prescribed subsidence tendencies
519       IF ( large_scale_subsidence .AND. use_subsidence_tendencies )  THEN
520
521          SELECT CASE ( prog_var )
522
523             CASE ( 'pt' )
524
525                DO  i = nxl, nxr
526                   DO  j = nys, nyn
527                      DO  k = nzb_u_inner(j,i)+1, nzt
528                         tend(k,j,i) = tend(k,j,i) + td_sub_lpt(k,nt) + fac *  &
529                                       ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) )
530                      ENDDO
531                   ENDDO
532                ENDDO
533 
534             CASE ( 'q' )
535
536                DO  i = nxl, nxr
537                   DO  j = nys, nyn
538                      DO  k = nzb_u_inner(j,i)+1, nzt
539                         tend(k,j,i) = tend(k,j,i) + td_sub_q(k,nt) + fac *    &
540                                       ( td_sub_q(k,nt+1) - td_sub_q(k,nt) )
541                      ENDDO
542                   ENDDO
543                ENDDO
544
545          END SELECT
546
547       ENDIF
548
549    END SUBROUTINE ls_advec
550
551
552!------------------------------------------------------------------------------!
553! Description:
554! ------------
555!> Call for grid point i,j
556!------------------------------------------------------------------------------!
557    SUBROUTINE ls_advec_ij ( i, j, time, prog_var )
558
559       USE arrays_3d,                                                          &
560           ONLY:  td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q, tend, time_vert       
561       
562       USE control_parameters,                                                 &
563           ONLY:  large_scale_subsidence, use_subsidence_tendencies
564       
565       USE indices
566       
567       USE kinds
568
569       IMPLICIT NONE
570
571       CHARACTER (LEN=*) ::  prog_var   !<
572
573       REAL(wp), INTENT(in)  :: time    !<
574       REAL(wp) :: fac                  !<
575
576       INTEGER(iwp) ::  i               !<
577       INTEGER(iwp) ::  j               !<
578       INTEGER(iwp) ::  k               !<
579       INTEGER(iwp) ::  nt               !<
580
581!
582!--    Interpolation in time of LSF_DATA
583       nt = 1
584       DO WHILE ( time > time_vert(nt) )
585          nt = nt + 1
586       ENDDO
587       IF ( time /= time_vert(nt) )  THEN
588         nt = nt - 1
589       ENDIF
590
591       fac = ( time-time_vert(nt) ) / ( time_vert(nt+1)-time_vert(nt) )
592
593!
594!--    Add horizontal large scale advection tendencies of pt and q
595       SELECT CASE ( prog_var )
596
597          CASE ( 'pt' )
598
599             DO  k = nzb_u_inner(j,i)+1, nzt
600                tend(k,j,i) = tend(k,j,i) + td_lsa_lpt(k,nt)                   &
601                              + fac * ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) )
602             ENDDO
603
604          CASE ( 'q' )
605
606             DO  k = nzb_u_inner(j,i)+1, nzt
607                tend(k,j,i) = tend(k,j,i) + td_lsa_q(k,nt)                     &
608                              + fac * ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) )
609             ENDDO
610
611       END SELECT
612
613!
614!--    Subsidence of pt and q with prescribed profiles
615       IF ( large_scale_subsidence .AND. use_subsidence_tendencies )  THEN
616
617          SELECT CASE ( prog_var )
618
619             CASE ( 'pt' )
620
621                DO  k = nzb_u_inner(j,i)+1, nzt
622                   tend(k,j,i) = tend(k,j,i) + td_sub_lpt(k,nt) + fac *        &
623                                 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) )
624                ENDDO
625 
626             CASE ( 'q' )
627
628                DO  k = nzb_u_inner(j,i)+1, nzt
629                   tend(k,j,i) = tend(k,j,i) + td_sub_q(k,nt) + fac *          &
630                                 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) )
631                ENDDO
632
633          END SELECT
634
635       ENDIF
636
637    END SUBROUTINE ls_advec_ij
638
639
640 END MODULE ls_forcing_mod
Note: See TracBrowser for help on using the repository browser.