Changeset 103 for palm/trunk/SOURCE

Timestamp:

Jul 30, 2007 10:49:46 AM (17 years ago)

Author:

raasch

Message:

further preliminary changes concerning coupling

Location:

palm/trunk/SOURCE

Files:

• check_parameters.f90 (modified) (1 diff)
• diffusion_u.f90 (modified) (2 diffs)
• diffusion_v.f90 (modified) (2 diffs)
• header.f90 (modified) (3 diffs)
• modules.f90 (modified) (3 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -1083 +1083 @@
     IF ( surface_heatflux == 9999999.9 )  constant_heatflux     = .FALSE.
     IF ( top_heatflux     == 9999999.9 )  constant_top_heatflux = .FALSE.
+    IF ( top_momentumflux_u /= 9999999.9  .AND.  &
+         top_momentumflux_v /= 9999999.9 )  THEN
+       constant_top_momentumflux = .TRUE.
+    ELSEIF (  .NOT. ( top_momentumflux_u == 9999999.9  .AND.  &
+           top_momentumflux_v == 9999999.9 ) )  THEN   
+       IF ( myid == 0 )  THEN
+          PRINT*, '+++ check_parameters:'
+          PRINT*, '    both, top_momentumflux_u AND top_momentumflux_v'
+          PRINT*, '    must be set'
+       ENDIF
+       CALL local_stop
+    ENDIF
 
 !

palm/trunk/SOURCE/diffusion_u.f90

@@ -204 +204 @@
 !--          Vertical diffusion at the first gridpoint below the top boundary,
 !--          if the momentum flux at the top is prescribed by the user
-             IF ( use_top_fluxes )  THEN
+             IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
                 k = nzt
 !
@@ -368 +368 @@
 !--    Vertical diffusion at the first gridpoint below the top boundary,
 !--    if the momentum flux at the top is prescribed by the user
-       IF ( use_top_fluxes )  THEN
+       IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
           k = nzt
 !

palm/trunk/SOURCE/diffusion_v.f90

@@ -202 +202 @@
 !--          Vertical diffusion at the first gridpoint below the top boundary,
 !--          if the momentum flux at the top is prescribed by the user
-             IF ( use_top_fluxes )  THEN
+             IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
                 k = nzt
 !
@@ -367 +367 @@
 !--    Vertical diffusion at the first gridpoint below the top boundary,
 !--    if the momentum flux at the top is prescribed by the user
-       IF ( use_top_fluxes )  THEN
+       IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
           k = nzt
 !

palm/trunk/SOURCE/header.f90

@@ -5 +5 @@
 ! -----------------
 ! Output of informations for coupled model runs (boundary conditions etc.)
+! + output of momentumfluxes at the top boundary
 !
 ! Former revisions:
@@ -462 +463 @@
        WRITE ( io, 304 )
        IF ( coupling_mode == 'uncoupled' )  THEN
+          WRITE ( io, 319 )  top_momentumflux_u, top_momentumflux_v
           IF ( constant_top_heatflux )  THEN
              WRITE ( io, 306 )  top_heatflux
@@ -1349 +1351 @@
 318 FORMAT (/'       outflow damping layer width: ',I3,' gridpoints with km_', &
                     'max =',F5.1,' m**2/s')
+319 FORMAT ('       Predefined constant momentumflux:  u: ',F9.6,' m**2/s**2'/ &
+            '                                          v: ',F9.6,' m**2/s**2')
 320 FORMAT (//' List output:'/ &
              ' -----------'//  &

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Actual revisions:
 ! -----------------
-! +comm_inter, coupling_char, coupling_mode, dt_coupling, ngp_xy, port_name,
-! time_coupling, type_xy, uswst*, vswst*
+! +comm_inter, constant_top_momentumflux, coupling_char, coupling_mode,
+! dt_coupling, ngp_xy, port_name, time_coupling, top_momentumflux_u|v,
+! type_xy, uswst*, vswst*
 !
 ! Former revisions:
@@ -311 +312 @@
                 conserve_volume_flow = .FALSE., constant_diffusion = .FALSE., &
                 constant_heatflux = .TRUE., constant_top_heatflux = .TRUE., &
+                constant_top_momentumflux = .FALSE., &
                 constant_top_salinityflux = .TRUE., &
                 constant_waterflux = .TRUE., create_disturbances = .TRUE., &
@@ -396 +398 @@
              time_do_sla = 0.0, time_dvrp = 0.0, time_prel = 0.0, &
              time_restart = 9999999.9, time_run_control = 0.0, &
-             top_heatflux = 9999999.9, top_momentumflux_u = 0.0, &
-             top_momentumflux_v = 0.0, top_salinityflux = 9999999.9, &
+             top_heatflux = 9999999.9, top_momentumflux_u = 9999999.9, &
+             top_momentumflux_v = 9999999.9, top_salinityflux = 9999999.9, &
              ug_surface = 0.0, u_gtrans = 0.0, &
              ups_limit_e = 0.0, ups_limit_pt = 0.0, ups_limit_u = 0.0, &