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

Timestamp:

Jan 20, 2017 3:54:21 PM (8 years ago)

Author:

raasch

Message:

for ocean simulations, density in SGS-TKE-production is replaced by potential density

File:

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

palm/trunk/SOURCE/production_e.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! density in ocean case replaced by potential density
 ! 
 ! Former revisions:
@@ -138 +138 @@
 
        USE arrays_3d,                                                          &
-           ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
+           ONLY:  ddzw, dd2zu, kh, km, prho, pt, q, ql, qsws, qswst, shf,      &
                   tend, tswst, u, v, vpt, w
 
@@ -537 +537 @@
                       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 / rho_reference *   &
-                                          ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
+                            tend(k,j,i) = tend(k,j,i) +                        &
+                                          kh(k,j,i) * g / rho_reference *      &
+                                          ( prho(k+1,j,i) - prho(k-1,j,i) ) *  &
                                           dd2zu(k)
                          ENDDO
@@ -577 +577 @@
                       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 / rho_ocean(k,j,i) *      &
-                                          ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
+                            tend(k,j,i) = tend(k,j,i) +                        &
+                                          kh(k,j,i) * g / prho(k,j,i) *        &
+                                          ( prho(k+1,j,i) - prho(k-1,j,i) ) *  &
                                           dd2zu(k)
                          ENDDO
@@ -744 +744 @@
 
        USE arrays_3d,                                                          &
-           ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
+           ONLY:  ddzw, dd2zu, kh, km, prho, pt, q, ql, qsws, qswst, shf,      &
                   tend, tswst, u, v, vpt, w
 
@@ -1100 +1100 @@
                       tend(k,j,i) = tend(k,j,i) +                   &
                                     kh(k,j,i) * g / rho_reference * &
-                                    ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
+                                    ( prho(k+1,j,i) - prho(k-1,j,i) ) * dd2zu(k)
                    ENDDO
 
@@ -1130 +1130 @@
 !--                the 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 / rho_ocean(k,j,i) * &
-                                    ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
+                      tend(k,j,i) = tend(k,j,i) +                              &
+                                    kh(k,j,i) * g / prho(k,j,i) *              &
+                                    ( prho(k+1,j,i) - prho(k-1,j,i) ) * dd2zu(k)
                    ENDDO