Changeset 57 for palm/trunk/SOURCE/production_e.f90

Timestamp:

Mar 9, 2007 12:05:41 PM (17 years ago)

Author:

raasch

Message:

preliminary update of further changes, advec_particles is not running!

File:

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

palm/trunk/SOURCE/production_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Wall functions now include diabatic conditions, call of routine wall_fluxes_e
+! Wall functions now include diabatic conditions, call of routine wall_fluxes_e,
+! reference temperature pt_reference can be used in buoyancy term
 !
 ! Former revisions:
@@ -359 +360 @@
           IF ( .NOT. moisture )  THEN
 
-             DO  j = nys, nyn
-
-                DO  k = nzb_diff_s_inner(j,i), nzt_diff
-                   tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / pt(k,j,i) * &
-                                    ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
+             IF ( use_pt_reference )  THEN
+
+                DO  j = nys, nyn
+                   DO  k = nzb_diff_s_inner(j,i), nzt_diff
+                      tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g/pt_reference * &
+                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
+                   ENDDO
+
+                   IF ( use_surface_fluxes )  THEN
+                      k = nzb_diff_s_inner(j,i)-1
+                      tend(k,j,i) = tend(k,j,i) + g / pt_reference * shf(j,i)
+                   ENDIF
+
+                   IF ( use_top_fluxes )  THEN
+                      k = nzt
+                      tend(k,j,i) = tend(k,j,i) + g / pt_reference * tswst(j,i)
+                   ENDIF
                 ENDDO
 
-                IF ( use_surface_fluxes )  THEN
-                   k = nzb_diff_s_inner(j,i)-1
-                   tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * shf(j,i)
-                ENDIF
-
-                IF ( use_top_fluxes )  THEN
-                   k = nzt
-                   tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * tswst(j,i)
-                ENDIF
-
-             ENDDO
+             ELSE
+
+                DO  j = nys, nyn
+                   DO  k = nzb_diff_s_inner(j,i), nzt_diff
+                      tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / pt(k,j,i) * &
+                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
+                   ENDDO
+
+                   IF ( use_surface_fluxes )  THEN
+                      k = nzb_diff_s_inner(j,i)-1
+                      tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * shf(j,i)
+                   ENDIF
+
+                   IF ( use_top_fluxes )  THEN
+                      k = nzt
+                      tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * tswst(j,i)
+                   ENDIF
+                ENDDO
+
+             ENDIF
 
           ELSE
@@ -736 +758 @@
        IF ( .NOT. moisture )  THEN
 
-          DO  k = nzb_diff_s_inner(j,i), nzt_diff
-             tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / pt(k,j,i) * &
-                                    ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
-          ENDDO
-
-          IF ( use_surface_fluxes )  THEN
-             k = nzb_diff_s_inner(j,i)-1
-             tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * shf(j,i)
-          ENDIF
-
-          IF ( use_top_fluxes )  THEN
-             k = nzt
-             tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * tswst(j,i)
-          ENDIF
+          IF ( use_pt_reference )  THEN
+
+             DO  k = nzb_diff_s_inner(j,i), nzt_diff
+                tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / pt_reference * &
+                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
+             ENDDO
+
+             IF ( use_surface_fluxes )  THEN
+                k = nzb_diff_s_inner(j,i)-1
+                tend(k,j,i) = tend(k,j,i) + g / pt_reference * shf(j,i)
+             ENDIF
+
+             IF ( use_top_fluxes )  THEN
+                k = nzt
+                tend(k,j,i) = tend(k,j,i) + g / pt_reference * tswst(j,i)
+             ENDIF
+
+          ELSE
+
+             DO  k = nzb_diff_s_inner(j,i), nzt_diff
+                tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / pt(k,j,i) * &
+                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
+             ENDDO
+
+             IF ( use_surface_fluxes )  THEN
+                k = nzb_diff_s_inner(j,i)-1
+                tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * shf(j,i)
+             ENDIF
+
+             IF ( use_top_fluxes )  THEN
+                k = nzt
+                tend(k,j,i) = tend(k,j,i) + g / pt(k,j,i) * tswst(j,i)
+             ENDIF
+
+         ENDIF
 
        ELSE