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

Timestamp:

Sep 23, 2009 9:40:33 AM (15 years ago)

Author:

raasch

Message:

in-situ AND potential density are calculated and used in the ocean version

File:

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

palm/trunk/SOURCE/production_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Bugfix: wrong sign in buoyancy production of ocean part in case of not using
+!         the reference density (only in 3D routine production_e)
 ! Bugfix to avoid zero division by km_neutral
 !
@@ -450 +452 @@
                    DO  j = nys, nyn
                       DO  k = nzb_s_inner(j,i)+1, nzt
-                         tend(k,j,i) = tend(k,j,i) +                    &
-                                       kh(k,j,i) * g / prho_reference * &
+                         tend(k,j,i) = tend(k,j,i) +                   &
+                                       kh(k,j,i) * g / rho_reference * &
                                        ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k)
                       ENDDO
@@ -487 +489 @@
                    DO  j = nys, nyn
                       DO  k = nzb_s_inner(j,i)+1, nzt
-                         tend(k,j,i) = tend(k,j,i) -                &
+                         tend(k,j,i) = tend(k,j,i) +                &
                                        kh(k,j,i) * g / rho(k,j,i) * &
                                        ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k)
@@ -946 +948 @@
 !--             bottom and top surface layer
                 DO  k = nzb_s_inner(j,i)+1, nzt
-                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
+                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / rho_reference * &
                                       ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
                 ENDDO