Changeset 4648 for palm/trunk/SOURCE/init_slope.f90

Timestamp:

Aug 25, 2020 7:52:08 AM (4 years ago)

Author:

raasch

Message:

files re-formatted to follow the PALM coding standard

File:

• palm/trunk/SOURCE/init_slope.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_slope.f90

@@ -1 +1 @@
 !> @file init_slope.f90
-!------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
 ! This file is part of the PALM model system.
 !
-! PALM is free software: you can redistribute it and/or modify it under the
-! terms of the GNU General Public License as published by the Free Software
-! Foundation, either version 3 of the License, or (at your option) any later
-! version.
+! PALM is free software: you can redistribute it and/or modify it under the terms of the GNU General
+! Public License as published by the Free Software Foundation, either version 3 of the License, or
+! (at your option) any later version.
 !
-! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
-! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
-! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+! PALM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the
+! implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
+! Public License for more details.
 !
-! You should have received a copy of the GNU General Public License along with
-! PALM. If not, see <http://www.gnu.org/licenses/>.
+! You should have received a copy of the GNU General Public License along with PALM. If not, see
+! <http://www.gnu.org/licenses/>.
 !
 ! Copyright 1997-2020 Leibniz Universitaet Hannover
-!------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! file re-formatted to follow the PALM coding standard
+!
+! 4360 2020-01-07 11:25:50Z suehring
 ! Corrected "Former revisions" section
-! 
+!
 ! 3655 2019-01-07 16:51:22Z knoop
 ! Modularization of all bulk cloud physics code components
@@ -36 +38 @@
 ! Description:
 ! ------------
-!> Initialization of the temperature field and other variables used in case
-!> of a sloping surface.
+!> Initialization of the temperature field and other variables used in case of a sloping surface.
 !> @note when a sloping surface is used, only one constant temperature
 !>       gradient is allowed!
-!------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
  SUBROUTINE init_slope
- 
 
-    USE arrays_3d,                                                             &
+
+    USE arrays_3d,                                                                                 &
         ONLY:  pt, pt_init, pt_slope_ref, zu
-        
-    USE basic_constants_and_equations_mod,                                     &
+
+    USE basic_constants_and_equations_mod,                                                         &
         ONLY:  pi
-                    
-    USE control_parameters,                                                    &
-        ONLY:  alpha_surface, initializing_actions, pt_slope_offset,           &
-               pt_surface, pt_vertical_gradient, sin_alpha_surface
-        
-    USE grid_variables,                                                        &
+
+    USE control_parameters,                                                                        &
+        ONLY:  alpha_surface, initializing_actions, pt_slope_offset, pt_surface,                   &
+               pt_vertical_gradient, sin_alpha_surface
+
+    USE grid_variables,                                                                            &
         ONLY:  dx
-        
-    USE indices,                                                               &
+
+    USE indices,                                                                                   &
         ONLY:  ngp_2dh, nx, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
-        
+
     USE kinds
 
@@ -70 +71 @@
     INTEGER(iwp) ::  j        !<
     INTEGER(iwp) ::  k        !<
-    
+
     REAL(wp)     ::  alpha    !<
     REAL(wp)     ::  height   !<
     REAL(wp)     ::  pt_value !<
     REAL(wp)     ::  radius   !<
-    
+
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt_init_local !<
 
@@ -86 +87 @@
 
 !
-!--       Compute height of grid-point relative to lower left corner of
-!--       the total domain.
-!--       First compute the distance between the actual grid point and the
-!--       lower left corner as well as the angle between the line connecting
-!--       these points and the bottom of the model.
+!--       Compute height of grid-point relative to lower left corner of the total domain.
+!--       First compute the distance between the actual grid point and the lower left corner as well
+!--       as the angle between the line connecting these points and the bottom of the model.
           IF ( k /= nzb )  THEN
              radius = SQRT( ( i * dx )**2 + zu(k)**2 )
@@ -106 +105 @@
 !--       Compute temperatures in the rotated coordinate system
           alpha    = alpha + alpha_surface / 180.0_wp * pi
-          pt_value = pt_surface + radius * SIN( alpha ) * &
-                                  pt_vertical_gradient(1) / 100.0_wp
+          pt_value = pt_surface + radius * SIN( alpha ) * pt_vertical_gradient(1) / 100.0_wp
           pt_slope_ref(k,i) = pt_value
-       ENDDO                
+       ENDDO
     ENDDO
 
 !
-!-- Temperature difference between left and right boundary of the total domain,
-!-- used for the cyclic boundary in x-direction
-    pt_slope_offset = (nx+1) * dx * sin_alpha_surface * &
-                      pt_vertical_gradient(1) / 100.0_wp
+!-- Temperature difference between left and right boundary of the total domain, used for the cyclic
+!-- boundary in x-direction
+    pt_slope_offset = (nx+1) * dx * sin_alpha_surface * pt_vertical_gradient(1) / 100.0_wp
 
 
@@ -129 +126 @@
 
 !
-!--    Recompute the mean initial temperature profile (mean along x-direction of
-!--    the rotated coordinate system)
+!--    Recompute the mean initial temperature profile (mean along x-direction of the rotated
+!--    coordinate system)
        ALLOCATE( pt_init_local(nzb:nzt+1) )
        pt_init_local = 0.0_wp
@@ -143 +140 @@
 #if defined( __parallel )
        IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-       CALL MPI_ALLREDUCE( pt_init_local, pt_init, nzt+2-nzb, MPI_REAL, &
-                            MPI_SUM, comm2d, ierr )
+       CALL MPI_ALLREDUCE( pt_init_local, pt_init, nzt+2-nzb, MPI_REAL, MPI_SUM, comm2d, ierr )
 #else
        pt_init = pt_init_local