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

Timestamp:

Aug 24, 2007 3:10:38 PM (17 years ago)

Author:

letzel

Message:

• Improved coupler: evaporation - salinity-flux coupling for humidity = .T.,

avoid MPI hangs when coupled runs terminate, add DOC/app/chapter_3.8;

• Optional calculation of km and kh from initial TKE e_init;
• Default initialization of km,kh = 0.00001 for ocean = .T.;
• Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T.;
• Bugfix: Rayleigh damping for ocean fixed.

File:

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

palm/trunk/SOURCE/production_e.f90

@@ -8 +8 @@
 ! u_0 and v_0 are calculated for nxr+1, nyn+1 also (otherwise these values are
 ! not available in case of non-cyclic boundary conditions)
+! Bugfix for ocean density flux at bottom
 !
 ! Former revisions:
@@ -810 +811 @@
 !--             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
+                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 * &
                                       ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)