Changeset 19 for palm/trunk/SOURCE/header.f90

Timestamp:

Feb 23, 2007 4:53:48 AM (17 years ago)

Author:

raasch

Message:

preliminary version of modified boundary conditions at top

File:

• palm/trunk/SOURCE/header.f90 (modified) (4 diffs)

palm/trunk/SOURCE/header.f90

@@ -345 +345 @@
     ENDIF
     IF ( ibc_pt_t == 0 )  THEN
-       roben  = TRIM( roben  ) // ' pt(nzt) = pt_top'
-    ELSE
-       roben  = TRIM( roben  ) // ' pt(nzt) = pt(nzt-1) + dpt/dz'
+       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt_top'
+    ELSEIF( ibc_pt_t == 1 )  THEN
+       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt)'
+    ELSEIF( ibc_pt_t == 2 )  THEN
+       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
     ENDIF
 
@@ -404 +406 @@
        IF ( passive_scalar  .AND.  constant_waterflux )  THEN
           WRITE ( io, 313 ) surface_waterflux
+       ENDIF
+    ENDIF
+
+    IF ( use_top_fluxes )  THEN
+       WRITE ( io, 304 )
+       IF ( constant_top_heatflux )  THEN
+          WRITE ( io, 306 )  top_heatflux
+       ENDIF
+       IF ( moisture  .OR.  passive_scalar )  THEN
+          WRITE ( io, 315 )
        ENDIF
     ENDIF
@@ -1220 +1232 @@
              ' B. bound.: ',A/ &
              ' T. bound.: ',A)
-303 FORMAT (/' Surface fluxes are used in diffusion terms at k=1')
-305 FORMAT (//'    Prandtl-Layer between surface and first computational ', &
-               'u,v-level:'// &
+303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
+304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
+305 FORMAT (//'    Prandtl-Layer between bottom surface and first ', &
+               'computational u,v-level:'// &
              '       zp = ',F6.2,' m   z0 = ',F6.4,' m   kappa = ',F4.2/ &
              '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
@@ -1234 +1247 @@
 313 FORMAT ('       Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
 314 FORMAT ('       Predefined scalar value at the surface')
+315 FORMAT ('       Humidity / scalar flux at top surface is 0.0')
 317 FORMAT (//' Lateral boundaries:'/ &
             '       left/right:  ',A/    &