Changeset 1174 for palm/trunk/SOURCE

Timestamp:

May 31, 2013 10:28:08 AM (11 years ago)

Author:

gryschka

Message:

Bugfix in computing initial profiles for ug, vg, lad, q in case of Atmosphere

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (6 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! Bugfix in computing initial profiles for ug, vg, lad, q in case of Atmosphere
 !
 ! Former revisions:
@@ -955 +955 @@
                    ug(k) = ug(k-1) + dzu(k) * gradient
                 ELSE
-                   ug(k) = ug_surface + 0.5 * dzu(k) * gradient
+                   ug(k) = ug_surface + dzu(k) * gradient
                 ENDIF
              ELSE
@@ -1019 +1019 @@
                    vg(k) = vg(k-1) + dzu(k) * gradient
                 ELSE
-                   vg(k) = vg_surface + 0.5 * dzu(k) * gradient
+                   vg(k) = vg_surface + dzu(k) * gradient
                 ENDIF
              ELSE
@@ -1133 +1133 @@
                       pt_init(k) = pt_init(k-1) + dzu(k) * gradient
                    ELSE
-                      pt_init(k) = pt_surface   + 0.5 * dzu(k) * gradient
+                      pt_init(k) = pt_surface   + dzu(k) * gradient
                    ENDIF
                 ELSE
@@ -1213 +1213 @@
                    q_init(k) = q_init(k-1) + dzu(k) * gradient
                 ELSE
-                   q_init(k) = q_init(k-1) + 0.5 * dzu(k) * gradient
+                   q_init(k) = q_init(k-1) + dzu(k) * gradient
                 ENDIF
              ELSE
@@ -1295 +1295 @@
                       lad(k) = lad(k-1) + dzu(k) * gradient
                    ELSE
-                      lad(k) = lad_surface + 0.5 * dzu(k) *gradient
+                      lad(k) = lad_surface + dzu(k) *gradient
                    ENDIF
                 ELSE