Changeset 411 for palm/trunk/SOURCE

Timestamp:

Dec 11, 2009 2:15:58 PM (15 years ago)

Author:

heinze

Message:

Large scale vertical motion (subsidence/ascent) can be applied to the prognostic equation for the potential temperature

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• Makefile (modified) (6 diffs)
• check_parameters.f90 (modified) (3 diffs)
• header.f90 (modified) (7 diffs)
• modules.f90 (modified) (5 diffs)
• parin.f90 (modified) (2 diffs)
• prognostic_equations.f90 (modified) (5 diffs)
• read_var_list.f90 (modified) (3 diffs)
• subsidence.f90 (added)
• write_var_list.f90 (modified) (3 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -5 +5 @@
 (check_parameters, init_3d_model, parin)
 
-check parameters, init_3d_model, parin
+Large scale vertical motion (subsidence/ascent) can be applied to the 
+prognostic equation for the potential temperature. (check_parameters, header, 
+Makefile, modules, parin, prognostic_equations, read_var_list, subsidence,
+write_var_list)
+
+check parameters, header, init_3d_model, Makefile, modules, parin prognostic equations, read_var_list, subsidence, write_var_list
 
 Changed:

palm/trunk/SOURCE/Makefile

@@ -5 +5 @@
 # -----------------
 # +set_particle_attributes, set_slicer_attributes_dvrp
+# +subsidence
 #
 # Branch revisions:
@@ -81 +82 @@
         set_particle_attributes.f90 set_slicer_attributes_dvrp.f90 \
         singleton.f90 sor.f90 spline_x.f90 \
-        spline_y.f90 spline_z.f90 sum_up_3d_data.f90 surface_coupler.f90 \
+        spline_y.f90 spline_z.f90 subsidence.f90 \
+        sum_up_3d_data.f90 surface_coupler.f90 \
         swap_timelevel.f90 temperton_fft.f90 time_integration.f90 \
         time_to_string.f90 timestep.f90 timestep_scheme_steering.f90 \
@@ -124 +126 @@
         set_particle_attributes.o set_slicer_attributes_dvrp.o \
         singleton.o sor.o spline_x.o spline_y.o \
-        spline_z.o sum_up_3d_data.o surface_coupler.o swap_timelevel.o \
+        spline_z.o subsidence.o sum_up_3d_data.o surface_coupler.o \
+        swap_timelevel.o \
         temperton_fft.o time_integration.o time_to_string.o \
         timestep.o timestep_scheme_steering.o transpose.o \
@@ -188 +191 @@
 check_for_restart.o: modules.o
 check_open.o: modules.o
-check_parameters.o: modules.o
+check_parameters.o: modules.o subsidence.o
 close_file.o: modules.o
 compute_vpt.o: modules.o
@@ -260 +263 @@
         eqn_state_seawater.o impact_of_latent_heat.o \
         plant_canopy_model.o production_e.o \
-        user_actions.o
+        subsidence.o user_actions.o
 random_gauss.o: random_function.o
 read_3d_binary.o: modules.o random_function.o
@@ -272 +275 @@
 spline_y.o: modules.o
 spline_z.o: modules.o
+subsidence.o: modules.o
 sum_up_3d_data.o: modules.o
 surface_coupler.o: modules.o

palm/trunk/SOURCE/check_parameters.f90

@@ -5 +5 @@
 ! -----------------
 ! Enabled passive scalar/humidity wall fluxes for non-flat topography
+! Initialization of large scale vertical motion (subsidence/ascent)
 !
 ! Branch revisions:
@@ -114 +115 @@
     USE pegrid
     USE profil_parameter
+    USE subsidence_mod
     USE statistics
     USE transpose_indices
@@ -983 +985 @@
          
     ENDIF
+
+!
+!-- Initialize large scale subsidence if required
+    IF ( ws_vertical_gradient_level(1) /= -9999999.9 )  THEN
+       large_scale_subsidence = .TRUE.
+       CALL init_w_subsidence
+    END IF
+ 
              
 !

palm/trunk/SOURCE/header.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! +large scale vertical motion (subsidence/ascent)
 ! 
 !
@@ -84 +85 @@
 ! end of the run, then writing additional information about CPU-usage on file
 ! header.
-!------------------------------------------------------------------------------!
+!-----------------------------------------------------------------------------!
 
     USE arrays_3d
@@ -96 +97 @@
     USE particle_attributes
     USE pegrid
+    USE subsidence_mod
     USE spectrum
 
@@ -323 +325 @@
           WRITE ( io, 152 )  dpdxy, dp_level_b, '.'
        ENDIF
+    ENDIF
+    IF ( large_scale_subsidence )  THEN
+        WRITE ( io, 153 )
+        WRITE ( io, 154 )
     ENDIF
     WRITE ( io, 99 )
@@ -1357 +1363 @@
 
        WRITE ( io, 425 )  TRIM( coordinates ), TRIM( temperatures ), &
+                          TRIM( gradients ), TRIM( slices )
+    ENDIF
+
+!
+!-- Profile for the large scale vertial velocity
+!-- Building output strings, starting with surface value
+    IF ( large_scale_subsidence )  THEN
+       temperatures = '   0.0'
+       gradients = '------'
+       slices = '     0'
+       coordinates = '   0.0'
+       i = 1
+       DO  WHILE ( ws_vertical_gradient_level_ind(i) /= -9999 )
+
+          WRITE (coor_chr,'(E10.2,7X)')  &
+                                w_subs(ws_vertical_gradient_level_ind(i))
+          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(E10.2,7X)')  ws_vertical_gradient(i)
+          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(I10,7X)')  ws_vertical_gradient_level_ind(i)
+          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(F10.2,7X)')  ws_vertical_gradient_level(i)
+          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
+
+          i = i + 1
+       ENDDO
+
+       WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
                           TRIM( gradients ), TRIM( slices )
     ENDIF
@@ -1568 +1605 @@
            /'     ',2(1X,E12.5),'Pa/m in x/y direction', &
            /'     starting from dp_level_b =', F8.3, 'm', A /)
