Changeset 2327

Timestamp:

Aug 2, 2017 7:40:57 AM (7 years ago)

Author:

maronga

Message:

temporary bugfix for calculation of air density

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Temporary bugfix
+! 
+! 2320 2017-07-21 12:47:43Z suehring
 ! Modularize large-scale forcing and nudging
 ! 
@@ -343 +346 @@
 !> or
 !> c) read values of a previous run
+!> @todo fix crashes for runs with topography and density /= 1
 !------------------------------------------------------------------------------!
  SUBROUTINE init_3d_model
@@ -733 +737 @@
        rho_air_zw(nzt+1)  = rho_air_zw(nzt)                                    &
                             + 2.0_wp * ( rho_air(nzt+1) - rho_air_zw(nzt)  )
-    ENDIF
+
+!
+!--    Temporary workaround as runs with topography do not work with density /= 1
+       rho_air    = 1.0_wp
+       rho_air_zw = 1.0_wp
+
+    ENDIF
+
+
+
 
 !-- compute the inverse density array in order to avoid expencive divisions