! !> @file palm.f90 !------------------------------------------------------------------------------! ! This file is part of the PALM model system. ! ! PALM is free software: you can redistribute it and/or modify it under the ! terms of the GNU General Public License as published by the Free Software ! Foundation, either version 3 of the License, or (at your option) any later ! version. ! ! PALM is distributed in the hope that it will be useful, but WITHOUT ANY ! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR ! A PARTICULAR PURPOSE. See the GNU General Public License for more details. ! ! You should have received a copy of the GNU General Public License along with ! PALM. If not, see . ! ! Copyright 1997-2019 Leibniz Universitaet Hannover !------------------------------------------------------------------------------! ! ! Current revisions: ! ----------------- ! ! ! Former revisions: ! ----------------- ! $Id: palm.f90 4017 2019-06-06 12:16:46Z raasch $ ! new module for calculation and output of diagnostic quantities added ! ! 3885 2019-04-11 11:29:34Z kanani ! Changes related to global restructuring of location messages and introduction ! of additional debug messages ! ! 3761 2019-02-25 15:31:42Z raasch ! unused variable removed ! ! 3719 2019-02-06 13:10:18Z kanani ! Included cpu measurement for wall/soil spinup ! ! 3703 2019-01-29 16:43:53Z knoop ! Some interface calls moved to module_interface + cleanup ! ! 3648 2019-01-02 16:35:46Z suehring ! Rename subroutines for surface-data output ! ! 3524 2018-11-14 13:36:44Z raasch ! unused variable removed ! ! 3494 2018-11-06 14:51:27Z suehring ! Last actions for surface output added ! ! 3487 2018-11-05 07:18:02Z maronga ! Updates version to 6.0 ! ! 3484 2018-11-02 14:41:25Z hellstea ! pmci_ensure_nest_mass_conservation removed permanently ! ! 3458 2018-10-30 14:51:23Z kanani ! from chemistry branch r3443, forkel: ! removed double do_emis check around CALL chem_init ! replaced call to calc_date_and_time to init_date_and_time ! ! 3337 2018-10-12 15:17:09Z kanani ! (from branch resler) ! Fix chemistry call ! ! 3298 2018-10-02 12:21:11Z kanani ! - Minor formatting (kanani) ! - Added Call of date_and_time_init (Russo) ! - Added Call of calc_date_and_time before call of init_3d where emissions ! are initialized: ! we have to know the time indices to initialize emission values (Russo) ! - Added Call of netcdf_data_input_chemistry_data (Russo) ! ! 3274 2018-09-24 15:42:55Z knoop ! Modularization of all bulk cloud physics code components ! ! 3258 2018-09-18 07:23:31Z Giersch ! current revision for user interface has been changed to 3240 ! ! 3241 2018-09-12 15:02:00Z raasch ! unused variables removed ! ! 3235 2018-09-07 14:06:15Z sward ! Added mas_last_actions call and multi_agent_system_mod dependency ! ! 3232 2018-09-07 12:21:44Z raasch ! references to mrun replaced by palmrun, and updated ! ! 3182 2018-07-27 13:36:03Z suehring ! Deduct spinup_time from RUN_CONTROL output of main 3d run ! (use time_since_reference_point instead of simulated_time) ! ! 2951 2018-04-06 09:05:08Z kanani ! Add log_point_s for pmci_init ! ! 2903 2018-03-16 08:17:06Z hellstea ! Nesting-related calls to pmci_ensure_nest_mass_conservation and pres after ! the nest initialization are removed as they may create unwanted initial ! perturbation in some cases. ! ! 2894 2018-03-15 09:17:58Z Giersch ! Modified todo list, _write_restart_data/_last_actions has been renamed to ! _wrd_local, unit 14 will be opened now for each io_group ! write_3d_binary is called wrd_local now, wrd_global moved from wrd_local to ! palm.f90, unit 14 is closed directly after the wrd_local call, Module related ! routines for writing restart data have been moved to wrd_local ! ! 2801 2018-02-14 16:01:55Z suehring ! Changed lpm from subroutine to module. ! Introduce particle transfer in nested models. ! ! 2766 2018-01-22 17:17:47Z kanani ! Removed preprocessor directive __chem ! ! 2720 2018-01-02 16:27:15Z kanani ! Version update to 5.0 ! ! 2718 2018-01-02 08:49:38Z maronga ! Corrected "Former revisions" section ! ! 2696 2017-12-14 17:12:51Z kanani ! Change in file header (GPL part) ! Implementation of chemistry module (FK) ! Introduce input-data standard ! Rename lsm_last_actions into lsm_write_restart_data ! Move usm_write_restart_data into io_blocks loop (MS) ! ! 2512 2017-10-04 08:26:59Z raasch ! user interface required revision updated ! ! 2320 2017-07-21 12:47:43Z suehring ! Modularize large-scale forcing and nudging ! ! 2298 2017-06-29 09:28:18Z raasch ! type of write_binary changed from CHARACTER to LOGICAL, ! user interface required revision updated, MPI2 related part removed ! ! 2296 2017-06-28 07:53:56Z maronga ! Added call to new spinup routine ! ! 2292 2017-06-20 09:51:42Z schwenkel ! Implementation of new microphysic scheme: cloud_scheme = 'morrison' ! includes two more prognostic equations for cloud drop concentration (nc) ! and cloud water content (qc). ! ! 2261 2017-06-08 14:25:57Z raasch ! output of run number for mrun to create unified cycle numbers ! ! 2233 2017-05-30 18:08:54Z suehring ! ! 2232 2017-05-30 17:47:52Z suehring ! Renamed wall_flags_0 and wall_flags_00 into advc_flags_1 and advc_flags_2, ! respectively, within copyin statement. Moreover, introduced further flag ! array wall_flags_0. ! Remove unused variables from ONLY list. ! ! 2178 2017-03-17 11:07:39Z hellstea ! Calls for pmci_ensure_nest_mass_conservation and pres are added after ! the nest initialization ! ! 2118 2017-01-17 16:38:49Z raasch ! OpenACC directives and related code removed ! ! 2011 2016-09-19 17:29:57Z kanani ! Flag urban_surface is now defined in module control_parameters. ! ! 2007 2016-08-24 15:47:17Z kanani ! Temporarily added CALL for writing of restart data for urban surface model ! ! 2000 2016-08-20 18:09:15Z knoop ! Forced header and separation lines into 80 columns ! ! 1976 2016-07-27 13:28:04Z maronga ! Added call to radiation_last_actions for binary output of land surface model ! data ! ! 1972 2016-07-26 07:52:02Z maronga ! Added call to lsm_last_actions for binary output of land surface model data ! ! 1960 2016-07-12 16:34:24Z suehring ! Separate humidity and passive scalar ! ! 1834 2016-04-07 14:34:20Z raasch ! Initial version of purely vertical nesting introduced. ! ! 1833 2016-04-07 14:23:03Z raasch ! required user interface version changed ! ! 1808 2016-04-05 19:44:00Z raasch ! routine local_flush replaced by FORTRAN statement ! ! 1783 2016-03-06 18:36:17Z raasch ! required user interface version changed ! ! 1781 2016-03-03 15:12:23Z raasch ! pmc initialization moved from time_integration to here ! ! 1779 2016-03-03 08:01:28Z raasch ! setting of nest_domain and coupling_char moved to the pmci ! ! 1764 2016-02-28 12:45:19Z raasch ! cpp-statements for nesting removed, communicator settings cleaned up ! ! 1762 2016-02-25 12:31:13Z hellstea ! Introduction of nested domain feature ! ! 1747 2016-02-08 12:25:53Z raasch ! OpenACC-adjustment for new surface layer parameterization ! ! 1682 2015-10-07 23:56:08Z knoop ! Code annotations made doxygen readable ! ! 1668 2015-09-23 13:45:36Z raasch ! warning replaced by abort in case of failed user interface check ! ! 1666 2015-09-23 07:31:10Z raasch ! check for user's interface version added ! ! 1482 2014-10-18 12:34:45Z raasch ! adjustments for using CUDA-aware OpenMPI ! ! 1468 2014-09-24 14:06:57Z maronga ! Adapted for use on up to 6-digit processor cores ! ! 1402 2014-05-09 14:25:13Z raasch ! location messages added ! ! 1374 2014-04-25 12:55:07Z raasch ! bugfix: various modules added ! ! 1320 2014-03-20 08:40:49Z raasch ! ONLY-attribute added to USE-statements, ! kind-parameters added to all INTEGER and REAL declaration statements, ! kinds are defined in new module kinds, ! old module precision_kind is removed, ! revision history before 2012 removed, ! comment fields (!:) to be used for variable explanations added to ! all variable declaration statements ! ! 1318 2014-03-17 13:35:16Z raasch ! module interfaces removed ! ! 1241 2013-10-30 11:36:58Z heinze ! initialization of nuding and large scale forcing from external file ! ! 1221 2013-09-10 08:59:13Z raasch ! +wall_flags_00, rflags_invers, rflags_s_inner in copyin statement ! ! 1212 2013-08-15 08:46:27Z raasch ! +tri in copyin statement ! ! 1179 2013-06-14 05:57:58Z raasch ! ref_state added to copyin-list ! ! 1113 2013-03-10 02:48:14Z raasch ! openACC statements modified ! ! 1111 2013-03-08 23:54:10Z raasch ! openACC statements updated ! ! 1092 2013-02-02 11:24:22Z raasch ! unused variables removed ! ! 1036 2012-10-22 13:43:42Z raasch ! code put under GPL (PALM 3.9) ! ! 1015 2012-09-27 09:23:24Z raasch ! Version number changed from 3.8 to 3.8a. ! OpenACC statements added + code changes required for GPU optimization ! ! 849 2012-03-15 10:35:09Z raasch ! write_particles renamed lpm_write_restart_file ! ! Revision 1.1 1997/07/24 11:23:35 raasch ! Initial revision ! ! ! Description: ! ------------ !> Large-Eddy Simulation (LES) model for atmospheric and oceanic boundary-layer !> flows !> see the PALM homepage https://palm-model.org for further information !------------------------------------------------------------------------------! PROGRAM palm USE arrays_3d USE bulk_cloud_model_mod, & ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert USE control_parameters, & ONLY: constant_diffusion, child_domain, & coupling_char, do2d_at_begin, do3d_at_begin, humidity, & initializing_actions, io_blocks, io_group, message_string, & neutral, passive_scalar, runnr, simulated_time_chr, spinup, & time_since_reference_point, user_interface_current_revision, & user_interface_required_revision, version, write_binary USE cpulog, & ONLY: cpu_log, log_point, log_point_s, cpu_statistics USE date_and_time_mod, & ONLY: calc_date_and_time, init_date_and_time USE diagnostic_output_quantities_mod, & ONLY: diagnostic_output_quantities_calculate USE indices, & ONLY: nbgp USE kinds USE module_interface, & ONLY: module_interface_last_actions USE multi_agent_system_mod, & ONLY: agents_active, mas_last_actions USE netcdf_data_input_mod, & ONLY: netcdf_data_input_inquire_file, netcdf_data_input_init, & netcdf_data_input_surface_data, netcdf_data_input_topo USE pegrid USE pmc_particle_interface, & ONLY: pmcp_g_alloc_win USE pmc_interface, & ONLY: nested_run, pmci_child_initialize, pmci_init, & pmci_modelconfiguration, pmci_parent_initialize USE surface_data_output_mod, & ONLY: surface_data_output_last_action USE write_restart_data_mod, & ONLY: wrd_global, wrd_local #if defined( __parallel) && defined( _OPENACC ) USE openacc #endif IMPLICIT NONE ! !-- Local variables CHARACTER(LEN=9) :: time_to_string !< INTEGER(iwp) :: i !< loop counter for blocked I/O #if defined( __parallel) && defined( _OPENACC ) INTEGER(iwp) :: local_comm !< local communicator (shared memory) INTEGER(iwp) :: local_num_procs !< local number of processes INTEGER(iwp) :: local_id !< local id INTEGER(acc_device_kind) :: device_type !< device type for OpenACC INTEGER(iwp) :: num_devices !< number of devices visible to OpenACC INTEGER(iwp) :: my_device !< device used by this process #endif version = 'PALM 6.0' user_interface_required_revision = 'r3703' #if defined( __parallel ) ! !-- MPI initialisation. comm2d is preliminary set, because !-- it will be defined in init_pegrid but is used before in cpu_log. CALL MPI_INIT( ierr ) ! !-- Initialize the coupling for nested-domain runs !-- comm_palm is the communicator which includes all PEs (MPI processes) !-- available for this (nested) model. If it is not a nested run, comm_palm !-- is returned as MPI_COMM_WORLD CALL cpu_log( log_point_s(70), 'pmci_init', 'start' ) CALL pmci_init( comm_palm ) CALL cpu_log( log_point_s(70), 'pmci_init', 'stop' ) comm2d = comm_palm ! !-- Get the (preliminary) number of MPI processes and the local PE-id (in case !-- of a further communicator splitting in init_coupling, these numbers will !-- be changed in init_pegrid). IF ( nested_run ) THEN CALL MPI_COMM_SIZE( comm_palm, numprocs, ierr ) CALL MPI_COMM_RANK( comm_palm, myid, ierr ) ELSE CALL MPI_COMM_SIZE( MPI_COMM_WORLD, numprocs, ierr ) CALL MPI_COMM_RANK( MPI_COMM_WORLD, myid, ierr ) ! !-- Initialize PE topology in case of coupled atmosphere-ocean runs (comm_palm !-- will be splitted in init_coupling) CALL init_coupling ENDIF #ifdef _OPENACC ! !-- Select OpenACC device to use in this process. For this find out how many !-- neighbors there are running on the same node and which id this process is. IF ( nested_run ) THEN CALL MPI_COMM_SPLIT_TYPE( comm_palm, MPI_COMM_TYPE_SHARED, 0, & MPI_INFO_NULL, local_comm, ierr ) ELSE CALL MPI_COMM_SPLIT_TYPE( MPI_COMM_WORLD, MPI_COMM_TYPE_SHARED, 0, & MPI_INFO_NULL, local_comm, ierr ) ENDIF CALL MPI_COMM_SIZE( local_comm, local_num_procs, ierr ) CALL MPI_COMM_RANK( local_comm, local_id, ierr ) ! !-- This loop including the barrier is a workaround for PGI compiler versions !-- up to and including 18.4. Later releases are able to select their GPUs in !-- parallel, without running into spurious errors. DO i = 0, local_num_procs-1 CALL MPI_BARRIER( local_comm, ierr ) IF ( i == local_id ) THEN device_type = acc_get_device_type() num_devices = acc_get_num_devices( device_type ) my_device = MOD( local_id, num_devices ) CALL acc_set_device_num( my_device, device_type ) ENDIF ENDDO CALL MPI_COMM_FREE( local_comm, ierr ) #endif #endif ! !-- Initialize measuring of the CPU-time remaining to the run CALL local_tremain_ini ! !-- Start of total CPU time measuring. CALL cpu_log( log_point(1), 'total', 'start' ) CALL cpu_log( log_point(2), 'initialisation', 'start' ) ! !-- Open a file for debug output WRITE (myid_char,'(''_'',I6.6)') myid OPEN( 9, FILE='DEBUG'//TRIM( coupling_char )//myid_char, FORM='FORMATTED' ) ! !-- Initialize dvrp logging. Also, one PE maybe split from the global !-- communicator for doing the dvrp output. In that case, the number of !-- PEs available for PALM is reduced by one and communicator comm_palm !-- is changed respectively. #if defined( __parallel ) CALL MPI_COMM_RANK( comm_palm, myid, ierr ) #endif ! !-- Read control parameters from NAMELIST files and read environment-variables CALL parin ! !-- Check for the user's interface version IF ( user_interface_current_revision /= user_interface_required_revision ) & THEN message_string = 'current user-interface revision "' // & TRIM( user_interface_current_revision ) // '" does ' // & 'not match the required revision ' // & TRIM( user_interface_required_revision ) CALL message( 'palm', 'PA0169', 1, 2, 0, 6, 0 ) ENDIF ! !-- Determine processor topology and local array indices CALL init_pegrid ! !-- Check if input file according to input-data standard exists CALL netcdf_data_input_inquire_file ! !-- Read topography input data if required. This is required before the !-- numerical grid is finally created in init_grid CALL netcdf_data_input_topo ! !-- Generate grid parameters, initialize generic topography and further process !-- topography information if required CALL init_grid ! !-- Read global attributes if available. CALL netcdf_data_input_init ! !-- Read surface classification data, e.g. vegetation and soil types, water !-- surfaces, etc., if available. Some of these data is required before !-- check parameters is invoked. CALL netcdf_data_input_surface_data ! !-- Check control parameters and deduce further quantities CALL check_parameters ! !-- Initial time for chem_emissions_mod CALL init_date_and_time CALL init_3d_model ! !-- Coupling protocol setup for nested-domain runs IF ( nested_run ) THEN CALL pmci_modelconfiguration ! !-- Receive and interpolate initial data on children. !-- Child initialization must be made first if the model is both child and !-- parent if necessary IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN CALL pmci_child_initialize ! !-- Send initial condition data from parent to children CALL pmci_parent_initialize ! !-- Exchange_horiz is needed after the nest initialization IF ( child_domain ) THEN CALL exchange_horiz( u, nbgp ) CALL exchange_horiz( v, nbgp ) CALL exchange_horiz( w, nbgp ) IF ( .NOT. neutral ) THEN CALL exchange_horiz( pt, nbgp ) ENDIF IF ( .NOT. constant_diffusion ) CALL exchange_horiz( e, nbgp ) IF ( humidity ) THEN CALL exchange_horiz( q, nbgp ) IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN CALL exchange_horiz( qc, nbgp ) CALL exchange_horiz( nc, nbgp ) ENDIF IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN CALL exchange_horiz( qr, nbgp ) CALL exchange_horiz( nr, nbgp ) ENDIF ENDIF IF ( passive_scalar ) CALL exchange_horiz( s, nbgp ) ENDIF ENDIF CALL pmcp_g_alloc_win ! Must be called after pmci_child_initialize and pmci_parent_initialize ENDIF ! !-- Output of program header IF ( myid == 0 ) CALL header CALL cpu_log( log_point(2), 'initialisation', 'stop' ) ! !-- Integration of the non-atmospheric equations (land surface model, urban !-- surface model) IF ( spinup ) THEN CALL cpu_log( log_point(41), 'wall/soil spinup', 'start' ) CALL time_integration_spinup CALL cpu_log( log_point(41), 'wall/soil spinup', 'stop' ) ENDIF ! !-- Set start time in format hh:mm:ss simulated_time_chr = time_to_string( time_since_reference_point ) ! !-- If required, output of initial arrays IF ( do2d_at_begin ) THEN CALL diagnostic_output_quantities_calculate !TODO, will be called twice CALL data_output_2d( 'xy', 0 ) CALL data_output_2d( 'xz', 0 ) CALL data_output_2d( 'yz', 0 ) ENDIF IF ( do3d_at_begin ) THEN CALL diagnostic_output_quantities_calculate !TODO, will be called twice CALL data_output_3d( 0 ) ENDIF ! !-- Integration of the model equations using timestep-scheme CALL time_integration ! !-- If required, write binary data for restart runs IF ( write_binary ) THEN CALL cpu_log( log_point(22), 'wrd_local', 'start' ) CALL location_message( 'writing restart data', 'start' ) DO i = 0, io_blocks-1 IF ( i == io_group ) THEN ! !-- Open binary file CALL check_open( 14 ) ! !-- Write control parameters and other global variables for restart. IF ( myid == 0 ) CALL wrd_global ! !-- Write processor specific flow field data for restart runs CALL wrd_local ! !-- Close binary file CALL close_file( 14 ) ENDIF #if defined( __parallel ) CALL MPI_BARRIER( comm2d, ierr ) #endif ENDDO CALL location_message( 'writing restart data', 'finished' ) CALL cpu_log( log_point(22), 'wrd_local', 'stop' ) ENDIF ! !-- Last actions for surface output, for instantaneous and time-averaged data CALL surface_data_output_last_action( 0 ) CALL surface_data_output_last_action( 1 ) ! !-- If required, repeat output of header including the required CPU-time IF ( myid == 0 ) CALL header ! !-- Perform module specific last actions CALL cpu_log( log_point(4), 'last actions', 'start' ) IF ( myid == 0 .AND. agents_active ) CALL mas_last_actions ! ToDo: move to module_interface CALL module_interface_last_actions CALL cpu_log( log_point(4), 'last actions', 'stop' ) ! !-- Close files CALL close_file( 0 ) ! !-- Write run number to file (used by palmrun to create unified cycle numbers !-- for output files IF ( myid == 0 .AND. runnr > 0 ) THEN OPEN( 90, FILE='RUN_NUMBER', FORM='FORMATTED' ) WRITE( 90, '(I4)' ) runnr CLOSE( 90 ) ENDIF ! !-- Take final CPU-time for CPU-time analysis CALL cpu_log( log_point(1), 'total', 'stop' ) CALL cpu_statistics #if defined( __parallel ) CALL MPI_FINALIZE( ierr ) #endif END PROGRAM palm