Changeset 3648


Ignore:
Timestamp:
Jan 2, 2019 4:35:46 PM (3 years ago)
Author:
suehring
Message:

Rename surface_output_mod into surface_data_output_mod, same with all corresponding subroutines

Location:
palm/trunk/SOURCE
Files:
6 edited
1 moved

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r3637 r3648  
    2525# -----------------
    2626# $Id$
     27# Rename surface_output_mod into surface_data_output_mod
     28#
     29# 3637 2018-12-20 01:51:36Z knoop
    2730# Implementation of the PALM module interface
    2831#
     
    672675        surface_layer_fluxes_mod.f90 \
    673676        surface_mod.f90 \
    674         surface_output_mod.f90 \
     677        surface_data_output_mod.f90 \
    675678        swap_timelevel.f90 \
    676679        synthetic_turbulence_generator_mod.f90 \
     
    11071110        surface_layer_fluxes_mod.o \
    11081111        surface_mod.o \
    1109         surface_output_mod.o \
     1112        surface_data_output_mod.o \
    11101113        synthetic_turbulence_generator_mod.o \
    11111114        turbulence_closure_mod.o \
     
    13371340        salsa_mod.o \
    13381341        spectra_mod.o \
    1339         surface_output_mod.o \
     1342        surface_data_output_mod.o \
    13401343        synthetic_turbulence_generator_mod.o \
    13411344        urban_surface_mod.o \
     
    14221425        salsa_mod.o \
    14231426        surface_layer_fluxes_mod.o \
    1424         surface_output_mod.o \
     1427        surface_data_output_mod.o \
    14251428        time_to_string.o \
    14261429        write_restart_data_mod.o
     
    16541657        mod_kinds.o \
    16551658        modules.o
    1656 surface_output_mod.o: \
     1659surface_data_output_mod.o: \
    16571660        mod_kinds.o \
    16581661        modules.o \
     
    17071710        surface_layer_fluxes_mod.o \
    17081711        surface_mod.o \
    1709         surface_output_mod.o \
     1712        surface_data_output_mod.o \
    17101713        synthetic_turbulence_generator_mod.o \
    17111714        time_to_string.o \
  • palm/trunk/SOURCE/init_3d_model.f90

    r3636 r3648  
    1515! PALM. If not, see <http://www.gnu.org/licenses/>.
    1616!
    17 ! Copyright 1997-2018 Leibniz Universitaet Hannover
     17! Copyright 1997-2019 Leibniz Universitaet Hannover
    1818!------------------------------------------------------------------------------!
    1919!
     
    2525! -----------------
    2626! $Id$
     27! Rename subroutines for surface-data output
     28!
     29! 3636 2018-12-19 13:48:34Z raasch
    2730! nopointer option removed
    2831!
     
    687690                surf_usm_h, get_topography_top_index_ji, vertical_surfaces_exist
    688691   
    689     USE surface_output_mod,                                                    &
    690         ONLY:  surface_output_init
     692    USE surface_data_output_mod,                                               &
     693        ONLY:  surface_data_output_init
    691694   
    692695    USE transpose_indices
     
    24462449!
    24472450!-- Initialize surface data output
    2448     IF ( surface_data_output )  THEN
    2449        CALL surface_output_init
     2451    IF ( surface_output )  THEN
     2452       CALL surface_data_output_init
    24502453    ENDIF
    24512454!
  • palm/trunk/SOURCE/module_interface.f90

    r3641 r3648  
    1515! PALM. If not, see <http://www.gnu.org/licenses/>.
    1616!
    17 ! Copyright 1997-2018 Leibniz Universitaet Hannover
     17! Copyright 1997-2019 Leibniz Universitaet Hannover
    1818!------------------------------------------------------------------------------!
    1919!
     
    2525! -----------------
    2626! $Id$
     27! Rename subroutines for surface-data output
     28!
     29! 3641 2018-12-23 22:10:01Z knoop
    2730! Initial implementation of the PALM module interface
    2831!
     
    5356              plant_canopy,                                                    &
    5457              salsa,                                                           &
    55               surface_data_output,                                             &
     58              surface_output,                                                  &
    5659              syn_turb_gen,                                                    &
    5760              urban_surface,                                                   &
     
    204207              spectra_header
    205208
    206    USE surface_output_mod,                                                     &
    207        ONLY:  surface_output_parin,                                            &
    208               surface_output_check_parameters
     209   USE surface_data_output_mod,                                                &
     210       ONLY:  surface_data_output_parin,                                       &
     211              surface_data_output_check_parameters
    209212
    210213   USE synthetic_turbulence_generator_mod,                                     &
     
    368371   CALL salsa_parin
    369372   CALL spectra_parin
    370    CALL surface_output_parin
     373   CALL surface_data_output_parin
    371374   CALL stg_parin
    372375   CALL user_parin ! ToDo: make user code a single Fortran module
     
    401404   IF ( salsa )                CALL salsa_check_parameters
    402405   IF ( calculate_spectra )    CALL spectra_check_parameters
    403    IF ( surface_data_output )  CALL surface_output_check_parameters
     406   IF ( surface_output )       CALL surface_data_output_check_parameters
    404407   IF ( syn_turb_gen )         CALL stg_check_parameters
    405408   IF ( urban_surface )        CALL usm_check_parameters
  • palm/trunk/SOURCE/modules.f90

    r3636 r3648  
    1515! PALM. If not, see <http://www.gnu.org/licenses/>.
    1616!
    17 ! Copyright 1997-2018 Leibniz Universitaet Hannover
     17! Copyright 1997-2019 Leibniz Universitaet Hannover
    1818!------------------------------------------------------------------------------!
    1919!
     
    2525! -----------------
    2626! $Id$
     27! -surface_data_output +surface_output
     28!
     29! 3636 2018-12-19 13:48:34Z raasch
    2730! nopointer option removed
    2831!
     
    13511354    LOGICAL ::  sloping_surface = .FALSE.                        !< use sloped surface? (namelist parameter alpha_surface)
    13521355    LOGICAL ::  spinup = .FALSE.                                 !< perform model spinup without atmosphere code?
    1353     LOGICAL ::  surface_data_output = .FALSE.                    !< output of surface data
     1356    LOGICAL ::  surface_output = .FALSE.                         !< output of surface data
    13541357    LOGICAL ::  stop_dt = .FALSE.                                !< internal switch to stop the time stepping
    13551358    LOGICAL ::  synchronous_exchange = .FALSE.                   !< namelist parameter
  • palm/trunk/SOURCE/palm.f90

    r3524 r3648  
    1515! PALM. If not, see <http://www.gnu.org/licenses/>.
    1616!
    17 ! Copyright 1997-2018 Leibniz Universitaet Hannover
     17! Copyright 1997-2019 Leibniz Universitaet Hannover
    1818!------------------------------------------------------------------------------!
    1919!
     
    2525! -----------------
    2626! $Id$
     27! Rename subroutines for surface-data output
     28!
     29! 3524 2018-11-14 13:36:44Z raasch
    2730! unused variable removed
    2831!
     
    321324               pmci_modelconfiguration, pmci_parent_initialize
    322325               
    323     USE surface_output_mod,                                                    &
    324         ONLY:  surface_output_last_action
     326    USE surface_data_output_mod,                                               &
     327        ONLY:  surface_data_output_last_action
    325328
    326329    USE write_restart_data_mod,                                                &
     
    578581!
    579582!-- Last actions for surface output, for instantaneous and time-averaged data
    580     CALL surface_output_last_action( 0 )
    581     CALL surface_output_last_action( 1 )
     583    CALL surface_data_output_last_action( 0 )
     584    CALL surface_data_output_last_action( 1 )
    582585
    583586!
  • palm/trunk/SOURCE/surface_data_output_mod.f90

    r3647 r3648  
    1 !> @file surface_output_mod.f90
     1!> @file surface_data_output_mod.f90
    22!------------------------------------------------------------------------------!
    33! This file is part of the PALM model system.
     
    1515! PALM. If not, see <http://www.gnu.org/licenses/>.
    1616!
    17 ! Copyright 1997-2018 Leibniz Universitaet Hannover
     17! Copyright 1997-2019 Leibniz Universitaet Hannover
    1818!------------------------------------------------------------------------------!
    1919!
     
    2525! -----------------
    2626! $Id$
     27! Rename module and subroutines
     28!
     29! 3646 2018-12-28 17:58:49Z kanani
    2730! Bugfix: use time_since_reference_point instead of simulated_time (relevant
    2831! when using wall/soil spinup)
     
    6366!------------------------------------------------------------------------------!
    6467
    65 MODULE surface_output_mod
     68MODULE surface_data_output_mod
    6669
    6770   USE kinds
     
    7174       
    7275   USE control_parameters,                                                     &
    73        ONLY:  surface_data_output
     76       ONLY:  surface_output
    7477       
    7578   USE grid_variables,                                                         &
     
    126129   PRIVATE
    127130
    128    INTERFACE  surface_output
    129       MODULE PROCEDURE surface_output
    130    END INTERFACE  surface_output
     131   INTERFACE  surface_data_output
     132      MODULE PROCEDURE surface_data_output
     133   END INTERFACE  surface_data_output
    131134   
    132    INTERFACE  surface_output_averaging
    133       MODULE PROCEDURE surface_output_averaging
    134    END INTERFACE  surface_output_averaging
     135   INTERFACE  surface_data_output_averaging
     136      MODULE PROCEDURE surface_data_output_averaging
     137   END INTERFACE  surface_data_output_averaging
    135138   
    136    INTERFACE  surface_output_check_parameters
    137       MODULE PROCEDURE surface_output_check_parameters
    138    END INTERFACE  surface_output_check_parameters
     139   INTERFACE  surface_data_output_check_parameters
     140      MODULE PROCEDURE surface_data_output_check_parameters
     141   END INTERFACE  surface_data_output_check_parameters
    139142   
    140    INTERFACE  surface_output_init
    141       MODULE PROCEDURE surface_output_init
    142    END INTERFACE  surface_output_init
     143   INTERFACE  surface_data_output_init
     144      MODULE PROCEDURE surface_data_output_init
     145   END INTERFACE  surface_data_output_init
    143146   
    144    INTERFACE  surface_output_last_action
    145       MODULE PROCEDURE surface_output_last_action
    146    END INTERFACE  surface_output_last_action
     147   INTERFACE  surface_data_output_last_action
     148      MODULE PROCEDURE surface_data_output_last_action
     149   END INTERFACE  surface_data_output_last_action
    147150     
    148    INTERFACE  surface_output_parin
    149       MODULE PROCEDURE surface_output_parin
    150    END INTERFACE  surface_output_parin
     151   INTERFACE  surface_data_output_parin
     152      MODULE PROCEDURE surface_data_output_parin
     153   END INTERFACE  surface_data_output_parin
    151154
    152155!
    153156!--Public subroutines
    154    PUBLIC surface_output, surface_output_averaging,                            &
    155           surface_output_check_parameters, surface_output_init,                &
    156           surface_output_last_action, surface_output_parin
     157   PUBLIC surface_data_output, surface_data_output_averaging,                  &
     158          surface_data_output_check_parameters, surface_data_output_init,      &
     159          surface_data_output_last_action, surface_data_output_parin
    157160!
    158161!--Public variables
     
    168171!> and polygon data for the surface elements, allocation of required arrays.
    169172!------------------------------------------------------------------------------!
    170    SUBROUTINE surface_output_init
     173   SUBROUTINE surface_data_output_init
    171174   
    172175      IMPLICIT NONE
     
    855858#endif
    856859
    857    END SUBROUTINE surface_output_init
     860   END SUBROUTINE surface_data_output_init
    858861   
    859862!------------------------------------------------------------------------------!
     
    865868!> are called to write the surface data in the respective NetCDF files.
    866869!------------------------------------------------------------------------------!   
    867    SUBROUTINE surface_output( av )
     870   SUBROUTINE surface_data_output( av )
    868871   
    869872      USE control_parameters,                                                  &
     
    922925!--            Output of instantaneous data
    923926               IF ( av == 0 )  THEN
    924                   CALL surface_output_collect( surf_def_h(0)%us,               &
     927                  CALL surface_data_output_collect( surf_def_h(0)%us,          &
    925928                                               surf_def_h(1)%us,               &
    926929                                               surf_lsm_h%us,                  &
     
    951954!--            Output of instantaneous data
    952955               IF ( av == 0 )  THEN
    953                   CALL surface_output_collect( surf_def_h(0)%ts,               &
     956                  CALL surface_data_output_collect( surf_def_h(0)%ts,          &
    954957                                               surf_def_h(1)%ts,               &
    955958                                               surf_lsm_h%ts,                  &
     
    980983!--            Output of instantaneous data
    981984               IF ( av == 0 )  THEN
    982                   CALL surface_output_collect( surf_def_h(0)%qs,               &
     985                  CALL surface_data_output_collect( surf_def_h(0)%qs,          &
    983986                                               surf_def_h(1)%qs,               &
    984987                                               surf_lsm_h%qs,                  &
     
    10091012!--            Output of instantaneous data
    10101013               IF ( av == 0 )  THEN
    1011                   CALL surface_output_collect( surf_def_h(0)%ss,               &
     1014                  CALL surface_data_output_collect( surf_def_h(0)%ss,          &
    10121015                                               surf_def_h(1)%ss,               &
    10131016                                               surf_lsm_h%ss,                  &
     
    10381041!--            Output of instantaneous data
    10391042               IF ( av == 0 )  THEN
    1040                   CALL surface_output_collect( surf_def_h(0)%qcs,              &
     1043                  CALL surface_data_output_collect( surf_def_h(0)%qcs,         &
    10411044                                               surf_def_h(1)%qcs,              &
    10421045                                               surf_lsm_h%qcs,                 &
     
    10671070!--            Output of instantaneous data
    10681071               IF ( av == 0 )  THEN
    1069                   CALL surface_output_collect( surf_def_h(0)%ncs,              &
     1072                  CALL surface_data_output_collect( surf_def_h(0)%ncs,         &
    10701073                                               surf_def_h(1)%ncs,              &
    10711074                                               surf_lsm_h%ncs,                 &
     
    10961099!--            Output of instantaneous data
    10971100               IF ( av == 0 )  THEN
    1098                   CALL surface_output_collect( surf_def_h(0)%qrs,              &
     1101                  CALL surface_data_output_collect( surf_def_h(0)%qrs,         &
    10991102                                               surf_def_h(1)%qrs,              &
    11001103                                               surf_lsm_h%qrs,                 &
     
    11251128!--            Output of instantaneous data
    11261129               IF ( av == 0 )  THEN
    1127                   CALL surface_output_collect( surf_def_h(0)%nrs,              &
     1130                  CALL surface_data_output_collect( surf_def_h(0)%nrs,         &
    11281131                                               surf_def_h(1)%nrs,              &
    11291132                                               surf_lsm_h%nrs,                 &
     
    11541157!--            Output of instantaneous data
    11551158               IF ( av == 0 )  THEN
    1156                   CALL surface_output_collect( surf_def_h(0)%ol,               &
     1159                  CALL surface_data_output_collect( surf_def_h(0)%ol,          &
    11571160                                               surf_def_h(1)%ol,               &
    11581161                                               surf_lsm_h%ol,                  &
     
    11831186!--            Output of instantaneous data
    11841187               IF ( av == 0 )  THEN
    1185                   CALL surface_output_collect( surf_def_h(0)%z0,               &
     1188                  CALL surface_data_output_collect( surf_def_h(0)%z0,          &
    11861189                                               surf_def_h(1)%z0,               &
    11871190                                               surf_lsm_h%z0,                  &
     
    12121215!--            Output of instantaneous data
    12131216               IF ( av == 0 )  THEN
    1214                   CALL surface_output_collect( surf_def_h(0)%z0h,              &
     1217                  CALL surface_data_output_collect( surf_def_h(0)%z0h,         &
    12151218                                               surf_def_h(1)%z0h,              &
    12161219                                               surf_lsm_h%z0h,                 &
     
    12411244!--            Output of instantaneous data
    12421245               IF ( av == 0 )  THEN
    1243                   CALL surface_output_collect( surf_def_h(0)%z0q,              &
     1246                  CALL surface_data_output_collect( surf_def_h(0)%z0q,         &
    12441247                                               surf_def_h(1)%z0q,              &
    12451248                                               surf_lsm_h%z0q,                 &
     
    12701273!--            Output of instantaneous data
    12711274               IF ( av == 0 )  THEN
    1272                   CALL surface_output_collect( surf_def_h(0)%pt1,              &
     1275                  CALL surface_data_output_collect( surf_def_h(0)%pt1,         &
    12731276                                               surf_def_h(1)%pt1,              &
    12741277                                               surf_lsm_h%pt1,                 &
     
    12991302!--            Output of instantaneous data
    13001303               IF ( av == 0 )  THEN
    1301                   CALL surface_output_collect( surf_def_h(0)%qv1,              &
     1304                  CALL surface_data_output_collect( surf_def_h(0)%qv1,         &
    13021305                                               surf_def_h(1)%qv1,              &
    13031306                                               surf_lsm_h%qv1,                 &
     
    13281331!--            Output of instantaneous data
    13291332               IF ( av == 0 )  THEN
    1330                   CALL surface_output_collect( surf_def_h(0)%vpt1,             &
     1333                  CALL surface_data_output_collect( surf_def_h(0)%vpt1,        &
    13311334                                               surf_def_h(1)%vpt1,             &
    13321335                                               surf_lsm_h%vpt1,                &
     
    13571360!--            Output of instantaneous data
    13581361               IF ( av == 0 )  THEN
    1359                   CALL surface_output_collect( surf_def_h(0)%usws,             &
     1362                  CALL surface_data_output_collect( surf_def_h(0)%usws,        &
    13601363                                               surf_def_h(1)%usws,             &
    13611364                                               surf_lsm_h%usws,                &
     
    13861389!--            Output of instantaneous data
    13871390               IF ( av == 0 )  THEN
    1388                   CALL surface_output_collect( surf_def_h(0)%vsws,             &
     1391                  CALL surface_data_output_collect( surf_def_h(0)%vsws,        &
    13891392                                               surf_def_h(1)%vsws,             &
    13901393                                               surf_lsm_h%vsws,                &
     
    14151418!--            Output of instantaneous data
    14161419               IF ( av == 0 )  THEN
    1417                   CALL surface_output_collect( surf_def_h(0)%shf,              &
     1420                  CALL surface_data_output_collect( surf_def_h(0)%shf,         &
    14181421                                               surf_def_h(1)%shf,              &
    14191422                                               surf_lsm_h%shf,                 &
     
    14431446!--            Output of instantaneous data
    14441447               IF ( av == 0 )  THEN
    1445                   CALL surface_output_collect( surf_def_h(0)%qsws,             &
     1448                  CALL surface_data_output_collect( surf_def_h(0)%qsws,        &
    14461449                                               surf_def_h(1)%qsws,             &
    14471450                                               surf_lsm_h%qsws,                &
     
    14721475!--            Output of instantaneous data
    14731476               IF ( av == 0 )  THEN
    1474                   CALL surface_output_collect( surf_def_h(0)%ssws,             &
     1477                  CALL surface_data_output_collect( surf_def_h(0)%ssws,        &
    14751478                                               surf_def_h(1)%ssws,             &
    14761479                                               surf_lsm_h%ssws,                &
     
    15011504!--            Output of instantaneous data
    15021505               IF ( av == 0 )  THEN
    1503                   CALL surface_output_collect( surf_def_h(0)%qcsws,            &
     1506                  CALL surface_data_output_collect( surf_def_h(0)%qcsws,       &
    15041507                                               surf_def_h(1)%qcsws,            &
    15051508                                               surf_lsm_h%qcsws,               &
     
    15301533!--            Output of instantaneous data
    15311534               IF ( av == 0 )  THEN
    1532                   CALL surface_output_collect( surf_def_h(0)%ncsws,            &
     1535                  CALL surface_data_output_collect( surf_def_h(0)%ncsws,       &
    15331536                                               surf_def_h(1)%ncsws,            &
    15341537                                               surf_lsm_h%ncsws,               &
     
    15591562!--            Output of instantaneous data
    15601563               IF ( av == 0 )  THEN
    1561                   CALL surface_output_collect( surf_def_h(0)%qrsws,            &
     1564                  CALL surface_data_output_collect( surf_def_h(0)%qrsws,       &
    15621565                                               surf_def_h(1)%qrsws,            &
    15631566                                               surf_lsm_h%qrsws,               &
     
    15881591!--            Output of instantaneous data
    15891592               IF ( av == 0 )  THEN
    1590                   CALL surface_output_collect( surf_def_h(0)%nrsws,            &
     1593                  CALL surface_data_output_collect( surf_def_h(0)%nrsws,       &
    15911594                                               surf_def_h(1)%nrsws,            &
    15921595                                               surf_lsm_h%nrsws,               &
     
    16171620!--            Output of instantaneous data
    16181621               IF ( av == 0 )  THEN
    1619                   CALL surface_output_collect( surf_def_h(0)%sasws,            &
     1622                  CALL surface_data_output_collect( surf_def_h(0)%sasws,       &
    16201623                                               surf_def_h(1)%sasws,            &
    16211624                                               surf_lsm_h%sasws,               &
     
    16461649!--            Output of instantaneous data
    16471650               IF ( av == 0 )  THEN
    1648                   CALL surface_output_collect( surf_def_h(0)%q_surface,        &
     1651                  CALL surface_data_output_collect( surf_def_h(0)%q_surface,   &
    16491652                                               surf_def_h(1)%q_surface,        &
    16501653                                               surf_lsm_h%q_surface,           &
     
    16751678!--            Output of instantaneous data
    16761679               IF ( av == 0 )  THEN
    1677                   CALL surface_output_collect( surf_def_h(0)%pt_surface,       &
     1680                  CALL surface_data_output_collect( surf_def_h(0)%pt_surface,  &
    16781681                                               surf_def_h(1)%pt_surface,       &
    16791682                                               surf_lsm_h%pt_surface,          &
     
    17041707!--            Output of instantaneous data
    17051708               IF ( av == 0 )  THEN
    1706                   CALL surface_output_collect( surf_def_h(0)%vpt_surface,      &
     1709                  CALL surface_data_output_collect( surf_def_h(0)%vpt_surface, &
    17071710                                               surf_def_h(1)%vpt_surface,      &
    17081711                                               surf_lsm_h%vpt_surface,         &
     
    17331736!--            Output of instantaneous data
    17341737               IF ( av == 0 )  THEN
    1735                   CALL surface_output_collect( surf_def_h(0)%rad_net,          &
     1738                  CALL surface_data_output_collect( surf_def_h(0)%rad_net,     &
    17361739                                               surf_def_h(1)%rad_net,          &
    17371740                                               surf_lsm_h%rad_net,             &
     
    17621765!--            Output of instantaneous data
    17631766               IF ( av == 0 )  THEN
    1764                   CALL surface_output_collect( surf_def_h(0)%rad_lw_in,        &
     1767                  CALL surface_data_output_collect( surf_def_h(0)%rad_lw_in,   &
    17651768                                               surf_def_h(1)%rad_lw_in,        &
    17661769                                               surf_lsm_h%rad_lw_in,           &
     
    17911794!--            Output of instantaneous data
    17921795               IF ( av == 0 )  THEN
    1793                   CALL surface_output_collect( surf_def_h(0)%rad_lw_out,       &
     1796                  CALL surface_data_output_collect( surf_def_h(0)%rad_lw_out,  &
    17941797                                               surf_def_h(1)%rad_lw_out,       &
    17951798                                               surf_lsm_h%rad_lw_out,          &
     
    18201823!--            Output of instantaneous data
    18211824               IF ( av == 0 )  THEN
    1822                   CALL surface_output_collect( surf_def_h(0)%rad_sw_in,        &
     1825                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_in,   &
    18231826                                               surf_def_h(1)%rad_sw_in,        &
    18241827                                               surf_lsm_h%rad_sw_in,           &
     
    18491852!--            Output of instantaneous data
    18501853               IF ( av == 0 )  THEN
    1851                   CALL surface_output_collect( surf_def_h(0)%rad_sw_out,       &
     1854                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_out,  &
    18521855                                               surf_def_h(1)%rad_sw_out,       &
    18531856                                               surf_lsm_h%rad_sw_out,          &
     
    18781881!--            Output of instantaneous data
    18791882               IF ( av == 0 )  THEN
    1880                   CALL surface_output_collect( surf_def_h(0)%ghf,              &
     1883                  CALL surface_data_output_collect( surf_def_h(0)%ghf,         &
    18811884                                               surf_def_h(1)%ghf,              &
    18821885                                               surf_lsm_h%ghf,                 &
     
    19071910!--            Output of instantaneous data
    19081911               IF ( av == 0 )  THEN
    1909                   CALL surface_output_collect( surf_def_h(0)%r_a,              &
     1912                  CALL surface_data_output_collect( surf_def_h(0)%r_a,         &
    19101913                                               surf_def_h(1)%r_a,              &
    19111914                                               surf_lsm_h%r_a,                 &
     
    19361939!--            Output of instantaneous data
    19371940               IF ( av == 0 )  THEN
    1938                   CALL surface_output_collect( surf_def_h(0)%r_soil,           &
     1941                  CALL surface_data_output_collect( surf_def_h(0)%r_soil,      &
    19391942                                               surf_def_h(1)%r_soil,           &
    19401943                                               surf_lsm_h%r_soil,              &
     
    19651968!--            Output of instantaneous data
    19661969               IF ( av == 0 )  THEN
    1967                   CALL surface_output_collect( surf_def_h(0)%r_canopy,         &
     1970                  CALL surface_data_output_collect( surf_def_h(0)%r_canopy,    &
    19681971                                               surf_def_h(1)%r_canopy,         &
    19691972                                               surf_lsm_h%r_canopy,            &
     
    19941997!--            Output of instantaneous data
    19951998               IF ( av == 0 )  THEN
    1996                   CALL surface_output_collect( surf_def_h(0)%r_s,              &
     1999                  CALL surface_data_output_collect( surf_def_h(0)%r_s,         &
    19972000                                               surf_def_h(1)%r_s,              &
    19982001                                               surf_lsm_h%r_s,                 &
     
    20232026!--            Output of instantaneous data
    20242027               IF ( av == 0 )  THEN
    2025                   CALL surface_output_collect( surf_def_h(0)%rad_sw_dir,       &
     2028                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_dir,  &
    20262029                                               surf_def_h(1)%rad_sw_dir,       &
    20272030                                               surf_lsm_h%rad_sw_dir,          &
     
    20522055!--            Output of instantaneous data
    20532056               IF ( av == 0 )  THEN
    2054                   CALL surface_output_collect( surf_def_h(0)%rad_sw_dif,       &
     2057                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_dif,  &
    20552058                                               surf_def_h(1)%rad_sw_dif,       &
    20562059                                               surf_lsm_h%rad_sw_dif,          &
     
    20812084!--            Output of instantaneous data
    20822085               IF ( av == 0 )  THEN
    2083                   CALL surface_output_collect( surf_def_h(0)%rad_sw_ref,       &
     2086                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_ref,  &
    20842087                                               surf_def_h(1)%rad_sw_ref,       &
    20852088                                               surf_lsm_h%rad_sw_ref,          &
     
    21102113!--            Output of instantaneous data
    21112114               IF ( av == 0 )  THEN
    2112                   CALL surface_output_collect( surf_def_h(0)%rad_sw_res,       &
     2115                  CALL surface_data_output_collect( surf_def_h(0)%rad_sw_res,  &
    21132116                                               surf_def_h(1)%rad_sw_res,       &
    21142117                                               surf_lsm_h%rad_sw_res,          &
     
    21392142!--            Output of instantaneous data
    21402143               IF ( av == 0 )  THEN
    2141                   CALL surface_output_collect( surf_def_h(0)%rad_lw_dif,       &
     2144                  CALL surface_data_output_collect( surf_def_h(0)%rad_lw_dif,  &
    21422145                                               surf_def_h(1)%rad_lw_dif,       &
    21432146                                               surf_lsm_h%rad_lw_dif,          &
     
    21682171!--            Output of instantaneous data
    21692172               IF ( av == 0 )  THEN
    2170                   CALL surface_output_collect( surf_def_h(0)%rad_lw_ref,       &
     2173                  CALL surface_data_output_collect( surf_def_h(0)%rad_lw_ref,  &
    21712174                                               surf_def_h(1)%rad_lw_ref,       &
    21722175                                               surf_lsm_h%rad_lw_ref,          &
     
    21972200!--            Output of instantaneous data
    21982201               IF ( av == 0 )  THEN
    2199                   CALL surface_output_collect( surf_def_h(0)%rad_lw_res,       &
     2202                  CALL surface_data_output_collect( surf_def_h(0)%rad_lw_res,  &
    22002203                                               surf_def_h(1)%rad_lw_res,       &
    22012204                                               surf_lsm_h%rad_lw_res,          &
     
    22542257      IF ( av == 1 )  average_count_surf = 0
    22552258             
    2256    END SUBROUTINE surface_output
     2259   END SUBROUTINE surface_data_output
    22572260   
    22582261!------------------------------------------------------------------------------!
     
    22612264!> Routine for controlling the data averaging.
    22622265!------------------------------------------------------------------------------!     
    2263    SUBROUTINE surface_output_averaging
     2266   SUBROUTINE surface_data_output_averaging
    22642267   
    22652268      IMPLICIT NONE
     
    22782281
    22792282            CASE ( 'us' )
    2280                CALL surface_output_sum_up( surf_def_h(0)%us,                   &
     2283               CALL surface_data_output_sum_up( surf_def_h(0)%us,              &
    22812284                                           surf_def_h(1)%us,                   &
    22822285                                           surf_lsm_h%us,                      &
     
    22962299
    22972300            CASE ( 'ts' )
    2298                CALL surface_output_sum_up( surf_def_h(0)%ts,                   &
     2301               CALL surface_data_output_sum_up( surf_def_h(0)%ts,              &
    22992302                                           surf_def_h(1)%ts,                   &
    23002303                                           surf_lsm_h%ts,                      &
     
    23142317
    23152318            CASE ( 'qs' )
    2316                CALL surface_output_sum_up( surf_def_h(0)%qs,                   &
     2319               CALL surface_data_output_sum_up( surf_def_h(0)%qs,              &
    23172320                                           surf_def_h(1)%qs,                   &
    23182321                                           surf_lsm_h%qs,                      &
     
    23322335
    23332336            CASE ( 'ss' )
    2334                CALL surface_output_sum_up( surf_def_h(0)%ss,                   &
     2337               CALL surface_data_output_sum_up( surf_def_h(0)%ss,              &
    23352338                                           surf_def_h(1)%ss,                   &
    23362339                                           surf_lsm_h%ss,                      &
     
    23502353
    23512354            CASE ( 'qcs' )
    2352                CALL surface_output_sum_up( surf_def_h(0)%qcs,                  &
     2355               CALL surface_data_output_sum_up( surf_def_h(0)%qcs,             &
    23532356                                           surf_def_h(1)%qcs,                  &
    23542357                                           surf_lsm_h%qcs,                     &
     
    23682371
    23692372            CASE ( 'ncs' )
    2370                CALL surface_output_sum_up( surf_def_h(0)%ncs,                  &
     2373               CALL surface_data_output_sum_up( surf_def_h(0)%ncs,             &
    23712374                                           surf_def_h(1)%ncs,                  &
    23722375                                           surf_lsm_h%ncs,                     &
     
    23862389
    23872390            CASE ( 'qrs' )
    2388                CALL surface_output_sum_up( surf_def_h(0)%qrs,                  &
     2391               CALL surface_data_output_sum_up( surf_def_h(0)%qrs,             &
    23892392                                           surf_def_h(1)%qrs,                  &
    23902393                                           surf_lsm_h%qrs,                     &
     
    24042407
    24052408            CASE ( 'nrs' )
    2406                CALL surface_output_sum_up( surf_def_h(0)%nrs,                  &
     2409               CALL surface_data_output_sum_up( surf_def_h(0)%nrs,             &
    24072410                                           surf_def_h(1)%nrs,                  &
    24082411                                           surf_lsm_h%nrs,                     &
     
    24222425
    24232426            CASE ( 'ol' )
    2424                CALL surface_output_sum_up( surf_def_h(0)%ol,                   &
     2427               CALL surface_data_output_sum_up( surf_def_h(0)%ol,              &
    24252428                                           surf_def_h(1)%ol,                   &
    24262429                                           surf_lsm_h%ol,                      &
     
    24402443
    24412444            CASE ( 'z0' )
    2442                CALL surface_output_sum_up( surf_def_h(0)%z0,                   &
     2445               CALL surface_data_output_sum_up( surf_def_h(0)%z0,              &
    24432446                                           surf_def_h(1)%z0,                   &
    24442447                                           surf_lsm_h%z0,                      &
     
    24582461
    24592462            CASE ( 'z0h' )
    2460                CALL surface_output_sum_up( surf_def_h(0)%z0h,                  &
     2463               CALL surface_data_output_sum_up( surf_def_h(0)%z0h,             &
    24612464                                           surf_def_h(1)%z0h,                  &
    24622465                                           surf_lsm_h%z0h,                     &
     
    24762479
    24772480            CASE ( 'z0q' )
    2478                CALL surface_output_sum_up( surf_def_h(0)%z0q,                  &
     2481               CALL surface_data_output_sum_up( surf_def_h(0)%z0q,             &
    24792482                                           surf_def_h(1)%z0q,                  &
    24802483                                           surf_lsm_h%z0q,                     &
     
    24942497
    24952498            CASE ( 'theta1' )                                                 
    2496                CALL surface_output_sum_up( surf_def_h(0)%pt1,                  &
     2499               CALL surface_data_output_sum_up( surf_def_h(0)%pt1,             &
    24972500                                           surf_def_h(1)%pt1,                  &
    24982501                                           surf_lsm_h%pt1,                     &
     
    25122515                                                                               
    25132516            CASE ( 'qv1' )                                                     
    2514                CALL surface_output_sum_up( surf_def_h(0)%qv1,                  &
     2517               CALL surface_data_output_sum_up( surf_def_h(0)%qv1,             &
    25152518                                           surf_def_h(1)%qv1,                  &
    25162519                                           surf_lsm_h%qv1,                     &
     
    25302533
    25312534            CASE ( 'thetav1' )
    2532                CALL surface_output_sum_up( surf_def_h(0)%vpt1,                 &
     2535               CALL surface_data_output_sum_up( surf_def_h(0)%vpt1,            &
    25332536                                           surf_def_h(1)%vpt1,                 &
    25342537                                           surf_lsm_h%vpt1,                    &
     
    25482551
    25492552            CASE ( 'usws' )
    2550                CALL surface_output_sum_up( surf_def_h(0)%usws,                 &
     2553               CALL surface_data_output_sum_up( surf_def_h(0)%usws,            &
    25512554                                           surf_def_h(1)%usws,                 &
    25522555                                           surf_lsm_h%usws,                    &
     
    25662569
    25672570            CASE ( 'vsws' )
    2568                CALL surface_output_sum_up( surf_def_h(0)%vsws,                 &
     2571               CALL surface_data_output_sum_up( surf_def_h(0)%vsws,            &
    25692572                                           surf_def_h(1)%vsws,                 &
    25702573                                           surf_lsm_h%vsws,                    &
     
    25842587
    25852588            CASE ( 'shf' )
    2586                CALL surface_output_sum_up( surf_def_h(0)%shf,                  &
     2589               CALL surface_data_output_sum_up( surf_def_h(0)%shf,             &
    25872590                                           surf_def_h(1)%shf,                  &
    25882591                                           surf_lsm_h%shf,                     &
     
    26022605
    26032606            CASE ( 'qsws' )
    2604                CALL surface_output_sum_up( surf_def_h(0)%qsws,                 &
     2607               CALL surface_data_output_sum_up( surf_def_h(0)%qsws,            &
    26052608                                           surf_def_h(1)%qsws,                 &
    26062609                                           surf_lsm_h%qsws,                    &
     
    26202623
    26212624            CASE ( 'ssws' )
    2622                CALL surface_output_sum_up( surf_def_h(0)%ssws,                 &
     2625               CALL surface_data_output_sum_up( surf_def_h(0)%ssws,            &
    26232626                                           surf_def_h(1)%ssws,                 &
    26242627                                           surf_lsm_h%ssws,                    &
     
    26382641
    26392642            CASE ( 'qcsws' )
    2640                CALL surface_output_sum_up( surf_def_h(0)%qcsws,                &
     2643               CALL surface_data_output_sum_up( surf_def_h(0)%qcsws,           &
    26412644                                           surf_def_h(1)%qcsws,                &
    26422645                                           surf_lsm_h%qcsws,                   &
     
    26562659
    26572660            CASE ( 'ncsws' )
    2658                CALL surface_output_sum_up( surf_def_h(0)%ncsws,                &
     2661               CALL surface_data_output_sum_up( surf_def_h(0)%ncsws,           &
    26592662                                           surf_def_h(1)%ncsws,                &
    26602663                                           surf_lsm_h%ncsws,                   &
     
    26742677
    26752678            CASE ( 'qrsws' )
    2676                CALL surface_output_sum_up( surf_def_h(0)%qrsws,                &
     2679               CALL surface_data_output_sum_up( surf_def_h(0)%qrsws,           &
    26772680                                           surf_def_h(1)%qrsws,                &
    26782681                                           surf_lsm_h%qrsws,                   &
     
    26922695
    26932696            CASE ( 'nrsws' )
    2694                CALL surface_output_sum_up( surf_def_h(0)%nrsws,                &
     2697               CALL surface_data_output_sum_up( surf_def_h(0)%nrsws,           &
    26952698                                           surf_def_h(1)%nrsws,                &
    26962699                                           surf_lsm_h%nrsws,                   &
     
    27102713
    27112714            CASE ( 'sasws' )
    2712                CALL surface_output_sum_up( surf_def_h(0)%sasws,                &
     2715               CALL surface_data_output_sum_up( surf_def_h(0)%sasws,           &
    27132716                                           surf_def_h(1)%sasws,                &
    27142717                                           surf_lsm_h%sasws,                   &
     
    27282731
    27292732            CASE ( 'q_surface' )
    2730                CALL surface_output_sum_up( surf_def_h(0)%q_surface,            &
     2733               CALL surface_data_output_sum_up( surf_def_h(0)%q_surface,       &
    27312734                                           surf_def_h(1)%q_surface,            &
    27322735                                           surf_lsm_h%q_surface,               &
     
    27462749
    27472750            CASE ( 'theta_surface' )
    2748                CALL surface_output_sum_up( surf_def_h(0)%pt_surface,           &
     2751               CALL surface_data_output_sum_up( surf_def_h(0)%pt_surface,      &
    27492752                                           surf_def_h(1)%pt_surface,           &
    27502753                                           surf_lsm_h%pt_surface,              &
     
    27642767
    27652768            CASE ( 'thetav_surface' )
    2766                CALL surface_output_sum_up( surf_def_h(0)%vpt_surface,          &
     2769               CALL surface_data_output_sum_up( surf_def_h(0)%vpt_surface,     &
    27672770                                           surf_def_h(1)%vpt_surface,          &
    27682771                                           surf_lsm_h%vpt_surface,             &
     
    27822785
    27832786            CASE ( 'rad_net' )
    2784                CALL surface_output_sum_up( surf_def_h(0)%rad_net,              &
     2787               CALL surface_data_output_sum_up( surf_def_h(0)%rad_net,         &
    27852788                                           surf_def_h(1)%rad_net,              &
    27862789                                           surf_lsm_h%rad_net,                 &
     
    28002803
    28012804            CASE ( 'rad_lw_in' )
    2802                CALL surface_output_sum_up( surf_def_h(0)%rad_lw_in,            &
     2805               CALL surface_data_output_sum_up( surf_def_h(0)%rad_lw_in,       &
    28032806                                           surf_def_h(1)%rad_lw_in,            &
    28042807                                           surf_lsm_h%rad_lw_in,               &
     
    28182821
    28192822            CASE ( 'rad_lw_out' )
    2820                CALL surface_output_sum_up( surf_def_h(0)%rad_lw_out,           &
     2823               CALL surface_data_output_sum_up( surf_def_h(0)%rad_lw_out,      &
    28212824                                           surf_def_h(1)%rad_lw_out,           &
    28222825                                           surf_lsm_h%rad_lw_out,              &
     
    28362839
    28372840            CASE ( 'rad_sw_in' )
    2838                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_in,            &
     2841               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_in,       &
    28392842                                           surf_def_h(1)%rad_sw_in,            &
    28402843                                           surf_lsm_h%rad_sw_in,               &
     
    28542857
    28552858            CASE ( 'rad_sw_out' )
    2856                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_out,           &
     2859               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_out,      &
    28572860                                           surf_def_h(1)%rad_sw_out,           &
    28582861                                           surf_lsm_h%rad_sw_out,              &
     
    28722875
    28732876            CASE ( 'ghf' )
    2874                CALL surface_output_sum_up( surf_def_h(0)%ghf,                  &
     2877               CALL surface_data_output_sum_up( surf_def_h(0)%ghf,             &
    28752878                                           surf_def_h(1)%ghf,                  &
    28762879                                           surf_lsm_h%ghf,                     &
     
    28902893                                                                               
    28912894            CASE ( 'r_a' )                                                     
    2892                CALL surface_output_sum_up( surf_def_h(0)%r_a,                  &
     2895               CALL surface_data_output_sum_up( surf_def_h(0)%r_a,             &
    28932896                                           surf_def_h(1)%r_a,                  &
    28942897                                           surf_lsm_h%r_a,                     &
     
    29082911                                                                               
    29092912            CASE ( 'r_soil' )                                                 
    2910                CALL surface_output_sum_up( surf_def_h(0)%r_soil,               &
     2913               CALL surface_data_output_sum_up( surf_def_h(0)%r_soil,          &
    29112914                                           surf_def_h(1)%r_soil,               &
    29122915                                           surf_lsm_h%r_soil,                  &
     
    29262929
    29272930            CASE ( 'r_canopy' )
    2928                CALL surface_output_sum_up( surf_def_h(0)%r_canopy,             &
     2931               CALL surface_data_output_sum_up( surf_def_h(0)%r_canopy,        &
    29292932                                           surf_def_h(1)%r_canopy,             &
    29302933                                           surf_lsm_h%r_canopy,                &
     
    29442947
    29452948            CASE ( 'r_s' )
    2946                CALL surface_output_sum_up( surf_def_h(0)%r_s,                  &
     2949               CALL surface_data_output_sum_up( surf_def_h(0)%r_s,             &
    29472950                                           surf_def_h(1)%r_s,                  &
    29482951                                           surf_lsm_h%r_s,                     &
     
    29632966
    29642967            CASE ( 'rad_sw_dir' )
    2965                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_dir,           &
     2968               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_dir,      &
    29662969                                           surf_def_h(1)%rad_sw_dir,           &
    29672970                                           surf_lsm_h%rad_sw_dir,              &
     
    29802983                                           surf_usm_v(3)%rad_sw_dir, n_out )
    29812984            CASE ( 'rad_sw_dif' )
    2982                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_dif,           &
     2985               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_dif,      &
    29832986                                           surf_def_h(1)%rad_sw_dif,           &
    29842987                                           surf_lsm_h%rad_sw_dif,              &
     
    29983001
    29993002            CASE ( 'rad_sw_ref' )
    3000                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_ref,           &
     3003               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_ref,      &
    30013004                                           surf_def_h(1)%rad_sw_ref,           &
    30023005                                           surf_lsm_h%rad_sw_ref,              &
     
    30163019               
    30173020            CASE ( 'rad_sw_res' )
    3018                CALL surface_output_sum_up( surf_def_h(0)%rad_sw_res,           &
     3021               CALL surface_data_output_sum_up( surf_def_h(0)%rad_sw_res,      &
    30193022                                           surf_def_h(1)%rad_sw_res,           &
    30203023                                           surf_lsm_h%rad_sw_res,              &
     
    30343037
    30353038            CASE ( 'rad_lw_dif' )
    3036                CALL surface_output_sum_up( surf_def_h(0)%rad_lw_dif,           &
     3039               CALL surface_data_output_sum_up( surf_def_h(0)%rad_lw_dif,      &
    30373040                                           surf_def_h(1)%rad_lw_dif,           &
    30383041                                           surf_lsm_h%rad_lw_dif,              &
     
    30523055
    30533056            CASE ( 'rad_lw_ref' )
    3054                CALL surface_output_sum_up( surf_def_h(0)%rad_lw_ref,           &
     3057               CALL surface_data_output_sum_up( surf_def_h(0)%rad_lw_ref,      &
    30553058                                           surf_def_h(1)%rad_lw_ref,           &
    30563059                                           surf_lsm_h%rad_lw_ref,              &
     
    30703073
    30713074            CASE ( 'rad_lw_res' )
    3072                CALL surface_output_sum_up( surf_def_h(0)%rad_lw_res,           &
     3075               CALL surface_data_output_sum_up( surf_def_h(0)%rad_lw_res,      &
    30733076                                           surf_def_h(1)%rad_lw_res,           &
    30743077                                           surf_lsm_h%rad_lw_res,              &
     
    30913094     
    30923095   
    3093    END SUBROUTINE surface_output_averaging
     3096   END SUBROUTINE surface_data_output_averaging
    30943097   
    30953098!------------------------------------------------------------------------------!
     
    30983101!> Sum-up the surface data for average output variables. 
    30993102!------------------------------------------------------------------------------!   
    3100    SUBROUTINE surface_output_sum_up( var_def_h0, var_def_h1,                   &
     3103   SUBROUTINE surface_data_output_sum_up( var_def_h0, var_def_h1,              &
    31013104                                     var_lsm_h,  var_usm_h,                    &
    31023105                                     var_def_v0, var_lsm_v0, var_usm_v0,       &
     
    32953298      ENDIF
    32963299   
    3297    END SUBROUTINE surface_output_sum_up
     3300   END SUBROUTINE surface_data_output_sum_up
    32983301   
    32993302!------------------------------------------------------------------------------!
     
    33023305!> Collect the surface data from different types and different orientation.
    33033306!------------------------------------------------------------------------------!
    3304    SUBROUTINE surface_output_collect( var_def_h0, var_def_h1,                  &
     3307   SUBROUTINE surface_data_output_collect( var_def_h0, var_def_h1,             &
    33053308                                      var_lsm_h,  var_usm_h,                   &
    33063309                                      var_def_v0, var_lsm_v0, var_usm_v0,      &
     
    34823485      ENDIF
    34833486   
    3484    END SUBROUTINE surface_output_collect
    3485    
    3486    
    3487    
    3488 !    SUBROUTINE data_output_surf_driver (time)
    3489 !       IMPLICIT NONE
    3490 !
    3491 !       REAL(kind=wp),INTENT(IN)    :: time
    3492 !       CHARACTER(LEN=64)           :: filename
    3493 !       INTEGER,SAVE                :: file_number=1
    3494 !
    3495 !       if(output_mode == "legacy_VTK")   then
    3496 !           CALL prepare_vtk_data
    3497 !
    3498 !           write(filename,'(a,i5.5,a)') 'surface_output_',file_number,'.vtk'
    3499 !           if(myid == surf_io_pe)   write(9,*) 'Legacy VTK surface output to file: ',trim(filename)
    3500 !           file_number = file_number+1
    3501 !
    3502 !           CALL write_ascii_vtk (filename, time)
    3503 !       END IF
    3504 !
    3505 !       RETURN
    3506 !    END SUBROUTINE data_output_surf_driver
    3507    
     3487   END SUBROUTINE surface_data_output_collect
    35083488   
    35093489!------------------------------------------------------------------------------!
     
    35123492!> Parin for output of surface parameters
    35133493!------------------------------------------------------------------------------!
    3514     SUBROUTINE surface_output_parin
     3494    SUBROUTINE surface_data_output_parin
    35153495
    35163496       IMPLICIT NONE
     
    35193499       
    35203500
    3521        NAMELIST /surface_output_parameters/                                    &
     3501       NAMELIST /surface_data_output_parameters/                               &
    35223502                                  averaging_interval_surf, data_output_surf,   &
    35233503                                  dt_dosurf, dt_dosurf_av,                     &
     
    35303510       REWIND ( 11 )
    35313511       line = ' '
    3532        DO WHILE ( INDEX( line, '&surface_output_parameters' ) == 0 )
     3512       DO WHILE ( INDEX( line, '&surface_data_output_parameters' ) == 0 )
    35333513          READ ( 11, '(A)', END=14 )  line
    35343514       ENDDO
     
    35373517!
    35383518!--    Read namelist
    3539        READ ( 11, surface_output_parameters, ERR = 10 )
     3519       READ ( 11, surface_data_output_parameters, ERR = 10 )
    35403520!
    35413521!--    Set flag that indicates that surface data output is switched on
    3542        surface_data_output = .TRUE.
     3522       surface_output = .TRUE.
    35433523       GOTO 14
    35443524
    35453525 10    BACKSPACE( 11 )
    35463526       READ( 11 , '(A)') line
    3547        CALL parin_fail_message( 'surface_output_parameters', line )     
     3527       CALL parin_fail_message( 'surface_data_output_parameters', line )     
    35483528       
    35493529 14    CONTINUE
    35503530       
    35513531
    3552     END SUBROUTINE surface_output_parin
     3532    END SUBROUTINE surface_data_output_parin
    35533533   
    35543534   
     
    35603540!> variables and map them onto internal output array.
    35613541!------------------------------------------------------------------------------!
    3562     SUBROUTINE surface_output_check_parameters
     3542    SUBROUTINE surface_data_output_check_parameters
    35633543
    35643544       USE control_parameters,                                                 &
     
    35873567                averaging_interval_surf, ' must be <= dt_dosurf_av = ',        &
    35883568                dt_dosurf_av
    3589           CALL message( 'surface_output_check_parameters',                     &
     3569          CALL message( 'surface_data_output_check_parameters',                &
    35903570                        'PA0087', 1, 2, 0, 6, 0 )
    35913571       ENDIF
     
    36213601               message_string = TRIM( trimvar ) //                             &
    36223602                             ' is not yet implemented in the surface output'
    3623                CALL message( 'surface_output_check_parameters',                &
     3603               CALL message( 'surface_data_output_check_parameters',           &
    36243604                             'PA0087', 1, 2, 0, 6, 0 )
    36253605
     
    36533633               message_string = TRIM( trimvar ) //                             &
    36543634                             ' is not part of the surface output'
    3655                CALL message( 'surface_output_check_parameters',                &
     3635               CALL message( 'surface_data_output_check_parameters',           &
    36563636                             'PA0087', 1, 2, 0, 6, 0 )
    36573637         END SELECT
     
    36593639       ENDDO
    36603640
    3661     END SUBROUTINE surface_output_check_parameters
     3641    END SUBROUTINE surface_data_output_check_parameters
    36623642   
    36633643   
     
    36673647!> Last action.
    36683648!------------------------------------------------------------------------------!
    3669     SUBROUTINE surface_output_last_action( av )
     3649    SUBROUTINE surface_data_output_last_action( av )
    36703650 
    36713651      USE control_parameters,                                                  &
     
    36983678      ENDDO   
    36993679
    3700     END SUBROUTINE surface_output_last_action   
    3701 
    3702 !--Private Subroutines (Only called from tinside his module)
    3703 
    3704 
    3705 
    3706 !    SUBROUTINE write_ascii_vtk (filename,time)
    3707 !       IMPLICIT NONE
    3708 !       CHARACTER(LEN=*),INTENT(IN)   :: filename
    3709 !       REAL(kind=wp),INTENT(IN)      :: time
    3710 !
    3711 !       INTEGER                       :: file_id=251
    3712 !       INTEGER                       :: n
    3713 !
    3714 !       IF(myid == surf_io_pe)  then
    3715 !          open(file_id,file=TRIM(filename))                                 ! Please use check_open in the final version
    3716 !
    3717 !          write(file_id,'(a)') "# vtk DataFile Version 3.0"
    3718 !          write(file_id,'(a,f8.2,a)') "legacy vtk File generated by PALM,  simulation time = ",time," sec"
    3719 !          write(file_id,'(a)') "ASCII"
    3720 !
    3721 !          write(file_id,'(a)') "DATASET POLYDATA"
    3722 !          write(file_id,'(a,i5,a)') "POINTS ",npo," float"
    3723 !          do n=1,npo
    3724 !             write (file_id,'(8f10.1)') points(1:3,n)
    3725 !          end do
    3726 !
    3727 !          write(file_id,'(a,8i10)') "POLYGONS ",npg, 5*npg
    3728 !          do n=1,npg
    3729 !             write(file_id,'(8i10)') polygons(:,n)
    3730 !          end do
    3731 !
    3732 !          write(file_id,'(a,i10)') "CELL_DATA ",npg
    3733 !          write(file_id,'(a,i10)') "SCALARS cell_scalars float 1 "
    3734 !          write(file_id,'(a,i10)') "LOOKUP_TABLE default "
    3735 !
    3736 !          do n=1,npg
    3737 !             write(file_id,'(f10.2)') cell_scalars(n)
    3738 !          end do
    3739 !
    3740 !          close (file_id)
    3741 !
    3742 !          DEALLOCATE (points,polygons,cell_scalars)
    3743 !       ENDIF
    3744 !
    3745 !       RETURN
    3746 !    END SUBROUTINE write_ascii_vtk
    3747 
    3748 END MODULE surface_output_mod
     3680    END SUBROUTINE surface_data_output_last_action   
     3681
     3682
     3683END MODULE surface_data_output_mod
  • palm/trunk/SOURCE/time_integration.f90

    r3647 r3648  
    2525! -----------------
    2626! $Id$
     27! Rename subroutines for surface-data output
     28!
     29! 3647 2019-01-02 14:10:44Z kanani
    2730! Bugfix: add time_since_reference_point to IF clause for data_output calls
    2831! (otherwise skip_time_* values don't come into affect with dt_do* = 0.0).
     
    488491               skip_time_do3d, skip_time_domask, skip_time_dopr,               &
    489492               skip_time_data_output_av, sloping_surface, stop_dt,             &
    490                surface_data_output, terminate_coupled, terminate_run,          &
     493               surface_output, terminate_coupled, terminate_run,               &
    491494               timestep_scheme,                                                &
    492495               time_coupling, time_do2d_xy, time_do2d_xz, time_do2d_yz,        &
     
    595598               enter_surface_arrays, exit_surface_arrays
    596599
    597     USE surface_output_mod,                                                    &
     600    USE surface_data_output_mod,                                               &
    598601        ONLY:  average_count_surf, averaging_interval_surf, dt_dosurf,         &
    599                dt_dosurf_av, surface_output, surface_output_averaging,         &
    600                skip_time_dosurf, skip_time_dosurf_av, time_dosurf,             &
    601                time_dosurf_av
     602               dt_dosurf_av, surface_data_output,                              &
     603               surface_data_output_averaging, skip_time_dosurf,                &
     604               skip_time_dosurf_av, time_dosurf, time_dosurf_av
    602605
    603606    USE turbulence_closure_mod,                                                &
     
    14661469!
    14671470!--    Increment time-counter for surface output
    1468        IF ( surface_data_output )  THEN
     1471       IF ( surface_output )  THEN
    14691472          IF ( time_since_reference_point >= skip_time_dosurf )  THEN
    14701473             time_dosurf    = time_dosurf + dt_3d
     
    15921595!
    15931596!--    Average surface data
    1594        IF ( surface_data_output )  THEN
     1597       IF ( surface_output )  THEN
    15951598          IF (         averaging_interval_surf /= 0.0_wp                               &
    15961599                .AND.  ( dt_dosurf_av - time_dosurf_av ) <= averaging_interval_surf    &
    15971600                .AND.  time_since_reference_point >= skip_time_dosurf_av )  THEN
    15981601             IF ( time_dosurf_av >= dt_averaging_input )  THEN       
    1599                 CALL surface_output_averaging
     1602                CALL surface_data_output_averaging
    16001603                average_count_surf = average_count_surf + 1
    16011604             ENDIF
     
    17141717!
    17151718!--    Output of surface data, instantaneous and averaged data
    1716        IF ( surface_data_output )  THEN
     1719       IF ( surface_output )  THEN
    17171720          IF (         time_dosurf >= dt_dosurf                                    &
    17181721                .AND.  time_since_reference_point >= skip_time_dosurf )  THEN
    1719              CALL surface_output( 0 )
     1722             CALL surface_data_output( 0 )
    17201723             time_dosurf = MOD( time_dosurf, MAX( dt_dosurf, dt_3d ) )
    17211724          ENDIF
    17221725          IF (         time_dosurf_av >= dt_dosurf_av                              &
    17231726                .AND.  time_since_reference_point >= skip_time_dosurf_av )  THEN
    1724              CALL surface_output( 1 )
     1727             CALL surface_data_output( 1 )
    17251728             time_dosurf_av = MOD( time_dosurf_av, MAX( dt_dosurf_av, dt_3d ) )
    17261729          ENDIF
Note: See TracChangeset for help on using the changeset viewer.