Index: palm/trunk/SOURCE/header.f90
===================================================================
--- palm/trunk/SOURCE/header.f90	(revision 4573)
+++ palm/trunk/SOURCE/header.f90	(revision 4586)
@@ -20,56 +20,60 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Renamed rif to Ri (gradient Richardson number, 1D model)
+! and zeta (= z_mo / ol, stability parameter, 3D model)
+!
+! 4573 2020-06-24 13:08:47Z oliver.maas
 ! added statement for pt_surface_heating_rate
-! 
+!
 ! 4536 2020-05-17 17:24:13Z raasch
 ! output of restart data format added
-! 
+!
 ! 4473 2020-03-25 21:04:07Z gronemeier
 ! revised message if wall_adjustment is used
-! 
+!
 ! 4444 2020-03-05 15:59:50Z raasch
 ! bugfix: cpp-directives for serial mode added
-! 
+!
 ! 4360 2020-01-07 11:25:50Z suehring
 ! Bugfix, character length too short, caused crash on NEC.
-! 
+!
 ! 4309 2019-11-26 18:49:59Z suehring
 ! replaced recycling_yshift by y_shift
-! 
+!
 ! 4301 2019-11-22 12:09:09Z oliver.maas
 !
 ! 4297 2019-11-21 10:37:50Z oliver.maas
 ! Adjusted format for simulated time and related quantities
-! 
+!
 ! 4297 2019-11-21 10:37:50Z oliver.maas
 ! adjusted message to the changed parameter recycling_yshift
-! 
+!
 ! 4227 2019-09-10 18:04:34Z gronemeier
 ! implement new palm_date_time_mod
-! 
+!
 ! 4223 2019-09-10 09:20:47Z gronemeier
 ! Write information about rotation angle
-! 
+!
 ! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
-! 
+!
 ! 4168 2019-08-16 13:50:17Z suehring
 ! Replace function get_topography_top_index by topo_top_ind
-! 
+!
 ! 4069 2019-07-01 14:05:51Z Giersch
-! Masked output running index mid has been introduced as a local variable to 
-! avoid runtime error (Loop variable has been modified) in time_integration 
-! 
+! Masked output running index mid has been introduced as a local variable to
+! avoid runtime error (Loop variable has been modified) in time_integration
+!
 ! 4023 2019-06-12 13:20:01Z maronga
 ! Renamed "coupling start time" to "spinup time"
-! 
+!
 ! 4017 2019-06-06 12:16:46Z schwenkel
 ! unused variable removed
-! 
+!
 ! 3655 2019-01-07 16:51:22Z knoop
 ! Implementation of the PALM module interface
@@ -82,5 +86,5 @@
 ! ------------
 !> Writing a header with all important information about the current run.
-!> This subroutine is called three times, two times at the beginning 
+!> This subroutine is called three times, two times at the beginning
 !> (writing information on files RUN_CONTROL and HEADER) and one time at the
 !> end of the run, then writing additional information about CPU-usage on file
@@ -88,5 +92,5 @@
 !-----------------------------------------------------------------------------!
  SUBROUTINE header
- 
+
 
     USE arrays_3d,                                                             &
@@ -148,22 +152,22 @@
     IMPLICIT NONE
 
-    
+
     CHARACTER (LEN=2)  ::  do2d_mode           !< mode of 2D data output (xy, xz, yz)
-    
+
     CHARACTER (LEN=5)  ::  section_chr         !< string indicating grid information where to output 2D slices
-    
+
     CHARACTER (LEN=10) ::  host_chr            !< string for hostname
-    
+
     CHARACTER (LEN=16) ::  begin_chr           !< string indication start time for the data output
     CHARACTER (LEN=16) ::  coor_chr            !< dummy string
-    
+
     CHARACTER (LEN=26) ::  ver_rev             !< string for run identification
 
 #if defined( __parallel )
-    CHARACTER (LEN=32) ::  cpl_name            !< name of child domain (nesting mode only) 
+    CHARACTER (LEN=32) ::  cpl_name            !< name of child domain (nesting mode only)
 #endif
-    
+
     CHARACTER (LEN=40) ::  output_format       !< netcdf format
-        
+
     CHARACTER (LEN=70) ::  char1               !< dummy varialbe used for various strings
     CHARACTER (LEN=70) ::  char2               !< string containing informating about the advected distance in case of Galilei transformation
@@ -175,9 +179,9 @@
     CHARACTER (LEN=70) ::  do3d_chr            !< string indicating 3D output variables
     CHARACTER (LEN=70) ::  domask_chr          !< string indicating masked output variables
-    CHARACTER (LEN=70) ::  run_classification  !< string classifying type of run, e.g. nested, coupled, etc. 
-    
+    CHARACTER (LEN=70) ::  run_classification  !< string classifying type of run, e.g. nested, coupled, etc.
+
     CHARACTER (LEN=85) ::  r_upper             !< string indicating model top boundary condition for various quantities
     CHARACTER (LEN=85) ::  r_lower             !< string indicating bottom boundary condition for various quantities
-    
+
     CHARACTER (LEN=86) ::  coordinates         !< string indicating height coordinates for profile-prescribed variables
     CHARACTER (LEN=86) ::  gradients           !< string indicating gradients of profile-prescribed variables between the prescribed height coordinates
@@ -218,5 +222,5 @@
     INTEGER(iwp) ::  npe_total      !< number of total PEs in a coupler (parent + child)
 #endif
-    
+
 
     REAL(wp) ::  cpuseconds_per_simulated_second  !< CPU time (in s) per simulated second