+153 FORMAT (' --> Large-scale vertical motion is used in the ', &
+                  'prognostic equation for')
+154 FORMAT ('     the potential temperature')
 200 FORMAT (//' Run time and time step information:'/ &
              ' ----------------------------------'/)
@@ -1814 +1854 @@
             '       Gradient:   ',A,'  psu/100m'/ &
             '       Gridpoint:  ',A)
+426 FORMAT (/'    Characteristic levels of the subsidence/ascent profile:'// &
+            '       Height:      ',A,'  m'/ &
+            '       w_subs:      ',A,'  m/s'/ &
+            '       Gradient:    ',A,'  (m/s)/100m'/ &
+            '       Gridpoint:   ',A)
 450 FORMAT (//' LES / Turbulence quantities:'/ &
               ' ---------------------------'/)

palm/trunk/SOURCE/modules.f90

@@ -8 +8 @@
 ! initial data assignments to some dvrp arrays changed due to error messages
 ! from gfortran compiler
+! +large_scale_subsidence, ws_vertical_gradient, ws_vertical_gradient_level,
+! ws_vertical_gradient_level_ind, w_subs
 !
 ! Branch revisions:
@@ -155 +157 @@
           ddzu, dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, km_damp_x,  &
           km_damp_y, lad, l_grid, pt_init, q_init, rdf, sa_init, ug, u_init,   &
-          u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, zu, zw
+          u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
 
     REAL, DIMENSION(:,:), ALLOCATABLE ::                                       &
@@ -398 +400 @@
                 section_xz(100) = -9999, section_yz(100) = -9999, &
                 ug_vertical_gradient_level_ind(10) = -9999, &
-                vg_vertical_gradient_level_ind(10) = -9999
+                vg_vertical_gradient_level_ind(10) = -9999, &
+                ws_vertical_gradient_level_ind(10) = -9999
 
     INTEGER, DIMENSION(:), ALLOCATABLE ::  grid_level_count
@@ -426 +429 @@
                 humidity = .FALSE., humidity_remote = .FALSE., &
                 inflow_l = .FALSE., inflow_n = .FALSE., &
-                inflow_r = .FALSE., inflow_s = .FALSE., iso2d_output = .FALSE.,&
+                inflow_r = .FALSE., inflow_s = .FALSE., &
+                iso2d_output = .FALSE., large_scale_subsidence = .FALSE., &
                 mg_switch_to_pe0 = .FALSE., &
                 netcdf_output = .FALSE., netcdf_64bit = .FALSE., &
@@ -545 +549 @@
              volume_flow_initial(1:2) = 0.0, wall_heatflux(0:4) = 0.0, &
              wall_humidityflux(0:4) = 0.0, wall_qflux(0:4) = 0.0, &
