Changeset 4765 for palm/trunk/TUTORIALS/cases/example_neutral/INPUT/example_neutral_p3d

Timestamp:

Nov 2, 2020 10:31:45 AM (4 years ago)

Author:

Giersch

Message:

Tutorials updated for r4761

File:

• palm/trunk/TUTORIALS/cases/example_neutral/INPUT/example_neutral_p3d (modified) (4 diffs)

palm/trunk/TUTORIALS/cases/example_neutral/INPUT/example_neutral_p3d

@@ -17 +17 @@
 !-- initialization
 !-------------------------------------------------------------------------------
-    initializing_actions       = 'set_constant_profiles' , ! initial conditions
+    initializing_actions       = 'set_constant_profiles', 
 
-    ug_surface                 = 5.0 , ! u-comp of geostrophic wind at surface
-    vg_surface                 = 0.0 , ! v-comp of geostrophic wind at surface
+    ug_surface                 = 5.0, 
+    vg_surface                 = 0.0,
 
-    pt_surface                 = 300.0, ! initial surface potential temp
-    pt_vertical_gradient       = 0.0,   ! piecewise temp gradients
-    pt_vertical_gradient_level = 0.0,   ! height level of temp gradients
+    pt_surface                 = 300.0,     
+    pt_vertical_gradient       = 0.0, ! Piecewise temp gradients
+    pt_vertical_gradient_level = 0.0, ! Height level of temp gradients
 !
 !-- boundary conditions
@@ -31 +31 @@
     bc_uv_t                    = 'dirichlet',
 
-    bc_p_b                     = 'neumann', 
-    bc_p_t                     = 'dirichlet',
+    bc_pt_b                    = 'dirichlet',        ! Bottom boundary condition of pt
+    bc_pt_t                    = 'initial_gradient', ! Top boundary condition of pt 
+
+    bc_p_b                     = 'neumann',   ! Bottom boundary condition of p.
+    bc_p_t                     = 'dirichlet', ! Top boundary condition of p.
 
     constant_flux_layer        = .TRUE.,
     
-    roughness_length           = 0.1 ,
+    roughness_length           = 0.1,
 !
 !-- numerics
 !-------------------------------------------------------------------------------
-    fft_method                 = 'temperton-algorithm',   ! build-in fft method     
+    fft_method                 = 'temperton-algorithm', ! Build-in fft method     
 
-    rayleigh_damping_height    = 1250.0,
-    rayleigh_damping_factor    = 0.01,
+    rayleigh_damping_height    = 1250.0, ! Height above which the Rayleigh damping starts
+    rayleigh_damping_factor    = 0.01,   ! Factor for Rayleigh damping.
 
 !
 !-- mode
 !-------------------------------------------------------------------------------
-    neutral                    = .TRUE., ! switch off calculation of temperature equation
-    galilei_transformation     = .TRUE., ! application of a Galilei-transformation
+    neutral                    = .TRUE., ! Switch off calculation of temp equation
+    galilei_transformation     = .TRUE.,      
     
-/ ! end of initialization parameter namelist
+/ ! End of initialization parameter namelist
 
 !-------------------------------------------------------------------------------
@@ -61 +64 @@
 !-- run steering
 !-------------------------------------------------------------------------------
-    end_time                   = 172800.0, ! simulation time of the 3D model
+    end_time                   = 172800.0, ! Simulation time of the 3D model
 
-    create_disturbances        = .TRUE., ! randomly perturbate horiz. velocity
-    dt_disturb                 = 150.0,  ! interval for random perturbations
-    disturbance_energy_limit   = 0.01,   ! upper limit for perturbation energy
+    create_disturbances        = .TRUE., ! Randomly perturbate horiz. velocity
+    dt_disturb                 = 150.0,  ! Interval for random perturbations
+    disturbance_energy_limit   = 0.01,   ! Upper limit for perturbation energy
 !
 !-- data output
 !-------------------------------------------------------------------------------
-    netcdf_data_format         = 2, ! use NetCDF3
+    netcdf_data_format         = 2, ! Data format for netCDF files
     
-    dt_run_control             = 60.0,
+    dt_run_control             = 60.0, ! Temporal interval of run control output
 
-    dt_dopr                    = 10800.0,
-         
-    averaging_interval_pr      = 10800.0, 
-    dt_averaging_input_pr      = 50.0,
-    data_output_pr             = '#u', '#v', 'w', '#theta', 'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv', 'w*2', 'theta*2',
+    dt_dopr                    = 10800.0, ! Temporal interval for vertical profile output         
+    averaging_interval_pr      = 10800.0, ! Averaging interval for vertical profile output 
+    dt_averaging_input_pr      = 50.0,    ! Temporal interval of data that is used for temporal averaging of vertical profiles and/or spectra
+
+    data_output_pr             = '#u', '#v', 'w', '#theta', 
+                                 'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv',
+                                 'w*2', 'theta*2', ! Quantities for vertical profile output
                         
-/! end of runtime parameter namelist
+/ ! End of runtime parameter namelist
 
 !-------------------------------------------------------------------------------
@@ -87 +92 @@
 &spectra_parameters
  
-    data_output_sp            =  'u', 'v', 'w',
-    dt_dosp                   = 10800.0,
-    comp_spectra_level        =  2, 5, 10, 15, 20,
-    spectra_direction         =  'x' , 'y' , 'x',
-    averaging_interval_sp     = 10800.0,
+    dt_dosp                   = 10800.0, ! Temporal interval for spectra output
+    averaging_interval_sp     = 10800.0, ! Averaging interval for vertical profile output
+
+    comp_spectra_level        =  2, 5, 10, 15, 20, ! Vertical level for horizontal spectra output
+    spectra_direction         =  'x' , 'y' , 'x',  ! Direction(s) along which spectra are calculated 
+
+    data_output_sp            =  'u', 'v', 'w',  ! Quantities for spectra output
          
-/! end of spectra parameter namelist
+/ ! End of spectra parameter namelist