Changeset 2478 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Sep 18, 2017 1:37:24 PM (7 years ago)

Author:

suehring

Message:

Bugfixes concerning top fluxes and TKE production

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfixes in initializing model top
+!
+! 2378 2017-08-31 13:57:27Z suehring
 ! Bugfix in write restart data
 ! 
@@ -1492 +1495 @@
 !--          Initialize top heat flux
              IF ( constant_top_heatflux )                                      &
-                surf%shf = top_heatflux * heatflux_input_conversion(nzt+1)
+                surf%shf(num_h) = top_heatflux * heatflux_input_conversion(nzt+1)
 !
 !--          Initialization in case of a coupled model run
              IF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
-                surf%shf = 0.0_wp
+                surf%shf(num_h) = 0.0_wp
              ENDIF
 !
 !--          Prescribe latent heat flux at the top      
              IF ( humidity )  THEN
-             surf%qsws = 0.0_wp
+             surf%qsws(num_h) = 0.0_wp
                 IF ( cloud_physics  .AND.  microphysics_morrison ) THEN
-                   surf%ncsws = 0.0_wp
-                   surf%qcsws = 0.0_wp
+                   surf%ncsws(num_h) = 0.0_wp
+                   surf%qcsws(num_h) = 0.0_wp
                 ENDIF
                 IF ( cloud_physics  .AND.  microphysics_seifert ) THEN
-                   surf%nrsws = 0.0_wp
-                   surf%qrsws = 0.0_wp
+                   surf%nrsws(num_h) = 0.0_wp
+                   surf%qrsws(num_h) = 0.0_wp
                 ENDIF
              ENDIF
@@ -1514 +1517 @@
 !--          Prescribe top scalar flux
              IF ( passive_scalar .AND. constant_top_scalarflux )               &
-                surf%ssws = top_scalarflux
+                surf%ssws(num_h) = top_scalarflux
 !
 !--          Prescribe top salinity flux
-             IF ( ocean .AND. constant_top_salinityflux)                          &
-                surf%sasws = top_salinityflux
+             IF ( ocean .AND. constant_top_salinityflux)                       &
+                surf%sasws(num_h) = top_salinityflux
 !
 !--          Initialization in case of a coupled model run
              IF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
-                surf%shf = 0.0_wp
+                surf%shf(num_h) = 0.0_wp
              ENDIF
 !
 !--          Top momentum fluxes
-             surf%usws = top_momentumflux_u * momentumflux_input_conversion(nzt+1)
-             surf%vsws = top_momentumflux_v * momentumflux_input_conversion(nzt+1)
+             surf%usws(num_h) = top_momentumflux_u *                           &
+                                            momentumflux_input_conversion(nzt+1)
+             surf%vsws(num_h) = top_momentumflux_v *                           &
+                                            momentumflux_input_conversion(nzt+1)
 !
 !--          Increment surface indices