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

Timestamp:

Feb 28, 2016 12:45:19 PM (8 years ago)

Author:

raasch

Message:

update of the nested domain system + some bugfixes

File:

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

palm/trunk/SOURCE/header.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! output of nesting informations
 !
 ! Former revisions:
@@ -205 +205 @@
 ! Description:
 ! ------------
-!> Writing a header with all important information about the actual run.
+!> Writing a header with all important information about the current run.
 !> This subroutine is called three times, two times at the beginning 
 !> (writing information on files RUN_CONTROL and HEADER) and one time at the
@@ -269 +269 @@
                plant_canopy
 
+    USE pmc_interface,                                                         &
+        ONLY:  cpl_id, cpl_parent_id, cpl_name, lower_left_coord_x,            &
+               lower_left_coord_y, nested_run, nesting_mode
+
     USE radiation_model_mod,                                                   &
         ONLY:  albedo, albedo_type, albedo_type_name, constant_albedo,         &
@@ -366 +370 @@
 !-- Determine kind of model run
     IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
-       run_classification = '3D - restart run'
+       run_classification = 'restart run'
     ELSEIF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
-       run_classification = '3D - run with cyclic fill of 3D - prerun data'
+       run_classification = 'run with cyclic fill of 3D - prerun data'
     ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 )  THEN
-       run_classification = '3D - run without 1D - prerun'
+       run_classification = 'run without 1D - prerun'
     ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
-       run_classification = '3D - run with 1D - prerun'
+       run_classification = 'run with 1D - prerun'
     ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
-       run_classification = '3D - run initialized by user'
+       run_classification = 'run initialized by user'
     ELSE
        message_string = ' unknown action(s): ' // TRIM( initializing_actions )
        CALL message( 'header', 'PA0191', 0, 0, 0, 6, 0 )
     ENDIF
+    IF ( nested_run )  run_classification = 'nested ' // run_classification
     IF ( ocean )  THEN
        run_classification = 'ocean - ' // run_classification
@@ -444 +449 @@
     IF ( num_acc_per_node /= 0 )  WRITE ( io, 120 )  num_acc_per_node
 #endif
+
+!
+!-- Nesting informations
+    IF ( nested_run )  THEN
+       WRITE ( io, 600 )  cpl_id, TRIM( cpl_name ), cpl_parent_id,             &
+                          nesting_mode, lower_left_coord_x, lower_left_coord_y
+    ENDIF
     WRITE ( io, 99 )
 
@@ -2399 +2411 @@
 513 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order ' // & 
             '+ monotonic adjustment')
-
+600 FORMAT (/' Nesting informations:'/                                        &
+            ' Nest id / name:                   ',I2.2,' / ',A,' (parent id: ',I2.2,')'/ &
+            ' Nesting mode:                     ',A/ &
+            ' Lower left corner coordinates:    ','x = ',F8.2,' m, y = ',F8.2,' m'/)
 
  END SUBROUTINE header