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

Timestamp:

Mar 6, 2014 1:15:21 PM (10 years ago)

Author:

heinze

Message:

enable usage of large_scale subsidence in combination with large_scale_forcing

File:

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

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! output for using large_scale subsidence in combination 
+! with large_scale_forcing
+! reformatting, more detailed explanations
 !
 ! Former revisions:
@@ -368 +370 @@
     IF ( momentum_advec == 'pw-scheme' )  THEN
        WRITE ( io, 113 )
-    ELSEIF (momentum_advec == 'ws-scheme' ) THEN
+    ELSEIF (momentum_advec == 'ws-scheme' )  THEN
        WRITE ( io, 503 )
     ENDIF
@@ -578 +580 @@
     ENDIF
 
-    IF ( plant_canopy ) THEN
+    IF ( plant_canopy )  THEN
 
        WRITE ( io, 280 ) canopy_mode, pch_index, drag_coefficient
-       IF ( passive_scalar ) THEN
+       IF ( passive_scalar )  THEN
           WRITE ( io, 281 ) scalar_exchange_coefficient,   &
                             leaf_surface_concentration
@@ -717 +719 @@
        WRITE ( io, 303 )
        IF ( constant_heatflux )  THEN
-          IF ( large_scale_forcing .AND. lsf_surf ) THEN
+          IF ( large_scale_forcing .AND. lsf_surf )  THEN
              WRITE ( io, 306 )  shf(0,0)
           ELSE
@@ -725 +727 @@
        ENDIF
        IF ( humidity  .AND.  constant_waterflux )  THEN
-          IF ( large_scale_forcing .AND. lsf_surf ) THEN
+          IF ( large_scale_forcing .AND. lsf_surf )  THEN
              WRITE ( io, 311 ) qsws(0,0)
           ELSE
@@ -1227 +1229 @@
 !
 !-- Spectra output
-    IF ( dt_dosp /= 9999999.9 ) THEN
+    IF ( dt_dosp /= 9999999.9 )  THEN
        WRITE ( io, 370 )
 
@@ -1269 +1271 @@
 !
 !-- Cloud physics parameters
-    IF ( cloud_physics ) THEN
+    IF ( cloud_physics )  THEN
        WRITE ( io, 415 )
        WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
@@ -1307 +1309 @@
     ENDDO
 
-    IF ( .NOT. large_scale_forcing ) THEN
+    IF ( .NOT. large_scale_forcing )  THEN
        WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
                           TRIM( gradients ), TRIM( slices )
@@ -1343 +1345 @@
     ENDDO
 
-    IF ( .NOT. large_scale_forcing ) THEN
+    IF ( .NOT. large_scale_forcing )  THEN
        WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
                           TRIM( gradients ), TRIM( slices )
@@ -1382 +1384 @@
     ENDDO
 
-    IF ( .NOT. nudging ) THEN
+    IF ( .NOT. nudging )  THEN
        WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
                           TRIM( gradients ), TRIM( slices )
@@ -1421 +1423 @@
 
        IF ( humidity )  THEN
-          IF ( .NOT. nudging ) THEN
+          IF ( .NOT. nudging )  THEN
              WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
                                 TRIM( gradients ), TRIM( slices )
@@ -1498 +1500 @@
        ENDDO
 
-       WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
-                          TRIM( gradients ), TRIM( slices )
+ 
+       IF ( .NOT. large_scale_forcing )  THEN
+          WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
+                             TRIM( gradients ), TRIM( slices )
+        ELSE
+           WRITE ( io, 460 ) 
+       ENDIF
+
+
     ENDIF
 
@@ -1734 +1743 @@
            /'     starting from dp_level_b =', F8.3, 'm', A /)
 153 FORMAT (' --> Large-scale vertical motion is used in the ', &
-                  'prognostic equation for')
-154 FORMAT ('     the potential temperature')
+                  'prognostic equation(s) for')
+154 FORMAT ('     the scalar(s) only')
 155 FORMAT (' --> Nudging is used - initial profiles for u, v, pt and q ',& 
                   'correspond to the')
@@ -1996 +2005 @@
              '    NUDGING_DATA')
 429 FORMAT (/'    Geostrophic wind profiles (ug, vg) are are taken from file '/ &
-             '    LSF_DATA_DATA')
+             '    LSF_DATA')
 430 FORMAT (//' Cloud physics quantities / methods:'/ &
               ' ----------------------------------'/)
@@ -2023 +2032 @@
 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')
+460 FORMAT (/'    Profiles for large scale vertical velocity are '/ &
+             '    taken from file LSF_DATA')
 470 FORMAT (//' Actions during the simulation:'/ &
               ' -----------------------------'/)