Ignore:
Timestamp:
Jul 1, 2020 4:16:43 PM (8 months ago)
Author:
gronemeier
Message:

renamed Richardson flux number into gradient Richardson number (model_1d_mod.f90, turbulence_closure_mod.f90, header.f90, surface_mod.f90) and zeta (header.f90);
do not add whitespaces at current-revision section (document_changes)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/header.f90

    r4573 r4586  
    2020! Current revisions:
    2121! -----------------
    22 ! 
    23 ! 
     22!
     23!
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! Renamed rif to Ri (gradient Richardson number, 1D model)
     28! and zeta (= z_mo / ol, stability parameter, 3D model)
     29!
     30! 4573 2020-06-24 13:08:47Z oliver.maas
    2731! added statement for pt_surface_heating_rate
    28 ! 
     32!
    2933! 4536 2020-05-17 17:24:13Z raasch
    3034! output of restart data format added
    31 ! 
     35!
    3236! 4473 2020-03-25 21:04:07Z gronemeier
    3337! revised message if wall_adjustment is used
    34 ! 
     38!
    3539! 4444 2020-03-05 15:59:50Z raasch
    3640! bugfix: cpp-directives for serial mode added
    37 ! 
     41!
    3842! 4360 2020-01-07 11:25:50Z suehring
    3943! Bugfix, character length too short, caused crash on NEC.
    40 ! 
     44!
    4145! 4309 2019-11-26 18:49:59Z suehring
    4246! replaced recycling_yshift by y_shift
    43 ! 
     47!
    4448! 4301 2019-11-22 12:09:09Z oliver.maas
    4549!
    4650! 4297 2019-11-21 10:37:50Z oliver.maas
    4751! Adjusted format for simulated time and related quantities
    48 ! 
     52!
    4953! 4297 2019-11-21 10:37:50Z oliver.maas
    5054! adjusted message to the changed parameter recycling_yshift
    51 ! 
     55!
    5256! 4227 2019-09-10 18:04:34Z gronemeier
    5357! implement new palm_date_time_mod
    54 ! 
     58!
    5559! 4223 2019-09-10 09:20:47Z gronemeier
    5660! Write information about rotation angle
    57 ! 
     61!
    5862! 4182 2019-08-22 15:20:23Z scharf
    5963! Corrected "Former revisions" section
    60 ! 
     64!
    6165! 4168 2019-08-16 13:50:17Z suehring
    6266! Replace function get_topography_top_index by topo_top_ind
    63 ! 
     67!
    6468! 4069 2019-07-01 14:05:51Z Giersch
    65 ! Masked output running index mid has been introduced as a local variable to 
    66 ! avoid runtime error (Loop variable has been modified) in time_integration 
    67 ! 
     69! Masked output running index mid has been introduced as a local variable to
     70! avoid runtime error (Loop variable has been modified) in time_integration
     71!
    6872! 4023 2019-06-12 13:20:01Z maronga
    6973! Renamed "coupling start time" to "spinup time"
    70 ! 
     74!
    7175! 4017 2019-06-06 12:16:46Z schwenkel
    7276! unused variable removed
    73 ! 
     77!
    7478! 3655 2019-01-07 16:51:22Z knoop
    7579! Implementation of the PALM module interface
     
    8286! ------------
    8387!> Writing a header with all important information about the current run.
    84 !> This subroutine is called three times, two times at the beginning 
     88!> This subroutine is called three times, two times at the beginning
    8589!> (writing information on files RUN_CONTROL and HEADER) and one time at the
    8690!> end of the run, then writing additional information about CPU-usage on file
     
    8892!-----------------------------------------------------------------------------!
    8993 SUBROUTINE header
    90  
     94
    9195
    9296    USE arrays_3d,                                                             &
     
    148152    IMPLICIT NONE
    149153
    150    
     154
    151155    CHARACTER (LEN=2)  ::  do2d_mode           !< mode of 2D data output (xy, xz, yz)
    152    
     156
    153157    CHARACTER (LEN=5)  ::  section_chr         !< string indicating grid information where to output 2D slices
    154    
     158
    155159    CHARACTER (LEN=10) ::  host_chr            !< string for hostname
    156    
     160
    157161    CHARACTER (LEN=16) ::  begin_chr           !< string indication start time for the data output
    158162    CHARACTER (LEN=16) ::  coor_chr            !< dummy string
    159    
     163
    160164    CHARACTER (LEN=26) ::  ver_rev             !< string for run identification
    161165
    162166#if defined( __parallel )
    163     CHARACTER (LEN=32) ::  cpl_name            !< name of child domain (nesting mode only) 
     167    CHARACTER (LEN=32) ::  cpl_name            !< name of child domain (nesting mode only)
    164168#endif
    165    
     169
    166170    CHARACTER (LEN=40) ::  output_format       !< netcdf format
    167        
     171
    168172    CHARACTER (LEN=70) ::  char1               !< dummy varialbe used for various strings
    169173    CHARACTER (LEN=70) ::  char2               !< string containing informating about the advected distance in case of Galilei transformation
     
    175179    CHARACTER (LEN=70) ::  do3d_chr            !< string indicating 3D output variables
    176180    CHARACTER (LEN=70) ::  domask_chr          !< string indicating masked output variables
    177     CHARACTER (LEN=70) ::  run_classification  !< string classifying type of run, e.g. nested, coupled, etc. 
    178    
     181    CHARACTER (LEN=70) ::  run_classification  !< string classifying type of run, e.g. nested, coupled, etc.
     182
    179183    CHARACTER (LEN=85) ::  r_upper             !< string indicating model top boundary condition for various quantities
    180184    CHARACTER (LEN=85) ::  r_lower             !< string indicating bottom boundary condition for various quantities
    181    
     185
    182186    CHARACTER (LEN=86) ::  coordinates         !< string indicating height coordinates for profile-prescribed variables
    183187    CHARACTER (LEN=86) ::  gradients           !< string indicating gradients of profile-prescribed variables between the prescribed height coordinates
     
    218222    INTEGER(iwp) ::  npe_total      !< number of total PEs in a coupler (parent + child)
    219223#endif
    220    
     224
    221225
    222226    REAL(wp) ::  cpuseconds_per_simulated_second  !< CPU time (in s) per simulated second
     
    506510    IF ( .NOT. ocean_mode )  THEN
    507511       WRITE ( io, 250 )  dx, dy
    508        
     512
    509513       DO i = 1, number_stretch_level_start+1
    510514          WRITE ( io, 253 )  i, dz(i)
    511515       ENDDO
    512        
     516
    513517       WRITE( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
    514        
     518
    515519       IF ( ANY( dz_stretch_level_start_index < nzt+1 ) )  THEN
    516520          WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
     
    535539          ENDDO
    536540       ENDIF
    537        
     541
    538542    ELSE
    539543       WRITE ( io, 250 )  dx, dy
     
    541545          WRITE ( io, 253 )  i, dz(i)
    542546       ENDDO
    543        
     547
    544548       WRITE ( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(0)
    545        
     549
    546550       IF ( ANY( dz_stretch_level_start_index > 0 ) )  THEN
    547551          WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
     
    603607       ENDDO
    604608
    605  
     609
    606610       IF ( .NOT. large_scale_forcing )  THEN
    607611          WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
     
    620624    i = 1
    621625    DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
    622      
     626
    623627       WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
    624628       ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
     
    672676          i = i + 1
    673677       ENDIF
    674  
     678
    675679    ENDDO
    676680
     
    782786!-- Complex terrain
    783787    IF ( complex_terrain )  THEN
    784        WRITE( io, 280 ) 
     788       WRITE( io, 280 )
    785789       IF ( turbulent_inflow )  THEN
    786790          WRITE( io, 281 )  zu(topo_top_ind(0,0,0))
     
    846850       r_upper = 'e(nzt+1) = e(nzt) = e(nzt-1)'
    847851
    848        WRITE ( io, 301 )  'e', r_lower, r_upper       
     852       WRITE ( io, 301 )  'e', r_lower, r_upper
    849853
    850854    ENDIF
     
    953957    WRITE ( io, 317 )  bc_lr, bc_ns
    954958    IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
    955        WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
     959       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor
    956960       IF ( turbulent_inflow )  THEN
    957961          IF ( y_shift == 0 ) THEN
     
    10091013                          TRIM( gradients ), TRIM( slices )
    10101014    ELSE
    1011        WRITE ( io, 428 ) 
     1015       WRITE ( io, 428 )
    10121016    ENDIF
    10131017
     
    10221026       i = 1
    10231027       DO  WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
    1024          
     1028
    10251029          WRITE (coor_chr,'(E8.1,4X)')  q_init(q_vertical_gradient_level_ind(i))
    10261030          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
     
    10281032          WRITE (coor_chr,'(E8.1,4X)')  q_vertical_gradient(i)
    10291033          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
    1030          
     1034
    10311035          WRITE (coor_chr,'(I8,4X)')  q_vertical_gradient_level_ind(i)
    10321036          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
    1033          
     1037
    10341038          WRITE (coor_chr,'(F8.1,4X)')  q_vertical_gradient_level(i)
    10351039          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
     
    10581062       i = 1
    10591063       DO  WHILE ( s_vertical_gradient_level_ind(i) /= -9999 )
    1060          
     1064
    10611065          WRITE (coor_chr,'(E8.1,4X)')  s_init(s_vertical_gradient_level_ind(i))
    10621066          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
     
    10641068          WRITE (coor_chr,'(E8.1,4X)')  s_vertical_gradient(i)
    10651069          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
    1066          
     1070
    10671071          WRITE (coor_chr,'(I8,4X)')  s_vertical_gradient_level_ind(i)
    10681072          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
    1069          
     1073
    10701074          WRITE (coor_chr,'(F8.1,4X)')  s_vertical_gradient_level(i)
    10711075          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
     
    10811085       WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ),           &
    10821086                          TRIM( gradients ), TRIM( slices )
    1083     ENDIF   
     1087    ENDIF
    10841088
    10851089!
     
    15861590    ENDIF
    15871591    IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
    1588        WRITE ( io, 477 )  q_surface_initial_change       
     1592       WRITE ( io, 477 )  q_surface_initial_change
    15891593    ENDIF
    15901594    IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
    1591        WRITE ( io, 478 )  q_surface_initial_change       
     1595       WRITE ( io, 478 )  q_surface_initial_change
    15921596    ENDIF
    15931597
     
    17271731              ' Tunnel width: ', F6.2 )
    17281732274 FORMAT (  ' Tunnel in ', A, ' direction.' / &
    1729               ' Tunnel height: ', F6.2, / &   
     1733              ' Tunnel height: ', F6.2, / &
    17301734              ' Tunnel-wall depth: ', F6.2      / &
    17311735              ' Tunnel width: ', F6.2, / &
     
    17561760             '       z_mo = ',F6.2,' m   z0 =',F7.4,' m   z0h =',F8.5,&
    17571761             ' m   kappa =',F5.2/ &
    1758              '       Rif value range:   ',F8.2,' <= rif <=',F6.2)
     1762             '       zeta value range:   ',F8.2,' <= zeta <=',F6.2)
    17591763306 FORMAT ('       Predefined constant heatflux:   ',F9.6,' K m/s')
    17601764307 FORMAT ('       Heatflux has a random normal distribution')
     
    17621766309 FORMAT ('       Predefined constant salinityflux:   ',F9.6,' psu m/s')
    17631767310 FORMAT (//'    1D-Model:'// &
    1764              '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
     1768             '       Ri value range:   ',F6.2,' <= Ri <=',F6.2)
    17651769311 FORMAT ('       Predefined constant humidity flux: ',E10.3,' kg/kg m/s')
    17661770312 FORMAT ('       Predefined surface humidity')
Note: See TracChangeset for help on using the changeset viewer.