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

Timestamp:

Oct 19, 2007 3:47:46 PM (17 years ago)

Author:

raasch

Message:

further bugfixes for the ocean part

File:

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

palm/trunk/SOURCE/production_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: calculation of density flux in the ocean now starts from nzb+1
 !
 ! Former revisions:
@@ -382 +382 @@
 !--                the bottom and top surface layer
                    DO  j = nys, nyn
-                      DO  k = nzb_s_inner(j,i), nzt
+                      DO  k = nzb_s_inner(j,i)+1, nzt
                          tend(k,j,i) = tend(k,j,i) +                    &
                                        kh(k,j,i) * g / prho_reference * &
@@ -419 +419 @@
 !--                the bottom and top surface layer
                    DO  j = nys, nyn
-                      DO  k = nzb_s_inner(j,i), nzt
+                      DO  k = nzb_s_inner(j,i)+1, nzt
                          tend(k,j,i) = tend(k,j,i) -                &
                                        kh(k,j,i) * g / rho(k,j,i) * &
@@ -814 +814 @@
 !--             So far in the ocean no special treatment of density flux in the
 !--             bottom and top surface layer
-                k= nzb_s_inner(j,i)
-                tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
-                     ( rho(k+1,j,i) - rho(k,j,i) ) * dzu(k+1)
                 DO  k = nzb_s_inner(j,i)+1, nzt
                    tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
@@ -847 +844 @@
 !--             So far in the ocean no special treatment of density flux in the
 !--             bottom and top surface layer
-                DO  k = nzb_s_inner(j,i), nzt
+                DO  k = nzb_s_inner(j,i)+1, nzt
                    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)