Changeset 108 for palm/trunk/SOURCE/init_3d_model.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/init_3d_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -9 +9 @@
 ! Flux initialization in case of coupled runs, +momentum fluxes at top boundary,
 ! +arrays for phase speed c_u, c_v, c_w, indices for u|v|w_m_l|r changed
+! +qswst_remote in case of atmosphere model with humidity coupled to ocean
+! Rayleigh damping for ocean
+! optionally calculate km and kh from initial TKE e_init
 !
 ! Former revisions:
@@ -204 +207 @@
        rho => rho_1  ! routine calc_mean_profile requires density to be a
                      ! pointer
+       IF ( humidity_remote )  THEN
+          ALLOCATE( qswst_remote(nys-1:nyn+1,nxl-1:nxr+1) )
+          qswst_remote = 0.0
+       ENDIF
     ENDIF
 
@@ -507 +514 @@
              km   = km_constant
              kh   = km / prandtl_number
+             e    = 0.0
+          ELSEIF ( e_init > 0.0 )  THEN
+             DO  k = nzb+1, nzt
+                km(k,:,:) = 0.1 * l_grid(k) * SQRT( e_init )
+             ENDDO
+             km(nzb,:,:)   = km(nzb+1,:,:)
+             km(nzt+1,:,:) = km(nzt,:,:)
+             kh   = km / prandtl_number
+             e    = e_init
           ELSE
-             kh   = 0.01   ! there must exist an initial diffusion, because 
-             km   = 0.01   ! otherwise no TKE would be produced by the
-                           ! production terms, as long as not yet
-                           ! e = (u*/cm)**2 at k=nzb+1
-          ENDIF
-          e     = 0.0
+             IF ( .NOT. ocean )  THEN
+                kh   = 0.01   ! there must exist an initial diffusion, because 
+                km   = 0.01   ! otherwise no TKE would be produced by the
+                              ! production terms, as long as not yet
+                              ! e = (u*/cm)**2 at k=nzb+1
+             ELSE
+                kh   = 0.00001
+                km   = 0.00001
+             ENDIF
+             e    = 0.0
+          ENDIF
           rif   = 0.0
           ts    = 0.0
@@ -892 +913 @@
     rdf = 0.0
     IF ( rayleigh_damping_factor /= 0.0 )  THEN
-       DO  k = nzb+1, nzt
-          IF ( zu(k) >= rayleigh_damping_height )  THEN
-             rdf(k) = rayleigh_damping_factor * &
+       IF ( .NOT. ocean )  THEN
+          DO  k = nzb+1, nzt
+             IF ( zu(k) >= rayleigh_damping_height )  THEN
+                rdf(k) = rayleigh_damping_factor * &
                       ( SIN( pi * 0.5 * ( zu(k) - rayleigh_damping_height )    &
                                       / ( zu(nzt) - rayleigh_damping_height ) )&
                       )**2
-          ENDIF
-       ENDDO
+             ENDIF
+          ENDDO
+       ELSE
+          DO  k = nzt, nzb+1, -1
+             IF ( zu(k) <= rayleigh_damping_height )  THEN
+                rdf(k) = rayleigh_damping_factor * &
+                      ( SIN( pi * 0.5 * ( rayleigh_damping_height - zu(k) )    &
+                                      / ( rayleigh_damping_height - zu(nzb+1)))&
+                      )**2
+             ENDIF
+          ENDDO
+       ENDIF
     ENDIF