Changeset 4648 for palm/trunk/SOURCE/init_vertical_profiles.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_vertical_profiles.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_vertical_profiles.f90

@@ -1 +1 @@
 !> @file ocean_mod.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 2017-2020 Leibniz Universitaet Hannover
-!--------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
-! split from check_parameters as separate file to avoid circular dependency
-! with ocean_mod
+! file re-formatted to follow the PALM coding standard
 !
-! 
+! 4481 2020-03-31 18:55:54Z maronga
+! split from check_parameters as separate file to avoid circular dependency with ocean_mod
+!
+!
 !
 !
@@ -38 +39 @@
 ! ------------
 !> Inititalizes the vertical profiles of scalar quantities.
-!------------------------------------------------------------------------------!
- SUBROUTINE init_vertical_profiles( vertical_gradient_level_ind,               &
-                                    vertical_gradient_level,                   &
-                                    vertical_gradient, initial_profile,        &
+!--------------------------------------------------------------------------------------------------!
+ SUBROUTINE init_vertical_profiles( vertical_gradient_level_ind,                                   &
+                                    vertical_gradient_level,                                       &
+                                    vertical_gradient, initial_profile,                            &
                                     surface_value, bc_top_gradient )
 
-    USE arrays_3d,                                                             &
+    USE arrays_3d,                                                                                 &
         ONLY:  dzu, zu
 
-    USE control_parameters,                                                    &
+    USE control_parameters,                                                                        &
         ONLY:  ocean_mode
 
-    USE indices,                                                               &
+    USE indices,                                                                                   &
         ONLY:  nz, nzt
 
@@ -59 +60 @@
     INTEGER(iwp) ::  i  !< loop counter
     INTEGER(iwp) ::  k  !< loop counter
+
     INTEGER(iwp), DIMENSION(1:10) ::  vertical_gradient_level_ind  !< vertical grid indices for gradient levels
 
@@ -77 +79 @@
        DO  k = 1, nzt+1
           IF ( i < 11 )  THEN
-             IF ( vertical_gradient_level(i) < zu(k)  .AND.            &
+             IF ( vertical_gradient_level(i) < zu(k)  .AND.                                        &
                   vertical_gradient_level(i) >= 0.0_wp )  THEN
                 gradient = vertical_gradient(i) / 100.0_wp
@@ -103 +105 @@
 
 !
-!--    In ocean mode, profiles are constructed starting from the ocean surface,
-!--    which is at the top of the model domain
+!--    In ocean mode, profiles are constructed starting from the ocean surface, which is at the top
+!--    of the model domain
        vertical_gradient_level_ind(1) = nzt+1
        DO  k = nzt, 0, -1
           IF ( i < 11 )  THEN
-             IF ( vertical_gradient_level(i) > zu(k)  .AND.            &
+             IF ( vertical_gradient_level(i) > zu(k)  .AND.                                        &
                   vertical_gradient_level(i) <= 0.0_wp )  THEN
                 gradient = vertical_gradient(i) / 100.0_wp
@@ -119 +121 @@
                 initial_profile(k) = initial_profile(k+1) - dzu(k+1) * gradient
              ELSE
-                initial_profile(k)   = surface_value - 0.5_wp * dzu(k+1) *     &
-                                                       gradient
-                initial_profile(k+1) = surface_value + 0.5_wp * dzu(k+1) *     &
-                                                       gradient
+                initial_profile(k)   = surface_value - 0.5_wp * dzu(k+1) * gradient
+                initial_profile(k+1) = surface_value + 0.5_wp * dzu(k+1) * gradient
              ENDIF
           ELSE
@@ -143 +143 @@
 !
 !-- Store gradient at the top boundary for possible Neumann boundary condition
-    bc_top_gradient  = ( initial_profile(nzt+1) - initial_profile(nzt) ) /     &
-                       dzu(nzt+1)
+    bc_top_gradient  = ( initial_profile(nzt+1) - initial_profile(nzt) ) / dzu(nzt+1)
 
  END SUBROUTINE init_vertical_profiles