@@ -506,11 +510,11 @@
     IF ( .NOT. ocean_mode )  THEN
        WRITE ( io, 250 )  dx, dy
-       
+
        DO i = 1, number_stretch_level_start+1
           WRITE ( io, 253 )  i, dz(i)
        ENDDO
-       
+
        WRITE( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
-       
+
        IF ( ANY( dz_stretch_level_start_index < nzt+1 ) )  THEN
           WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
@@ -535,5 +539,5 @@
           ENDDO
        ENDIF
-       
+
     ELSE
        WRITE ( io, 250 )  dx, dy
@@ -541,7 +545,7 @@
           WRITE ( io, 253 )  i, dz(i)
        ENDDO
-       
+
        WRITE ( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(0)
-       
+
        IF ( ANY( dz_stretch_level_start_index > 0 ) )  THEN
           WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
@@ -603,5 +607,5 @@
        ENDDO
 
- 
+
        IF ( .NOT. large_scale_forcing )  THEN
           WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
@@ -620,5 +624,5 @@
     i = 1
     DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
-      
+
        WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
        ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
@@ -672,5 +676,5 @@
           i = i + 1
        ENDIF
- 
+
     ENDDO
 
@@ -782,5 +786,5 @@
 !-- Complex terrain
     IF ( complex_terrain )  THEN
-       WRITE( io, 280 ) 
+       WRITE( io, 280 )
        IF ( turbulent_inflow )  THEN
           WRITE( io, 281 )  zu(topo_top_ind(0,0,0))
@@ -846,5 +850,5 @@
        r_upper = 'e(nzt+1) = e(nzt) = e(nzt-1)'
 
-       WRITE ( io, 301 )  'e', r_lower, r_upper       
+       WRITE ( io, 301 )  'e', r_lower, r_upper
 
     ENDIF
@@ -953,5 +957,5 @@
     WRITE ( io, 317 )  bc_lr, bc_ns
     IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
-       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
+       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor
        IF ( turbulent_inflow )  THEN
           IF ( y_shift == 0 ) THEN
@@ -1009,5 +1013,5 @@
                           TRIM( gradients ), TRIM( slices )
     ELSE
-       WRITE ( io, 428 ) 
+       WRITE ( io, 428 )
     ENDIF
 
@@ -1022,5 +1026,5 @@
        i = 1
        DO  WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
-          
+
           WRITE (coor_chr,'(E8.1,4X)')  q_init(q_vertical_gradient_level_ind(i))
           temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
@@ -1028,8 +1032,8 @@
           WRITE (coor_chr,'(E8.1,4X)')  q_vertical_gradient(i)
           gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
-          
+
           WRITE (coor_chr,'(I8,4X)')  q_vertical_gradient_level_ind(i)
           slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
-          
+
           WRITE (coor_chr,'(F8.1,4X)')  q_vertical_gradient_level(i)
           coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
@@ -1058,5 +1062,5 @@
        i = 1
        DO  WHILE ( s_vertical_gradient_level_ind(i) /= -9999 )
-          
+
           WRITE (coor_chr,'(E8.1,4X)')  s_init(s_vertical_gradient_level_ind(i))
           temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
@@ -1064,8 +1068,8 @@
           WRITE (coor_chr,'(E8.1,4X)')  s_vertical_gradient(i)
           gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
-          
+
           WRITE (coor_chr,'(I8,4X)')  s_vertical_gradient_level_ind(i)
           slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
-          
+
           WRITE (coor_chr,'(F8.1,4X)')  s_vertical_gradient_level(i)
           coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
@@ -1081,5 +1085,5 @@
        WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ),           &
                           TRIM( gradients ), TRIM( slices )
-    ENDIF    
+    ENDIF
 
 !
@@ -1586,8 +1590,8 @@
     ENDIF
     IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
-       WRITE ( io, 477 )  q_surface_initial_change       
+       WRITE ( io, 477 )  q_surface_initial_change
     ENDIF
     IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
-       WRITE ( io, 478 )  q_surface_initial_change       
+       WRITE ( io, 478 )  q_surface_initial_change
     ENDIF
 
@@ -1727,5 +1731,5 @@
               ' Tunnel width: ', F6.2 )
 274 FORMAT (  ' Tunnel in ', A, ' direction.' / &
-              ' Tunnel height: ', F6.2, / &    
+              ' Tunnel height: ', F6.2, / &
               ' Tunnel-wall depth: ', F6.2      / &
               ' Tunnel width: ', F6.2, / &
@@ -1756,5 +1760,5 @@
              '       z_mo = ',F6.2,' m   z0 =',F7.4,' m   z0h =',F8.5,&
              ' m   kappa =',F5.2/ &
-             '       Rif value range:   ',F8.2,' <= rif <=',F6.2)
+             '       zeta value range:   ',F8.2,' <= zeta <=',F6.2)
 306 FORMAT ('       Predefined constant heatflux:   ',F9.6,' K m/s')
 307 FORMAT ('       Heatflux has a random normal distribution')
@@ -1762,5 +1766,5 @@
 309 FORMAT ('       Predefined constant salinityflux:   ',F9.6,' psu m/s')
 310 FORMAT (//'    1D-Model:'// &
-             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
+             '       Ri value range:   ',F6.2,' <= Ri <=',F6.2)
 311 FORMAT ('       Predefined constant humidity flux: ',E10.3,' kg/kg m/s')
 312 FORMAT ('       Predefined surface humidity')
