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

palm/trunk/SOURCE/header.f90

@@ -6 +6 @@
 ! Output of informations for coupled model runs (boundary conditions etc.)
 ! + output of momentumfluxes at the top boundary
+! Rayleigh damping for ocean, e_init
 !
 ! Former revisions:
@@ -241 +242 @@
     IF ( use_upstream_for_tke )  WRITE ( io, 143 )
     IF ( rayleigh_damping_factor /= 0.0 )  THEN
-       WRITE ( io, 123 )  rayleigh_damping_height, rayleigh_damping_factor
+       IF ( .NOT. ocean )  THEN
+          WRITE ( io, 123 )  'above', rayleigh_damping_height, &
+               rayleigh_damping_factor
+       ELSE
+          WRITE ( io, 123 )  'below', rayleigh_damping_height, &
+               rayleigh_damping_factor
+       ENDIF
     ENDIF
     IF ( humidity )  THEN
@@ -1091 +1098 @@
     ENDIF
     IF ( .NOT. constant_diffusion)  THEN
+       IF ( e_init > 0.0 )  WRITE ( io, 455 )  e_init
        IF ( e_min > 0.0 )  WRITE ( io, 454 )  e_min
        IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
@@ -1240 +1248 @@
                   ' timestep changes')
 122 FORMAT (' --> Time differencing scheme: ',A)
-123 FORMAT (' --> Rayleigh-Damping active, starts above z = ',F8.2,' m'/ &
+123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
             '     maximum damping coefficient: ',F5.3, ' 1/s')
 124 FORMAT ('     Spline-overshoots are being suppressed')
@@ -1468 +1476 @@
 453 FORMAT ('   Mixing length is limited to ',F4.2,' * z')
 454 FORMAT ('   TKE is not allowed to fall below ',E9.2,' (m/s)**2')
+455 FORMAT ('   initial TKE is prescribed as ',E9.2,' (m/s)**2')
 470 FORMAT (//' Actions during the simulation:'/ &
               ' -----------------------------'/)