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

Timestamp:

Oct 30, 2013 11:36:58 AM (10 years ago)

Author:

heinze

Message:

Nudging and large scale forcing from external file implemented

File:

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

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! output for nudging + large scale forcing from external file
 !
 ! Former revisions:
@@ -428 +428 @@
         WRITE ( io, 154 )
     ENDIF
+    IF ( nudging )  THEN
+        WRITE ( io, 155 )
+        WRITE ( io, 156 )
+    ENDIF
+    IF ( large_scale_forcing )  THEN
+        WRITE ( io, 157 )
+        WRITE ( io, 158 )
+    ENDIF
     WRITE ( io, 99 )
 
@@ -706 +714 @@
        WRITE ( io, 303 )
        IF ( constant_heatflux )  THEN
-          WRITE ( io, 306 )  surface_heatflux
+          IF ( large_scale_forcing .AND. lsf_surf ) THEN
+             WRITE ( io, 306 )  shf(0,0)
+          ELSE
+             WRITE ( io, 306 )  surface_heatflux
+          ENDIF
           IF ( random_heatflux )  WRITE ( io, 307 )
        ENDIF
        IF ( humidity  .AND.  constant_waterflux )  THEN
-          WRITE ( io, 311 ) surface_waterflux
+          IF ( large_scale_forcing .AND. lsf_surf ) THEN
+             WRITE ( io, 311 ) qsws(0,0)
+          ELSE
+             WRITE ( io, 311 ) surface_waterflux
+          ENDIF
        ENDIF
        IF ( passive_scalar  .AND.  constant_waterflux )  THEN
@@ -1288 +1304 @@
     ENDDO
 
-    WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
-                       TRIM( gradients ), TRIM( slices )
+    IF ( .NOT. large_scale_forcing ) THEN
+       WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
+                          TRIM( gradients ), TRIM( slices )
+    ELSE
+       WRITE ( io, 429 ) 
+    ENDIF
 
 !-- Profile of the geostrophic wind (component vg)
@@ -1320 +1340 @@
     ENDDO
 
-    WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
-                       TRIM( gradients ), TRIM( slices )
+    IF ( .NOT. large_scale_forcing ) THEN
+       WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
+                          TRIM( gradients ), TRIM( slices )
+    ENDIF
 
 !
@@ -1357 +1379 @@
     ENDDO
 
-    WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
-                       TRIM( gradients ), TRIM( slices )
+    IF ( .NOT. nudging ) THEN
+       WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
+                          TRIM( gradients ), TRIM( slices )
+    ELSE
+       WRITE ( io, 428 ) 
+    ENDIF
 
 !
@@ -1392 +1418 @@
 
        IF ( humidity )  THEN
-          WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
-                             TRIM( gradients ), TRIM( slices )
+          IF ( .NOT. nudging ) THEN
+             WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
+                                TRIM( gradients ), TRIM( slices )
+          ENDIF
        ELSE
           WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ), &
@@ -1705 +1733 @@
                   'prognostic equation for')
 154 FORMAT ('     the potential temperature')
+155 FORMAT (' --> Nudging is used - initial profiles for u, v, pt and q ',& 
+                  'correspond to the')
+156 FORMAT ('     first profiles in NUDGING_DATA')
+157 FORMAT (' --> Large scale forcing from external file (LSF_DATA) is used: ')
+158 FORMAT ('     prescribed surfaces fluxes and geostrophic wind components')
 200 FORMAT (//' Run time and time step information:'/ &
              ' ----------------------------------'/)
@@ -1957 +1990 @@
 427 FORMAT (/'    Initial wind profiles (u,v) are interpolated from given'// &
                   ' profiles')
+428 FORMAT (/'    Initial profiles (u, v, pt, q) are taken from file '/ &
+             '    NUDGING_DATA')
+429 FORMAT (/'    Geostrophic wind profiles (ug, vg) are are taken from file '/ &
+             '    LSF_DATA_DATA')
 430 FORMAT (//' Cloud physics quantities / methods:'/ &
               ' ----------------------------------'/)