Changeset 4328 for palm/trunk/SOURCE/advec_ws.f90

Timestamp:

Dec 9, 2019 6:53:04 PM (4 years ago)

Author:

suehring

Message:

Minor change in nzb_max computation; commentation added; minor formatting adjustments in advec_ws

File:

• palm/trunk/SOURCE/advec_ws.f90 (modified) (9 diffs)

palm/trunk/SOURCE/advec_ws.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor formatting adjustments
+! 
+! 4327 2019-12-06 14:48:31Z Giersch
 ! Setting of advection flags for vertical fluxes of w revised, air density for
 ! vertical flux calculation of w at k=1 is considered now
@@ -2088 +2091 @@
                                            )                                  
                                                                               
-              diss_l_u(k,j,tn) = - ABS( u_comp_l ) * (                         &
+             diss_l_u(k,j,tn) = - ABS( u_comp_l ) * (                          &
                               ( 10.0_wp * ibit2 * adv_mom_5                    &
                            +     3.0_wp * ibit1 * adv_mom_3                    &
@@ -2101 +2104 @@
                               ) *                                              &
                                         ( u(k,j,i+2) - u(k,j,i-3) )            &
-                                                     )
+                                                    )
 
           ENDDO
@@ -2387 +2390 @@
                 - ( u(k,j,i)   + u(k,j,i-1)   )                                &
                                     * (                                        &
-                     REAL( IBITS(advc_flags_m(k,j,i-1),0,1),  KIND = wp )      &
+                     REAL( IBITS(advc_flags_m(k,j,i-1),0,1), KIND = wp )       &
                    + REAL( IBITS(advc_flags_m(k,j,i-1),1,1), KIND = wp )       &
                    + REAL( IBITS(advc_flags_m(k,j,i-1),2,1), KIND = wp )       &
@@ -2591 +2594 @@
                                             )
 
-              diss_l_v(k,j,tn) = - ABS( u_comp(k) ) * (                       &
+             diss_l_v(k,j,tn) = - ABS( u_comp(k) ) * (                        &
                               ( 10.0_wp * ibit11 * adv_mom_5                  &
                            +     3.0_wp * ibit10 * adv_mom_3                  &
@@ -2604 +2607 @@
                               ) *                                             &
                                         ( v(k,j,i+2) - v(k,j,i-3) )           &
-                                                      )
+                                                     )
 
           ENDDO
@@ -2935 +2938 @@
           ibit11 = REAL( IBITS(advc_flags_m(k,j,i),11,1), KIND = wp )
           ibit10 = REAL( IBITS(advc_flags_m(k,j,i),10,1), KIND = wp )
-          ibit9  = REAL( IBITS(advc_flags_m(k,j,i),9,1), KIND = wp )
+          ibit9  = REAL( IBITS(advc_flags_m(k,j,i),9,1),  KIND = wp )
           
           ibit14 = REAL( IBITS(advc_flags_m(k,j,i),14,1), KIND = wp )
@@ -3213 +3216 @@
                                             )
 
-               diss_l_w(k,j,tn) = - ABS( u_comp(k) ) * (                      &
+             diss_l_w(k,j,tn) = - ABS( u_comp(k) ) * (                        &
                              ( 10.0_wp * ibit20 * adv_mom_5                   &
                           +     3.0_wp * ibit19 * adv_mom_3                   &
@@ -3226 +3229 @@
                              ) *                                              &
                                         ( w(k,j,i+2) - w(k,j,i-3) )           &
-                                                       )
+                                                     )
 
           ENDDO