Changeset 4877 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Feb 17, 2021 4:17:35 PM (3 years ago)

Author:

suehring

Message:

Bugfix in initialization of vertical surfaces with roughness and surface heat fluxes

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Bugfix in initialization of vertical surfaces with roughness and surface heat fluxes
+!
+! 4848 2021-01-21 15:51:51Z gronemeier
 ! replaced use_syn_turb_gen by syn_turb_gen
 !
@@ -166 +169 @@
 
     USE control_parameters
+
+    USE exchange_horiz_mod,                                                                        &
+        ONLY:  exchange_horiz_2d
 
     USE grid_variables,                                                                            &
@@ -207 +213 @@
                pids_id,                                                                            &
                real_2d,                                                                            &
+               resize_array,                                                                       &
                vars_pids
 
@@ -1508 +1515 @@
 !--       Read variable
           CALL get_variable( pids_id, 'z0', tmp_2d%var, nxl, nxr, nys, nyn )
+          CALL resize_array( tmp_2d%var, nys, nyn, nxl, nxr )
+          CALL exchange_horiz_2d( tmp_2d%var )
 !
 !--       Initialize roughness length. Note, z0 will be only initialized at default-type surfaces.
@@ -1534 +1543 @@
 !--       Read variable
           CALL get_variable( pids_id, 'shf', tmp_2d%var, nxl, nxr, nys, nyn )
+          CALL resize_array( tmp_2d%var, nys, nyn, nxl, nxr )
+          CALL exchange_horiz_2d( tmp_2d%var )
 !
 !--       Initialize heat flux. Note, shf will be only initialized at default-type surfaces. At
@@ -1561 +1572 @@
 !--          Read variable
              CALL get_variable( pids_id, 'qsws', tmp_2d%var, nxl, nxr, nys, nyn )
+             CALL resize_array( tmp_2d%var, nys, nyn, nxl, nxr )
+             CALL exchange_horiz_2d( tmp_2d%var )
 !
 !--          Initialize latent heat flux. Note, qsws will be only initialized at default-type surfaces.
@@ -1590 +1603 @@
 !--          Read variable
              CALL get_variable( pids_id, 'ssws', tmp_2d%var, nxl, nxr, nys, nyn )
+             CALL resize_array( tmp_2d%var, nys, nyn, nxl, nxr )
+             CALL exchange_horiz_2d( tmp_2d%var )
 !
 !--          Initialize passive scalar flux. Initialize horizontal surface elements.