Changeset 2669 for palm/trunk/UTIL

Timestamp:

Dec 6, 2017 4:03:27 PM (7 years ago)

Author:

raasch

Message:

file attributes and activation strings in .palm.iofiles revised, file extensions for nesting, masked output, wind turbine data, etc. changed

File:

• palm/trunk/UTIL/combine_plot_fields.f90 (modified) (10 diffs)

palm/trunk/UTIL/combine_plot_fields.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! data of 3d-nest runs are completely processed now
+! 
+! 2523 2017-10-05 14:42:47Z kanani
 ! Increased LEN for CHARACTER variable var_name, equal to the value in PALM
 ! 
@@ -105 +108 @@
 !
 !-- Local variables
-    CHARACTER (LEN=2)    ::  modus, model_string
+    CHARACTER (LEN=2)    ::  modus
+    CHARACTER (LEN=4)    ::  model_string
     CHARACTER (LEN=6)    ::  id_string
     CHARACTER (LEN=30)   ::  dimname, var_name
@@ -116 +120 @@
     INTEGER(iwp) ::  av, danz, i, id, j, k, model, models, nc_stat,            &
                      nxa, nxag, nxe, nxeg, nya, nyag, nye, nyeg,               &
-                     nza, nze, pos, time_step,                                 &
-                     xa, xe, xxa, xxe, ya, ya_do, ya_tot, ye, ye_do, ye_tot, yya, yye, za,   &
-                     ze, zza, zze
+                     nza, nze, pos, time_step, xa, xe, xxa, xxe, ya, ya_do,    &
+                     ya_tot, ye, ye_do, ye_tot, yya, yye, za, ze, zza, zze
 
     INTEGER(8)                        ::  count, count_rate
@@ -128 +131 @@
     INTEGER(iwp), DIMENSION(0:1,1000) ::  id_var, levels
 
-    LOGICAL  ::  found, netcdf_output, netcdf_parallel, netcdf_0, netcdf_1
-    LOGICAL  ::  vnest
+    LOGICAL  ::  found, nest3d, netcdf_output, netcdf_parallel, netcdf_0,      &
+                 netcdf_1, vnest
 
     REAL(wp) ::  cpu_start_time, cpu_end_time, dx, simulated_time
@@ -142 +145 @@
 
 !
-!-- Find out if a coupled run has been carried out
+!-- Find out if a coupled or nested run has been carried out
     INQUIRE( FILE='COUPLING_PORT_OPENED', EXIST=found )
     INQUIRE( FILE='VNESTING_PORT_OPENED', EXIST=vnest )
+    INQUIRE( FILE='3DNESTING', EXIST=nest3d )
     IF ( found )  THEN
        models = 2
@@ -151 +155 @@
        models = 2
        PRINT*, '    Vertically nested grid coupling'
+    ELSEIF ( nest3d )  THEN
+       OPEN( 90, FILE='3DNESTING', FORM='FORMATTED' )
+       READ ( 90, '(I2)' )  models
+       CLOSE ( 90 )
+       PRINT*, '    3d-nest run'
+       PRINT*, '    number of nest domains = ', models
     ELSE
        models = 1
@@ -171 +181 @@
 !
 !--    Set the model string used to identify the filenames
-       IF ( models == 2 )  THEN
+       IF ( found  .OR.  vnest )  THEN
           PRINT*, ''
-          PRINT*, '*** combine_plot_fields ***'
           IF ( model == 2 )  THEN
              IF ( vnest )  THEN
@@ -193 +202 @@
              ENDIF
           ENDIF
+       ELSEIF ( nest3d )  THEN
+          PRINT*, ''
+          PRINT*, '--> processing nest id = ', model
+          IF ( model == 1 )  THEN
+             model_string = ''
+          ELSE
+             WRITE( model_string, '(A2,I2.2)' )  '_N', model
+          ENDIF
        ENDIF
 !
@@ -260 +277 @@
 !--       Info-output
           PRINT*, ''
-          PRINT*, '*** combine_plot_fields ***'
 #if defined( __netcdf )
           IF ( netcdf_output )  THEN
@@ -568 +584 @@
 !--    Info-output
        PRINT*, ' '
-       PRINT*, '*** combine_plot_fields ***'
 
 !