Changeset 1359 for palm/trunk/SOURCE/lpm_init_sgs_tke.f90

Timestamp:

Apr 11, 2014 5:15:14 PM (10 years ago)

Author:

hoffmann

Message:

new Lagrangian particle structure integrated

File:

• palm/trunk/SOURCE/lpm_init_sgs_tke.f90 (modified) (6 diffs)

palm/trunk/SOURCE/lpm_init_sgs_tke.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! New particle structure integrated.
+! Kind definition added to all floating point numbers.
 !
 ! Former revisions:
@@ -84 +85 @@
                   k  > nzb_s_inner(j,i+1) )                                    &
              THEN
-                de_dx(k,j,i) = 2.0 * sgs_wfu_part * ( e(k,j,i+1) - e(k,j,i) )  &
-                               * ddx
+                de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
+                               ( e(k,j,i+1) - e(k,j,i) ) * ddx
              ELSEIF ( k  > nzb_s_inner(j,i-1)  .AND.  k > nzb_s_inner(j,i)     &
                       .AND.  k <= nzb_s_inner(j,i+1) )                         &
              THEN
-                de_dx(k,j,i) = 2.0 * sgs_wfu_part * ( e(k,j,i) - e(k,j,i-1) )  &
-                               * ddx
+                de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
+                               ( e(k,j,i) - e(k,j,i-1) ) * ddx
              ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j,i+1) )    &
              THEN
-                de_dx(k,j,i) = 0.0
+                de_dx(k,j,i) = 0.0_wp
              ELSEIF ( k < nzb_s_inner(j,i-1)  .AND.  k < nzb_s_inner(j,i) )    &
              THEN
-                de_dx(k,j,i) = 0.0
+                de_dx(k,j,i) = 0.0_wp
              ELSE
                 de_dx(k,j,i) = sgs_wfu_part * ( e(k,j,i+1) - e(k,j,i-1) ) * ddx
@@ -104 +105 @@
                   k  > nzb_s_inner(j+1,i) )                                    &
              THEN
-                de_dy(k,j,i) = 2.0 * sgs_wfv_part * ( e(k,j+1,i) - e(k,j,i) )  &
-                               * ddy
+                de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
+                               ( e(k,j+1,i) - e(k,j,i) ) * ddy
              ELSEIF ( k  > nzb_s_inner(j-1,i)  .AND.  k  > nzb_s_inner(j,i)    &
                       .AND.  k <= nzb_s_inner(j+1,i) )                         &
              THEN
-                de_dy(k,j,i) = 2.0 * sgs_wfv_part * ( e(k,j,i) - e(k,j-1,i) )  &
-                               * ddy
+                de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
+                               ( e(k,j,i) - e(k,j-1,i) ) * ddy
              ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j+1,i) )    &
              THEN
-                de_dy(k,j,i) = 0.0
+                de_dy(k,j,i) = 0.0_wp
              ELSEIF ( k < nzb_s_inner(j-1,i)  .AND.  k < nzb_s_inner(j,i) )    &
              THEN
-                de_dy(k,j,i) = 0.0
+                de_dy(k,j,i) = 0.0_wp
              ELSE
                 de_dy(k,j,i) = sgs_wfv_part * ( e(k,j+1,i) - e(k,j-1,i) ) * ddy
@@ -131 +132 @@
 
           DO  k = nzb_s_inner(j,i)+2, nzt-1
-             de_dz(k,j,i)  = 2.0 * sgs_wfw_part * &
+             de_dz(k,j,i)  = 2.0_wp * sgs_wfw_part *                           &
                              ( e(k+1,j,i) - e(k-1,j,i) ) / ( zu(k+1)-zu(k-1) )
           ENDDO
 
           k = nzb_s_inner(j,i)
-          de_dz(nzb:k,j,i) = 0.0
-          de_dz(k+1,j,i)   = 2.0 * sgs_wfw_part * ( e(k+2,j,i) - e(k+1,j,i) ) &
-                                                / ( zu(k+2) - zu(k+1) )
-          de_dz(nzt,j,i)   = 0.0
-          de_dz(nzt+1,j,i) = 0.0
+          de_dz(nzb:k,j,i) = 0.0_wp
+          de_dz(k+1,j,i)   = 2.0_wp * sgs_wfw_part *                           &
+                             ( e(k+2,j,i) - e(k+1,j,i) ) / ( zu(k+2) - zu(k+1) )
+          de_dz(nzt,j,i)   = 0.0_wp
+          de_dz(nzt+1,j,i) = 0.0_wp
        ENDDO
     ENDDO
@@ -162 +163 @@
 !--    First calculate horizontally averaged profiles of the horizontal
 !--    velocities.
-       sums_l(:,1,0) = 0.0
-       sums_l(:,2,0) = 0.0
+       sums_l(:,1,0) = 0.0_wp
+       sums_l(:,2,0) = 0.0_wp
 
        DO  i = nxl, nxr
@@ -197 +198 @@
 !--    Now calculate the profiles of SGS TKE and the resolved-scale
 !--    velocity variances
-       sums_l(:,8,0)  = 0.0
-       sums_l(:,30,0) = 0.0
-       sums_l(:,31,0) = 0.0
-       sums_l(:,32,0) = 0.0
+       sums_l(:,8,0)  = 0.0_wp
+       sums_l(:,30,0) = 0.0_wp
+       sums_l(:,31,0) = 0.0_wp
+       sums_l(:,32,0) = 0.0_wp
        DO  i = nxl, nxr
           DO  j = nys, nyn