-             wall_salinityflux(0:4) = 0.0, wall_scalarflux(0:4) = 0.0
+             wall_salinityflux(0:4) = 0.0, wall_scalarflux(0:4) = 0.0, &
+             ws_vertical_gradient(10) = 0.0, &
+             ws_vertical_gradient_level(10) = -9999999.9
 
     REAL, DIMENSION(:), ALLOCATABLE ::  dp_smooth_factor

palm/trunk/SOURCE/parin.f90

@@ -5 +5 @@
 ! -----------------
 ! +wall_humidityflux, wall_scalarflux
+! +ws_vertical_gradient, ws_vertical_gradient_level
 !
 ! Branch revisions:
@@ -137 +138 @@
              use_upstream_for_tke, vg_surface, vg_vertical_gradient, &
              vg_vertical_gradient_level, v_bulk, wall_adjustment, &
-             wall_heatflux, wall_humidityflux, wall_scalarflux
+             wall_heatflux, wall_humidityflux, wall_scalarflux, &
+             ws_vertical_gradient, ws_vertical_gradient_level
 
 

palm/trunk/SOURCE/prognostic_equations.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! add call of subsidence in the prognostic equation for
+! the potential temperature
 ! 
 !
@@ -94 +96 @@
     USE plant_canopy_model_mod
     USE production_e_mod
+    USE subsidence_mod
     USE user_actions_mod
 
@@ -452 +455 @@
           IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
              CALL plant_canopy_model( i, j, 4 )
+          ENDIF
+
+!
+!--       If required compute influence of large-scale subsidence/ascent
+          IF ( large_scale_subsidence ) THEN
+             CALL subsidence ( i, j, tend, pt, pt_init )
           ENDIF
 
@@ -1094 +1103 @@
           ENDIF
 
+
+!--       If required compute influence of large-scale subsidence/ascent
+          IF ( large_scale_subsidence ) THEN
+             CALL subsidence ( i, j, tend, pt, pt_init )
+          ENDIF
+
+
           CALL user_actions( i, j, 'pt-tendency' )
 
@@ -1680 +1696 @@
     ENDIF
 
+!--If required compute influence of large-scale subsidence/ascent
+   IF ( large_scale_subsidence ) THEN
+      CALL subsidence ( tend, pt, pt_init )
+   ENDIF
 
     CALL user_actions( 'pt-tendency' )

palm/trunk/SOURCE/read_var_list.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! +large_scale_subsidence, ws_vertical_gradient, ws_vertical_gradient_level,
+! ws_vertical_gradient_level_ind
 !
 ! Former revisions:
@@ -346 +347 @@
           CASE ( 'last_dt_change' )
              READ ( 13 )  last_dt_change
+           CASE ( 'large_scale_subsidence' )
+             READ ( 13 )  large_scale_subsidence
           CASE ( 'leaf_surface_concentration' )
              READ ( 13 )  leaf_surface_concentration
@@ -593 +596 @@
           CASE ( 'wall_adjustment' )
              READ ( 13 )  wall_adjustment
+          CASE ( 'ws_vertical_gradient' )
+             READ ( 13 )  ws_vertical_gradient
+          CASE ( 'ws_vertical_gradient_level' )
+             READ ( 13 )  ws_vertical_gradient_level
+          CASE ( 'ws_vertical_gradient_level_ind' )
+             READ ( 13 )  ws_vertical_gradient_level_ind
           CASE ( 'w_max' )
              READ ( 13 )  w_max

palm/trunk/SOURCE/write_var_list.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! +large_scale_subsidence, ws_vertical_gradient, ws_vertical_gradient_level,
+! ws_vertical_gradient_level_ind
 !
 ! Former revisions:
@@ -267 +268 @@
     WRITE ( 14 )  'last_dt_change                '
     WRITE ( 14 )  last_dt_change
+    WRITE ( 14 )  'large_scale_subsidence        '
+    WRITE ( 14 )  large_scale_subsidence
     WRITE ( 14 )  'leaf_surface_concentration    '
     WRITE ( 14 )  leaf_surface_concentration
@@ -511 +514 @@
     WRITE ( 14 )  'wall_adjustment               '
     WRITE ( 14 )  wall_adjustment
+    WRITE ( 14 )  'ws_vertical_gradient          '
+    WRITE ( 14 )  ws_vertical_gradient
+    WRITE ( 14 )  'ws_vertical_gradient_level    '
+    WRITE ( 14 )  ws_vertical_gradient_level
+    WRITE ( 14 )  'ws_vertical_gradient_level_ind'
+    WRITE ( 14 )  ws_vertical_gradient_level_ind
     WRITE ( 14 )  'w_max                         '
     WRITE ( 14 )  w_max