Changeset 1781 for palm/trunk/SOURCE/pmc_interface.f90

Timestamp:

Mar 3, 2016 3:12:23 PM (8 years ago)

Author:

raasch

Message:

some pmc calls completely removed, some moved from time_integration to main program, t_..m arrays removed from pmc

File:

• palm/trunk/SOURCE/pmc_interface.f90 (modified) (4 diffs)

palm/trunk/SOURCE/pmc_interface.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! _p arrays are set zero within buildings too, t.._m arrays and respective
+! settings within buildings completely removed
 !
 ! Former revisions:
@@ -53 +54 @@
 
 #if defined( __nopointer )
-!-- TO_DO: is it really required to use arrays te_m, tu_m, etc. and to set them
-!--        to zero withi buildings?
     USE arrays_3d,                                                             &
-        ONLY:  dzu, dzw, e, e_p, pt, pt_p, q, q_p, te_m, tu_m, tv_m, tw_m, u,  &
-               u_p, v, v_p, w, w_p, zu, zw, z0
+        ONLY:  dzu, dzw, e, e_p, pt, pt_p, q, q_p, u, u_p, v, v_p, w, w_p, zu, &
+               zw, z0
 #else
    USE arrays_3d,                                                              &
-        ONLY:  dzu, dzw, e, e_p, e_1, e_2, pt, pt_p, pt_1, pt_2, q, q_p, q_1, q_2, te_m, tu_m, &
-               tv_m, tw_m, u, u_p, u_1, u_2, v, v_p, v_1, v_2, w, w_p, w_1, w_2, zu, zw, z0
+        ONLY:  dzu, dzw, e, e_p, e_1, e_2, pt, pt_p, pt_1, pt_2, q, q_p, q_1,  &
+               q_2, u, u_p, u_1, u_2, v, v_p, v_1, v_2, w, w_p, w_1, w_2, zu,  &
+               zw, z0
 #endif
 
@@ -2337 +2337 @@
              DO   i = nxlg, nxrg
                 DO   j = nysg, nyng
-                   u(nzb:nzb_u_inner(j,i),j,i)     = 0.0_wp
-                   v(nzb:nzb_v_inner(j,i),j,i)     = 0.0_wp
-                   w(nzb:nzb_w_inner(j,i),j,i)     = 0.0_wp
-                   e(nzb:nzb_w_inner(j,i),j,i)     = 0.0_wp
-                   tu_m(nzb:nzb_u_inner(j,i),j,i)  = 0.0_wp
-                   tv_m(nzb:nzb_v_inner(j,i),j,i)  = 0.0_wp
-                   tw_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0_wp
-                   te_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0_wp
+                   u(nzb:nzb_u_inner(j,i),j,i) = 0.0_wp
+                   v(nzb:nzb_v_inner(j,i),j,i) = 0.0_wp
+                   w(nzb:nzb_w_inner(j,i),j,i) = 0.0_wp
+                   e(nzb:nzb_w_inner(j,i),j,i) = 0.0_wp
                 ENDDO
              ENDDO
@@ -3621 +3617 @@
           DO   i = nxlg, nxrg
              DO   j = nysg, nyng
-                u(nzb:nzb_u_inner(j,i),j,i)     = 0.0_wp
-                v(nzb:nzb_v_inner(j,i),j,i)     = 0.0_wp
-                w(nzb:nzb_w_inner(j,i),j,i)     = 0.0_wp
-                e(nzb:nzb_w_inner(j,i),j,i)     = 0.0_wp              
-                tu_m(nzb:nzb_u_inner(j,i),j,i)  = 0.0_wp
-                tv_m(nzb:nzb_v_inner(j,i),j,i)  = 0.0_wp
-                tw_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0_wp
-                te_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0_wp
+                u(nzb:nzb_u_inner(j,i),j,i)   = 0.0_wp
+                v(nzb:nzb_v_inner(j,i),j,i)   = 0.0_wp
+                w(nzb:nzb_w_inner(j,i),j,i)   = 0.0_wp
+                e(nzb:nzb_s_inner(j,i),j,i)   = 0.0_wp
+                u_p(nzb:nzb_u_inner(j,i),j,i) = 0.0_wp
+                v_p(nzb:nzb_v_inner(j,i),j,i) = 0.0_wp
+                w_p(nzb:nzb_w_inner(j,i),j,i) = 0.0_wp
+                e_p(nzb:nzb_s_inner(j,i),j,i) = 0.0_wp
              ENDDO
           ENDDO