Changeset 3 for palm/trunk/SOURCE


Ignore:
Timestamp:
Feb 13, 2007 11:30:58 AM (17 years ago)
Author:
raasch
Message:

RCS Log replace by Id keyword, revision history cleaned up

Location:
palm/trunk/SOURCE
Files:
106 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r1 r3  
    88# Former revisions:
    99# -----------------
    10 # $Log: Makefile,v $
     10# $Id $
     11# RCS Log replace by Id keyword, revision history cleaned up
     12#
    1113# Revision 1.13  2006/08/04 14:45:31  raasch
    1214# +data_output_ptseries
    13 #
    14 # Revision 1.12  2006/04/26 15:10:49  raasch
    15 # modules.o added to dependencies of user_interface.o
    16 #
    17 # Revision 1.11  2006/02/23 13:17:11  raasch
    18 # +average_3d_data, data_log, sum_up_3d_data,
    19 # plot_.. renamed data_output_..
    20 #
    21 # Revision 1.10  2005/06/26 20:00:15  raasch
    22 # +interaction_droplets_ptq
    23 #
    24 # Revision 1.9  2005/05/18 14:52:20  raasch
    25 # +netcdf
    26 #
    27 # Revision 1.8  2004/04/30 13:15:46  raasch
    28 # dependencies for calc_spectra changed, +read_3d_binary,
    29 # leap_frog changed to time_integration, fft_for_1d_decomp to fft_xy,
    30 # init_pegrid and pres now depend on poisfft, all transpose routines joined
    31 # to one file transpose, local_getenv depends on modules
    32 #
    33 # Revision 1.7  2004/01/28 14:58:55  raasch
    34 # +timestep_scheme_steering
    35 #
    36 # Revision 1.6  2003/10/29 08:21:40  raasch
    37 # advec_particles, disturb_field, init_particles, init_3d_model, random_gauss
    38 # and write_3d_binary are using module random_function_mod
    39 #
    40 # Revision 1.5  2003/05/09 14:04:39  raasch
    41 # poisfft depends on fft_for_1d_decomp and no more on singleton,
    42 # init_pegrid depends on fft_for_1d_decomp, temperton_fft depends on modules
    43 #
    44 # Revision 1.4  2003/03/13 08:46:52  raasch
    45 # fft_for_1d_decomp has to depend on temperton_fft
    46 #
    47 # Revision 1.3  2003/03/12 15:57:37  raasch
    48 # modules and user_interface now depend on the f90 files, dependency of
    49 # singleton added, LDFLAGS moved to the end of PROG rule, +temperton_fft
    50 #
    51 # Revision 1.2  2003/01/06 08:48:19  raasch
    52 # +local_stop, check_cpu_time renamed to check_for_restart
    5315#
    5416# Revision 1.1  2002/06/11 13:42:20  raasch
  • palm/trunk/SOURCE/advec_particles.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_particles.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.32  2007/02/11 12:48:20  raasch
    1214! Allways the lower level k is used for interpolation
  • palm/trunk/SOURCE/advec_s_bc.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_s_bc.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.22  2006/02/23 09:42:08  raasch
    1214! anz renamed ngp
    13 !
    14 ! Revision 1.21  2003/03/16 09:22:13  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.20  2002/12/19 13:41:50  raasch
    18 ! STOP statement replaced by call of subroutine local_stop
    19 !
    20 ! Revision 1.19  2001/03/29 17:29:52  raasch
    21 ! Result is stored as tendency on tend, sk_p is now a local array, rho is
    22 ! replaced by d
    23 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    24 !
    25 ! Revision 1.18  2001/01/22 04:25:52  raasch
    26 ! Module test_variables removed
    27 !
    28 ! Revision 1.17  2000/07/03 12:54:35  raasch
    29 ! Arguments rho, sk and sk_p now declared as pointers
    30 !
    31 ! Revision 1.16  2000/04/27 07:03:07  raasch
    32 ! Treatment at cyclic boundaries is simplified by adding one additional
    33 ! ghost layer to tend at each cyclic boundary. Result is that most of
    34 ! mpi_sendrecv-calls (for polynom-coefficients, switches and fluxes) are no
    35 ! more necessary so that much cpu-time is spared.
    36 ! Adjustment of cyclic boundaries in x-direction in case of a sloping surface
    37 ! (temperature only).
    38 ! Vertical boundary conditions for humidity added.
    39 !
    40 ! Revision 1.15  2000/01/20 09:07:01  letzel
    41 ! All comments translated into English
    42 !
    43 ! Revision 1.14  1999/02/17 09:01:38  raasch
    44 ! Verfahren erweitert fuer TKE-Advektion
    45 ! Dichteberechnung (Korrekturverfahren) korrigiert (Galilei-Abzug war falsch)
    46 !
    47 ! Revision 1.10  1998/04/21 15:51:34  raasch
    48 ! Galilei-Transformation eingebaut
    49 !
    50 ! Revision 1.4  1997/09/08 10:42:35  raasch
    51 ! Optimierungen im Exponentialteil, MIN/MAX-Bildung aus progn. Gleichung
    52 ! entfernt, sw und m1 auf t3e einfach genau definiert, ffmax aus progn.
    53 ! Gleichung entfernt (Skalarfeld wird glatter im gradientschwachen Bereich)
    54 ! Geschwindigkeitsgewinn: ca 32%
    5515!
    5616! Revision 1.1  1997/08/29 08:53:46  raasch
  • palm/trunk/SOURCE/advec_s_pw.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_s_pw.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.12  2006/02/23 09:42:55  raasch
    1214! nzb_2d replaced by nzb_s_inner
    13 !
    14 ! Revision 1.11  2004/01/30 10:10:03  raasch
    15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    16 !
    17 ! Revision 1.10  2003/03/12 16:15:32  raasch
    18 ! Full code replaced in the call for all gridpoints instead of calling the
    19 ! _ij version (required by NEC, because otherwise no vectorization)
    20 !
    21 ! Revision 1.9  2002/06/11 12:22:14  raasch
    22 ! Former subroutine changed to a module which allows to be called for all grid
    23 ! points of a single vertical column with index i,j or for all grid points by
    24 ! using function overloading.
    25 ! Array sk now declared as a pointer.
    26 !
    27 ! Revision 1.8  2001/03/29 17:30:14  raasch
    28 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    29 !
    30 ! Revision 1.7  2001/01/22 04:41:39  raasch
    31 ! Module test_variables removed
    32 !
    33 ! Revision 1.6  2000/01/20 09:13:33  letzel
    34 ! All comments translated into English
    35 !
    36 ! Revision 1.5  1998/07/06 12:03:09  raasch
    37 ! + USE test_variables
    38 !
    39 ! Revision 1.4  1998/04/21 15:51:53  raasch
    40 ! Galilei-Transformation eingebaut
    41 !
    42 ! Revision 1.3  1998/02/10 14:57:53  raasch
    43 ! Es wird jetzt grundsaetzlich die Form C3 verwendet
    44 !
    45 ! Revision 1.2  1997/09/09 08:27:10  raasch
    46 ! Kehrwerte der Gitterweiten implementiert
    4715!
    4816! Revision 1.1  1997/08/29 08:54:20  raasch
  • palm/trunk/SOURCE/advec_s_up.f90

    r1 r3  
    88! Forner revisions:
    99! -----------------
    10 ! $Log: advec_s_up.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2006/02/23 09:43:44  raasch
    1214! nzb_2d replaced by nzb_s_inner
    13 !
    14 ! Revision 1.10  2004/01/30 10:10:39  raasch
    15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    16 !
    17 ! Revision 1.9  2003/03/12 16:16:25  raasch
    18 ! Full code replaced in the call for all gridpoints instead of calling the
    19 ! _ij version (required by NEC, because otherwise no vectorization)
    20 !
    21 ! Revision 1.8  2002/06/11 12:23:31  raasch
    22 ! Former subroutine changed to a module which allows to be called for all grid
    23 ! points of a single vertical column with index i,j or for all grid points by
    24 ! using function overloading.
    25 ! Array sk now declared as a pointer.
    26 !
    27 ! Revision 1.7  2001/03/29 17:30:31  raasch
    28 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    29 !
    30 ! Revision 1.6  2001/01/22 04:43:56  raasch
    31 ! Module test_variables removed
    32 !
    33 ! Revision 1.5  2000/01/20 09:16:59  letzel
    34 ! All comments translated into English
    35 !
    36 ! Revision 1.4  1998/07/06 12:03:46  raasch
    37 ! + USE test_variables
    38 !
    39 ! Revision 1.3  1998/04/21 15:52:14  raasch
    40 ! Galilei-Transformation eingebaut
    41 !
    42 ! Revision 1.2  1997/09/09 08:27:24  raasch
    43 ! Kehrwerte der Gitterweiten implementiert
    4415!
    4516! Revision 1.1  1997/08/29 08:54:33  raasch
  • palm/trunk/SOURCE/advec_s_ups.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_s_ups.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.6  2004/04/30 08:02:43  raasch
    1214! Enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.5  2003/03/16 09:23:56  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.4  2001/03/29 17:34:05  raasch
    18 ! v_ad is now allocated locally, several temporary arrays removed from
    19 ! argument lists of spline_* and long_filter_*
    20 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    21 !
    22 ! Revision 1.3  2001/01/22 04:46:10  raasch
    23 ! Module test_variables removed
    24 !
    25 ! Revision 1.2  2000/01/20 09:28:52  letzel
    26 ! All comments translated into English
    2715!
    2816! Revision 1.1  1999/02/05 08:44:47  raasch
  • palm/trunk/SOURCE/advec_u_pw.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_u_pw.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 09:44:21  raasch
    1214! nzb_2d replaced by nzb_u_inner
    13 !
    14 ! Revision 1.14  2005/03/26 14:49:31  raasch
    15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions
    16 !
    17 ! Revision 1.13  2004/01/30 10:11:21  raasch
    18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    19 !
    20 ! Revision 1.12  2003/03/12 16:17:06  raasch
    21 ! Full code replaced in the call for all gridpoints instead of calling the
    22 ! _ij version (required by NEC, because otherwise no vectorization)
    23 !
    24 ! Revision 1.11  2002/06/11 12:25:00  raasch
    25 ! Former subroutine changed to a module which allows to be called for all grid
    26 ! points of a single vertical column with index i,j or for all grid points by
    27 ! using function overloading.
    28 !
    29 ! Revision 1.10  2001/03/29 17:34:21  raasch
    30 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    31 !
    32 ! Revision 1.9  2001/01/22 04:55:36  raasch
    33 ! Module test_variables removed
    34 !
    35 ! Revision 1.8  2000/01/20 09:34:26  letzel
    36 ! All comments translated into English
    37 !
    38 ! Revision 1.7  1998/08/05 06:47:04  raasch
    39 ! Alle Gitterpunkte gleich behandelt (keine zentralen Differenzenquotienten
    40 ! mehr bei Vertikaladvektion bei k=1)
    4115!
    4216! Revision 1.1  1997/08/11 06:09:21  raasch
  • palm/trunk/SOURCE/advec_u_up.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_u_up.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.12  2006/02/23 09:45:04  raasch
    1214! nzb_2d replaced by nzb_u_inner
    13 !
    14 ! Revision 1.11  2005/03/26 14:50:43  raasch
    15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions
    16 !
    17 ! Revision 1.10  2004/01/30 10:11:47  raasch
    18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    19 !
    20 ! Revision 1.9  2003/03/12 16:17:42  raasch
    21 ! Full code replaced in the call for all gridpoints instead of calling the
    22 ! _ij version (required by NEC, because otherwise no vectorization)
    23 !
    24 ! Revision 1.8  2002/06/11 12:25:43  raasch
    25 ! Former subroutine changed to a module which allows to be called for all grid
    26 ! points of a single vertical column with index i,j or for all grid points by
    27 ! using function overloading.
    28 !
    29 ! Revision 1.7  2001/03/29 17:35:08  raasch
    30 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    31 !
    32 ! Revision 1.6  2001/01/22 04:57:23  raasch
    33 ! Module test_variables removed
    34 !
    35 ! Revision 1.5  2000/01/20 09:40:00  letzel
    36 ! All comments translated into English
    37 !
    38 ! Revision 1.4  1998/07/06 12:04:43  raasch
    39 ! + USE test_variables
    40 !
    41 ! Revision 1.3  1998/04/21 15:52:51  raasch
    42 ! Galilei-Transformation eingebaut
    43 !
    44 ! Revision 1.2  1997/09/09 08:27:49  raasch
    45 ! Kehrwerte der Gitterweiten implementiert
    4615!
    4716! Revision 1.1  1997/08/29 08:55:25  raasch
  • palm/trunk/SOURCE/advec_u_ups.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_u_ups.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.6  2004/04/30 08:03:19  raasch
    1214! Enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.5  2003/03/16 09:25:07  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.4  2001/03/29 17:35:32  raasch
    18 ! v_ad is now allocated locally, several temporary arrays removed from
    19 ! argument lists of spline_* and long_filter_*
    20 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    21 !
    22 ! Revision 1.3  2001/01/22 04:59:29  raasch
    23 ! Module test_variables removed
    24 !
    25 ! Revision 1.2  2000/01/20 09:46:52  letzel
    26 ! All comments translated into English
    2715!
    2816! Revision 1.1  1999/02/05 08:49:08  raasch
  • palm/trunk/SOURCE/advec_v_pw.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_v_pw.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 09:46:08  raasch
    1214! nzb_2d replaced by nzb_v_inner
    13 !
    14 ! Revision 1.14  2005/03/26 14:52:51  raasch
    15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions
    16 !
    17 ! Revision 1.13  2004/01/30 10:12:25  raasch
    18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    19 !
    20 ! Revision 1.12  2003/03/12 16:18:19  raasch
    21 ! Full code replaced in the call for all gridpoints instead of calling the
    22 ! _ij version (required by NEC, because otherwise no vectorization)
    23 !
    24 ! Revision 1.11  2002/06/11 12:26:18  raasch
    25 ! Former subroutine changed to a module which allows to be called for all grid
    26 ! points of a single vertical column with index i,j or for all grid points by
    27 ! using function overloading.
    28 !
    29 ! Revision 1.10  2001/03/29 17:36:06  raasch
    30 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    31 !
    32 ! Revision 1.9  2001/01/22 05:14:48  raasch
    33 ! Module test_variables removed
    34 !
    35 ! Revision 1.8  2000/01/20 10:16:37  letzel
    36 ! All comments translated into English
    37 !
    38 ! Revision 1.7  1998/08/05 06:49:43  raasch
    39 ! Alle Gitterpunkte gleich behandelt (keine zentralen Differenzenquotienten mehr
    40 ! bei Vertikaladvektion bei k=1)
    4115!
    4216! Revision 1.1  1997/08/11 06:09:57  raasch
  • palm/trunk/SOURCE/advec_v_up.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_v_up.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.12  2006/02/23 09:46:37  raasch
    1214! nzb_2d replaced by nzb_v_inner
    13 !
    14 ! Revision 1.11  2005/03/26 14:55:59  raasch
    15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions
    16 !
    17 ! Revision 1.10  2004/01/30 10:12:45  raasch
    18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    19 !
    20 ! Revision 1.9  2003/03/12 16:19:00  raasch
    21 ! Full code replaced in the call for all gridpoints instead of calling the
    22 ! _ij version (required by NEC, because otherwise no vectorization)
    23 !
    24 ! Revision 1.8  2002/06/11 12:26:48  raasch
    25 ! Former subroutine changed to a module which allows to be called for all grid
    26 ! points of a single vertical column with index i,j or for all grid points by
    27 ! using function overloading.
    28 !
    29 ! Revision 1.7  2001/03/29 17:36:22  raasch
    30 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    31 !
    32 ! Revision 1.6  2001/01/22 05:16:19  raasch
    33 ! Module test_variables removed
    34 !
    35 ! Revision 1.5  2000/01/20 10:20:31  letzel
    36 ! All comments translated into English
    37 !
    38 ! Revision 1.4  1998/07/06 12:05:37  raasch
    39 ! + USE test_variables
    40 !
    41 ! Revision 1.3  1998/04/21 15:53:36  raasch
    42 ! Galilei-Transformation eingebaut
    43 !
    44 ! Revision 1.2  1997/09/09 08:28:21  raasch
    45 ! Kehrwerte der Gitterweiten implementiert
    4615!
    4716! Revision 1.1  1997/08/29 08:56:05  raasch
  • palm/trunk/SOURCE/advec_v_ups.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_v_ups.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.7  2004/04/30 08:03:52  raasch
    1214! Enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.6  2003/03/16 09:25:59  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.5  2001/03/29 17:36:43  raasch
    18 ! v_ad is now allocated locally, several temporary arrays removed from
    19 ! argument lists of spline_* and long_filter_*
    20 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    21 !
    22 ! Revision 1.4  2001/01/22 08:31:10  raasch
    23 ! Module test_variables removed
    24 !
    25 !
    26 ! Revision 1.2  2000/01/20 09:52:54  letzel
    27 ! All comments translated into English
    2815!
    2916! Revision 1.1  1999/02/05 08:50:32  raasch
  • palm/trunk/SOURCE/advec_w_pw.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_w_pw.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 09:47:01  raasch
    1214! nzb_2d replaced by nzb_w_inner
    13 !
    14 ! Revision 1.14  2004/01/30 10:13:07  raasch
    15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    16 !
    17 ! Revision 1.13  2003/03/12 16:19:29  raasch
    18 ! Full code replaced in the call for all gridpoints instead of calling the
    19 ! _ij version (required by NEC, because otherwise no vectorization)
    20 !
    21 ! Revision 1.12  2002/06/11 12:27:33  raasch
    22 ! Former subroutine changed to a module which allows to be called for all grid
    23 ! points of a single vertical column with index i,j or for all grid points by
    24 ! using function overloading.
    25 !
    26 ! Revision 1.11  2001/03/29 17:37:10  raasch
    27 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    28 !
    29 ! Revision 1.10  2001/01/22 05:19:42  raasch
    30 ! Module test_variables removed
    31 !
    32 ! Revision 1.9  2000/01/20 10:18:31  letzel
    33 ! All comments translated into English
    34 !
    35 ! Revision 1.8  1998/08/05 06:51:10  raasch
    36 ! Alle Gitterpunkte werden gleich behandelt (keine zentralen Differenzenquotien-
    37 ! ten mehr bei Vertikaladvektion bei k=1)
    3815!
    3916! Revision 1.1  1997/08/11 06:10:29  raasch
  • palm/trunk/SOURCE/advec_w_up.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_w_up.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2006/02/23 09:47:23  raasch
    1214! *** empty log message ***
    13 !
    14 ! Revision 1.10  2004/01/30 10:13:31  raasch
    15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    16 !
    17 ! Revision 1.9  2003/03/12 16:20:07  raasch
    18 ! Full code replaced in the call for all gridpoints instead of calling the
    19 ! _ij version (required by NEC, because otherwise no vectorization)
    20 !
    21 ! Revision 1.8  2002/06/11 12:28:33  raasch
    22 ! Former subroutine changed to a module which allows to be called for all grid
    23 ! points of a single vertical column with index i,j or for all grid points by
    24 ! using function overloading.
    25 !
    26 ! Revision 1.7  2001/03/29 17:37:31  raasch
    27 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    28 !
    29 ! Revision 1.6  2001/01/22 05:21:01  raasch
    30 ! Module test_variables removed
    31 !
    32 ! Revision 1.5  2000/01/20 10:22:08  letzel
    33 ! All comments translated into English
    34 !
    35 ! Revision 1.4  1998/07/06 12:06:25  raasch
    36 ! + USE test_variables
    37 !
    38 ! Revision 1.3  1998/04/21 15:54:20  raasch
    39 ! Galilei-Transformation eingebaut
    40 !
    41 ! Revision 1.2  1997/09/09 08:28:55  raasch
    42 ! Kehrwerte der Gitterweiten implementiert
    4315!
    4416! Revision 1.1  1997/08/29 08:56:33  raasch
  • palm/trunk/SOURCE/advec_w_ups.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: advec_w_ups.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.6  2004/04/30 08:05:05  raasch
    1214! Enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.5  2003/03/16 09:26:46  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.4  2001/03/29 17:37:53  raasch
    18 ! v_ad is now allocated locally, several temporary arrays removed from
    19 ! argument lists of spline_* and long_filter_*
    20 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    21 !
    22 ! Revision 1.3  2001/01/22 05:23:31  raasch
    23 ! Module test_variables removed
    24 !
    25 ! Revision 1.2  2000/01/20 09:57:15  letzel
    26 ! All comments translated into English
    2715!
    2816! Revision 1.1  1999/02/05 08:52:09  raasch
  • palm/trunk/SOURCE/asselin_filter.f90

    r1 r3  
    11 SUBROUTINE asselin_filter
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! ---------------------
    10 ! $Log: asselin_filter.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.8  2004/01/30 10:14:02  raasch
    1214! Scalar lower k index nzb replaced by 2d-array nzb_2d
    13 !
    14 ! Revision 1.7  2003/03/16 09:26:58  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.6  2003/03/14 13:37:38  raasch
    18 ! Do loop constructs only for IBM, otherwise array notation re-introduced
    19 !
    20 ! Revision 1.5  2002/12/19 13:45:57  raasch
    21 ! Array notation changed to do-loop constructs due to better performance,
    22 ! OpenMP directives added
    23 !
    24 ! Revision 1.4  2002/05/15 10:00:18  raasch
    25 ! Nothing changed
    26 !
    27 ! Revision 1.3  2002/05/02  19:39:07  19:39:07  raasch (Siegfried Raasch)
    28 ! Use of module interfaces included
    29 !
    30 ! Revision 1.2  2002/05/02 14:06:02  raasch
    31 ! *** empty log message ***
    3215!
    3316! Revision 1.1  2002/05/02  13:43:53  13:43:53  raasch (Siegfried Raasch)
     
    3821! -------------
    3922! Time filter needed for the leap-frog method
    40 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    4124
    4225    USE arrays_3d
  • palm/trunk/SOURCE/average_3d_data.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: average_3d_data.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.1  2006/02/23 09:48:58  raasch
    1214! Initial revision
  • palm/trunk/SOURCE/boundary_conds.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: boundary_conds.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 09:54:55  raasch
    1214! Surface boundary conditions in case of topography: nzb replaced by
     
    1416! unchanged (still using nzb) because a non-flat topography must use a
    1517! Prandtl-layer, which don't requires explicit setting of the surface values.
    16 !
    17 ! Revision 1.14  2005/06/29 09:49:40  steinfeld
    18 ! Baroclinicity considered in the Dirichlet boundary condition for u and v at
    19 ! the top boundary
    20 !
    21 ! Revision 1.13  2005/04/23 08:35:30  raasch
    22 ! Error removed in Dirichlet bottom boundary conditions for pt and q in case
    23 ! of Runge-Kutta schemes. So far, new surface values have been taken from
    24 ! timelevel t-dt, which does not exist in case of Runge-Kutta schemes.
    25 !
    26 ! Revision 1.12  2005/03/26 14:58:24  raasch
    27 ! Non-cyclic boundary conditions included, argument range added
    28 !
    29 ! Revision 1.11  2001/03/30 06:54:18  raasch
    30 ! Timelevel t+dt replaced by timelevel t,
    31 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    32 !
    33 ! Revision 1.10  2001/01/29 12:19:27  raasch
    34 ! Passive scalar is considered
    35 !
    36 ! Revision 1.9  2001/01/22 05:25:56  raasch
    37 ! Module test_variables removed
    38 !
    39 ! Revision 1.8  2000/07/04 14:07:59  raasch
    40 ! Missing diriclet boundary conditions for temperature and total water
    41 ! content added
    42 !
    43 ! Revision 1.7  2000/04/13 13:56:05  schroeter
    44 ! Boundaray conditions for total water content
    45 !
    46 ! Revision 1.6  2000/01/20  10:44:48  10:44:48  letzel (Marcus Letzel)
    47 ! All comments translated into English
    48 !
    49 ! Revision 1.5  2000/01/10 09:28:53  raasch
    50 ! Variablenuebergabe jetzt per Modul
    51 ! Randbedingungen fuer w im Rahmen der pointer-Einfuehrung
    52 !
    53 ! Revision 1.4  1998/07/06 12:07:05  raasch
    54 ! + USE test_variables
    55 !
    56 ! Revision 1.3  1998/03/10 07:19:37  raasch
    57 ! Beschreibung ergaenzt um Randbedingung fuer TKE
    58 !
    59 ! Revision 1.2  1997/09/19 07:39:03  raasch
    60 ! Randbedingungen fuer TKE
    6118!
    6219! Revision 1.1  1997/09/12 06:21:34  raasch
  • palm/trunk/SOURCE/buoyancy.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: buoyancy.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.19  2006/04/26 12:09:56  raasch
    1214! OpenMP optimization (one dimension added to sums_l)
    13 !
    14 ! Revision 1.18  2006/02/23 09:55:33  raasch
    15 ! nzb_2d replaced by nzb_s_inner/outer, nanz_2dh(0) replaced by
    16 ! ngp_2dh_outer(:,0)
    17 !
    18 ! Revision 1.17  2005/12/06 15:29:16  raasch
    19 ! Mean pt profile is calculated only in case of the first intermediate timestep
    20 !
    21 ! Revision 1.16  2004/01/30 10:14:31  raasch
    22 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    23 !
    24 ! Revision 1.15  2003/10/29 08:41:47  raasch
    25 ! Horizontal mean temperature is now taken from array hom instead of array sums
    26 !
    27 ! Revision 1.14  2003/03/16 09:27:22  raasch
    28 ! Two underscores (_) are placed in front of all define-strings
    29 !
    30 ! Revision 1.13  2003/03/12 16:20:47  raasch
    31 ! Full code replaced in the call for all gridpoints instead of calling the
    32 ! _ij version (required by NEC, because otherwise no vectorization)
    33 !
    34 ! Revision 1.12  2002/12/19 13:47:38  raasch
    35 ! STOP statement replaced by call of subroutine local_stop
    36 !
    37 ! Revision 1.11  2002/06/11 12:33:07  raasch
    38 ! Former subroutine changed to a module which allows to be called for all grid
    39 ! points of a single vertical column with index i,j or for all grid points by
    40 ! using function overloading.
    41 ! Calculation of the horizontally averaged temperature profile moved to new
    42 ! subroutine calc_mean_pt_profile, which is part of this module.
    43 !
    44 ! Revision 1.10  2001/03/30 06:54:56  raasch
    45 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    46 !
    47 ! Revision 1.9  2001/01/22 05:27:03  raasch
    48 ! Module test_variables removed
    49 !
    50 ! Revision 1.8  2000/07/03 12:55:56  raasch
    51 ! Argument theta declared as pointer
    52 !
    53 ! Revision 1.7  2000/04/27 07:04:19  raasch
    54 ! in case of a sloping surface the temperature difference is computed relative
    55 ! to the initial (2D) temperature field
    56 !
    57 ! Revision 1.6  2000/04/13 14:27:34  schroeter
    58 ! considering the influence of humidity to the buoyancy term
    59 ! (only for sloping_surface = false),
    60 !
    61 ! Revision 1.5  2000/01/21  16:17:23  16:17:23  letzel (Marcus Letzel)
    62 ! All comments translated into English
    63 !
    64 ! Revision 1.4  1998/09/22 17:16:22  raasch
    65 ! Auftriebsterm erweitert fuer Rechnungen mit geneigter Ebene
    66 !
    67 ! Revision 1.3  1998/07/06 12:07:27  raasch
    68 ! + USE test_variables
    69 !
    70 ! Revision 1.2  1998/03/30 11:33:29  raasch
    71 ! nanz_2dh in nanz_2dh(0) geaendert
    7215!
    7316! Revision 1.1  1997/08/29 08:56:48  raasch
  • palm/trunk/SOURCE/calc_liquid_water_content.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: calc_liquid_water_content.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2005/03/26 15:22:06  raasch
    1214! Arguments for non-cyclic boundary conditions added to argument list of
    1315! routine exchange_horiz,
    1416! ql calculated for the ghost points, exchange of ghost points removed
    15 !
    16 ! Revision 1.4  2004/01/30 10:15:13  raasch
    17 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    18 !
    19 ! Revision 1.3  2001/01/22 05:30:59  raasch
    20 ! Module test_variables removed
    21 !
    22 ! Revision 1.2  2000/06/13 07:49:32  schroeter
    23 ! Horizontal boundary conditions (cyclic) of ql will be set now
    2417!
    2518! Revision 1.1  2000/04/13 14:50:45  schroeter
  • palm/trunk/SOURCE/calc_precipitation.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: calc_precipitation.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2004/01/30 10:15:57  raasch
    1214! Scalar lower k index nzb replaced by 2d-array nzb_2d
    13 !
    14 ! Revision 1.4  2003/03/12 16:21:23  raasch
    15 ! Full code replaced in the call for all gridpoints instead of calling the
    16 ! _ij version (required by NEC, because otherwise no vectorization)
    17 !
    18 ! Revision 1.3  2002/06/11 12:38:58  raasch
    19 ! Former subroutine changed to a module which allows to be called for all grid
    20 ! points of a single vertical column with index i,j or for all grid points by
    21 ! using function overloading.
    22 !
    23 ! Revision 1.2  2001/01/22 05:32:25  raasch
    24 ! Module test_variables removed
    2515!
    2616! Revision 1.1  2000/04/13 14:45:22  schroeter
  • palm/trunk/SOURCE/calc_radiation.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: calc_radiation.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.6  2004/01/30 10:17:03  raasch
    1214! Scalar lower k index nzb replaced by 2d-array nzb_2d
    13 !
    14 ! Revision 1.5  2003/03/12 16:21:56  raasch
    15 ! Full code replaced in the call for all gridpoints instead of calling the
    16 ! _ij version (required by NEC, because otherwise no vectorization)
    17 !
    18 ! Revision 1.4  2002/06/11 12:40:33  raasch
    19 ! Former subroutine changed to a module which allows to be called for all grid
    20 ! points of a single vertical column with index i,j or for all grid points by
    21 ! using function overloading.
    22 ! 1D-arrays lwp_groud,... are only allocated once in the first call.
    23 !
    24 ! Revision 1.3  2001/03/30 06:55:27  raasch
    25 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    26 !
    27 ! Revision 1.2  2001/01/22 05:35:05  raasch
    28 ! Module test_variables removed
    2915!
    3016! Revision 1.1  2000/04/13 14:42:45  schroeter
  • palm/trunk/SOURCE/calc_spectra.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: calc_spectra.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.9  2006/04/11 14:56:00  raasch
    1214! pl_spectra renamed data_output_sp
    13 !
    14 ! Revision 1.8  2005/05/18 15:08:18  raasch
    15 ! n is assigned the number of spectra (formerly, the number of spectra was
    16 ! given by n-1)
    17 !
    18 ! Revision 1.7  2005/04/23 08:41:39  raasch
    19 ! Abort in case of non-cyclic lateral boundaries
    20 !
    21 ! Revision 1.6  2004/04/30 08:23:12  raasch
    22 ! Usage of all available FFT-methods is possible, direct transpositions from
    23 ! z to x and from z to y are used, FFT routines renamed, module
    24 ! fft_for_1d_decomp renamed fft_xy, extensions for enlarged subdomain added
    25 !
    26 ! Revision 1.5  2003/03/16 09:27:35  raasch
    27 ! Two underscores (_) are placed in front of all define-strings
    28 !
    29 ! Revision 1.4  2002/06/11 12:44:15  raasch
    30 ! TRIM function added (needed on IBM-systems)
    31 !
    32 ! Revision 1.3  2001/03/30  06:55:58  06:55:58  raasch (Siegfried Raasch)
    33 ! Work array replaced by tend,
    34 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    35 !
    36 ! Revision 1.2  2001/01/25 06:52:21  raasch
    37 ! Module singleton moved to modules.f90
    3815!
    3916! Revision 1.1  2001/01/05 15:08:07  raasch
  • palm/trunk/SOURCE/check_for_restart.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: check_for_restart.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2007/02/11 12:55:13  raasch
    1214! Informative output to the job protocol
    13 !
    14 ! Revision 1.10  2003/03/16 09:27:49  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.9  2002/12/19 13:48:59  raasch
    18 ! Stop flag is also set in case of user forced restarts, routine name
    19 ! changed from check_cpu_time to check_for_restart, the run description
    20 ! header will be written on file CONTINUE_RUN
    21 !
    22 ! Revision 1.8  2001/08/21 08:22:23  raasch
    23 ! MPI_ALLREDUCE moved out of IF-clause
    24 !
    25 ! Revision 1.7  2001/03/30 06:56:28  raasch
    26 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    27 !
    28 ! Revision 1.6  2001/01/22 05:36:33  raasch
    29 ! Module test_variables removed
    30 !
    31 ! Revision 1.5  2001/01/02 17:19:27  raasch
    32 ! Initial value of terminate_run_l is set to .FALSE. (undefined before).
    33 ! Flag file will only be created by PE0. Unit 90 is used instead of 91.
    34 !
    35 ! Revision 1.4  2000/01/21 16:27:59  letzel
    36 ! All comments translated into English
    37 !
    38 ! Revision 1.3  1998/07/06 12:07:47  raasch
    39 ! + USE test_variables
    40 !
    41 ! Revision 1.2  1998/03/24 15:24:54  raasch
    42 ! Lauf wird nur dann beendet und Fortsetzungslauf gestartet, falls auch eine
    43 ! Abspeicherung der entsprechenden Binaerdateien vorgesehen ist.
    44 ! Beendigung der Rechnung auf mehreren Prozessoren richtig implementiert.
    4515!
    4616! Revision 1.1  1998/03/18 20:06:51  raasch
  • palm/trunk/SOURCE/check_open.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: check_open.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.44  2006/08/22 13:48:34  raasch
    1214! xz and yz cross sections now up to nzt+1
    13 !
    14 ! Revision 1.43  2006/08/04 14:05:52  raasch
    15 ! Opening of unit 109 (particle time series)
    16 !
    17 ! Revision 1.42  2006/03/01 16:39:47  raasch
    18 ! In parallel mode, all PEs are opening the same copy of file PARIN
    19 !
    20 ! Revision 1.41  2006/02/23 10:05:05  raasch
    21 ! New unit 20 for data logging,
    22 ! new units (111:113,116) for 2d/3d data averaging
    23 ! NetCDF files are allowed a 64 bit offset format,
    24 ! error number argument for handle_netcdf_error
    25 !
    26 ! Revision 1.40  2005/06/30 10:39:04  raasch
    27 ! No opening of units 101-108 in case of missing netcdf support
    28 !
    29 ! Revision 1.39  2005/06/26 19:51:36  raasch
    30 ! Data format version number for unit 85 changed
    31 !
    32 ! Revision 1.38  2005/05/18 15:09:01  raasch
    33 ! Extensions for NetCDF output
    34 !
    35 ! Revision 1.37  2005/04/23 08:50:18  raasch
    36 ! batch_scp calls is given the remote username as an additional argument
    37 !
    38 ! Revision 1.36  2004/01/28 15:03:21  raasch
    39 ! In case of opening unit 80, a barrier is set only for the first call of
    40 ! check_open (from routine init_particles), in order to avoid the possibility
    41 ! of hanging jobs, which may occur if unit 80 is opened within routine
    42 ! allocate_prt_memory.
    43 !
    44 ! Revision 1.35  2003/10/29 08:47:03  raasch
    45 ! Call of batch_scp on IBM corrected, output of data format version and
    46 ! particle group data on unit 85
    47 !
    48 ! Revision 1.34  2003/05/09 14:06:58  raasch
    49 ! Batch_scp paths for IBM, NEC and linux added
    50 !
    51 ! Revision 1.33  2003/03/16 09:28:00  raasch
    52 ! Two underscores (_) are placed in front of all define-strings
    53 !
    54 ! Revision 1.32  2002/12/19 14:03:15  raasch
    55 ! Unit 18 (CPU_MEASURES) will be opened on PE0 only, unit 14 is opened
    56 ! with special process-id-string myid_char_14 (see routine parin),
    57 ! STOP statement replaced by call of subroutine local_stop
    58 !
    59 ! Revision 1.31  2002/04/24 19:01:15  raasch
    60 ! Unit 85 is allowed to be re-opened. In case that a re-open of a file is not
    61 ! allowed, the subroutine will be exit.
    62 !
    63 ! Revision 1.30  2002/04/16 08:00:27  raasch
    64 ! ftpcopy replaced by batch_scp, retry of opening CPU_MEASURES/ in case of
    65 ! errors forced by adding ERR-parameter in OPEN statement.
    66 ! New unit 85 for particle data
    67 !
    68 ! Revision 1.29  2001/07/12 12:07:21  raasch
    69 ! Opening of unit 33 moved out of ELSE-block
    70 !
    71 ! Revision 1.28  2001/03/30 06:56:49  raasch
    72 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    73 !
    74 ! Revision 1.27  2001/01/29 12:21:52  raasch
    75 ! In case of AVS-output: if no path for ftpcopy is given, filenames are set
    76 ! to "unknown"
    77 !
    78 ! Revision 1.26  2001/01/25 06:54:28  raasch
    79 ! Error concerning initialization of unit 23 removed (initial informations
    80 ! were written also when reopened)
    81 !
    82 ! Revision 1.25  2001/01/05 15:09:29  raasch
    83 ! New units 81-84 for spectra
    84 !
    85 ! Revision 1.24  2001/01/02 17:21:11  raasch
    86 ! Re-opening of files is allowed, variable openfile now is of type file_status,
    87 ! new unit 80 for particle informations
    88 ! old revison remarks deleted
    89 !
    90 ! Revision 1.23  2000/04/27 06:51:12  raasch
    91 ! remote_addres changed to return_addres, old revision remarks deleted
    92 !
    93 ! Revision 1.22  2000/01/21 17:07:18  letzel
    94 ! All comments translated into English
    9515!
    9616! Revision 1.1  1997/08/11 06:10:55  raasch
  • palm/trunk/SOURCE/check_parameters.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: check_parameters.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.61  2006/08/04 14:20:25  raasch
    1214! do2d_unit and do3d_unit now defined as 2d-arrays, check of
     
    1416! generation of file header moved from routines palm and header to here
    1517!
    16 ! Revision 1.60  2006/06/02 15:09:15  raasch
    17 ! data_output_xy/xz/yz are set
    18 !
    19 ! Revision 1.59  2006/03/03 15:55:34  letzel
    20 ! Spelling inconsistency removed
    21 !
    22 ! Revision 1.58  2006/03/03 08:50:08  raasch
    23 ! Default setting of constant_heatflux now only depends on the value of
    24 ! surface_heatflux again
    25 !
    26 ! Revision 1.57  2006/02/27 10:16:44  letzel
    27 ! Error removed in initialization of dt_averaging_input_pr
    28 !
    29 ! Revision 1.56  2006/02/23 10:06:55  raasch
    30 ! Volume flow conservation is not allowed with non-cyclic lateral boundaries,
    31 ! calculation of nzb_diff moved to init_grid, the default setting of
    32 ! constant_heatflux = .F. additionally requires a surface initial temperature
    33 ! change /= 0, check of topography, dissipation_1d, mixing_length_1d,
    34 ! check and setting of new data output steering parameters and user defined
    35 ! output quantities,
    36 ! ebene renamed section, pl1d_anz, plts_anz renamed dopr_n, dots_n,
    37 ! plts, pl1d renamed data_output_ts/1d
    38 !
    39 ! Revision 1.55  2005/12/06 16:39:34  raasch
    40 ! Output of ql profile is allowed in case of using cloud droplets
    41 !
    42 ! Revision 1.54  2005/06/30 11:27:33  steinfeld
    43 ! Initialization of the profiles of the geostrophic wind components
    44 ! ug and vg;
    45 ! the usage of the geostrophic wind in a galilei transformation is
    46 ! not allowed with a simultaneous specification of a baroclinicity
    47 !
    48 ! Revision 1.53  2005/05/18 15:12:36  raasch
    49 ! Check of data_output_format and netcdf precision,
    50 ! abort if poisfft_hybrid is called in a non-parallel environment
    51 !
    52 ! Revision 1.52  2005/04/23 09:04:10  raasch
    53 ! fcl_factor renamed cfl_factor, default setting of outflow_damping_width
    54 ! corrected, additional check for non-cyclic lateral boundary conditions
    55 !
    56 ! Revision 1.51  2005/03/26 19:47:32  raasch
    57 ! Check of bc_lr, bc_ns, km_damp_max, inflow_disturbance_begin,
    58 ! inflow_disturbance_end, outflow_damping_width
    59 !
    60 ! Revision 1.50  2004/04/30 11:20:27  raasch
    61 ! Check of grid_matching, check of time series profiles modified,
    62 ! impulse_advec renamed momentum_advec
    63 !
    64 ! Revision 1.49  2004/01/28 15:04:56  raasch
    65 ! Additional checks and settings for timestep schemes because of the
    66 ! newly implemented Runge-Kutta schemes
    67 !
    68 ! Revision 1.48  2003/04/17 07:07:43  raasch
    69 ! Hybrid solver check in parallel environment only
    70 !
    71 ! Revision 1.47  2003/04/16 09:43:20  raasch
    72 ! Additional checks for hybrid solver and host, Monin Obukhov length added to
    73 ! time series
    74 !
    75 ! Revision 1.46  2003/03/31 15:43:19  raasch
    76 ! Check for horizontal temperature flux profiles cleaned up
    77 !
    78 ! Revision 1.45  2003/03/14 13:38:02  raasch
    79 ! Random generator check included
    80 !
    81 ! Revision 1.44  2003/03/12 16:23:10  raasch
    82 ! Small checks for NEC and temperton-algorithm
    83 !
    84 ! Revision 1.43  2002/12/19 14:04:26  raasch
    85 ! Setting default values for termination_time_needed and check the
    86 ! value for IBM-regatta, STOP statement replaced by call of subroutine
    87 ! local_stop
    88 !
    89 ! Revision 1.42  2002/09/12 13:02:25  raasch
    90 ! Error in calculating initial pt and q profiles removed (array bound 10 in
    91 ! *_ind arrays could be exceeded)
    92 !
    93 ! Revision 1.41  2002/06/11 12:49:33  raasch
    94 ! Hybrid solver allowed for solving the Poisson equation
    95 !
    96 ! Revision 1.40  2002/05/02  18:48:56  18:48:56  raasch (Siegfried Raasch)
    97 ! SELECT case of cross_normalized_x and ..y requires usage of TRIM-function
    98 ! on IBM
    99 !
    100 ! Revision 1.39  2002/04/16 08:03:27  raasch
    101 ! Parameters for scalar transport included, output of scalar profiles changed
    102 !
    103 ! Revision 1.38  2001/11/28 09:48:32  schroeter
    104 ! Profiles for horizontal heat fluxes and
    105 ! flux divergences allowed (58-72)
    106 !
    107 ! Revision 1.37  2001/09/04  11:58:55  11:58:55  raasch (Siegfried Raasch)
    108 ! Profiles for vertical flux of resolved/subgrid scale energy and pressure
    109 ! fluctuations allowed (#55-57)
    110 !
    111 ! Revision 1.36  2001/07/20 12:52:32  raasch
    112 ! Multigrid method allowed for solving the Poisson equation
    113 !
    114 ! Revision 1.35  2001/03/30 06:57:38  raasch
    115 ! Asselin filter factor is set to 0, if Euler scheme is used,
    116 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    117 !
    118 ! Revision 1.34  2001/02/09 07:04:45  raasch
    119 ! Checking of boundary conditions for pt and q and their combination with
    120 ! prescribed surface fluxes revised
    121 !
    122 ! Revision 1.33  2001/01/29 12:23:31  raasch
    123 ! New parameter passive_scalar is checked
    124 !
    125 ! Revision 1.32  2001/01/25 06:55:37  raasch
    126 ! New control-parameters fft_method and use_surface_fluxes are checked
    127 !
    128 ! Revision 1.31  2000/04/18 11:19:43  schroeter
    129 ! Bott-Chlond Advektionsschema nun auch fuer Rechnungen
    130 ! mit Feuchte/Wolkenphysik zugelassen
    131 !
    132 ! Revision 1.30  2000/04/13 13:28:56  schroeter
    133 ! + initialize profiles of total water content
    134 ! + check boundary conditions of total water content
    135 ! + initialize output of cloud-physics-profiles
    136 !
    137 ! Revision 1.29  2000/01/25  14:38:37  14:38:37  letzel (Marcus Letzel)
    138 ! All comments translated into English
    139 !
    14018! Revision 1.1  1997/08/26 06:29:23  raasch
    14119! Initial revision
  • palm/trunk/SOURCE/close_file.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: close_file.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.10  2006/08/22 13:50:01  raasch
    1214! xz and yz cross sections now up to nzt+1
    13 !
    14 ! Revision 1.9  2006/08/04 14:22:53  raasch
    15 ! Closing unit 109 (particle timeseries)
    16 !
    17 ! Revision 1.8  2006/02/23 10:08:23  raasch
    18 ! Extensions for 2d/3d data averaging
    19 ! pl2d_.._anz, pl3d_anz renamed do2d_.._n, do3d_avs_n, respectively
    20 ! error number argument for handle_netcdf_error
    21 !
    22 ! Revision 1.7  2005/05/18 15:13:12  raasch
    23 ! Extensions for NetCDF output
    24 !
    25 ! Revision 1.6  2005/03/26 19:54:23  raasch
    26 ! +datform in global-namelist (iso2d)
    27 !
    28 ! Revision 1.5  2003/03/04 11:24:46  raasch
    29 ! AVS data format changed from float to xdr_float (needed on linux machines
    30 ! due to the little/big endian problem)
    31 !
    32 ! Revision 1.4  2002/12/19 14:06:34  raasch
    33 ! ISO2D parameter dp set to true for ibm hosts
    34 !
    35 ! Revision 1.3  2001/03/30 06:58:36  raasch
    36 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    37 !
    38 ! Revision 1.2  2001/01/22 05:39:06  raasch
    39 ! Module test_variables removed
    4015!
    4116! Revision 1.1  2001/01/02 17:23:41  raasch
  • palm/trunk/SOURCE/compute_vpt.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: compute_vpt.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2001/03/30 06:58:52  raasch
    1214! Translation of remaining German identifiers (variables, subroutines, etc.)
    13 !
    14 ! Revision 1.4  2001/01/22 05:42:55  raasch
    15 ! Minor corrections within file header
    16 !
    17 ! Revision 1.2  2000/09/07 09:00:36  schroeter
    18 ! Now, vpt will be computed for k=nzb and nzt+1
    19 ! additionally to keep boundary values
    2015!
    2116! Revision 1.1  2000/04/13 14:40:53  schroeter
  • palm/trunk/SOURCE/coriolis.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: coriolis.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.12  2006/02/23 10:08:57  raasch
    1214! nzb_2d replaced by nzb_u/v/w_inner
    13 !
    14 ! Revision 1.11  2005/06/29 10:30:17  steinfeld
    15 ! A potential dependency of the geostrophic wind components ug and vg on height
    16 ! is considered
    17 !
    18 ! Revision 1.10  2005/03/26 20:02:04  raasch
    19 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions
    20 !
    21 ! Revision 1.9  2004/01/30 10:17:34  raasch
    22 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    23 !
    24 ! Revision 1.8  2003/03/12 16:23:47  raasch
    25 ! Full code replaced in the call for all gridpoints instead of calling the
    26 ! _ij version (required by NEC, because otherwise no vectorization)
    27 !
    28 ! Revision 1.7  2002/12/19 14:08:14  raasch
    29 ! STOP statement replaced by call of subroutine local_stop
    30 !
    31 ! Revision 1.6  2002/06/11 12:51:09  raasch
    32 ! Former subroutine changed to a module which allows to be called for all grid
    33 ! points of a single vertical column with index i,j or for all grid points by
    34 ! using function overloading.
    35 !
    36 ! Revision 1.5  2001/03/30 06:59:09  raasch
    37 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    38 !
    39 ! Revision 1.4  2001/01/22 05:44:16  raasch
    40 ! Module test_variables removed
    41 !
    42 ! Revision 1.3  2000/12/20 09:50:27  letzel
    43 ! All comments translated into English.
    44 !
    45 ! Revision 1.2  1998/07/06 12:09:54  raasch
    46 ! + USE test_variables
    4715!
    4816! Revision 1.1  1997/08/29 08:57:38  raasch
  • palm/trunk/SOURCE/cpu_log.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: cpu_log.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.24  2006/06/02 15:12:17  raasch
    1214! cpp-directives extended for lctit
    13 !
    14 ! Revision 1.23  2004/01/28 15:08:52  raasch
    15 ! Type log changed to logpoint due to name conflict with intrinsic log
    16 !
    17 ! Revision 1.22  2003/05/09 14:15:18  raasch
    18 ! Time measurement on linux machines included, measurements on IBM are now using
    19 ! function irtc, which allows correct measurements for jobs running over
    20 ! the 24:00 timeline
    21 !
    22 ! Revision 1.21  2003/03/16 09:30:04  raasch
    23 ! Two underscores (_) are placed in front of all define-strings
    24 !
    25 ! Revision 1.20  2003/03/12 16:24:18  raasch
    26 ! Time measurement on NEC implemented
    27 !
    28 ! Revision 1.19  2002/12/20 09:38:43  raasch
    29 ! variable first is defined for all hosts
    30 !
    31 ! Revision 1.18  2002/12/19 14:09:36  raasch
    32 ! Output of warnings in case of negative cpu-times, STOP statement replaced by
    33 ! call of subroutine local_stop
    34 !
    35 ! Revision 1.17  2002/05/02 18:49:37  raasch
    36 ! Time measurement on IBM implemented
    37 !
    38 ! Revision 1.16  2001/03/30 07:00:43  raasch
    39 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    40 !
    41 ! Revision 1.14  2001/01/25 06:56:49  raasch
    42 ! +cpp-directives for dec-alpha-workstations
    43 !
    44 ! Revision 1.13  2001/01/22 08:32:46  raasch
    45 ! Module test_variables removed
    46 !
    47 ! Revision 1.11  2000/12/20 10:18:58  letzel
    48 ! All comments translated into English.
    49 !
    50 ! Revision 1.9  1999/11/25 16:19:07  raasch
    51 ! TIMEF jetzt als Function-Aufruf
    52 !
    53 ! Revision 1.8  1999/03/03 09:26:14  raasch
    54 ! Zeitmessung auf T3E grundsaetzlich auf TIMEF umgestellt
    55 !
    56 ! Revision 1.7  1998/12/15 07:32:39  raasch
    57 ! Messung der wahren I/O-Zeiten auf T3E mit TIMEF moeglich
    58 !
    59 ! Revision 1.4  1997/09/16 06:37:24  raasch
    60 ! Zeitmessungen auf t3e umgestellt auf tsecnd (schneller als mpi_wtime)
    6115!
    6216! Revision 1.1  1997/07/24 11:12:29  raasch
     
    6721! ------------
    6822! Cpu-time measurements for any program part whatever.
    69 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    7024
    7125    USE cpulog
  • palm/trunk/SOURCE/cpu_statistics.f90

    r1 r3  
    11 SUBROUTINE cpu_statistics
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: cpu_statistics.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.13  2006/04/26 12:10:51  raasch
    1214! Output of number of threads per task, max = min in case of 1 PE
    13 !
    14 ! Revision 1.12  2003/03/16 09:30:22  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.11  2002/12/19 14:12:08  raasch
    18 ! PE results are collected on PE0 in order to calculate mean cpu statistics.
    19 ! Output format changed.
    20 !
    21 ! Revision 1.10  2001/03/30 07:01:02  raasch
    22 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    23 !
    24 ! Revision 1.9  2001/01/22 05:49:49  raasch
    25 ! Module test_variables removed
    26 !
    27 ! Revision 1.8  2001/01/02 17:25:23  raasch
    28 ! Unit 18 closed at the end of the subroutine
    29 !
    30 ! Revision 1.7  2000/12/20 10:10:26  letzel
    31 ! All comments translated into English.
    3215!
    3316! Revision 1.1  1997/07/24 11:11:11  raasch
     
    4124! statistics. The output is sorted according to the amount of cpu-time consumed
    4225! and output on PE0.
    43 !-------------------------------------------------------------------------------!
     26!------------------------------------------------------------------------------!
    4427
    4528    USE cpulog
  • palm/trunk/SOURCE/data_log.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_log.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.1  2006/02/23 10:09:29  raasch
    1214! Initial revision
  • palm/trunk/SOURCE/data_output_2d.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_2d.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2006/08/22 13:50:29  raasch
    1214! xz and yz cross sections now up to nzt+1
    13 !
    14 ! Revision 1.4  2006/06/06 07:55:18  raasch
    15 ! Error removed: p is assigned to to_be_resorted instead of w
    16 !
    17 ! Revision 1.3  2006/06/02 15:18:21  raasch
    18 ! +argument "found" in call of routine user_data_output_2d,
    19 ! immediate return, if no output is requested
    2015!
    2116! Revision 1.2  2006/02/23 10:19:22  raasch
     
    2621! Former subroutine plot_2d renamed data_output_2d, pl2d.. renamed do2d..,
    2722! anz renamed ngp, ebene renamed section, pl2d_.._anz renamed do2d_.._n
    28 !
    29 ! Revision 1.1  2006/02/23 10:12:31  raasch
    30 ! Initial revision
    31 !
    32 ! Revision 1.25  2005/10/20 15:41:41  raasch
    33 ! Error in output of netcdf yz-cross-sections removed
    34 !
    35 ! Revision 1.24  2005/06/26 20:07:25  raasch
    36 ! +ql_c, ql_v, ql_vp
    37 !
    38 ! Revision 1.23  2005/05/18 15:46:56  raasch
    39 ! Extensions for NetCDF output, output of vpt, qv and lwp
    40 !
    41 ! Revision 1.22  2005/03/26 20:47:45  raasch
    42 ! local_2d is allocated with upper bound nzt (former nzt+1 was wrong),
    43 ! additional barrier in case that PE0 gathers data from the other
    44 ! PEs, this barrier had to be set at the end of the gathering because
    45 ! otherwise PE0 may receive wrong data on tag 0.
    46 !
    47 ! Revision 1.21  2004/04/30 12:32:50  raasch
    48 ! Adjustments for enlarged grids (unequal subdomain sizes)
    49 !
    50 ! Revision 1.20  2003/03/16 09:42:27  raasch
    51 ! Two underscores (_) are placed in front of all define-strings
    52 !
    53 ! Revision 1.19  2002/12/19 15:56:24  raasch
    54 ! STOP statement replaced by call of subroutine local_stop
    55 !
    56 ! Revision 1.18  2002/06/11 13:13:00  raasch
    57 ! TRIM function added (needed on IBM-systems)
    58 !
    59 ! Revision 1.17  2002/04/16  08:11:16  08:11:16  raasch (Siegfried Raasch)
    60 ! Output of scalar field allowed with strings s_xy, s_xz, s_yz
    61 !
    62 ! Revision 1.16  2001/03/30 07:44:29  raasch
    63 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    64 !
    65 ! Revision 1.15  2001/03/09 21:36:29  letzel
    66 ! All comments translated into English.
    67 !
    68 ! Revision 1.14  2001/01/29 12:32:07  raasch
    69 ! Passive scalar is considered, xy-plot of q and ql corrected (level_z was
    70 ! missing)
    71 !
    72 ! Revision 1.13  2001/01/22 07:38:20  raasch
    73 ! Module test_variables removed
    74 !
    75 ! Revision 1.12  2001/01/02 17:33:39  raasch
    76 ! Plot output file is closed at the end of the subroutine
    77 !
    78 ! Revision 1.11  2000/04/13 14:01:44  schroeter
    79 ! Plotausgaben fuer Gesamtwasser- und Fluessigwassergehalt
    80 !
    81 ! Revision 1.10  99/02/05  09:14:06  09:14:06  raasch (Siegfried Raasch)
    82 ! u=v=0 am unteren Rand bei Vertikalschnitten wird jetzt ausschliesslich
    83 ! bei Spiegelungsrandbedingungen verwendet
    84 !
    85 ! Revision 1.9  1998/09/22 17:25:58  raasch
    86 ! u*- und theta*-Ausgabe korrigiert (muessen auf nzb+1 des Hilfsfeldes local_2d
    87 ! umgespeichert werden.
    88 !
    89 ! Revision 1.8  1998/07/06 12:26:06  raasch
    90 ! Exakte Ausgabe der Lage der Schnittebenen fuer die verschiedenen Felder
    91 ! bei Horizontalschnitten, + USE test_variables
    92 !
    93 ! Revision 1.7  1998/03/09 16:24:19  raasch
    94 ! TRIM( simulated_time_chr )
    95 !
    96 ! Revision 1.6  1998/02/19 07:10:30  raasch
    97 ! Ausgabe der Teilfelder jetzt auch wieder wahlweise durch die einzelnen PEs
    98 !
    99 ! Revision 1.5  1998/02/10 15:08:25  raasch
    100 ! Spiegelungsrandbedingung fuer u und v bei Vertikalschnitten durch
    101 ! u=v=0 ersetzt, Ausgabe der Schnittlevel,
    102 ! Empfang der Teilfelder durch PE0 erfolgt in beliebiger Reihenfolge durch
    103 ! Verwendung von MPI_ANY_TAG, Barrieren werden gesetzt, um MPI-Fehler zu umgehen
    104 !
    105 ! Revision 1.4  1998/01/23 10:05:42  raasch
    106 ! Ausgabe von 2D-Feldern (u*_xy, t*_xy)
    107 !
    108 ! Revision 1.3  1997/09/08 10:47:23  raasch
    109 ! Prozessor 0 gibt jetzt gesamte Schnittebene des Modells aus
    110 !
    111 ! Revision 1.2  1997/08/29 09:00:36  raasch
    112 ! Ausgabe der potentiellen Temperatur
    11323!
    11424! Revision 1.1  1997/08/11 06:24:09  raasch
  • palm/trunk/SOURCE/data_output_3d.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_3d.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.3  2006/06/02 15:18:59  raasch
    1214! +argument "found", -argument grid in call of routine user_data_output_3d
     
    1618! .._anz renamed .._n,
    1719! output extended to (almost) all quantities, output of user-defined quantities
    18 !
    19 ! Revision 1.1  2006/02/23 10:21:34  raasch
    20 ! Initial revision
    21 !
    22 ! Revision 1.13  2005/06/26 20:08:36  raasch
    23 ! +ql_c, ql_v, ql_vp
    24 !
    25 ! Revision 1.12  2005/05/18 15:47:34  raasch
    26 ! Extensions for NetCDF output
    27 !
    28 ! Revision 1.11  2003/03/16 09:42:31  raasch
    29 ! Two underscores (_) are placed in front of all define-strings
    30 !
    31 ! Revision 1.10  2002/06/11 13:13:27  raasch
    32 ! TRIM function added (needed on IBM-systems)
    33 !
    34 ! Revision 1.9  2002/04/16  08:11:50  08:11:50  raasch (Siegfried Raasch)
    35 ! Output of scalar field s allowed
    36 !
    37 ! Revision 1.8  2001/03/30 07:44:45  raasch
    38 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    39 !
    40 ! Revision 1.7  2001/03/09 21:46:23  letzel
    41 ! All comments translated into English.
    42 !
    43 ! Revision 1.6  2001/01/22 07:40:10  raasch
    44 ! Module test_variables removed
    45 !
    46 ! Revision 1.5  2000/04/13 14:02:39  schroeter
    47 ! Plotausgabe fuer Gesamtwasser- und Fluessigwassergehalt sowie potentielle
    48 ! Fluessigwassertemperatur
    49 !
    50 ! Revision 1.4  99/03/02  09:24:51  09:24:51  raasch (Siegfried Raasch)
    51 ! Datenkompression implementiert
    52 !
    53 ! Revision 1.3  1998/07/06 12:26:35  raasch
    54 ! + USE test_variables
    55 !
    56 ! Revision 1.2  1998/03/09 16:24:33  raasch
    57 ! simulated_time_mod verlaengert
    5820!
    5921! Revision 1.1  1997/09/03 06:29:36  raasch
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r1 r3  
    3232! Actual revisions:
    3333! -----------------
    34 !
    3534! TEST: write statements
    3635!
    3736! Former revisions:
    3837! -----------------
    39 ! $Log: data_output_dvrp.f90,v $
     38! $Id$
     39! RCS Log replace by Id keyword, revision history cleaned up
     40!
    4041! Revision 1.13  2006/02/23 10:25:12  raasch
    4142! Former routine plot_dvrp renamed data_output_dvrp,
    4243! Only a fraction of the particles may have a tail,
    4344! pl.. replaced by do.., %size renamed %dvrp_psize
    44 !
    45 ! Revision 1.12  2005/03/31 07:00:46  raasch
    46 ! Testversion with print statements checked in
    47 !
    48 ! Revision 1.11  2003/08/01 08:44:50  raasch
    49 ! Routine defined as recursive, because it is called from DVRP_STEERING_UPDATE
    50 ! in case that the user has suspended the run by steering control.
    51 ! Routine color_dvrp placed at the beginning of this file. This routine is used
    52 ! for slicer colortable settings.
    53 !
    54 ! Revision 1.10  2003/03/16 09:42:34  raasch
    55 ! Two underscores (_) are placed in front of all define-strings
    56 !
    57 ! Revision 1.9  2003/03/04 11:35:07  raasch
    58 ! Update for new version of dvrp software (e.g. using module dvrp is now
    59 ! mandatory)
    60 !
    61 ! Revision 1.8  2001/11/09 13:41:21  raasch
    62 ! Colour can change along the particle tails
    63 !
    64 ! Revision 1.7  2001/08/21 09:55:20  raasch
    65 ! Superelevation for particles implemented, plotting of particle tails added,
    66 ! colourtable definitions moved to subroutine user_dvrp_coltab
    67 !
    68 ! Revision 1.6  2001/07/12 12:44:12  raasch
    69 ! Test: color definitions for particles changed
    70 !
    71 ! Revision 1.5  2001/03/30 07:45:00  raasch
    72 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    73 !
    74 ! Revision 1.4  2001/01/29 12:32:51  raasch
    75 ! Passive scalar considered
    76 !
    77 ! Revision 1.3  2001/01/22 07:41:12  raasch
    78 ! Module test_variables removed
    79 !
    80 ! Revision 1.2  2000/12/28 13:37:19  raasch
    81 ! Total revision due to new features of dvrp-software,
    82 ! plot of slicers included,
    83 ! code is used only optionally (cpp-directives are added)
    8445!
    8546! Revision 1.1  2000/04/27 06:27:17  raasch
  • palm/trunk/SOURCE/data_output_profiles.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_profiles.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.18  2006/08/16 14:27:04  raasch
    1214! PRINT* statements for testing removed
    13 !
    14 ! Revision 1.17  2006/08/08 12:45:47  raasch
    15 ! Immediate return in case of average_count_pr = 0.
    16 !
    17 ! Revision 1.16  2006/02/23 10:26:53  raasch
    18 ! Former routine plot_1d renamed data_output_profiles
    19 ! pl.. renamed do.., .._anz renamed .._no
    20 ! error number argument for handle_netcdf_error
    21 !
    22 ! Revision 1.15  2005/05/18 15:44:44  raasch
    23 ! Extensions for NetCDF output
    24 !
    25 ! Revision 1.14  2001/03/30 07:44:12  raasch
    26 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    27 !
    28 ! Revision 1.13  2001/03/02 14:21:17  letzel
    29 ! All comments translated into English.
    30 !
    31 ! Revision 1.12  2001/01/22 07:36:55  raasch
    32 ! Module test_variables removed
    33 !
    34 ! Revision 1.11  1998/07/06 12:25:36  raasch
    35 ! + USE test_variables
    36 !
    37 ! Revision 1.10  1998/04/15 11:22:02  raasch
    38 ! zeitlich gemittelte Profile werden jetzt auch mit zeitlich gemittelten
    39 ! Groessen normiert
    40 !
    41 ! Revision 1.9  1998/04/06 14:33:53  raasch
    42 ! Implementierung weiterer Normierungsmethoden
    43 !
    44 ! Revision 1.8  1998/03/30 11:37:30  raasch
    45 ! Erweiterungen fuer Gebietsstatistiken korrigiert
    46 !
    47 ! Revision 1.7  1998/03/25 13:55:45  raasch
    48 ! Erweiterungen fuer Gebietsstatistiken
    49 !
    50 ! Revision 1.6  1998/03/09 16:24:04  raasch
    51 ! Ueberpruefung der Normierungsfaktoren auf 0.0 und evtl.
    52 ! Abschaltung der Normalisierung
    53 !
    54 ! Revision 1.5  1998/03/03 07:58:11  raasch
    55 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier
    56 ! Ausgabe von zeitlich gemittelten Profilen, Abspeichern von Normierungsgroessen
    57 !
    58 ! Revision 1.4  1998/02/19 07:10:04  raasch
    59 ! x-Wertebereichsermittelung nur bezueglich des auszugebenden Hoehenbereiches,
    60 ! x-Wertebereichsermittelung auch fuer Anfangsprofile
    61 !
    62 ! Revision 1.3  1998/02/04 16:09:11  raasch
    63 ! x-Wertebereiche fuer die einzelnen Koordinatenkreuze ermitteln
    64 !
    65 ! Revision 1.2  1997/09/16 06:38:20  raasch
    66 ! style- und color-Indexberechnung korrigiert
    6715!
    6816! Revision 1.1  1997/09/12 06:28:48  raasch
  • palm/trunk/SOURCE/data_output_ptseries.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_ptseries.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2006/08/22 13:51:13  raasch
    1214! Seperate output for particle groups
  • palm/trunk/SOURCE/data_output_spectra.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_spectra.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.7  2006/04/11 14:56:38  raasch
    1214! pl_spectra renamed data_output_sp
    1315!
    14 ! Revision 1.6  2006/02/23 10:29:38  raasch
    15 ! Former routine plot_spectra renamed data_output_spectra
    16 ! error number argument for handle_netcdf_error
    17 !
    18 ! Revision 1.5  2005/05/18 15:48:29  raasch
    19 ! Extensions for NetCDF output
    20 !
    21 ! Revision 1.4  2003/03/16 09:42:42  raasch
    22 ! Two underscores (_) are placed in front of all define-strings
    23 !
    24 ! Revision 1.3  2002/06/11 13:14:02  raasch
    25 ! TRIM function added (needed on IBM-systems), additional translations into
    26 ! English and further formatting
    27 !
    28 ! Revision 1.2  2001/03/30  07:45:18  07:45:18  raasch (Siegfried Raasch)
    29 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    30 !
    3116! Revision 1.1  2001/01/05 15:14:20  raasch
    3217! Initial revision
  • palm/trunk/SOURCE/data_output_tseries.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: data_output_tseries.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.13  2006/03/14 12:42:51  raasch
    1214! Error removed: NetCDF output only if switched on
    13 !
    14 ! Revision 1.12  2006/02/23 10:30:54  raasch
    15 ! Former routine plot_ts renamed data_output_tseries
    16 ! error number argument for handle_netcdf_error
    17 ! plts_anz renamed dots_n
    18 !
    19 ! Revision 1.11  2005/05/18 15:50:54  raasch
    20 ! Extensions for NetCDF output
    21 !
    22 ! Revision 1.10  2003/04/16 13:10:40  raasch
    23 ! Output of Monin Obukhov length included
    24 !
    25 ! Revision 1.9  2001/03/30 07:45:34  raasch
    26 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    27 !
    28 ! Revision 1.8  2001/01/22 20:36:59  letzel
    29 ! All comments translated into English.
    30 !
    31 ! Revision 1.7  2001/01/22 07:42:50  raasch
    32 ! Module test_variables removed
    33 !
    34 ! Revision 1.6  1999/02/05 09:14:35  raasch
    35 ! Erweiterung fuer Ausgabe von Upstream-Anteilen, Erweiterung von ts_value
    36 ! auf 30 Elemente
    37 !
    38 ! Revision 1.5  1998/07/06 12:26:59  raasch
    39 ! + USE test_variables
    40 !
    41 ! Revision 1.4  1998/04/15 11:22:22  raasch
    42 ! Erweiterung der Zeitreihenausgabe auf oberflaechennahe Temperaturen und
    43 ! Waermestroeme
    44 !
    45 ! Revision 1.3  1998/03/30 11:48:03  raasch
    46 ! Divergenzen stehen jetzt in hom
    47 !
    48 ! Revision 1.2  1998/03/25 13:56:01  raasch
    49 ! Erweiterungen fuer Gebietsstatistiken,
    50 ! dt in dt_3d umbenannt
    5115!
    5216! Revision 1.1  1998/03/03 08:00:13  raasch
  • palm/trunk/SOURCE/diffusion_e.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusion_e.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.18  2006/08/04 14:29:43  raasch
    1214! dissipation is stored in extra array diss if needed later on for calculating
    1315! the sgs particle velocities
    14 !
    15 ! Revision 1.17  2006/02/23 10:31:46  raasch
    16 ! nzb_2d replaced by nzb_s_inner
    17 !
    18 ! Revision 1.16  2004/01/30 10:18:18  raasch
    19 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    20 !
    21 ! Revision 1.15  2003/03/14 13:39:33  raasch
    22 ! Loop optimization for diffusion_e, l and ll are now automatic arrays
    23 !
    24 ! Revision 1.14  2003/03/12 16:25:03  raasch
    25 ! Full code replaced in the call for all gridpoints instead of calling the
    26 ! _ij version (required by NEC, because otherwise no vectorization)
    27 !
    28 ! Revision 1.13  2002/12/19 14:25:28  raasch
    29 ! Correction of mixing length term (l(k)/ll(k)). The condition kh=3*km in
    30 ! the unstable case is now also exactly met in the wall adjustment region.
    31 ! Factor 0.7 in wall adjustment part replaced by variable
    32 ! wall_adjustment_factor, which is set to ... in modules.f90.
    33 !
    34 ! Revision 1.12  2002/06/11 12:51:59  raasch
    35 ! Former subroutine changed to a module which allows to be called for all grid
    36 ! points of a single vertical column with index i,j or for all grid points by
    37 ! using function overloading.
    38 ! 1D-array l is allocated only once in the first call.
    39 !
    40 ! Revision 1.11  2001/08/21 08:24:34  raasch
    41 ! Wall adjustment of mixing length to 0.7 z can be switched off
    42 !
    43 ! Revision 1.10  2001/03/30 07:06:30  raasch
    44 ! Near surface mixing length is limited to 0.7*zu,
    45 ! e**1.5 replaced by e*SQRT(e) (more than 10% total increase in performance
    46 ! of this routine),
    47 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    48 !
    49 ! Revision 1.9  2001/01/22 06:05:28  raasch
    50 ! Module test_variables removed
    51 !
    52 ! Revision 1.8  2001/01/02 17:27:00  raasch
    53 ! -dpt_dz_d, dpt_dz_u
    54 !
    55 ! Revision 1.7  2000/07/03 12:56:34  raasch
    56 ! array l changed from dummy argument to local allocatable array,
    57 ! dummy arguments, whose corresponding actual arguments are pointers,
    58 ! are now also defined as pointers
    59 ! all comments translated into English
    60 !
    61 ! Revision 1.6  2000/04/18 08:10:12  schroeter
    62 ! Revision 1.4 wieder rueckgaengig gemacht, das Stabilitaets-
    63 ! kriterium basiert nun wieder auf zentralen Differenzen
    64 !
    65 ! Revision 1.5  2000/04/13 14:33:08  schroeter
    66 ! je nach Initialisierungsmodus (trocken/feucht) fliesst in die
    67 ! Berechnung des Mischungsweges pt oder vpt ein, wird durch
    68 ! entsprechende Variablenuebergabe geregelt
    69 !
    70 ! Revision 1.4  99/02/17  09:15:52  09:15:52  raasch (Siegfried Raasch)
    71 ! Dissipation jetzt gemaess dem originalen Deardorff-Ansatz
    72 ! Kriterium fuer reduzierten Mischungsweg im stabil geschichteten Fall enger
    73 ! gefasst (Schichtung muss sowohl oberhalb als auch unterhalb des betrachteten
    74 ! Gitterpunkts stabil sein)
    75 !
    76 ! Revision 1.3  1998/07/06 12:10:56  raasch
    77 ! + USE test_variables
    78 !
    79 ! Revision 1.2  1998/03/11 11:48:59  raasch
    80 ! Anpassung des Mischungsweges an den Prandtlschen Mischungsweg moeglich
    8116!
    8217! Revision 1.1  1997/09/19 07:40:24  raasch
  • palm/trunk/SOURCE/diffusion_s.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusion_s.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.8  2006/02/23 10:34:17  raasch
    1214! nzb_2d replaced by nzb_s_outer in horizontal diffusion and by nzb_s_inner
    1315! or nzb_diff_s_inner, respectively, in vertical diffusion, prescribed surface
    1416! fluxes at vertically oriented topography
    15 !
    16 ! Revision 1.7  2004/01/30 10:20:56  raasch
    17 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    18 !
    19 ! Revision 1.6  2003/03/12 16:25:32  raasch
    20 ! Full code replaced in the call for all gridpoints instead of calling the
    21 ! _ij version (required by NEC, because otherwise no vectorization)
    22 !
    23 ! Revision 1.5  2002/06/11 12:52:41  raasch
    24 ! Former subroutine changed to a module which allows to be called for all grid
    25 ! points of a single vertical column with index i,j or for all grid points by
    26 ! using function overloading.
    27 !
    28 ! Revision 1.4  2001/03/30 07:11:44  raasch
    29 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    30 !
    31 ! Revision 1.3  2001/01/25 06:58:14  raasch
    32 ! Variable "prandtl_layer replaced by "use_surface_fluxes"
    33 !
    34 ! Revision 1.2  2000/07/03 12:57:13  raasch
    35 ! dummy arguments, whose corresponding actual arguments are pointers,
    36 ! are now also defined as pointers,
    37 ! all comments translated into English
    3817!
    3918! Revision 1.1  2000/04/13 14:54:02  schroeter
  • palm/trunk/SOURCE/diffusion_u.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusion_u.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 10:35:35  raasch
    1214! nzb_2d replaced by nzb_u_outer in horizontal diffusion and by nzb_u_inner
     
    1618! because they cause errors at the edges of topography
    1719! WARNING: loops containing the MAX function are still not properly vectorized!
    18 !
    19 ! Revision 1.14  2005/03/26 20:08:21  raasch
    20 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,
    21 ! additional damping layer at the outflow in case of non-cyclic lateral
    22 ! boundaries, additional argument km_damp_y
    23 !
    24 ! Revision 1.13  2004/01/30 10:21:32  raasch
    25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    26 !
    27 ! Revision 1.12  2003/03/12 16:26:14  raasch
    28 ! Full code replaced in the call for all gridpoints instead of calling the
    29 ! _ij version (required by NEC, because otherwise no vectorization)
    30 !
    31 ! Revision 1.11  2002/06/11 12:53:12  raasch
    32 ! Former subroutine changed to a module which allows to be called for all grid
    33 ! points of a single vertical column with index i,j or for all grid points by
    34 ! using function overloading.
    35 !
    36 ! Revision 1.10  2001/08/21 08:25:55  raasch
    37 ! Special treatment at k=1 generally if momentum flux is prescribed (not only in
    38 ! case of a Prandtl layer)
    39 !
    40 ! Revision 1.9  2001/03/30 07:12:02  raasch
    41 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    42 !
    43 ! Revision 1.8  2001/01/22 06:14:00  raasch
    44 ! Module test_variables removed
    45 !
    46 ! Revision 1.7  2000/07/03 12:57:46  raasch
    47 ! dummy arguments, whose corresponding actual arguments are pointers,
    48 ! are now also defined as pointers,
    49 ! all comments translated into English
    50 !
    51 ! Revision 1.6  2000/03/14 13:55:33  schroeter
    52 ! vertical impulse fluxes are computed over dz at the
    53 ! top of the Prandtl-layer
    54 !
    55 ! Revision 1.5  98/07/06  12:11:42  12:11:42  raasch (Siegfried Raasch)
    56 ! + USE test_variables
    57 !
    58 ! Revision 1.4  1998/01/23 09:56:44  raasch
    59 ! Einbau der Prandtl-Schicht
    60 !
    61 ! Revision 1.3  1997/09/12 07:22:50  raasch
    62 ! Leerzeilen mussten entfernt werden
    63 !
    64 ! Revision 1.2  1997/09/12 06:43:27  raasch
    65 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen
    6620!
    6721! Revision 1.1  1997/09/12 06:23:51  raasch
  • palm/trunk/SOURCE/diffusion_v.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusion_v.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.15  2006/02/23 10:36:00  raasch
    1214! nzb_2d replaced by nzb_v_outer in horizontal diffusion and by nzb_v_inner
     
    1618! because they cause errors at the edges of topography
    1719! WARNING: loops containing the MAX function are still not properly vectorized!
    18 !
    19 ! Revision 1.14  2005/03/26 20:09:56  raasch
    20 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,
    21 ! additional damping layer at the outflow in case of non-cyclic lateral
    22 ! boundaries, additional argument km_damp_x
    23 !
    24 ! Revision 1.13  2004/01/30 10:21:56  raasch
    25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    26 !
    27 ! Revision 1.12  2003/03/12 16:26:31  raasch
    28 ! Full code replaced in the call for all gridpoints instead of calling the
    29 ! _ij version (required by NEC, because otherwise no vectorization)
    30 !
    31 ! Revision 1.11  2002/06/11 12:53:46  raasch
    32 ! Former subroutine changed to a module which allows to be called for all grid
    33 ! points of a single vertical column with index i,j or for all grid points by
    34 ! using function overloading.
    35 !
    36 ! Revision 1.10  2001/08/21 08:29:31  raasch
    37 ! Special treatment at k=1 generally if momentum flux is prescribed (not only in
    38 ! case of a Prandtl layer)
    39 !
    40 ! Revision 1.9  2001/03/30 07:12:19  raasch
    41 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    42 !
    43 ! Revision 1.8  2001/01/22 06:18:41  raasch
    44 ! Module test_variables removed
    45 !
    46 ! Revision 1.7  2000/07/03 12:58:14  raasch
    47 ! dummy arguments, whose corresponding actual arguments are pointers,
    48 ! are now also defined as pointers,
    49 ! all comments translated into English
    50 !
    51 ! Revision 1.6  2000/03/14 13:58:18  schroeter
    52 ! vertical impulse fluxes are computed over dz at the
    53 ! top of the Prandtl-layer
    54 !
    55 ! Revision 1.5  98/07/06  12:12:05  12:12:05  raasch (Siegfried Raasch)
    56 ! + USE test_variables
    57 !
    58 ! Revision 1.4  1998/01/23 09:57:18  raasch
    59 ! Einbau der Prandtl-Schicht
    60 !
    61 ! Revision 1.3  1997/09/12 07:24:20  raasch
    62 ! Leerzeilen mussten entfernt werden
    63 !
    64 ! Revision 1.2  1997/09/12 06:43:46  raasch
    65 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen
    6620!
    6721! Revision 1.1  1997/09/12 06:24:01  raasch
  • palm/trunk/SOURCE/diffusion_w.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusion_w.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.12  2006/02/23 10:38:03  raasch
    1214! nzb_2d replaced by nzb_w_outer, wall functions added for all vertical walls,
    1315! +z0 in argument list
    1416! WARNING: loops containing the MAX function are still not properly vectorized!
    15 !
    16 ! Revision 1.11  2005/03/26 20:11:05  raasch
    17 ! Additional damping layer at the outflow in case of non-cyclic lateral
    18 ! boundaries, additional arguments km_damp_x, km_damp_y
    19 !
    20 ! Revision 1.10  2004/01/30 10:22:32  raasch
    21 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    22 !
    23 ! Revision 1.9  2003/03/12 16:26:49  raasch
    24 ! Full code replaced in the call for all gridpoints instead of calling the
    25 ! _ij version (required by NEC, because otherwise no vectorization)
    26 !
    27 ! Revision 1.8  2002/06/11 12:54:16  raasch
    28 ! Former subroutine changed to a module which allows to be called for all grid
    29 ! points of a single vertical column with index i,j or for all grid points by
    30 ! using function overloading.
    31 !
    32 ! Revision 1.7  2001/03/30 07:12:35  raasch
    33 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    34 !
    35 ! Revision 1.6  2001/01/22 06:20:12  raasch
    36 ! Module test_variables removed
    37 !
    38 ! Revision 1.5  2000/07/03 12:58:40  raasch
    39 ! dummy arguments, whose corresponding actual arguments are pointers,
    40 ! are now also defined as pointers,
    41 ! all comments translated into English
    42 !
    43 ! Revision 1.4  1998/07/06 12:12:24  raasch
    44 ! + USE test_variables
    45 !
    46 ! Revision 1.3  1997/09/12 07:24:58  raasch
    47 ! Leerzeilen mussten entfernt werden
    48 !
    49 ! Revision 1.2  1997/09/12 06:44:00  raasch
    50 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen
    5117!
    5218! Revision 1.1  1997/09/12 06:24:11  raasch
  • palm/trunk/SOURCE/diffusivities.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: diffusivities.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.24  2006/04/26 12:16:26  raasch
    1214! OpenMP optimization (+sums_l_l_t), sqrt_e must be private
    13 !
    14 ! Revision 1.23  2006/02/23 12:13:29  raasch
    15 ! nzb_2d replaced by nzb_s_inner/outer, option of a minimum tke
    16 ! quasi boundary conditions is now set for sums_l_l instead of sums_l
    17 !
    18 ! Revision 1.22  2005/03/26 20:16:26  raasch
    19 ! Eddy diffusivities are also calculated for the ghost points in order to spare
    20 ! the communication otherwise needed for the exchange of these points.
    21 ! Neumann boundary conditions set at the outflow boundaries in case of
    22 ! non-cyclic lateral boundaries
    23 !
    24 ! Revision 1.21  2004/01/30 10:23:12  raasch
    25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    26 !
    27 ! Revision 1.20  2003/03/14 13:40:09  raasch
    28 ! Loop optimization (SQRT(e) solved in a seperate loop, joining of loops)
    29 !
    30 ! Revision 1.19  2003/03/12 16:28:56  raasch
    31 ! phi_m is initialized in declaration statement in order to avoid run time
    32 ! errors due to compiler problems on ibm and nec
    33 !
    34 ! Revision 1.18  2002/12/19 14:32:23  raasch
    35 ! Correction of mixing length term (l(k)/ll(k)). The condition kh=3*km in
    36 ! the unstable case is now also exactly met in the wall adjustment region.
    37 ! Factor 0.7 in wall adjustment part replaced by variable
    38 ! wall_adjustment_factor, which is set to 1.8 in modules.f90.
    39 ! Missing OMP DO statement added.
    40 !
    41 ! Revision 1.17  2002/06/11 12:56:49  raasch
    42 ! OpenMP directives added,
    43 ! array l changed from allocatable to automatic, because it must be on the stack
    44 ! in case of using OpenMP
    45 !
    46 ! Revision 1.16  2001/08/21  08:30:34  08:30:34  raasch (Siegfried Raasch)
    47 ! Wall adjustment of mixing length to 0.7 z can be switched off
    48 ! Lower boundary condition for km and kh changed to Neumann
    49 !
    50 ! Revision 1.15  2001/03/30 07:22:09  raasch
    51 ! Near surface mixing length is limited to 0.7*zu,
    52 ! if adjust_mixing_length=T, l is modified at all gridpoints (previously at
    53 ! k=1 only),
    54 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    55 !
    56 ! Revision 1.14  2001/01/22 06:31:36  raasch
    57 ! Module test_variables removed
    58 !
    59 ! Revision 1.13  2001/01/02 17:27:41  raasch
    60 ! -dpt_dz_d, dpt_dz_u
    61 !
    62 ! Revision 1.12  2000/12/20 10:35:47  letzel
    63 ! All comments translated into English.
    64 !
    65 ! Revision 1.11  2000/04/18 08:10:42  schroeter
    66 ! Revision 1.4 wieder rueckgaengig gemacht, das Stabilitaets-
    67 ! kriterium basiert nun wieder auf zentralen Differenzen
    68 !
    69 ! Revision 1.9  2000/04/13 13:42:14  schroeter
    70 ! je nach Initialisierungsmodus (trocken/feucht) fliesst in
    71 ! Mischungswegberechnung pt oder vpt ein, wird durch entsprechende
    72 ! Variablenuebergabe geregelt
    73 !
    74 ! Revision 1.8  99/02/17  09:17:39  09:17:39  raasch (Siegfried Raasch)
    75 ! Kriterium fuer reduzierten Mischungsweg im stabil geschichteten Fall enger
    76 ! gefasst (Schichtung muss sowohl oberhalb als auch unterhalb des betrachteten
    77 ! Gitterpunkts stabil sein)
    78 !
    79 ! Revision 1.7  1998/09/22 17:18:41  raasch
    80 ! Testweise cm = 0.4 bei k=1, aber vorerst auskommentiert
    81 !
    82 ! Revision 1.6  1998/08/05 06:52:44  raasch
    83 ! Mischungsweganpassung an Prandtl nur bei k=1
    84 !
    85 ! Revision 1.2  1998/03/11 11:49:26  raasch
    86 ! Anpassung des Mischungsweges an den Prandtlschen Mischungsweg moeglich
    8715!
    8816! Revision 1.1  1997/09/19 07:41:10  raasch
  • palm/trunk/SOURCE/disturb_field.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: disturb_field.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2006/08/04 14:31:59  raasch
    1214! izuf renamed iran
    13 !
    14 ! Revision 1.10  2006/02/23 12:14:06  raasch
    15 ! Additional argument nzb_uv_inner (no perturbations below the topography)
    16 !
    17 ! Revision 1.9  2005/03/26 20:18:30  raasch
    18 ! Arguments for non-cyclic boundary conditions added to argument list and to
    19 ! argument list of routine exchange_horiz.
    20 ! Perturbations can be imposed for different ranges (steered by dist_range).
    21 !
    22 ! Revision 1.8  2003/10/29 08:48:51  raasch
    23 ! Module random_function_mod is used
    24 !
    25 ! Revision 1.7  2003/03/14 13:41:48  raasch
    26 ! Different random number generators available
    27 !
    28 ! Revision 1.6  2001/03/30 07:22:26  raasch
    29 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    30 !
    31 ! Revision 1.5  2001/01/25 06:59:38  raasch
    32 ! Module interfaces added (missing this interface caused errors on machines
    33 ! other than Cray)
    34 !
    35 ! Revision 1.4  2001/01/22 06:32:49  raasch
    36 ! Module test_variables removed
    37 !
    38 ! Revision 1.3  2000/07/03 12:59:03  raasch
    39 ! All comments translated into English
    40 !
    41 ! Revision 1.2  1998/07/06 12:13:01  raasch
    42 ! + USE test_variables
    4315!
    4416! Revision 1.1  1998/02/04 15:40:45  raasch
  • palm/trunk/SOURCE/disturb_heatflux.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: disturb_heatflux.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.7  2006/08/04 14:35:07  raasch
    1214! Additional parameter in function random_gauss which limits the range of the
    1315! created random numbers, izuf renamed iran
    14 !
    15 ! Revision 1.6  2006/02/23 12:14:38  raasch
    16 ! Over topography surface_heatflux is replaced by wall_heatflux(0)
    17 !
    18 ! Revision 1.5  2001/03/30 07:22:44  raasch
    19 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    20 !
    21 ! Revision 1.4  2001/01/22 06:34:34  raasch
    22 ! Module test_variables removed
    23 !
    24 ! Revision 1.3  2000/12/20 11:57:38  letzel
    25 ! All comments translated into English.
    26 !
    27 ! Revision 1.2  1998/07/06 12:13:19  raasch
    28 ! + USE test_variables
    2916!
    3017! Revision 1.1  1998/03/25 20:03:47  raasch
  • palm/trunk/SOURCE/exchange_horiz.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: exchange_horiz.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.16  2006/02/23 12:19:08  raasch
    1214! anz_yz renamed ngp_yz
    13 !
    14 ! Revision 1.15  2005/03/26 20:20:39  raasch
    15 ! Extensions for non-cyclic boundary conditions along x or y
    16 !
    17 ! Revision 1.14  2003/10/29 08:51:10  raasch
    18 ! Exchange of boundaries within the PE memory in case that multigrid method
    19 ! has switched to PE0 only
    20 !
    21 ! Revision 1.13  2003/03/16 09:30:38  raasch
    22 ! Two underscores (_) are placed in front of all define-strings
    23 !
    24 ! Revision 1.12  2002/12/20 10:40:08  raasch
    25 ! Integer variables will be defined in the parallel case only
    26 !
    27 ! Revision 1.11  2002/12/19 14:46:43  raasch
    28 ! Cyclic boundary conditions along y used instead of sendrecv in case of
    29 ! pdims(2)=1. SENDRECV replaced by nonblocking routines ISEND and IRECV.
    30 ! Array "feld" is renamed to "ar".
    31 !
    32 ! Revision 1.9  2002/04/16  08:16:50  08:16:50  raasch (Siegfried Raasch)
    33 ! Cyclic boundary conditions along x used instead of sendrecv in case of
    34 ! pdims(1)=1
    35 !
    36 ! Revision 1.8  2001/07/20 13:04:59  raasch
    37 ! anz_yz and type_xz changed to arrays to allow for the boundary value exchange
    38 ! for the different grid levels used in the multigrid method,
    39 ! +module control_parameters
    40 !
    41 ! Revision 1.7  2001/03/30 07:23:00  raasch
    42 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    43 !
    44 ! Revision 1.6  2001/01/22 06:36:07  raasch
    45 ! Module test_variables removed
    46 !
    47 ! Revision 1.5  2000/12/20 12:06:44  letzel
    48 ! All comments translated into English.
    49 !
    50 ! Revision 1.4  1998/07/06 12:13:36  raasch
    51 ! + USE test_variables
    52 !
    53 ! Revision 1.3  1998/01/23 09:57:45  raasch
    54 ! Umbenennung von anz_yz_p in anz_yz sowie von type_xz_p in type_xz
    55 !
    56 ! Revision 1.2  1997/08/11 06:14:55  raasch
    57 ! Felder werden jetzt ueber Formalparameter uebergeben.
    5815!
    5916! Revision 1.1  1997/07/24 11:13:29  raasch
  • palm/trunk/SOURCE/exchange_horiz_2d.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: exchange_horiz_2d.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.9  2006/05/12 19:15:52  letzel
    1214! MPI_REAL replaced by MPI_INTEGER in exchange_horiz_2d_int
    13 !
    14 ! Revision 1.8  2006/02/23 12:18:32  raasch
    15 ! Additional subroutine exchange_horiz_2d_int for 2D integer arrays,
    16 ! extensions for non-cyclic boundary conditions along x or y for non-parallel
    17 ! case, anz_y renamed ngp_y
    18 !
    19 ! Revision 1.7  2003/03/16 09:30:43  raasch
    20 ! Two underscores (_) are placed in front of all define-strings
    21 !
    22 ! Revision 1.6  2002/06/11 12:59:35  raasch
    23 ! Cyclic boundary conditions are used instead of sendrecv in case of
    24 ! pdims(..)=1. Array "feld" is renamed to "ar".
    25 !
    26 ! Revision 1.5  2001/03/30  07:23:16  07:23:16  raasch (Siegfried Raasch)
    27 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    28 !
    29 ! Revision 1.4  2001/01/22 06:43:50  raasch
    30 ! Module test_variables removed
    31 !
    32 ! Revision 1.3  2000/12/20 12:09:27  letzel
    33 ! All comments translated into English.
    34 !
    35 ! Revision 1.2  1998/07/06 12:13:53  raasch
    36 ! + USE test_variables
    3715!
    3816! Revision 1.1  1998/01/23 09:58:21  raasch
  • palm/trunk/SOURCE/fft_xy.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: fft_xy.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.4  2006/03/28 12:27:09  raasch
    1214! Stop when system-specific fft is selected on NEC. For unknown reasons this
    1315! causes a program abort during first allocation in init_grid.
    14 !
    15 ! Revision 1.3  2006/02/23 12:22:41  raasch
    16 ! System specfic routines are not allowed to be used on ibmy
    1716!
    1817! Revision 1.2  2004/04/30 11:44:27  raasch
     
    2524! decalpha),
    2625! non vectorized FFT for nec included
    27 !
    28 ! Revision 1.1  2004/04/30 11:37:14  raasch
    29 ! Initial revision
    30 !
    31 ! Revision 1.6  2003/08/01 08:15:09  raasch
    32 ! No abort in fft_1dd_init on t3e-systems if system specific routines are used
    33 !
    34 ! Revision 1.5  2003/04/16 12:51:47  raasch
    35 ! Temperton-fft implemented in routines fft_x_1dd and fft_y_1dd
    36 !
    37 ! Revision 1.4  2003/03/16 09:37:26  raasch
    38 ! Two underscores (_) are placed in front of all define-strings
    39 !
    40 ! Revision 1.3  2003/03/12 16:29:39  raasch
    41 ! Routines fft_x_1dd_m and fft_y_1dd_m added (suitable for multiple 1d-fft on
    42 ! vector processors)
    43 !
    44 ! Revision 1.2  2002/12/19 14:49:18  raasch
    45 ! STOP statement replaced by call of subroutine local_stop
    4626!
    4727! Revision 1.1  2002/06/11 13:00:49  raasch
  • palm/trunk/SOURCE/flow_statistics.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: flow_statistics.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.41  2006/08/04 14:37:50  raasch
    1214! Error removed in non-parallel part (sums_l)
    13 !
    14 ! Revision 1.40  2006/04/26 14:06:01  raasch
    15 ! Call of omp_get_thread_num written as comment line. Since lcmuk does not
    16 ! understand this comment, a cpp-directive is required.
    17 !
    18 ! Revision 1.39  2006/04/26 12:19:16  raasch
    19 ! OpenMP optimization, thread number added as third dimension of sums_l,
    20 ! nzb_s_inner(j,i) replaced by nzb in determination of z_i
    21 !
    22 ! Revision 1.38  2006/02/23 12:23:15  raasch
    23 ! nzb and nzb_diff replaced by nzb_s_outer and nzb_diff_s_outer resrectively,
    24 ! nanz_sums, nanz_2dh, nanz_3d renamed ngp_sums, ngp_2dh and ngp_3d,
    25 ! respectively,
    26 ! ngp_2dh and ngp_3d partially replaced by ngp_2dh_outer and ngp_3d_inner,
    27 ! respectively
    28 !
    29 ! Revision 1.37  2005/12/06 16:40:21  raasch
    30 ! Calculation of ql (54) profile in case of using cloud droplets
    31 !
    32 ! Revision 1.36  2005/04/23 09:10:27  raasch
    33 ! +local array sums_ll instead of temporarily using sums_l
    34 !
    35 ! Revision 1.35  2004/04/30 11:46:56  raasch
    36 ! leap_frog changed to time_integration
    37 !
    38 ! Revision 1.34  2003/05/09 14:25:12  raasch
    39 ! Parameters in MPI_ALLREDUCE call changed
    40 !
    41 ! Revision 1.33  2003/03/31 15:42:53  raasch
    42 ! Calculation of horizontal temperature fluxes modified and cleaned up
    43 !
    44 ! Revision 1.32  2003/03/16 09:39:07  raasch
    45 ! Two underscores (_) are placed in front of all define-strings
    46 !
    47 ! Revision 1.31  2003/03/12 16:30:07  raasch
    48 ! Calculation of perturbation and total energy of the total domain
    49 ! re-arranged in order to allow vectorization of the respective loop
    50 !
    51 ! Revision 1.30  2003/03/04 11:26:35  raasch
    52 ! Error in calculation of the vertical flux of resolved scale energy
    53 ! (profile 57) removed
    54 !
    55 ! Revision 1.29  2002/12/19 14:50:02  raasch
    56 ! STOP statement replaced by call of subroutine local_stop
    57 !
    58 ! Revision 1.28  2002/04/16 08:04:19  raasch
    59 ! Error in passive scalar profile (hom(..41..)) removed
    60 !
    61 ! Revision 1.27  2002/04/11 16:31:03  schroeter
    62 ! No changes.
    63 !
    64 ! Revision 1.26  2001/09/04  12:00:12  12:00:12  raasch (Siegfried Raasch)
    65 ! Calculation of the vertical flux of resolved scale energy and pressure
    66 ! fluctuations
    67 !
    68 ! Revision 1.25  2001/08/21 08:32:24  raasch
    69 ! Special calculation of surface fluxes whenever they are prescribed (not only
    70 ! in case of a Prandtl layer), call of user_interface
    71 !
    72 ! Revision 1.24  2001/03/30 07:23:34  raasch
    73 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    74 !
    75 ! Revision 1.23  2001/02/07 06:48:25  raasch
    76 ! Error in non-parallel passive scalar code removed
    77 !
    78 ! Revision 1.22  2001/01/29 12:25:02  raasch
    79 ! Passive scalar is considered
    80 !
    81 ! Revision 1.21  2001/01/22 06:49:27  raasch
    82 ! Module test_variables removed
    83 !
    84 ! Revision 1.20  2000/04/13 13:54:15  schroeter
    85 ! + computation of new profiles for calculations with
    86 !   cloud_physics
    87 !
    88 ! Revision 1.19  2000/01/26  09:32:58  09:32:58  letzel (Marcus Letzel)
    89 ! All comments translated into English
    90 !
    91 ! Revision 1.16  1998/04/21 15:54:54  raasch
    92 ! Berechnung von w* beruecksichtigt jetzt den Fall z_i=0,
    93 ! Faktor 0.5 bei Energietransport hinzugefuegt
    94 !
    95 ! Revision 1.6  1998/02/10 15:06:02  raasch
    96 ! Berechnung des Gesamtwaermestroms, in aufgeloesten Waermestroemen wird das
    97 ! horizontale Temperaturmittel nicht mehr abgezogen, Berechnung der
    98 ! Gesamtimpulsfluesse
    9915!
    10016! Revision 1.1  1997/08/11 06:15:17  raasch
  • palm/trunk/SOURCE/global_min_max.f90

    r1 r3  
    99! Former revisions:
    1010! -----------------
    11 ! $Log: global_min_max.f90,v $
     11! $Id$
     12! RCS Log replace by Id keyword, revision history cleaned up
     13!
    1214! Revision 1.11  2003/04/16 12:56:58  raasch
    1315! Index values of the extrema are limited to the range 0..nx, 0..ny
    14 !
    15 ! Revision 1.10  2003/03/16 09:39:21  raasch
    16 ! Two underscores (_) are placed in front of all define-strings
    17 !
    18 ! Revision 1.9  2003/03/14 13:42:50  raasch
    19 ! IF clause extracted from loop which determines the absolut array maximum,
    20 ! because otherwise it prevents vectorization of that loop
    21 !
    22 ! Revision 1.8  2002/12/19 14:51:03  raasch
    23 ! Speed optimization by removing MINVAL/MAXVAL calls and by handling
    24 ! the "abs" case in a different way than the min/max cases,
    25 ! translation of remaining variables
    26 !
    27 ! Revision 1.7  2002/06/11 13:08:16  raasch
    28 ! Last change limited to IBM since it does not work on 64-bit machines
    29 !
    30 ! Revision 1.6  2002/05/02  18:50:36  18:50:36  raasch (Siegfried Raasch)
    31 ! Type of fmax, fmax_l, fmin and fmin_l changed to REAL (KIND=4).
    32 !
    33 ! Revision 1.5  2001/01/22 06:51:02  raasch
    34 ! Module test_variables removed
    35 !
    36 ! Revision 1.4  2000/12/20 12:19:57  letzel
    37 ! All comments translated into English.
    38 !
    39 ! Revision 1.3  1998/07/06 12:14:37  raasch
    40 ! + USE test_variables
    41 !
    42 ! Revision 1.2  1997/08/11 06:16:54  raasch
    43 ! Indices des ermittelten Minumums/Maximums werden an PEs gesendet
    4416!
    4517! Revision 1.1  1997/07/24 11:14:03  raasch
  • palm/trunk/SOURCE/header.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: header.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.63  2006/08/22 13:53:13  raasch
    1214! Output of dz_max
    13 !
    14 ! Revision 1.62  2006/08/04 14:38:41  raasch
    15 ! Generation of the run description header as well as run date and run time
    16 ! moved to routine check_parameters,
    17 ! output for additional particle quantities
    18 !
    19 ! Revision 1.61  2006/04/26 12:19:51  raasch
    20 ! Output of threads per task
    21 !
    22 ! Revision 1.60  2006/04/11 14:57:01  raasch
    23 ! pl_spectra renamed data_output_sp
    24 !
    25 ! Revision 1.59  2006/03/14 12:44:10  raasch
    26 ! Output of the initial geostrophic wind profile corrected
    27 !
    28 ! Revision 1.58  2006/02/23 12:24:22  raasch
    29 ! Output of psl, psr, pdx, etc. for each particle group, output of e_min,
    30 ! information about volume flow conservation and usage of NetCDF 64 bit offset
    31 ! format, pl.. renamed do.., skip_time control parameters for do..,
    32 ! ebene renamed section, .._anz renamed .._no, dt_average renamed
    33 ! dt_averaging_input_pr
    34 ! output of 2d/3d data averaging informations, buffer is forced to be written
    35 ! on file at the end of the routine
    36 !
    37 ! Revision 1.57  2005/10/20 14:06:44  raasch
    38 ! Error in output of yz-slice information removed
    39 !
    40 ! Revision 1.56  2005/07/01 07:48:14  steinfeld
    41 ! Information on the initial profiles of ug and vg added;
    42 ! dependency of ug and vg on height considered in the output of
    43 ! the boundary conditions for u and v
    44 !
    45 ! Revision 1.55  2005/06/26 19:53:15  raasch
    46 ! gas_constant renamed r_d, latent_heat renamed l_v, radius is used instead
    47 ! of diameter, output formats 335+337 changed
    48 !
    49 ! Revision 1.54  2005/05/19 08:32:27  raasch
    50 ! Error in output of particle inertia information removed
    51 !
    52 ! Revision 1.53  2005/05/18 15:31:59  raasch
    53 ! Informations about data output format (netcdf, iso2d, etc.) are given
    54 !
    55 ! Revision 1.52  2005/04/23 09:16:10  raasch
    56 ! fcl_factor renamed cfl_factor
    57 !
    58 ! Revision 1.51  2005/03/26 20:23:14  raasch
    59 ! Output of horizontal boundary conditions, output of additional call of
    60 ! pressure solver in case of Runge-Kutta schemes,
    61 ! calculate the number of particle groups (formerly unknown = 0 )
    62 !
    63 ! Revision 1.50  2004/04/30 11:48:02  raasch
    64 ! Forcing a 1d-decomposition is output on linux-cluster and decalpha,
    65 ! output of subdomain sizes, impulse_advec renamed momentum_advec
    66 !
    67 ! Revision 1.49  2004/01/28 15:09:15  raasch
    68 ! Output of timestep schemes modified
    69 !
    70 ! Revision 1.48  2003/10/29 08:52:04  raasch
    71 ! Additional output for multigrid method (lowest levels gathered on PE0),
    72 ! particle output modified for particle groups
    73 !
    74 ! Revision 1.47  2003/04/16 12:58:18  raasch
    75 ! Output for mixing length limitations revised
    76 !
    77 ! Revision 1.46  2003/03/16 09:39:31  raasch
    78 ! Two underscores (_) are placed in front of all define-strings
    79 !
    80 ! Revision 1.45  2003/03/14 13:43:21  raasch
    81 ! Informations about random generator
    82 !
    83 ! Revision 1.44  2003/03/12 16:31:45  raasch
    84 ! Small change for NEC system
    85 !
    86 ! Revision 1.43  2002/12/19 14:53:47  raasch
    87 ! Informations about user defined restart times added
    88 !
    89 ! Revision 1.42  2002/09/12 13:03:24  raasch
    90 ! Informations about particle inertia added
    91 !
    92 ! Revision 1.41  2002/06/11 13:09:11  raasch
    93 ! Information about hybrid solver, usage of 1d-decomposition and fast version
    94 ! of prognostic_equations included,
    95 ! no output of mixing length adjustments in case of constant diffusion
    96 !
    97 ! Revision 1.40  2002/04/16  08:06:13  08:06:13  raasch (Siegfried Raasch)
    98 ! Informations about particle boundary conditions and output of particle
    99 ! data added, direct output of informations about scalar transport (instead
    100 ! of using humidity)
    101 !
    102 ! Revision 1.39  2001/08/21 08:50:54  raasch
    103 ! Informations about mixing length adjustment, particle tails and dvrp PE usage
    104 ! included
    105 !
    106 ! Revision 1.38  2001/07/20 13:05:39  raasch
    107 ! Information about multigrid method included
    108 !
    109 ! Revision 1.37  2001/07/12 12:08:14  raasch
    110 ! Additional informations about dvrp-output (host and filenames) and
    111 ! particles (maximum age, random start positions), + module particle_attributes
    112 !
    113 ! Revision 1.36  2001/03/30 07:23:53  raasch
    114 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    115 !
    116 ! Revision 1.35  2001/01/29 12:25:55  raasch
    117 ! Informations about using passive scalar
    118 !
    119 ! Revision 1.34  2001/01/25 07:00:48  raasch
    120 ! Information about using surface fluxes and fft-method added
    121 !
    122 ! Revision 1.33  2001/01/05 15:11:41  raasch
    123 ! Spectra informations added. Old revision remarks deleted.
    124 !
    125 ! Revision 1.32  2001/01/03 13:08:47  letzel
    126 ! File output translated into English.
    127 !
    128 ! Revision 1.31  2000/12/28 13:04:54  raasch
    129 ! Slicer information added to dvrp-output.
    130 ! Informations about optionally used packages are put in cpp-define-brackets.
    131 !
    132 ! Revision 1.30  2000/04/27 07:08:14  raasch
    133 ! notice if processor topology is predefined by user,
    134 ! general notice about dvrp-output (including isosurface and particles),
    135 ! all comments translated into English
    13615!
    13716! Revision 1.1  1997/08/11 06:17:20  raasch
  • palm/trunk/SOURCE/impact_of_latent_heat.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: impact_of_latent_heat.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2004/01/30 10:25:59  raasch
    1214! Scalar lower k index nzb replaced by 2d-array nzb_2d
    13 !
    14 ! Revision 1.4  2003/03/12 16:32:14  raasch
    15 ! Full code replaced in the call for all gridpoints instead of calling the
    16 ! _ij version (required by NEC, because otherwise no vectorization)
    17 !
    18 ! Revision 1.3  2002/06/11 13:10:06  raasch
    19 ! Former subroutine changed to a module which allows to be called for all grid
    20 ! points of a single vertical column with index i,j or for all grid points by
    21 ! using function overloading.
    22 !
    23 ! Revision 1.2  2001/01/22 06:52:34  raasch
    24 ! Module test_variables removed
    2515!
    2616! Revision 1.1  2000/04/13 14:48:40  schroeter
  • palm/trunk/SOURCE/init_1d_model.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_1d_model.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.21  2006/06/02 15:19:57  raasch
    1214! cpp-directives extended for lctit
    1315!
    14 ! Revision 1.20  2006/04/20 13:16:24  letzel
    15 ! Errors removed in the computation of the diabatic mixing length
    16 !
    17 ! Revision 1.19  2006/03/14 12:48:52  raasch
    18 ! Error in initializing l_black removed,
    19 ! maximum possible timestep for runge schemes reduced, output for units 15 and
    20 ! 17 is flushed
    21 !
    22 ! Revision 1.18  2006/02/23 12:27:38  raasch
    23 ! Switches for calculating the mixing length and the dissipation,
    24 ! remaining e1d_m**1.5 replaced by e1d_m*SQRT(e1d_m),
    25 ! nt_anz_1d renamed current_timestep_number_1d
    26 !
    27 ! Revision 1.17  2005/06/29 09:58:25  steinfeld
    28 ! A potential dependency of ug and vg on height is considered in the
    29 ! tendency terms of u and v; calculation of lambda is commented out
    30 !
    31 ! Revision 1.16  2005/04/23 09:22:04  raasch
    32 ! Initial horizontal velocities at the lowest grid levels are set to a very
    33 ! small value in order to avoid too small time steps caused by the diffusion
    34 ! criterion
    35 !
    36 ! Revision 1.15  2005/03/26 20:27:33  raasch
    37 ! Output format of timestep changed, e1d_m**1.5 replaced by e1d_m*SQRT(e1d_m)
    38 ! because of wrong results on NEC-SX6
    39 !
    40 ! Revision 1.14  2004/04/30 11:55:02  raasch
    41 ! Runge-Kutta scheme implemented, routine leap_frog_1d renamed to
    42 ! time_integration_1d
    43 !
    44 ! Revision 1.13  2001/03/30 07:24:10  raasch
    45 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    46 !
    47 ! Revision 1.12  2001/01/29 12:26:57  raasch
    48 ! Passive scalar is considered
    49 !
    50 ! Revision 1.11  2001/01/22 06:54:13  raasch
    51 ! Module test_variables removed
    52 !
    53 ! Revision 1.10  2001/01/02 17:28:35  raasch
    54 ! Adjustment of computations of energy dissipation, mixing length and
    55 ! eddy diffusivities to methods used in the 3D-model:
    56 ! cm changed from 0.4 to 0.1, grid dependend mixing length instead of
    57 ! Blackadar mixing length is used, Deardorff modification of mxixing length
    58 ! when stratification is stable
    59 !
    60 ! Revision 1.9  2000/12/21 14:16:15  letzel
    61 ! All comments translated into English.
    62 !
    63 ! Revision 1.8  2000/04/13 13:23:47  schroeter
    64 ! + 1D-Modell an Rechnungen mit Feuchte angepasst
    65 !
    66 ! Revision 1.7  2000/03/14  12:56:02  12:56:02  schroeter (Michael Schroeter)
    67 ! Asselin-filter is called before swapping timelevels
    68 !
    6916! Revision 1.1  1998/03/09 16:22:10  raasch
    7017! Initial revision
  • palm/trunk/SOURCE/init_3d_model.f90

    r1 r3  
    1111! Former revisions:
    1212! -----------------
    13 ! $Log: init_3d_model.f90,v $
     13! $Id$
     14! RCS Log replace by Id keyword, revision history cleaned up
     15!
    1416! Revision 1.49  2006/08/22 15:59:07  raasch
    1517! No optimization of this file on the ibmy (Yonsei Univ.)
    16 !
    17 ! Revision 1.48  2006/08/22 13:58:30  raasch
    18 ! Initialisation of all tendency arrays (t.._m) needed for the Runge-Kutta
    19 ! schemes
    20 !
    21 ! Revision 1.47  2006/08/04 14:39:20  raasch
    22 ! Allocation of diss, needed for calculating the sgs particle velocities
    23 !
    24 ! Revision 1.46  2006/05/15 08:14:56  letzel
    25 ! Initialize volume_flow_initial_l and volume_flow_area_l with 0.0
    26 !
    27 ! Revision 1.45  2006/04/26 13:32:44  raasch
    28 ! Provisional correction for Piacsek & Williams advection scheme: keep u and
    29 ! v zero one layer below the topography (in case of set_1d-model_profiles).
    30 !
    31 ! Revision 1.44  2006/04/26 12:22:09  raasch
    32 ! sums_l and sums_l_l changed to three dimensions
    33 !
    34 ! Revision 1.43  2006/02/23 12:29:11  raasch
    35 ! Allocation and default setting of k index arrays nzb_2d and nzb_diff_2d
    36 ! moved to routine init_grid, nanz_sums, nanz_2dh and nanz_3d renamed ngp_sums,
    37 ! ngp_2dh and ngp_3d, respectively,
    38 ! allocation and setting of ngp_2dh and ngp_3d moved from parin and init_pegrid
    39 ! to here,
    40 ! read_fields_from_prior_run renamed read_restart_data,
    41 ! new grid point counters ngp_2dh_outer and ngp_3d_inner are calculated,
    42 ! additional argument nzb_.._inner added to call of disturb_field
    43 ! optional calculation of the initial volume flow through the right and north
    44 ! boundary
    45 !
    46 ! Revision 1.42  2005/06/26 19:55:01  raasch
    47 ! ql is now a pointer, init_cloud_physics is called anyway
    48 !
    49 ! Revision 1.41  2005/03/26 20:29:41  raasch
    50 ! Arguments for non-cyclic boundary conditions added in call of routine
    51 ! disturb_field, diffusivities km_damp_x and km_damp_y used within the outflow
    52 ! damping layer are allocated and initialized, uvmean_outflow is allocated
    53 !
    54 ! Revision 1.40  2004/04/30 11:59:05  raasch
    55 ! Reading of restart data is moved to new routine read_3d_binary,
    56 ! diffusion quantity arrays are allocated for previous timelevels only
    57 ! in case of leapfrog scheme, enlarged transposition array d
    58 !
    59 ! Revision 1.39  2004/01/30 10:27:20  raasch
    60 ! Allocation and default setting of index arrays nzb_2d and nzb_diff_2d
    61 !
    62 ! Revision 1.38  2004/01/28 15:09:56  raasch
    63 ! Time levels for Runge-Kutta scheme initialized, setting of bt before
    64 ! calling pres removed
    65 !
    66 ! Revision 1.37  2003/10/29 08:54:05  raasch
    67 ! Random numbers are initialized or random seed is read from file of
    68 ! previous run, module random_function_mod is used, new binary version 2.2
    69 !
    70 ! Revision 1.36  2003/08/01 08:25:08  raasch
    71 ! Array sums is initialized (must be defined in the first call of pres),
    72 ! call of init_cloud_physics moved before init_particles
    73 !
    74 ! Revision 1.35  2003/04/09 11:50:50  raasch
    75 ! In case of initialisation with constant profiles, the horizontal
    76 ! velocities at k=0,1 are explicitly set to zero
    77 !
    78 ! Revision 1.34  2003/03/16 09:39:39  raasch
    79 ! Two underscores (_) are placed in front of all define-strings
    80 !
    81 ! Revision 1.33  2002/12/19 15:32:42  raasch
    82 ! Subroutine local_tremain_ini (initializing the remaining cpu-time)
    83 ! moved to main program (palm), reading of array hom is moved to routine
    84 ! read_var_list, STOP statement replaced by call of subroutine local_stop
    85 !
    86 ! Revision 1.32  2002/05/02 18:51:11  raasch
    87 ! Time level t+dt re-introduced and initialized
    88 !
    89 ! Revision 1.31  2001/11/05  09:31:38  09:31:38  weinbrec (Sonja Weinbrecht)
    90 ! Initialization of spectrum_x and spectrum_y for continuation runs
    91 !
    92 ! Revision 1.30  2001/03/30 07:25:25  raasch
    93 ! Arguments removed from pres, arrays work and work1 removed,
    94 ! initialization of timelevel t+dt removed,
    95 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    96 !
    97 ! Revision 1.29  2001/01/29 12:27:49  raasch
    98 ! Passive scalar is considered
    99 !
    100 ! Revision 1.28  2001/01/25 07:02:19  raasch
    101 ! Initialization of surface fluxes changed due to new variable
    102 ! use_surface_fluxes
    103 !
    104 ! Revision 1.27  2001/01/05 15:12:23  raasch
    105 ! Allocate spectrum_x and spectrum_y, new binary version 2.1
    106 !
    107 ! Revision 1.26  2001/01/02 17:29:43  raasch
    108 ! Unit 13 closed after reading binary input (after call of user_init, where
    109 ! read from unit 13 may occur)
    110 !
    111 ! Revision 1.25  2000/12/28 13:05:53  raasch
    112 ! Init_particles is called unconditionally
    113 !
    114 ! Revision 1.24  2000/04/27 07:05:19  raasch
    115 ! Call of subroutine init_slope added to compute the initial temperature field
    116 ! and other constants used in case of a sloping surface,
    117 ! call of subroutine init_dvrp added to initialize dvrp-software,
    118 ! old revision remarks deleted
    119 !
    120 ! Revision 1.23  2000/04/13 13:22:09  schroeter
    121 ! Allokieren und Initialisierung der Feuchtefelder und Wolkenphysikparameter
    122 !
    123 ! Revision 1.22  2000/01/26  10:30:49  10:30:49  letzel (Marcus Letzel)
    124 ! All comments translated into English
    125 !
    126 ! Revision 1.21  2000/01/10 09:54:10  raasch
    127 ! Trennung der zeitabhaengigen Variablen in pointer und targets,
    128 ! Anfangszuordnung der pointer
    129 !
    130 ! Revision 1.17  1999/03/03 07:22:35  raasch
    131 ! Beim binaeren Einlesen werden alle Variablen mit exakter Laenge von Datei
    132 ! gelesen, da sich bei schneller I/O Records nicht einfach ueberspringen lassen
    133 !
    134 ! Revision 1.14  1998/12/15 07:33:50  raasch
    135 ! Zusaetzliche Ermittelung der tatsaechlichen Einlesezeiten auf T3E
    136 ! mit messpunkt_e(14) unter Verwendung von TIMEF (s. cpu_zeitmessung)
    13718!
    13819! Revision 1.1  1998/03/09 16:22:22  raasch
  • palm/trunk/SOURCE/init_advec.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_advec.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.6  2004/04/30 11:59:31  raasch
    1214! impulse_advec renamed momentum_advec
    13 !
    14 ! Revision 1.5  2001/03/30 07:25:53  raasch
    15 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    16 !
    17 ! Revision 1.4  2001/01/22 06:57:43  raasch
    18 ! Module test_variables removed
    19 !
    20 ! Revision 1.3  2000/12/21 16:12:52  letzel
    21 ! All comments translated into English.
    22 !
    23 ! Revision 1.2  1999/02/17 09:18:45  raasch
    24 ! Fehlerkorrektur, Verfahren muessen unabhaengig voneinander initialisierbar
    25 ! sein
    2615!
    2716! Revision 1.1  1999/02/05 09:07:38  raasch
  • palm/trunk/SOURCE/init_cloud_physics.f90

    r1 r3  
    88! Former revisions:
    99! ------------------
    10 ! $Log: init_cloud_physics.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2005/06/26 19:55:58  raasch
    1214! Initialization of cloud droplet constants, gas_constant renamed r_d,
    1315! latent_heat renamed l_v
    14 !
    15 ! Revision 1.4  2001/03/30 07:26:30  raasch
    16 ! Translation of remaining German identifiers (variables, subroutines, etc.),
    17 ! surface_pressure keeps unit hPa and is not converted to Pa
    18 !
    19 ! Revision 1.3  2001/01/25 07:03:10  raasch
    20 ! Module test_variables removed
    21 !
    22 ! Revision 1.2  2001/01/22 09:16:57  schroeter
    23 ! To calculate verticle pressure profile use actual surface-temperature in
    24 ! place of potential temperature as reference value.
    2516!
    2617! Revision 1.1  2000/04/13 14:37:22  schroeter
  • palm/trunk/SOURCE/init_dvrp.f90

    r1 r3  
    44! Actual revisions:
    55! -----------------
    6 !
    76! Test output of isosurface on camera file
    87! TEST: print* statements
     
    1110! Former revisions:
    1211! -----------------
    13 ! $Log: init_dvrp.f90,v $
     12! $Id$
     13! RCS Log replace by Id keyword, revision history cleaned up
     14!
    1415! Revision 1.12  2006/02/23 12:30:22  raasch
    1516! ebene renamed section, pl.. replaced by do..,
    16 !
    17 ! Revision 1.11  2005/03/31 07:01:20  raasch
    18 ! Test version with print statements checked in
    19 !
    20 ! Revision 1.10  2003/08/01 08:36:28  raasch
    21 ! Module dvrp added to routine close_dvrp, set the maximum time the program
    22 ! can be suspended on user request, every steering variable is defined for
    23 ! every stream
    24 !
    25 ! Revision 1.9  2003/03/16 09:39:44  raasch
    26 ! Two underscores (_) are placed in front of all define-strings
    27 !
    28 ! Revision 1.8  2003/03/04 11:28:53  raasch
    29 ! Update for new version of dvrp software (e.g. using module dvrp is now
    30 ! mandatory). Initialization of dvrp steering variables included.
    31 !
    32 ! Revision 1.7  2002/12/19 15:37:57  raasch
    33 ! STOP statement replaced by call of subroutine local_stop
    34 !
    35 ! Revision 1.6  2001/08/21 09:45:54  raasch
    36 ! Superelevation implemented, additional PE can be used for handling the
    37 ! dvrp output (see init_dvrp_logging and close_dvrp)
    38 !
    39 ! Revision 1.5  2001/07/12 12:10:16  raasch
    40 ! Adjustment of camera position, output of boundingbox moved from the scene
    41 ! files to camera.dvr, destination for dvrp_output is now controlled via
    42 ! parameters given in the parameter file.
    43 !
    44 ! Revision 1.4  2001/03/30 07:31:03  raasch
    45 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    46 !
    47 ! Revision 1.3  2001/01/25 07:04:23  raasch
    48 ! Subroutines init_dvrp_logging and close_dvrp attached.
    49 !
    50 ! Revision 1.2  2000/12/28 13:11:43  raasch
    51 ! Complete revision due to new features of dvrp-software
    5217!
    5318! Revision 1.1  2000/04/27 06:24:39  raasch
  • palm/trunk/SOURCE/init_grid.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_grid.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.17  2006/08/22 14:00:05  raasch
    1214! +dz_max to limit vertical stretching,
    1315! bugfix in index array initialization for line- or point-like topography
    1416! structures
    15 !
    16 ! Revision 1.16  2006/03/03 14:28:31  letzel
    17 ! call of user_init_grid now requires the actual parameter nzb_local
    18 !
    19 ! Revision 1.15  2006/02/24 16:05:59  raasch
    20 ! Error removed: l_wall modification (multiplication with wall adjustment
    21 ! factor and ghost point exchange) must be taken out of the if-clause (see end
    22 ! of this routine)
    23 !
    24 ! Revision 1.14  2006/02/23 12:30:59  raasch
    25 ! Definition of index and wall arrays describing topography
    26 !
    27 ! Revision 1.13  2002/12/19 15:38:38  raasch
    28 ! Calculation of dl_grid removed, STOP statement replaced by call of subroutine
    29 ! local_stop
    30 !
    31 ! Revision 1.12  2001/07/20 13:06:33  raasch
    32 ! +dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg
    33 !
    34 ! Revision 1.11  2001/03/30 07:31:20  raasch
    35 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    36 !
    37 ! Revision 1.10  2001/01/22 19:35:19  letzel
    38 ! All comments translated into English.
    39 !
    40 ! Revision 1.9  2001/01/22 07:06:16  raasch
    41 ! Module test_variables removed
    42 !
    43 ! Revision 1.8  1999/02/05 09:08:32  raasch
    44 ! Erweiterung von dzw um einen Gitterpunkt nach oben fuer Upstream-Spline-
    45 ! Verfahren
    46 !
    47 ! Revision 1.7  1998/07/06 12:16:43  raasch
    48 ! + USE test_variables
    49 !
    50 ! Revision 1.6  1998/03/18 20:13:16  raasch
    51 ! dz wird unveraendert gelassen (hatte bisher den Wert der groessten,
    52 ! gestreckten) Gitterweite
    53 !
    54 ! Revision 1.5  1997/09/19 07:43:27  raasch
    55 ! +dd2zu, dl_grid, l_grid
    56 !
    57 ! Revision 1.4  1997/09/12 06:26:02  raasch
    58 ! Kehrwerte der Gitterweitenquadrate werden berechnet
    59 !
    60 ! Revision 1.3  1997/09/09 08:29:12  raasch
    61 ! Kehrwerte der Gitterweiten werden berechnet
    62 !
    63 ! Revision 1.2  1997/08/26 06:31:11  raasch
    64 ! Berechnung der Vertikalniveaus und der Gitterweiten
    6517!
    6618! Revision 1.1  1997/08/11 06:17:45  raasch
  • palm/trunk/SOURCE/init_particles.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_particles.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.24  2007/02/11 13:00:17  raasch
    1214! Bugfix: allocation of tail_mask and new_tail_id in case of restart-runs
    1315! Bugfix: __ was missing in a cpp-directive
    14 !
    15 ! Revision 1.23  2006/08/22 14:11:08  raasch
    16 ! Default settings of particle start positions changed
    17 !
    18 ! Revision 1.22  2006/08/04 14:44:58  raasch
    19 ! New particles attributes initialized (e_m, dt_sum, speed_x/y/z_sgs),
    20 ! total length of mpi_particle_type adjusted, more than one particle can be
    21 ! released per point, izuf renamed iran, set a special seed value iran_part
    22 ! to be used for the particle advection code, determination of particle group
    23 ! number removed from header to here
    24 !
    25 ! Revision 1.21  2006/03/14 12:55:24  raasch
    26 ! Determination of the number of particle groups removed
    27 !
    28 ! Revision 1.20  2006/02/23 12:34:00  raasch
    29 ! Allocation of prt_start_index and prt_count + calling subroutine
    30 ! sort_particles,
    31 ! particle source parameters psl, psr, pss, psn, psb, pst, pdx, pdy, pdz are now
    32 ! 1D arrays (1:max_number_of_particle_groups) and allow up to specify
    33 ! max_number_of_particle_groups different sources by setting
    34 ! number_of_particle_groups >= 1 even if density_ratio is not specified,
    35 ! only a fraction of the particles may have tails, therefore input from unit 90
    36 ! has changed,
    37 ! 0.5 replaced by 0.4999999999 in setting the random fluctuations,
    38 ! Defaults of psl, psr, pss, psn adjusted.
    39 ! If the number of particle groups is set by the user, density_ratio is reset to
    40 ! default.
    41 ! In case of user-specified pdx, pdy, pdz = 0.0 reset them to defaults,
    42 ! improve particle release at PE boundaries, nt_anz renamed
    43 ! current_timestep_number
    44 !
    45 ! Revision 1.19  2005/06/26 19:56:47  raasch
    46 ! particle_groups%radius is used instead of diameter, particle%radius is
    47 ! initialized, new initial informations about radius and weighting factor
    48 ! stored on array particles, total length of mpi_particle_type adjusted
    49 !
    50 ! Revision 1.18  2004/04/30 12:03:28  raasch
    51 ! MPI data type for particles modified due to modification of FORTRAN particle
    52 ! type
    53 !
    54 ! Revision 1.17  2003/10/29 08:55:35  raasch
    55 ! Module random_function_mod is used, version check of the binary file,
    56 ! modifications for new particle group feature
    57 !
    58 ! Revision 1.16  2003/03/16 09:39:51  raasch
    59 ! Two underscores (_) are placed in front of all define-strings
    60 !
    61 ! Revision 1.15  2003/03/04 11:30:48  raasch
    62 ! Displacement for integers in mpi_particle_type reduced from 16 to 8 on
    63 ! ibm, since integers seem to be only 4 byte long. Accordingly, total length
    64 ! reduced from 104 to 96.
    65 !
    66 ! Revision 1.14  2002/12/19 15:39:48  raasch
    67 ! STOP statement replaced by call of subroutine local_stop
    68 !
    69 ! Revision 1.13  2002/09/12 13:04:13  raasch
    70 ! Read density_ratio from restart file
    71 !
    72 ! Revision 1.12  2002/04/16 08:08:34  raasch
    73 ! New initial informations about speed components and start positions stored on
    74 ! array particles. Check boundary conditions and read them from restart file.
    75 !
    76 ! Revision 1.11  2001/11/12 16:06:06  raasch
    77 ! Array particle_tail_coordinates extended to contain colour informations,
    78 ! reading particle data from restart file can be prohibited
    79 !
    80 ! Revision 1.9  2001/08/21 09:49:05  raasch
    81 ! Particle tail is initialized
    82 !
    83 ! Revision 1.8  2001/07/12 12:10:58  raasch
    84 ! Particle sources (starting positions) are now controlled by parameters
    85 !
    86 ! Revision 1.7  2001/03/30 07:31:37  raasch
    87 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    88 !
    89 ! Revision 1.6  2001/01/25 07:05:36  raasch
    90 ! Writing of particle informations is optional now
    91 !
    92 ! Revision 1.5  2001/01/22 07:16:17  raasch
    93 ! Module test_particles removed
    94 !
    95 ! Revision 1.4  2001/01/02 17:30:16  raasch
    96 ! Opening of unit 80 moved to check_open. Unit 90 instead of 81 is used
    97 ! for binary input.
    98 !
    99 ! Revision 1.3  2000/12/28 13:15:58  raasch
    100 ! Complete revision due to new features of dvrp-software,
    101 ! all comments translated into English,
    102 ! code is used only optionally (cpp-directives are added)
    103 !
    104 ! Revision 1.2  2000/04/27 06:38:48  raasch
    105 ! compute position of camera with new indepent routine vtk_camera and
    106 ! store result on a separate file
    10716!
    10817! Revision 1.1  1999/11/25 16:22:38  raasch
  • palm/trunk/SOURCE/init_pegrid.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_pegrid.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.28  2006/04/26 13:23:32  raasch
    1214! lcmuk does not understand the !$ comment so a cpp-directive is required
    13 !
    14 ! Revision 1.27  2006/04/26 12:50:32  raasch
    15 ! Get the number of OpenMP threads
    16 !
    17 ! Revision 1.26  2006/02/23 12:34:57  raasch
    18 ! Calculation of nanz_sums, nanz_2dh and nanz_3d moved to init_3d_model,
    19 ! anz_y.. renamed ngp_y.., + MPI type type_x_int
    20 !
    21 ! Revision 1.25  2005/10/20 15:37:36  raasch
    22 ! error removed: the last PE in a row is not allowed to have more grid points
    23 ! than the other PES (only less)
    24 !
    25 ! Revision 1.24  2005/03/26 20:35:06  raasch
    26 ! Extensions for non-cyclic boundary conditions
    27 ! MIN function must be used instead of MAX for deciding if a switch to PE0 is
    28 ! useful (multigrid-method), maximuma_grid_level_l is defined in case that a
    29 ! pre-defined value for mg_switch_to_pe0_level is given (caused errors on NEC)
    30 !
    31 ! Revision 1.23  2004/04/30 12:27:13  raasch
    32 ! Setting 1d-decomposition as default on decalpha, indices for direct
    33 ! transpositions z --> y and y --> x added, +module poisfft_mod,
    34 ! data types for mpi_alltoallv removed (not neeeded any more),
    35 ! poisfft_hybrid is no longer enforced in case of a 1d-decomposition along x,
    36 ! an enlarged grid which will fit all transposition demands is calculated,
    37 ! impulse_advec renamed momentum_advec
    38 !
    39 ! Revision 1.22  2003/10/29 08:56:20  raasch
    40 ! Number of possible grid levels for the mulitgrid-method are calculated
    41 ! based on the total domain grid
    42 !
    43 ! Revision 1.21  2003/08/01 08:40:03  raasch
    44 ! Error removed: for a 1D-decomposition along x, the number PEs must also
    45 ! be an integral divisor of the number of gridpoints along y
    46 !
    47 ! Revision 1.20  2003/05/09 14:27:26  raasch
    48 ! Setting 1d-decomposition as default on linux clusters
    49 !
    50 ! Revision 1.19  2003/04/16 13:00:11  raasch
    51 ! Call of fft_1dd_init, +module fft_for_1d_decomp
    52 !
    53 ! Revision 1.18  2003/03/16 09:39:57  raasch
    54 ! Two underscores (_) are placed in front of all define-strings
    55 !
    56 ! Revision 1.17  2003/03/12 16:32:47  raasch
    57 ! A 1d-decomposition along x is the default on NEC hosts
    58 !
    59 ! Revision 1.16  2002/12/19 15:40:39  raasch
    60 ! In case of 1d-decomposition along x only a part of the integral divisor
    61 ! conditions is checked, on IBM hosts a 1d-decomposition along x is the
    62 ! default, a switch from Poisson-FFT-solver to the hybrid-solver is made
    63 ! in case of 1d-decomposition along x, the hybrid-solver does not force a
    64 ! 1d-decomposition any more, STOP replaced by call of local_stop
    65 !
    66 ! Revision 1.15  2002/06/11 13:10:58  raasch
    67 ! Usage of the Poisson-hybrid-solver forces a 1d-decomposition,
    68 ! initialization of hybrid-solver
    69 !
    70 ! Revision 1.14  2001/08/21  09:49:45  09:49:45  raasch (Siegfried Raasch)
    71 ! MPI_COMM_WORLD replaced by comm_palm
    72 !
    73 ! Revision 1.13  2001/07/20 13:07:26  raasch
    74 ! In case of using multigrid method, grid levels and additional MPI derived data
    75 ! types are calculated. anz_yz and type_xz are now arrays.
    76 !
    77 ! Revision 1.12  2001/03/30 07:32:01  raasch
    78 ! Calculations of nw and nw1 removed,
    79 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    80 !
    81 ! Revision 1.11  2001/01/22 08:33:48  raasch
    82 ! Module test_variables removed
    83 !
    84 ! Revision 1.9  2000/12/28 13:16:37  raasch
    85 ! definition of mpi_particle_type moved to subroutine init_particles
    86 !
    87 ! Revision 1.8  2000/04/27 06:54:03  raasch
    88 ! Processor topology can now be prescribed by the user using npex and npey,
    89 ! all comments translated into English
    90 !
    91 ! Revision 1.7  1999/11/25 16:23:36  raasch
    92 ! MPI-Datentyp fuer Partikeltransfer vereinbart
    93 !
    94 ! Revision 1.6  1998/07/16 06:50:23  raasch
    95 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
    96 !
    97 ! Revision 1.5  1998/07/06 12:17:06  raasch
    98 ! + USE test_variables
    99 !
    100 ! Revision 1.4  1998/03/18 20:13:42  raasch
    101 ! myid_char auf 5 Zeichen erweitert
    102 !
    103 ! Revision 1.3  1998/01/23 10:02:08  raasch
    104 ! Gitterpunktzahlen und Datentyp fuer 2D-Randaustausch
    105 ! anz_yz_p, type_xz_p in anz_yz und type_xz umbenannt
    106 !
    107 ! Revision 1.2  1997/08/11 06:18:28  raasch
    108 ! Gitterpunktanzahl von globalen Feldern werden berechnet
    10915!
    11016! Revision 1.1  1997/07/24 11:15:09  raasch
  • palm/trunk/SOURCE/init_pt_anomaly.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_pt_anomaly.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.7  2005/03/26 20:36:55  raasch
    1214! Arguments for non-cyclic boundary conditions added to argument list of
    1315! routine exchange_horiz
    14 !
    15 ! Revision 1.6  2001/03/30 07:32:28  raasch
    16 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    17 !
    18 ! Revision 1.5  2001/01/22 19:40:09  letzel
    19 ! All comments translated into English.
    20 !
    21 ! Revision 1.4  2001/01/22 07:19:14  raasch
    22 ! Module test_variables removed
    23 !
    24 ! Revision 1.3  1998/07/06 12:17:25  raasch
    25 ! + USE test_variables
    26 !
    27 ! Revision 1.2  1997/09/12 06:26:16  raasch
    28 ! Stoerungsradius jetzt proportional zur Gitterweite
    2916!
    3017! Revision 1.1  1997/08/29 08:58:56  raasch
  • palm/trunk/SOURCE/init_rankine.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_rankine.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2005/03/26 20:38:49  raasch
    1214! Arguments for non-cyclic boundary conditions added to argument list of
    1315! routine exchange_horiz
    14 !
    15 ! Revision 1.10  2004/01/28 15:11:27  raasch
    16 ! Setting of bt before calling pres removed
    17 !
    18 ! Revision 1.9  2001/03/30 07:32:49  raasch
    19 ! Arguments removed from pres,
    20 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    21 !
    22 ! Revision 1.8  2001/01/22 19:44:52  letzel
    23 ! All comments translated into English.
    24 !
    25 ! Revision 1.7  2001/01/22 07:20:30  raasch
    26 ! Module test_variables removed
    27 !
    28 ! Revision 1.6  1998/07/06 12:17:44  raasch
    29 ! + USE test_variables
    30 !
    31 ! Revision 1.5  1998/04/15 11:20:38  raasch
    32 ! pt, usws und vsws werden an pres uebergeben
    33 !
    34 ! Revision 1.4  1997/09/12 06:26:30  raasch
    35 ! Wirbelradius jetzt proportional zur Gitterweite
    36 !
    37 ! Revision 1.3  1997/09/09 08:29:24  raasch
    38 ! Kehrwerte der Gitterweiten implementiert
    39 !
    40 ! Revision 1.2  1997/08/26 06:31:33  raasch
    41 ! +tend in pres-Aufruf
    4216!
    4317! Revision 1.1  1997/08/11 06:18:43  raasch
  • palm/trunk/SOURCE/init_slope.f90

    r1 r3  
    11 SUBROUTINE init_slope
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: init_slope.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2006/02/23 12:35:34  raasch
    1214! nanz_2dh renamed ngp_2dh
    13 !
    14 ! Revision 1.4  2003/03/16 09:40:02  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.3  2001/03/30 07:33:59  raasch
    18 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    19 !
    20 ! Revision 1.2  2001/01/22 07:21:35  raasch
    21 ! Module test_variables removed
    2215!
    2316! Revision 1.1  2000/04/27 07:06:24  raasch
     
    3124! Remember: when a sloping surface is used, only one constant temperature
    3225!           gradient is allowed!
    33 !-------------------------------------------------------------------------------!
     26!------------------------------------------------------------------------------!
    3427
    3528    USE arrays_3d
  • palm/trunk/SOURCE/interaction_droplets_ptq.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: interaction_droplets_ptq.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.1  2005/06/26 19:57:47  raasch
    1214! Initial revision
  • palm/trunk/SOURCE/local_getenv.f90

    r1 r3  
    11 SUBROUTINE local_getenv( var, ivar, value, ivalue )
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: local_getenv.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.5  2003/05/09 14:37:07  raasch
    1214! On the MUK cluster, only PE0 is able to read environment variables.
    1315! Therefore, they have to be communicated via broadcast to the other PEs.
    14 !
    15 ! Revision 1.4  2003/03/16 09:40:22  raasch
    16 ! Two underscores (_) are placed in front of all define-strings
    17 !
    18 ! Revision 1.3  2001/01/22 07:23:55  raasch
    19 ! Comments translated into English
    20 !
    21 ! Revision 1.2  1998/07/16 06:50:50  raasch
    22 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
    2316!
    2417! Revision 1.1  1997/08/11 06:21:01  raasch
     
    2922! ------------
    3023! Getting the values of environment-variabls (for different operating-systems)
    31 !-------------------------------------------------------------------------------!
     24!------------------------------------------------------------------------------!
    3225
    3326#if defined( __lcmuk )
  • palm/trunk/SOURCE/local_stop.f90

    r1 r3  
    11 SUBROUTINE local_stop
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: local_stop.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2003/03/16 09:40:28  raasch
    1214! Two underscores (_) are placed in front of all define-strings
     
    1921! ------------
    2022! Stop program execution
    21 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    2224
    2325    USE pegrid
  • palm/trunk/SOURCE/local_system.f90

    r1 r3  
    11 SUBROUTINE local_system( command )
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: local_system.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.4  2003/03/16 09:40:33  raasch
    1214! Two underscores (_) are placed in front of all define-strings
    13 !
    14 ! Revision 1.3  2001/01/22 07:25:35  raasch
    15 ! Comments translated into English
    16 !
    17 ! Revision 1.2  1998/07/16 06:51:19  raasch
    18 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
    1915!
    2016! Revision 1.1  1997/09/03 06:27:27  raasch
     
    2521! ------------
    2622! System calls for different operating systems
    27 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    2824
    2925    CHARACTER (LEN=*) ::  command
  • palm/trunk/SOURCE/local_tremain.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: local_tremain.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.14  2006/06/02 15:20:33  raasch
    1214! Extended to TIT Sun Fire X4600 System (lctit)
    13 !
    14 ! Revision 1.13  2003/05/09 14:39:55  raasch
    15 ! Time measurement on IBM is now done using function irtc, which allows correct
    16 ! measurements for jobs running over the 24:00 timeline
    17 !
    18 ! Revision 1.12  2003/03/16 09:41:49  raasch
    19 ! Two underscores (_) are placed in front of all define-strings
    20 !
    21 ! Revision 1.11  2003/03/14 13:44:20  raasch
    22 ! Print statement for testing removed
    23 !
    24 ! Revision 1.10  2003/03/12 16:33:55  raasch
    25 ! Extended to NEC-SX6 systems
    26 !
    27 ! Revision 1.9  2002/12/19 15:49:17  raasch
    28 ! Extended to IBM-Regatta systems
    29 !
    30 ! Revision 1.8  2001/03/30 07:35:00  raasch
    31 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    32 !
    33 ! Revision 1.7  2001/01/22 08:34:53  raasch
    34 ! Module test_variables removed
    35 !
    36 ! Revision 1.5  1998/07/16 06:51:44  raasch
    37 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
    38 !
    39 ! Revision 1.4  1998/07/06 12:19:02  raasch
    40 ! + USE test_variables
    41 !
    42 ! Revision 1.3  1998/03/24 15:27:23  raasch
    43 ! Zeitbestimmung auf t3eh aufgenommen
    44 !
    45 ! Revision 1.2  1998/03/23 08:39:53  raasch
    46 ! Bestimmung der verbleibenden Zeit auf vpp ergaenzt
    4715!
    4816! Revision 1.1  1998/03/18 20:14:47  raasch
  • palm/trunk/SOURCE/local_tremain_ini.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: local_tremain_ini.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.13  2007/02/11 13:07:03  raasch
    1214! Allowed cpu limit is now read from file instead of reading the value from
    1315! environment variable (see routine parin)
    14 !
    15 ! Revision 1.12  2006/06/02 15:21:09  raasch
    16 ! Extended to TIT Sun Fire X4600 System (lctit)
    17 !
    18 ! Revision 1.11  2003/05/09 14:40:20  raasch
    19 ! Time measurement on IBM is now done using function irtc, which allows correct
    20 ! measurements for jobs running over the 24:00 timeline
    21 !
    22 ! Revision 1.10  2003/03/16 09:42:01  raasch
    23 ! Two underscores (_) are placed in front of all define-strings
    24 !
    25 ! Revision 1.9  2002/12/19 15:50:31  raasch
    26 ! Measurements extended for IBM
    27 !
    28 ! Revision 1.8  2001/03/30 07:35:16  raasch
    29 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    30 !
    31 ! Revision 1.7  2001/01/22 08:35:38  raasch
    32 ! Module test_variables removed
    33 !
    34 ! Revision 1.5  1998/07/16 06:52:11  raasch
    35 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
    36 !
    37 ! Revision 1.4  1998/07/06 12:19:21  raasch
    38 ! + USE test_variables
    39 !
    40 ! Revision 1.3  1998/03/24 15:27:38  raasch
    41 ! Initialisierung auf t3eh ergaenzt
    42 !
    43 ! Revision 1.2  1998/03/23 08:40:23  raasch
    44 ! Initialisierung auf vpp ergaenzt
    4516!
    4617! Revision 1.1  1998/03/18 20:15:05  raasch
  • palm/trunk/SOURCE/modules.f90

    r1 r3  
    99! Former revisions:
    1010! -----------------
    11 ! $Log: modules.f90,v $
     11! $Id$
     12! RCS Log replace by Id keyword, revision history cleaned up
     13!
    1214! Revision 1.95  2007/02/11 13:18:30  raasch
    1315! version 3.1b (last under RCS control)
    14 !
    15 ! Revision 1.94  2006/08/22 14:12:18  raasch
    16 ! +dz_max,
    17 ! id_var_dopts extended for second dimension
    18 !
    19 ! Revision 1.93  2006/08/04 14:50:58  raasch
    20 ! +c_0, diss, dt_min_part, end_time_prel, particles_per_point, sgs_wfu/v/w_part,
    21 ! use_sgs_for_particles, use_upstream_for_tke needed for the particle package,
    22 ! +dt_sum, e_m, speed_x/y/z/_sgs in particle_type,
    23 ! +dt_dopts, time_dopts, netcdf ids for output of particle time series,
    24 ! do2d_unit and do3d_unit now defined as 2d-arrays, izuf renamed iran
    25 !
    26 ! Revision 1.92  2006/06/02 15:21:45  raasch
    27 ! +NetCDF ids for the shifted u- and v-grid, +data_output_xy/xz/yz
    28 !
    29 ! Revision 1.91  2006/04/26 12:22:52  raasch
    30 ! +threads_per_task,
    31 ! sums_l, sums_l_l arrays changed to three dimensions
    32 !
    33 ! Revision 1.90  2006/04/11 14:58:27  raasch
    34 ! pl_spectra renamed data_output_sp, length of run_identifier increased to 80
    35 !
    36 ! Revision 1.89  2006/03/14 12:55:52  raasch
    37 ! Default for number_of_particle_groups changed to 1
    38 !
    39 ! Revision 1.88  2006/03/03 14:26:53  letzel
    40 ! control parameter topography now has length 40
    41 !
    42 ! Revision 1.87  2006/03/03 08:50:25  raasch
    43 ! Default value for surface_heatflux changed to 9999999.9
    44 !
    45 ! Revision 1.86  2006/02/23 12:39:14  raasch
    46 ! Version 3.0 with topography
    47 ! index arrays nzb_2d replaced by several new variable dependend arrays
    48 ! (e.g. nzb_u_inner), nzb_diff_2d removed, +new grid point
    49 ! counters ngp_2dh_outer and ngp_3d_inner,
    50 ! +skip time control parameters skip_time_..., +log_message,
    51 ! particle parameters psl, psr, pss, psn, psb, pst, pdx, pdy, pdz
    52 ! are now 1D arrays (1:max_number_of_particle_groups), default of
    53 ! number_of_particle_groups changed to 999999, idum in particle_type renamed
    54 ! tail_id,
    55 ! +data_output, data_output_user, do2d, do3d, doav for new steering of data
    56 ! output,
    57 ! +type_x_int for ghost exchange of integer arrays,
    58 ! +netcdf_64bit for supporting NetCDF 64-bit offset format,
    59 ! +mixing_length_1d, dissipation_1d for steering the 1d-model,
    60 ! +e_min (control parameter for minimum tke), +parameters for assuring a
    61 ! constant volume flow, +prt_count, prt_start_index, particle_advection,
    62 ! ebene renamed section, nt_anz renamed current_timestep_number, anz_y..
    63 ! renamed ngp_y.., nanz_sums, nanz_2dh, nanz_3d renamed ngp_sums, ngp_2dh,
    64 ! ngp_3d, respectively, pl1d_anz, pl2d_.._anz, plts_anz, pl3d_anz renamed
    65 ! dopr_n, pl2d_.._no, dots_n, do3d_n, respectively,
    66 ! ql_v and ql_vp now have the target attribute,
    67 ! average_period_pr1d removed
    68 !
    69 ! Revision 1.85  2005/12/06 16:40:56  raasch
    70 ! Version 2.11c
    71 !
    72 ! Revision 1.84  2005/10/20 15:38:23  raasch
    73 ! +id_var_rnop_prt, some values of prt_var_names and prt_var_units changed
    74 !
    75 ! Revision 1.83  2005/06/29 10:40:30  steinfeld
    76 ! Scalars ug and vg changed into one-dimensional arrays;
    77 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
    78 !  ug_vertical_gradient_level_ind,
    79 !  vg_surface, vg_vertical_gradient, vg_vertical_gradient_level,
    80 !  vg_vertical_gradient_level_ind
    81 !
    82 ! Revision 1.82  2005/06/26 20:03:33  raasch
    83 ! +ql_1, ql_2, ql_m, ql is now a pointer, +variables for cloud droplets (arrays
    84 ! ql_v,ql_vp, pointer ql_c), gas_constant renamed r_d, latent_heat renamed l_v,
    85 ! radius instead of diameter, radius + weighting factor included in
    86 ! particle_type, pl2d and pl3d string length increased from 6 to 10,
    87 ! default of call_psolver_at_all_substeps is TRUE
    88 !
    89 ! Revision 1.81  2005/05/18 15:36:46  raasch
    90 ! New module netcdf_control, +data_output_format, pl1d_callcount renamed
    91 ! dopr_time_count,
    92 ! default values of pl_spectra and spectra_direction changed
    93 !
    94 ! Revision 1.80  2005/04/23 09:31:14  raasch
    95 ! Implicit counters i renamed i9 due to declaration conflicts with i in other
    96 ! subroutines (reported by Portland compiler),
    97 ! crmax increased to 100, fcl_factor renamed cfl_factor
    98 !
    99 ! Revision 1.79  2005/03/26 20:40:48  raasch
    100 ! +variables for non-cyclic boundary conditions, call_psolver_at_all_substeps,
    101 ! default value for residual limit increased from 1E-6 to 1E-4
    102 !
    103 ! Revision 1.78  2004/04/30 12:30:01  raasch
    104 ! idum added in particle_type for better alignment (decalpha),
    105 ! DATA statements replaced by (/ .../) because of internal compiler errors
    106 ! on decalpha, + indices for direct transpositions z --> y,
    107 ! dt_max_1d increased to 300.0 (from 60.0), +grid_matching,
    108 ! arrays cross_ts_numbers and cross_ts_number_count enlarged,
    109 ! impulse_advec renamed momentum_advec
    110 !
    111 ! Revision 1.77  2004/01/30 10:30:26  raasch
    112 ! +index arrays nzb_2d, nzb_diff_2d
    113 !
    114 ! Revision 1.76  2004/01/28 15:25:15  raasch
    115 ! +array tsc for timestep steering, +tu_m, tv_m etc. pointer for additional
    116 ! tendencies needed for Runge-Kutta scheme, +intermediate_timestep_count(_max)
    117 ! type log changed to logpoint due to name conflict with intrinsic log
    118 !
    119 ! Revision 1.75  2003/10/29 08:58:54  raasch
    120 ! +mg_cycles, mg_switch_to_pe0, mg_switch_to_pe0_level in module
    121 ! control_parameters,
    122 ! +mg_loc_ind in module indices, variables related to particle groups added
    123 ! to module particle_attributes
    124 !
    125 ! Revision 1.74  2003/08/01 08:40:59  raasch
    126 ! +isclice_dvrp, lock_steering_update, slicer_position_dvrp,
    127 ! slicer_range_limits_dvrp in module dvrp_variables
    128 !
    129 ! Revision 1.73  2003/05/09 14:41:11  raasch
    130 ! New version number 2.7
    131 !
    132 ! Revision 1.72  2003/04/25 18:54:39  raasch
    133 ! Module mpi added in module pegrid in case of linux cluster
    134 !
    135 ! Revision 1.71  2003/04/16 13:27:33  raasch
    136 ! Arrays cross_ts_uy... enlarged
    137 !
    138 ! Revision 1.70  2003/04/16 13:07:43  raasch
    139 ! Array cross_ts_profiles enlarged
    140 !
    141 ! Revision 1.69  2003/03/16 09:42:07  raasch
    142 ! Two underscores (_) are placed in front of all define-strings
    143 !
    144 ! Revision 1.68  2003/03/14 13:44:52  raasch
    145 ! +random_generator in control_parameters
    146 !
    147 ! Revision 1.67  2003/03/12 16:34:19  raasch
    148 ! +tasks_per_node in pegrid
    149 !
    150 ! Revision 1.66  2003/03/04 11:32:50  raasch
    151 ! +type steering, +steering_dvrp in dvrp_variables
    152 !
    153 ! Revision 1.65  2002/12/19 15:52:37  raasch
    154 ! +initial_wallclock_time in cpulog, +myid_char_14 in pegrid,
    155 ! default value for termination_time_needed now set in check_parameters,
    156 ! array dl_grid removed, +wall_adjustment_factor, dt_restart, restart_time,
    157 ! time_restart
    158 !
    159 ! Revision 1.64  2002/09/12 13:05:44  raasch
    160 ! Default values of *_vertical_gradient_level changed from 100000.0 to -1.0,
    161 ! +molecular_viscosity, density_ratio
    162 !
    163 ! Revision 1.63  2002/06/11 13:12:29  raasch
    164 ! Entrys buoyancy, diffusion_e, diffusion_s, diffusion_u, diffusion_v and
    165 ! diffusion_w removed from module pointer_interfaces,
    166 ! length of character variables host + psolver increased to 64 and 16
    167 !
    168 ! Revision 1.62  2002/05/02 18:53:18  raasch
    169 ! Module singleton moved to a separate file. Arrays e_3, pt_3, q_3, u_3, v_3,
    170 ! w_3 re-introduced. Pointer e_p, pt_p, q_p, u_p, v_p, w_p re-introduced.
    171 !
    172 ! Revision 1.61  2002/04/16 08:09:32  raasch
    173 ! Information about speed components and start positions included in
    174 ! particle_type, +dt_write_particle_data, particle_advection_start,
    175 ! time_write_particle_data and parameters for particle boundary conditions,
    176 ! +vertical_particle_advection,
    177 ! + variables for scalar boundary conditions and initialization
    178 !
    179 ! Revision 1.60  2001/11/12 16:05:01  raasch
    180 ! +read_particles_from_restartfile
    181 !
    182 ! Revision 1.58  2001/09/04 12:00:52  raasch
    183 ! New version number 2.3b
    184 !
    185 ! Revision 1.57  2001/08/21 09:52:08  raasch
    186 ! +comm_palm, maximum_number_of_tailpoints, maximum_tailpoint_age,
    187 ! minimum_tailpoint_distance, particle_tail_coordinates, superelevation_x,
    188 ! superelevation_y, wall_adjustment
    189 !
    190 ! Revision 1.56  2001/07/20 13:09:07  raasch
    191 ! +grid_level, maximum_grid_level, cycle_mg, gamma_mg, ngsrb,
    192 ! residual_limit in control_parameters, +dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg
    193 ! in arrays_3d, +superelevation in dvrp_variables, +nxl_mg .. nzt_mg in
    194 ! indices, +ddx2_mg, ddy2_mg in grid_variables
    195 !
    196 ! Revision 1.55  2001/07/12 12:13:24  raasch
    197 ! default of bc_e_b changed from (u*)**2+neumann to neumann,
    198 ! +particle_maximum_age, psl-pst, pdx-pdz,
    199 ! default values defined for dvrp_directory, dvrp_file, dvrp_host, dvrp_password
    200 ! SAVE-statements added to all modules
    201 !
    202 ! Revision 1.54  2001/03/30 07:35:57  raasch
    203 ! Interface for advec_s_bc changed, interface for production_e eliminated,
    204 ! arrays work, work1, e_3, pt_3, q_3, u_3, v_3, w_3 removed,
    205 ! pointer e_p, pt_p, q_p, u_p, v_p, w_p removed,
    206 ! -nw, nw1,
    207 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    208 !
    209 ! Revision 1.53  2001/01/30 01:57:16  raasch
    210 ! +passive_scalar
    211 !
    212 ! Revision 1.51  2001/01/26 00:31:49  raasch
    213 ! +fft_method, use_surface_fluxes in steuer, +module singleton,
    214 ! +write_particle_informations in particle_attributes
    215 !
    216 ! Revision 1.49  2001/01/05 15:13:16  raasch
    217 ! new module spectrum for computing spectra, +average_count_sp,
    218 ! averaging_interval_sp, dt_dosp, time_average_sp, time_dosp in module steuer,
    219 ! +spectrum_x, spectrum_y in module statistik
    220 !
    221 ! Revision 1.48  2001/01/02 17:31:36  raasch
    222 ! openfile has the new type file_status, version updated to 2.1a
    223 !
    224 ! Revision 1.47  2000/12/30 18:48:53  raasch
    225 ! Default value of mode_dvrp changed
    226 !
    227 ! Revision 1.46  2000/12/28 13:21:35  raasch
    228 ! version updated to 2.1
    229 !
    230 ! Revision 1.45  2000/12/28 13:19:35  raasch
    231 ! +mode_dvrp, +interval_..._dvrp, REAL*4 changed to REAL in particle_attributes,
    232 ! +dt_prel, time_prel, -particle_transport
    233 !
    234 ! Revision 1.44  2000/07/03 13:02:32  raasch
    235 ! version updated to 2.0b
    236 !
    237 ! Revision 1.43  2000/07/03 12:59:53  raasch
    238 ! +dvrp_directory, dvrp_host, dvrp_password, dvrp_username,
    239 ! new module pointer_interfaces which includes interfaces for subroutines with
    240 ! pointer arguments called in prognostic_equations
    241 !
    242 ! Revision 1.42  2000/04/27 07:35:56  raasch
    243 ! npex and npey must be defined in general (not only for parallel case)
    244 !
    245 ! Revision 1.41  2000/04/27 07:07:14  raasch
    246 ! +pt_slope_offset, pt_slope_ref, npex, npey, dt_dvrp, time_dvrp,
    247 ! new module dvrp_variables
    248 ! -dt_plisos, dt_plpart, isosurface_filecount, time_plisos, time_plpart
    249 ! Version 2.0a, all comments translated into English, old revision remarks
    250 ! deleted
    251 !
    252 ! Revision 1.40  2000/04/13 13:17:25  schroeter
    253 ! Felder fuer Wolkenphysik zu Modul arrays3d hinzugefuegt, neues Modul
    254 ! (cloud_parameters) fuer Wolkenphysikparameter eingefuegt
    255 !
    256 ! Revision 1.39  2000/01/10  10:05:11  10:05:11  raasch (Siegfried Raasch)
    257 ! +use_ug_for_galilei_tr, Aufteilung der zeitabhaengigen Groessen in pointer
    258 ! und targets, + timestep_count
    259 !
    260 ! Revision 1.38  1999/12/28 10:06:52  raasch
    261 ! +Variablen fuer Plot von Isooberflaechen
    262 !
    263 ! Revision 1.37  1999/11/25 16:26:38  raasch
    264 ! Neues Modul particle_attributes,
    265 ! +Variablen fuer Partikelverlagerung: particle_transport, dt_plpart,
    266 ! Version 1.1e
    26716!
    26817! Revision 1.1  1997/07/24 11:21:26  raasch
  • palm/trunk/SOURCE/netcdf.f90

    r1 r3  
    1111! Former revisions:
    1212! -----------------
    13 ! $Log: netcdf.f90,v $
     13! $Id$
     14! RCS Log replace by Id keyword, revision history cleaned up
     15!
    1416! Revision 1.12  2006/09/26 19:35:16  raasch
    1517! Bugfix yv coordinates for yz cross sections
    16 !
    17 ! Revision 1.11  2006/08/22 14:13:11  raasch
    18 ! Extension for timeseries information of particle groups,
    19 ! xz and yz cross sections now up to nzt+1
    20 !
    21 ! Revision 1.10  2006/08/04 14:51:42  raasch
    22 ! Extensions for output of particle time series
    23 ! do2d_unit and do3d_unit now defined as 2d-arrays
    24 !
    25 ! Revision 1.9  2006/06/13 13:01:08  raasch
    26 ! Some errors of revision 1.8 removed
    27 !
    28 ! Revision 1.8  2006/06/02 15:23:13  raasch
    29 ! x- and y-coordinates of 2d- and 3d-variables are now precisely defined
    30 ! depending on which grid (u-, v- or scalar) the respective variable is defined
    31 !
    32 ! Revision 1.7  2006/04/21 14:23:47  letzel
    33 ! Error removed in extend mode when checking whether the selected cross sections
    34 ! match those in the already existing NetCDF file.
    35 !
    36 ! Revision 1.6  2006/04/11 14:58:47  raasch
    37 ! pl_spectra renamed data_output_sp
    38 !
    39 ! Revision 1.5  2006/02/23 12:42:07  raasch
    40 ! Extensions for output of 2d/3d time averaged data, definition of grids for
    41 ! user-defined output quantities,
    42 ! error removed getting the variable ids for _ext for profiles and time series,
    43 ! ebene renamed section, pl1d_anz renamed dopr_n,
    44 ! error number argument for handle_netcdf_error,
    45 ! no optimization on ibmy due to compiler error
    46 !
    47 ! Revision 1.4  2005/12/06 16:42:57  raasch
    48 ! Error in netcdf variable declaration for statistic_regions /= 0 removed
    49 !
    50 ! Revision 1.3  2005/10/20 15:40:17  raasch
    51 ! real number of particles is additionally contained as a variable on the
    52 ! particle data file,
    53 ! pr_xz added to vertical cross sections,
    54 ! error removed in calculating y axis data for yz cross sections
    55 !
    56 ! Revision 1.2  2005/06/26 20:04:31  raasch
    57 ! Particle output extended, 2d- and 3d-files are not opened if there is no
    58 ! output for these files (pl2d=pl3d=' ')
    5918!
    6019! Revision 1.1  2005/05/18 15:37:16  raasch
  • palm/trunk/SOURCE/package_parin.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: package_parin.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.18  2006/08/04 14:52:23  raasch
    1214! +dt_dopts, dt_min_part, end_time_prel, particles_per_point,
    1315! use_sgs_for_particles in particles_par
    14 !
    15 ! Revision 1.17  2006/04/11 14:59:05  raasch
    16 ! pl_spectra renamed data_output_sp
    17 !
    18 ! Revision 1.16  2006/02/23 12:43:24  raasch
    19 ! +number_of_particle_groups, skip_particles_for_tail, use_particle_tails in
    20 ! particles_par,
    21 ! +skip_time_dosp in spectra_par,
    22 ! dt_plsp renamed dt_dosp, default setting of dt_dops,
    23 ! set flag particle_advection
    24 !
    25 ! Revision 1.15  2005/06/26 20:05:42  raasch
    26 ! diameter renamed radius
    27 !
    28 ! Revision 1.14  2003/10/29 09:01:29  raasch
    29 ! +diameter in particles_par, uniform_psize renamed dvrp_psize
    30 !
    31 ! Revision 1.13  2003/08/01 08:43:33  raasch
    32 ! +slicer_range_limits_dvrp in dvrp_graphics_par
    33 !
    34 ! Revision 1.12  2003/03/16 09:42:11  raasch
    35 ! Two underscores (_) are placed in front of all define-strings
    36 !
    37 ! Revision 1.11  2002/09/12 13:06:50  raasch
    38 ! +density_ratio in particles_par
    39 !
    40 ! Revision 1.10  2002/04/16 08:10:08  raasch
    41 ! +variables for particle boundary conditions, +vertical_particle_advection
    42 ! in particles_par
    43 !
    44 ! Revision 1.9  2001/11/12 16:07:04  raasch
    45 ! +read_particles_from_restartfile
    46 !
    47 ! Revision 1.7  2001/08/21 09:52:55  raasch
    48 ! +superelevation, superelevation_x, superelevation_y in dvrp_graphics_par,
    49 ! +maximum_number_of_particles, maximum_number_of_tailpoints,
    50 ! maximum_tailpoint_age, minimum_tailpoint_distance in particles_par
    51 !
    52 ! Revision 1.6  2001/07/12 12:42:49  raasch
    53 ! +several new parameters in dvrp_graphics_par and particles_par
    54 !
    55 ! Revision 1.5  2001/03/30 07:37:56  raasch
    56 ! Namelist names changed,
    57 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    58 !
    59 ! Revision 1.4  2001/01/25 07:10:48  raasch
    60 ! +write_particle_informations
    61 !
    62 ! Revision 1.2  2001/01/05 15:13:34  raasch
    63 ! Spectra package added
    6416!
    6517! Revision 1.1  2000/12/28 13:21:57  raasch
  • palm/trunk/SOURCE/palm.f90

    r1 r3  
    1212! Former revisions:
    1313! -----------------
    14 ! $Log: palm.f90,v $
     14! $Id$
     15! RCS Log replace by Id keyword, revision history cleaned up
     16!
    1517! Revision 1.10  2006/08/04 14:53:12  raasch
    1618! Distibution of run description header removed, call of header moved behind
    1719! init_3d_model
    18 !
    19 ! Revision 1.9  2006/02/23 12:44:43  raasch
    20 ! plot_.. renamed data_output_..
    21 !
    22 ! Revision 1.8  2004/04/30 12:31:07  raasch
    23 ! Missing cpp-instructions added, leap_frog changed to time_integration
    24 !
    25 ! Revision 1.7  2003/05/09 14:46:49  raasch
    26 ! String converted to integer before send via broadcast due to Linux MPICH
    27 ! problem
    28 !
    29 ! Revision 1.6  2003/03/16 09:42:17  raasch
    30 ! Two underscores (_) are placed in front of all define-strings
    31 !
    32 ! Revision 1.5  2002/12/19 15:53:59  raasch
    33 ! Routine local_tremain_ini moved from init_3d_model to here
    34 !
    35 ! Revision 1.4  2001/08/21 09:53:35  raasch
    36 ! comm_palm is now used as the global communicator since one PE maybe split from
    37 ! MPI_COMM_WORLD for usage in dvrp-graphics
    38 !
    39 ! Revision 1.3  2001/03/30 07:39:41  raasch
    40 ! Translation of remaining German identifiers (variables, subroutines, etc.),
    41 ! closing of unit 14 moved from routine write_3d_binary to this routine
    4220!
    4321! Revision 1.2  2001/01/25 07:15:06  raasch
     
    4523! Initialization of dvrp logging as well as exit of dvrp moved to new
    4624! subroutines init_dvrp_logging and close_dvrp (file init_dvrp.f90)
    47 !
    48 ! Revision 1.23  2001/01/05 15:13:51  raasch
    49 ! +module spectrum
    50 !
    51 ! Revision 1.22  2001/01/02 17:32:55  raasch
    52 ! Subroutine close_files renamed to close_file, now called with an argument
    53 !
    54 ! Revision 1.21  2000/12/28 13:32:16  raasch
    55 ! Call of MPI_FINALIZE moved to the end of the program. At the beginning
    56 ! comm2d is set to MPI_COMM_WORLD (it was undefined before).
    57 ! Call of write_particles at the end. Open and close of DVRP-Software at the
    58 ! beginning and at the end (optionally, by cpp-directives)
    59 !
    60 ! Revision 1.20  2000/04/27 06:42:43  raasch
    61 ! closing call to dvrp-software (if used) at the end of the run,
    62 ! old revision remarks deleted
    63 !
    64 ! Revision 1.19  2000/01/25 15:29:23  letzel
    65 ! All comments translated into English
    66 !
    67 ! Revision 1.18  1999/11/25 16:27:49  raasch
    68 ! Bei Benutzung von VTK-Grafik wird Modell als UP in C++ - Programm eingebunden
    6925!
    7026! Revision 1.1  1997/07/24 11:23:35  raasch
  • palm/trunk/SOURCE/parin.f90

    r2 r3  
    44! Actual revisions:
    55! -----------------
    6 ! +netcdf_64bit_3d in d3par, revision history cleaned up
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11! +netcdf_64bit_3d in d3par,
     12! RCS Log replace by Id keyword, revision history cleaned up
    1113!
    1214! Revision 1.57  2007/02/11 13:11:22  raasch
  • palm/trunk/SOURCE/poisfft.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: poisfft.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.24  2006/08/04 15:00:24  raasch
    1214! Default setting of the thread number tn in case of not using OpenMP
     
    1517! Additional compiler directive in routine tridia_1dd for preventing loop
    1618! exchange on NEC-SX6
    17 !
    18 ! Revision 1.22  2005/05/18 15:53:00  raasch
    19 ! Missing argument in ffty
    20 !
    21 ! Revision 1.21  2005/03/26 20:51:59  raasch
    22 ! Use of module indices removed from routine split_1dd since it caused errors
    23 ! in case of nx /= ny
    2419!
    2520! Revision 1.20  2004/04/30 12:38:09  raasch
     
    3025! is now a 1d-array, variables passed by modules instead of using parameter
    3126! lists, enlarged transposition arrays introduced
    32 !
    33 ! Revision 1.19  2003/04/16 13:11:16  raasch
    34 ! Internal fft method calls substituted by calls of fft_x_1dd and fft_y_1dd
    35 !
    36 ! Revision 1.18  2003/03/16 09:44:45  raasch
    37 ! Two underscores (_) are placed in front of all define-strings
    38 !
    39 ! Revision 1.17  2003/03/12 16:36:34  raasch
    40 ! System specific version for NEC included (equivalent to CRAY scilib
    41 ! calls, but different subroutine names; the CRAY names are replaced with
    42 ! NEC names by cpp), size of array work1 is expanded because NEC routines
    43 ! need more memory
    44 !
    45 ! Revision 1.15  2002/06/11 13:14:31  raasch
    46 ! System specific version for IBM included (parallel case only)
    47 !
    48 ! Revision 1.13  2001/01/25 22:19:18  letzel
    49 ! All comments translated into English.
    50 !
    51 ! Revision 1.12  2001/01/23 01:53:04  raasch
    52 ! Total revision of parallel and non-paralle fft-routines. Additional
    53 ! implementation of the singleton algorithm.
    54 !
    55 ! Revision 1.9  1998/07/27 08:12:56  schroeter
    56 ! Berechnung von l in maketri geschieht nun ausserhalb
    57 ! der k-Schleife, da unabhaengig von k, dadurch
    58 ! Rechenzeitersparnis
    59 !
    60 ! Revision 1.8  1998/07/24 07:52:09  schroeter
    61 ! Anpassung des Druckloesealgorithmusses an Vorgehensweise
    62 ! von Schmidt, Schumann und Volkert
    63 ! (DFVLR-Mitteilung 84-15, S. 20-26)
    6427!
    6528! Revision 1.1  1997/07/24 11:24:14  raasch
  • palm/trunk/SOURCE/poisfft_hybrid.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: poisfft_hybrid.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2004/04/30 12:43:14  raasch
    1214! Renaming of fft routines, additional argument in calls of fft_y_m
    13 !
    14 ! Revision 1.10  2003/08/01 08:49:30  raasch
    15 ! Error removed: reqa=req after calling routine MPI_IRECV (cascade)
    16 !
    17 ! Revision 1.9  2003/05/09 14:50:47  raasch
    18 ! Parameters in calls of MPI_ALLTOALL, MPI_WAITALL changed in order to fulfill
    19 ! f95 type and rank requirements
    20 !
    21 ! Revision 1.8  2003/03/16 11:17:04  raasch
    22 ! +module pegrid in poisfft_hybrid_ini
    23 !
    24 ! Revision 1.7  2003/03/16 09:45:52  raasch
    25 ! Two underscores (_) are placed in front of all define-strings
    26 !
    27 ! Revision 1.6  2003/03/16 09:11:42  raasch
    28 ! Fortran error removed
    29 !
    30 ! Revision 1.5  2003/03/14 13:45:53  raasch
    31 ! Loop optimization in maketri_hybrid
    32 !
    33 ! Revision 1.4  2003/03/12 16:37:22  raasch
    34 ! Vectorized version (poisfft_hybrid_omp_vec) added, arguments ar and tri
    35 ! passed to routine substi_hybrid (otherwise no vectorization and very
    36 ! bad performance on NEC), reading of environment variable tasks_per_node
    37 ! moved to routine parin
    38 !
    39 ! Revision 1.3  2002/12/20 15:18:02  raasch
    40 ! Variable statement restricted for parallel environment in subroutine
    41 ! cascade
    4215!
    4316! Revision 1.2  2002/12/19 16:08:31  raasch
  • palm/trunk/SOURCE/poismg.f90

    r1 r3  
    1212! Former revisions:
    1313! -----------------
    14 ! $Log: poismg.f90,v $
     14! $Id$
     15! RCS Log replace by Id keyword, revision history cleaned up
     16!
    1517! Revision 1.6  2005/03/26 20:55:54  raasch
    1618! Implementation of non-cyclic (Neumann) horizontal boundary conditions,
    1719! routine prolong simplified (one call of exchange_horiz spared)
    18 !
    19 ! Revision 1.5  2003/10/29 09:04:34  raasch
    20 ! On a defined level, data are gathered on PE0 and further calculations are
    21 ! carried out only on this PE.
    22 !
    23 ! Revision 1.4  2003/03/16 09:45:57  raasch
    24 ! Two underscores (_) are placed in front of all define-strings
    25 !
    26 ! Revision 1.3  2001/07/23 12:00:31  raasch
    27 ! Neumann boundary conditions for the upper boundary and Diriclet conditions
    28 ! for the lower boundary implemented
    29 !
    30 ! Revision 1.2  2001/07/20 15:21:55  raasch
    31 ! Praeprocessor directives included for allreduce
    3220!
    3321! Revision 1.1  2001/07/20 13:10:51  raasch
  • palm/trunk/SOURCE/prandtl_fluxes.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: prandtl_fluxes.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.19  2006/04/26 12:24:35  raasch
    1214! +OpenMP directives and optimization (array assignments replaced by DO loops)
    13 !
    14 ! Revision 1.18  2006/02/23 12:49:32  raasch
    15 ! shf, ts, rif, us, usws, vsws, qs and qsws are now defined at the actual
    16 ! surface of the model domain, i.e. either at the bottom or at the height
    17 ! of the topography. Thus, zu(1) is replaced
    18 ! by zp = ( zu(nzb_[variable]_inner(j,i)+1 - zw(nzb_[variable]_inner(j,i) ),
    19 ! [variable](1,j,i) by [variable](nzb_[variable]_inner(j,i)+1,j,i) and
    20 ! [variable](0,j,i) by [variable](nzb_[variable]_inner(j,i)  ,j,i).
    21 !
    22 ! Revision 1.17  2004/04/30 12:43:50  raasch
    23 ! rif_m replaced by rif (they are the same when used here)
    24 !
    25 ! Revision 1.16  2003/11/20 15:13:26  raasch
    26 ! Variable name (B) changed from capital to small
    27 !
    28 ! Revision 1.15  2001/03/30 07:46:47  raasch
    29 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    30 !
    31 ! Revision 1.14  2001/01/29 12:33:38  raasch
    32 ! Passive scalar is considered
    33 !
    34 ! Revision 1.13  2001/01/25 07:23:22  raasch
    35 ! Range of ts is limited, since otherwise overflow may occur on REAL*4
    36 ! machines
    37 !
    38 ! Revision 1.11  2001/01/22 07:50:34  raasch
    39 ! Module test_variables removed
    40 !
    41 ! Revision 1.10  2000/04/13 13:44:41  schroeter
    42 ! + implentation of the Prandtl layer for the total water content
    43 !
    44 ! Revision 1.9  2000/01/26  14:48:48  14:48:48  letzel (Marcus Letzel)
    45 ! correct wrong time level of rif in computation of theta*,
    46 ! all comments translated into English
    47 !
    48 ! Revision 1.8  1998/09/22 17:27:15  raasch
    49 ! Testweise Randbedingung fuer TKE mit cm = 0.4, aber vorerst auskommentiert
    50 !
    51 ! Revision 1.7  1998/08/05 06:54:23  raasch
    52 ! Begrenzung von rif ist jetzt variabel (rif_max, rif_min)
    53 !
    54 ! Revision 1.6  1998/07/06 12:30:01  raasch
    55 ! + USE test_variables
    56 !
    57 ! Revision 1.5  1998/04/15 11:22:58  raasch
    58 ! Berechnung von theta* ueber Oberflaechentemperatur moeglich
    59 !
    60 ! Revision 1.4  1998/03/11 11:53:22  raasch
    61 ! Zusaetzliche untere Randbedingung fuer TKE ( (u*/0.1)**2 )
    62 !
    63 ! Revision 1.3  1998/02/10 15:08:53  raasch
    64 ! Grenzfall bei labiler Schichtung ( a=1, b=1 ) wird stabil gerechnet
    65 ! Begrenzung von rif aktiviert
    66 !
    67 ! Revision 1.2  1998/02/04 16:09:29  raasch
    68 ! Berechnung der Waermefluesse
    6915!
    7016! Revision 1.1  1998/01/23 10:06:06  raasch
  • palm/trunk/SOURCE/pres.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: pres.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.25  2006/04/26 13:26:12  raasch
    1214! OpenMP optimization (+localsum, threadsum)
    13 !
    14 ! Revision 1.24  2006/02/23 12:50:08  raasch
    15 ! nzb replaced by nzb_.._inner, optional volume flow control, calculation of
    16 ! divergence sum reordered due to optimization
    17 !
    18 ! Revision 1.23  2005/05/18 15:53:44  raasch
    19 ! call conditions for poisfft_hybrid modified
    20 !
    21 ! Revision 1.22  2005/03/26 20:57:24  raasch
    22 ! Adjustments for non-cyclic boundary conditions: arguments added to routine
    23 ! exchange_horiz, mean horizontal wind parallel to the outflow is calculated.
    24 !
    25 ! Revision 1.21  2004/04/30 12:45:00  raasch
    26 ! Most of the arguments removed from poisfft call, +module poisfft_mod,
    27 ! poisfft is also used for 1d-decompositions, more memory has to be
    28 ! allocated for tend, when an enlarged domain is used for transposition
    29 !
    30 ! Revision 1.20  2003/03/16 09:46:02  raasch
    31 ! Two underscores (_) are placed in front of all define-strings
    32 !
    33 ! Revision 1.19  2003/03/14 13:46:32  raasch
    34 ! Optimization of loops (IBM-version and vectorizable version)
    35 !
    36 ! Revision 1.18  2003/03/12 16:38:20  raasch
    37 ! Simulated_time has been replaced by sums(nzb+1,4) in first if clause
    38 ! in order to avoid run time errors due to compiler bugs on ibm and nec
    39 ! machines, due to optimization, dividing d by dt_3d is now done in the
    40 ! first loop (not in a seperated loop),
    41 ! error removed: tend=p added after calling sor method
    42 !
    43 ! Revision 1.17  2002/06/11 13:17:03  raasch
    44 ! Hybrid solver included, array operations on d removed due to bad performance
    45 ! on IBM, loops including the array d are combined, OpenMP directives added.
    46 ! Total perturbation pressure is no more obtained iteratively (p=p+tend) from
    47 ! the current perturbation pressure.
    48 !
    49 ! Revision 1.16  2001/09/04  12:05:52  12:05:52  raasch (Siegfried Raasch)
    50 ! Value of perturbation pressure is no more obtained iteratively (p = p + tend,
    51 ! FFT-solver only)
    52 !
    53 ! Revision 1.15  2001/07/20 13:12:54  raasch
    54 ! Multigrid method included
    55 !
    56 ! Revision 1.14  2001/03/30 07:47:29  raasch
    57 ! All arguments removed, dp replaced by tend, array work removed from
    58 ! sor argument list,
    59 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    60 !
    61 ! Revision 1.13  2001/01/22 07:53:53  raasch
    62 ! Module test_variables removed
    63 !
    64 ! Revision 1.12  2000/01/26 13:23:32  letzel
    65 ! All comments translated into English
    66 !
    67 ! Revision 1.11  1998/09/22 17:28:14  raasch
    68 ! dp war bei SOR-Aufrufen bisher nicht initialisiert, entsprechend geaendert
    69 !
    70 ! Revision 1.10  1998/07/06 12:30:22  raasch
    71 ! + USE test_variables
    72 !
    73 ! Revision 1.9  1998/04/15 11:23:20  raasch
    74 ! Testweise Zusatzbedingung fuer den Druck am unteren Rand bei inhomogener
    75 ! Oberflaechentemperatur (pt wird uebergeben)
    76 !
    77 ! Revision 1.8  1998/03/30 11:37:59  raasch
    78 ! Divergenzsummenberechnung (PEs) nach flow_statistics verlagert
    79 !
    80 ! Revision 1.7  1998/03/25 13:56:16  raasch
    81 ! dt in dt_3d umbenannt
    82 !
    83 ! Revision 1.6  1997/09/16 06:38:37  raasch
    84 ! Ueberfluessige Initialisierung von d entfernt
    85 !
    86 ! Revision 1.5  1997/09/12 06:29:56  raasch
    87 ! Randbedingungen umbenannt
    88 !
    89 ! Revision 1.4  1997/09/09 08:30:27  raasch
    90 ! Kehrwerte der Gitterweiten implementiert
    91 !
    92 ! Revision 1.3  1997/08/26 06:34:16  raasch
    93 ! Gesamtstoerdruck p wird in Bew.gleichungen mitgerechnet und ergibt sich
    94 ! durch Summation des aktuellen Stoerdrucks dp, mit dem hier korrigiert wird
    95 !
    96 ! Revision 1.2  1997/08/11 06:25:24  raasch
    97 ! Felder werden ueber Formalparameterliste uebergeben.
    9815!
    9916! Revision 1.1  1997/07/24 11:24:44  raasch
  • palm/trunk/SOURCE/print_1d.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: print_1d.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.11  2006/02/23 12:50:43  raasch
    1214! Preliminary no output of time-averaged data
    13 !
    14 ! Revision 1.10  2001/03/30 07:48:24  raasch
    15 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    16 !
    17 ! Revision 1.9  2001/01/22 20:51:40  letzel
    18 ! All comments translated into English.
    19 !
    20 ! Revision 1.8  2001/01/22 07:55:14  raasch
    21 ! Module test_variables removed
    22 !
    23 ! Revision 1.7  1998/07/06 12:31:17  raasch
    24 ! Mittelungszeitraum wird direkt durch average_period_pr1d angegeben,
    25 ! + USE test_variables
    26 !
    27 ! Revision 1.6  1998/03/30 11:38:16  raasch
    28 ! Ausgaben fuer die (Teil-) Gebiete
    29 !
    30 ! Revision 1.5  1998/03/25 13:56:31  raasch
    31 ! Erweiterungen fuer Gebietsstatistiken
    32 !
    33 ! Revision 1.4  1998/03/18 20:18:33  raasch
    34 ! Zusaetzliche Header werden ausgegeben
    35 !
    36 ! Revision 1.3  1998/03/03 08:01:28  raasch
    37 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier
    38 ! Ausgabe von zeitlich gemittelten Profilen
    39 !
    40 ! Revision 1.2  1998/02/10 15:09:16  raasch
    41 ! Ausdruck von Werten im w-Niveau
    4215!
    4316! Revision 1.1  1997/09/19 07:45:22  raasch
  • palm/trunk/SOURCE/production_e.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: production_e.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.21  2006/04/26 12:45:35  raasch
    1214! OpenMP parallelization of production_e_init
    13 !
    14 ! Revision 1.20  2006/02/23 12:51:25  raasch
    15 ! nzb_2d and nzb_diff_2d replaced by nzb_s_outer and nzb_diff_s_outer
    16 ! respectively, momentum flux calculation at vertical walls using profile
    17 ! functions, loop rearrangement for better vectorization
    18 !
    19 ! Revision 1.19  2005/06/29 08:15:33  steinfeld
    20 ! Error in calculating u_0 and v_0 removed (multiplication instead of
    21 ! division by 2dz)
    22 !
    23 ! Revision 1.18  2004/04/30 12:45:54  raasch
    24 ! dptdz eliminated
    25 !
    26 ! Revision 1.17  2004/01/30 10:36:00  raasch
    27 ! Velocity gradients at the surface limited (see u_0, v_0),
    28 ! scalar lower k index nzb_diff replaced by 2d-array nzb_diff_2d
    29 !
    30 ! Revision 1.16  2003/03/12 16:40:22  raasch
    31 ! Full code replaced in the call for all gridpoints instead of calling the
    32 ! _ij version (required by NEC, because otherwise no vectorization)
    33 !
    34 ! Revision 1.15  2002/12/19 16:16:34  raasch
    35 ! Calculation of deformation tensor re-designed (most of the finite
    36 ! differences are now formed over two grid spacings, errors in derivatives
    37 ! along y removed)
    38 !
    39 ! Revision 1.14  2002/09/12 13:09:44  raasch
    40 ! Error in calculating v_0 removed
    41 !
    42 ! Revision 1.13  2002/06/11 13:19:04  raasch
    43 ! Former subroutine changed to a module which allows to be called for all grid
    44 ! points of a single vertical column with index i,j or for all grid points by
    45 ! using function overloading.
    46 ! Calculation of u_0 and v_0 moved to the new subroutine production_e_init
    47 ! due to the global communication necessary. These arrays are allocated only
    48 ! once during the first call.
    49 !
    50 ! Revision 1.12  2001/08/21 09:59:47  raasch
    51 ! Special treatment at k=1 generally if surface fluxes are prescribed (not only
    52 ! in case of a Prandtl layer). In these cases, the boundary conditions are
    53 ! also included in the shear production term.
    54 !
    55 ! Revision 1.11  2001/03/30 07:48:54  raasch
    56 ! Calculation of shear production and buoyancy production is split into
    57 ! two loops, arguments k1 and k2 are eliminated and now used as scalars,
    58 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    59 !
    60 ! Revision 1.10  2001/01/22 07:56:19  raasch
    61 ! Module test_variables removed
    62 !
    63 ! Revision 1.9  2001/01/02 17:34:12  raasch
    64 ! -dpt_dz_d, dpt_dz_u
    65 !
    66 ! Revision 1.8  2000/07/03 13:00:24  raasch
    67 ! Arguments k1 and k2 declared as pointers,
    68 ! all comments tranlated into English
    69 !
    70 ! Revision 1.7  2000/04/18 08:13:00  schroeter
    71 ! Revision 1.5 rueckgaengig gemacht, Temperaturgrandient im
    72 ! Auftriebtrem wird wieder durch zentrale Differenzen gebildet
    73 !
    74 ! Revision 1.6  2000/04/13 14:24:02  schroeter
    75 ! condsidering the influence of humidity to TKE-production
    76 !
    77 ! Revision 1.5  99/02/17  09:29:18  09:29:18  raasch (Siegfried Raasch)
    78 ! Vertikalscherung des Windes exakter formuliert
    79 ! Temperaturgradient im Auftriebsterm enger gefasst, im stabilen Fall wird
    80 ! jetzt das Minimum vom oberen und unteren Differenzenquotienten verwendet
    81 !
    82 ! Revision 1.4  1998/07/06 12:31:52  raasch
    83 ! + USE test_variables
    84 !
    85 ! Revision 1.3  1998/04/15 11:23:49  raasch
    86 ! Im Fall einer Prandtl-Schicht wird die Energieproduktion bei nzb+1
    87 ! jetzt immer mittels shf und sonst ueber dpt/dz berechnet (bisher nur bei
    88 ! vorgegebenem Waermestrom)
    89 !
    90 ! Revision 1.2  1998/02/19 07:10:51  raasch
    91 ! vorgegebener Waermestrom wird gegebenenfalls im Auftriebsterm bei
    92 ! k=nzb+1 direkt angegeben
    9315!
    9416! Revision 1.1  1997/09/19 07:45:35  raasch
  • palm/trunk/SOURCE/prognostic_equations.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: prognostic_equations.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.21  2006/08/04 15:01:07  raasch
    1214! upstream scheme can be forced to be used for tke (use_upstream_for_tke)
    1315! regardless of the timestep scheme used for the other quantities,
    1416! new argument diss in call of diffusion_e
    15 !
    16 ! Revision 1.20  2006/02/23 12:52:08  raasch
    17 ! nzb_2d replaced by nzb_u/v/w/s_inner, +z0 in argument list of diffusion_u/v/w
    18 !
    19 ! Revision 1.19  2005/06/29 10:33:41  steinfeld
    20 ! Scalars ug and vg are replaced by homonymous arrays in order to allow for
    21 ! the consideration of a potential dependency of the geostrophic wind
    22 ! components on height
    23 !
    24 ! Revision 1.18  2005/03/26 20:58:45  raasch
    25 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,
    26 ! additional arguments km_damp_x/y for diffusion_u/v/w
    27 !
    28 ! Revision 1.17  2004/04/30 12:46:40  raasch
    29 ! impulse_advec renamed momentum_advec
    30 !
    31 ! Revision 1.16  2004/01/30 10:36:57  raasch
    32 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
    33 !
    34 ! Revision 1.15  2004/01/28 15:24:35  raasch
    35 ! Runge-Kutta schemes available, steering variables at and bt in equations
    36 ! replaced by array sct
    37 !
    38 ! Revision 1.14  2003/03/12 16:40:48  raasch
    39 ! New routine prognostic_equations_vec, optimized for vector processors
    40 !
    41 ! Revision 1.13  2002/06/11 13:20:21  raasch
    42 ! Single node optimization: loops i and j extracted from all those
    43 ! tendency-subroutines which do not contain communication. Some communication
    44 ! parts moved before the i,j loops. Call of particle advection moved to
    45 ! subroutine leap_frog.
    46 ! Former subroutine changed to a module which contains two versions: one
    47 ! with loop optimization for the single equations and one version with one
    48 ! big optimized loop containing all prognostic equations.
    49 !
    50 ! Revision 1.12  2002/05/02 18:54:12  raasch
    51 ! Timelevel t+dt re-introduced, Asselin filter and exchange of ghost points
    52 ! moved to routine leap_frog
    53 !
    54 ! Revision 1.11  2002/04/16  08:12:40  08:12:40  raasch (Siegfried Raasch)
    55 ! Particle advection will be performed only after the simulated time has
    56 ! exceeded a user defined limit
    57 !
    58 ! Revision 1.10  2001/09/04 12:10:28  raasch
    59 ! Exchange of ghost points added for the time filtered arrays
    60 !
    61 ! Revision 1.9  2001/08/21 09:59:03  raasch
    62 ! Calls to user-interface for tendency terms added
    63 !
    64 ! Revision 1.8  2001/03/30 07:50:26  raasch
    65 ! Timelevel t+dt eliminated, Asselin filter included in prognostic equations,
    66 ! calling arguments eliminated or cut down from advec_s_bc, production_e,
    67 ! advec_*_ups,
    68 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    69 !
    70 ! Revision 1.7  2001/01/29 12:34:20  raasch
    71 ! Passive scalar is considered
    72 !
    73 ! Revision 1.6  2001/01/22 07:57:20  raasch
    74 ! Module test_variables removed
    75 !
    76 ! Revision 1.5  2000/12/28 13:38:00  raasch
    77 ! Advec_particles is called unconditionally
    78 !
    79 ! Revision 1.4  2000/07/03 13:01:02  raasch
    80 ! module pointer_interfaces added to reduce cpu-time and memory demands
    81 ! of the compilation process, actual argument "work" eliminated from
    82 ! parameter list of diffusion_e
    83 !
    84 ! Revision 1.3  2000/04/27 07:09:57  raasch
    85 ! Temperature offset is added at cyclic boundaries (x-direction) when using
    86 ! a sloping surface
    87 !
    88 ! Revision 1.2  2000/04/13 15:04:28  schroeter
    89 ! prognostic equation for the total water content added, new routines:
    90 ! calc_precipitation, impact_of_latent_heat, calc_radiation,
    91 ! diffusion_pt renamed to diffusion_s, changes in the parameter list of
    92 ! subroutines diffusion_e, buyoancy, production_e
    9317!
    9418! Revision 1.1  2000/04/13 14:56:27  schroeter
  • palm/trunk/SOURCE/random_function.f90

    r1 r3  
    11 MODULE random_function_mod
    2 !-------------------------------------------------------------------------------!
     2!------------------------------------------------------------------------------!
    33! Actual revisions:
    44! -----------------
     
    77! Former revisions:
    88! -----------------
    9 ! $Log: random_function.f90,v $
     9! $Id$
     10! RCS Log replace by Id keyword, revision history cleaned up
     11!
    1012! Revision 1.3  2003/10/29 09:06:57  raasch
    1113! Former function changed to a module.
    12 !
    13 ! Revision 1.2  2001/01/22 08:00:23  raasch
    14 ! Comments translated into English
    1514!
    1615! Revision 1.1  1998/02/04 16:09:45  raasch
     
    2221! Random number generator, produces numbers equally distributed in interval [0,1]
    2322! This routine is taken from the "numerical recipies"
    24 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    2524
    2625    IMPLICIT NONE
  • palm/trunk/SOURCE/random_gauss.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: random_gauss.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.4  2006/08/04 15:01:48  raasch
    1214! Range of random number is limited by an upper limit (new second parameter)
    13 !
    14 ! Revision 1.3  2003/10/29 09:08:04  raasch
    15 ! Module random_function_mod is used
    16 !
    17 ! Revision 1.2  2001/01/22 08:03:31  raasch
    18 ! Comments translated into English
    1915!
    2016! Revision 1.1  1998/03/25 20:09:47  raasch
  • palm/trunk/SOURCE/read_3d_binary.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: read_3d_binary.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.4  2006/08/04 15:02:32  raasch
    1214! +iran, iran_part
    13 !
    14 ! Revision 1.3  2006/02/23 12:52:45  raasch
    15 ! +volume_flow_area, volume_flow_initial, -hom_sum_pr,
    16 ! reading arrays used for time-averaging
    17 !
    18 ! Revision 1.2  2005/05/18 15:54:33  raasch
    19 ! pl1d_callcount renamed dopr_time_count
    2015!
    2116! Revision 1.1  2004/04/30 12:47:27  raasch
  • palm/trunk/SOURCE/read_var_list.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: read_var_list.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.34  2006/08/22 14:14:27  raasch
    1214! +dz_max
    1315!
    14 ! Revision 1.33  2006/08/04 15:04:08  raasch
    15 ! +first_call_advec_particles, time_dopts, use_upstream_for_tke
    16 !
    17 ! Revision 1.32  2006/02/23 12:53:31  raasch
    18 ! +building parameters, conserve_volume_flow, dissipation_1d, e_min,
    19 ! mixing_length_1d, topography,
    20 ! -average_count_pl, time_average_pr,
    21 ! nt_anz renamed current_timestep_number, ..pl.. renamed ..do..
    22 !
    23 ! Revision 1.31  2005/06/29 09:12:49  steinfeld
    24 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
    25 ! ug_vertical_gradient_level_ind, vg_surface, vg_vertical_gradient,
    26 ! vg_vertical_gradient_level, vg_vertical_gradient_level_ind;
    27 ! allocation of the new arrays ug and vg
    28 !
    29 ! Revision 1.30  2005/06/26 20:09:42  raasch
    30 ! +cloud_droplets
    31 !
    32 ! Revision 1.29  2005/05/18 15:55:12  raasch
    33 ! +netcdf_precision
    34 !
    35 ! Revision 1.28  2005/03/26 21:01:00  raasch
    36 ! +bc_lr, bc_ns, inflow_disturbance_begin, inflow_disturbance_end, km_damp_max,
    37 ! outflow_damping_width
    38 !
    39 ! Revision 1.27  2004/04/30 12:49:56  raasch
    40 ! +grid_matching, impulse_advec renamed momentum_advec
    41 !
    42 ! Revision 1.26  2004/01/28 15:27:33  raasch
    43 ! at and bt replaced by array tsc
    44 !
    45 ! Revision 1.25  2003/10/29 09:10:41  raasch
    46 ! +last_dt_change
    47 !
    48 ! Revision 1.24  2003/03/14 13:47:15  raasch
    49 ! +random_generator
    50 !
    51 ! Revision 1.23  2002/12/19 16:18:27  raasch
    52 ! Array hom, which was formerly read in init_3d_model, will now be allocated
    53 ! and read here, binary version incremented to 2.2, +time_restart,
    54 ! STOP statement replaced by call of subroutine local_stop
    55 !
    56 ! Revision 1.22  2002/04/16 08:21:06  raasch
    57 ! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,
    58 ! s_vertical_gradient, s_vertical_gradient_level
    59 !
    60 ! Revision 1.21  2001/08/21 10:01:23  raasch
    61 ! +wall_adjustment
    62 !
    63 ! Revision 1.20  2001/03/30 07:51:25  raasch
    64 ! Translation of remaining German identifiers (variables, subroutines, etc.),
    65 ! runtime and package parameters are not read from unit 13 any more
    66 !
    67 ! Revision 1.19  2001/01/29 12:35:03  raasch
    68 ! +passive_scalar
    69 !
    70 ! Revision 1.18  2001/01/25 07:24:28  raasch
    71 ! +fft_method, use_surface_fluxes
    72 !
    73 ! Revision 1.17  2001/01/05 15:16:18  raasch
    74 ! +average_count_sp, time_dosp
    75 !
    76 ! Revision 1.16  2000/12/28 13:38:46  raasch
    77 ! Module dvrp_variables removed, -particle_transport
    78 !
    79 ! Revision 1.15  2000/04/27 06:57:04  raasch
    80 ! +dt_dvrp, dvrp_filecount, npex, npey, threshold, time_dvrp
    81 ! -dt_plisos, dt_plpart, time_plisos, time_plpart,
    82 ! all comments translated into English, old revision remarks deleted
    83 !
    84 ! Revision 1.14  2000/04/13 13:06:16  schroeter
    85 ! Einlesen der Initialisierungsparameter fuer Rechnungen
    86 ! mit Feuchte/Wolkenphysik
    87 !
    88 ! Revision 1.13  2000/01/10  10:06:32  10:06:32  raasch (Siegfried Raasch)
    89 ! +use_ug_for_galilei_tr
    90 !
    9116! Revision 1.1  1998/03/18 20:18:48  raasch
    9217! Initial revision
  • palm/trunk/SOURCE/run_control.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: run_control.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.20  2006/06/02 15:23:47  raasch
    1214! cpp-directives extended for lctit
    13 !
    14 ! Revision 1.19  2006/02/23 12:54:21  raasch
    15 ! Output buffer is flushed also on NEC, nt_anz renamed current_timestep_number
    16 !
    17 ! Revision 1.18  2004/01/30 10:38:26  raasch
    18 ! Euler-step information removed in case of Runge-Kutta
    19 !
    20 ! Revision 1.17  2003/10/29 09:13:18  raasch
    21 ! Data are forced to be immediately written on disc on IBM and Linux clusters
    22 !
    23 ! Revision 1.16  2003/04/16 13:12:39  raasch
    24 ! Output format of iteration count enlarged
    25 !
    26 ! Revision 1.15  2001/07/20 13:13:45  raasch
    27 ! Additional output of the number of multigrid cycles used
    28 !
    29 ! Revision 1.14  2001/03/30 07:52:24  raasch
    30 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    31 !
    32 ! Revision 1.13  2001/01/22 08:04:54  raasch
    33 ! Module test_variables removed
    34 !
    35 ! Revision 1.12  2001/01/03 13:46:12  letzel
    36 ! File output and all comments translated into English.
    37 !
    38 ! Revision 1.11  1998/09/22 17:29:49  raasch
    39 ! Ausgabeformate t.w. um einige Kommastellen erweitert
    40 !
    41 ! Revision 1.10  1998/07/06 12:32:56  raasch
    42 ! + USE test_variables
    43 !
    44 ! Revision 1.9  1998/04/21 16:23:48  raasch
    45 ! Ausgabe von advected_distance
    46 !
    47 ! Revision 1.8  1998/03/30 11:50:16  raasch
    48 ! Divergenzen stehen jetzt in hom
    49 !
    50 ! Revision 1.7  1998/03/25 13:57:06  raasch
    51 ! Erweiterungen fuer Gebietsstatistiken,
    52 ! dt in dt_3d umbenannt
    53 !
    54 ! Revision 1.6  1998/03/18 20:19:50  raasch
    55 ! Pruefen, ob Einheit 15 bereits geoeffnet ist
    56 !
    57 ! Revision 1.5  1998/03/03 08:01:49  raasch
    58 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier,
    59 ! Ausgabe von w* und z_i
    60 !
    61 ! Revision 1.4  1998/02/04 16:10:46  raasch
    62 ! e_total und e_pertu stehen jetzt im Feld hom
    63 !
    64 ! Revision 1.3  1998/01/23 10:08:01  raasch
    65 ! Ausgabe von u* und theta*
    66 !
    67 ! Revision 1.2  1997/08/26 06:34:36  raasch
    68 ! Schoenheitskorrekturen
    6915!
    7016! Revision 1.1  1997/08/11 06:25:38  raasch
  • palm/trunk/SOURCE/singleton.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: singleton.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2004/04/30 12:52:09  raasch
    1214! Shape of arrays is explicitly stored in ishape and handled to the
  • palm/trunk/SOURCE/sor.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: sor.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.9  2005/03/26 21:02:23  raasch
    1214! Implementation of non-cyclic (Neumann) horizontal boundary conditions,
    1315! dx2,dy2 replaced by ddx2,ddy2
    14 !
    15 ! Revision 1.8  2001/03/30 07:52:44  raasch
    16 ! Arrays f1-f3 changed from dummy arguments to allocatable arrays,
    17 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    18 !
    19 ! Revision 1.7  2001/01/30 21:41:25  letzel
    20 ! All comments translated into English.
    21 !
    22 ! Revision 1.6  2001/01/22 08:06:14  raasch
    23 ! Module test_variables removed
    24 !
    25 ! Revision 1.5  1998/07/06 12:33:13  raasch
    26 ! + USE test_variables
    27 !
    28 ! Revision 1.4  1997/09/12 06:30:12  raasch
    29 ! Randbedingungen umbenannt
    30 !
    31 ! Revision 1.3  1997/09/09 08:30:44  raasch
    32 ! Kehrwerte der Gitterweiten implementiert
    33 !
    34 ! Revision 1.2  1997/08/29 09:00:58  raasch
    35 ! omega --> omega_sor
    3616!
    3717! Revision 1.1  1997/08/11 06:25:56  raasch
     
    4222! ------------
    4323! Solve the Poisson-equation with the SOR-Red/Black-scheme.
    44 !-------------------------------------------------------------------------------!
     24!------------------------------------------------------------------------------!
    4525
    4626    USE grid_variables
  • palm/trunk/SOURCE/spline_x.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: spline_x.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.8  2004/04/30 12:54:20  raasch
    1214! Names of transpose indices changed, enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.7  2003/03/16 09:47:17  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.6  2001/03/30 07:53:19  raasch
    18 ! Arrays r, vad and wrk_spline changed from 3D to 2D and removed from argument
    19 ! list. Several loops over k/j (parallel/non-parallel) combined to one loop.
    20 ! Application of long filter moved to this routine. All comments and identifiers
    21 ! translated into English.
    22 !
    23 ! Revision 1.5  2001/01/22 08:07:41  raasch
    24 ! Module test_variables removed
    25 !
    26 ! Revision 1.4  1999/03/25 07:32:31  raasch
    27 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,
    28 ! Ueberschwinger werden in gewissen Grenzen erlaubt
    29 !
    30 ! Revision 1.3  1999/02/26 17:53:26  schroeter
    31 ! - Gradientenkontrolle fuer den nicht-parallelen Teil
    32 ! - statistische Auswertung ueber den prozentualen Anteil des
    33 !   Upstream-Verfahrens an der Gesamtadvektion fuer nicht-
    34 !   parallelen Teil
    35 !
    36 ! Revision 1.2  1999/02/17 09:30:19  raasch
    37 ! Wertebegrenzung zur Verhinderung von Ueberschwingern
    3815!
    3916! Revision 1.1  1999/02/05 09:15:59  raasch
  • palm/trunk/SOURCE/spline_y.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: spline_y.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.9  2004/04/30 12:54:37  raasch
    1214! Names of transpose indices changed, enlarged transposition arrays introduced
    13 !
    14 ! Revision 1.8  2003/03/16 09:48:41  raasch
    15 ! Two underscores (_) are placed in front of all define-strings
    16 !
    17 ! Revision 1.7  2001/03/30 07:53:40  raasch
    18 ! Arrays r and wrk_spline changed from 3D to 2D and removed from argument
    19 ! list. Several loops over k/i (parallel/non-parallel) combined to one loop.
    20 ! Application of long filter moved to this routine. All comments and identifiers
    21 ! translated into English.
    22 !
    23 ! Revision 1.6  2001/01/22 08:09:04  raasch
    24 ! Module test_variables removed
    25 !
    26 ! Revision 1.5  1999/11/25 16:30:42  raasch
    27 ! Laufindex-Fehler korrigiert
    28 !
    29 ! Revision 1.4  1999/03/25 07:33:16  raasch
    30 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,
    31 ! Ueberschwinger werden in gewissen Grenzen erlaubt
    32 !
    33 ! Revision 1.3  1999/02/26 17:54:28  schroeter
    34 ! - Gradientenkontrolle fuer den nicht-parallelen Teil
    35 ! - statistische Auswertung ueber den prozentualen Anteil des
    36 !   Upstream-Verfahrens an der Gesamtadvektion fuer nicht-
    37 !   parallelen Teil
    38 !
    39 ! Revision 1.2  1999/02/17 09:31:34  raasch
    40 ! Test einer Wertebegrenzung zur Verhinderung von Ueberschwingern
    4115!
    4216! Revision 1.1  1999/02/05 09:16:31  raasch
  • palm/trunk/SOURCE/spline_z.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: spline_z.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.9  2005/06/29 08:22:56  steinfeld
    1214! Dependency of ug and vg on height considered in the determination of the
    1315! upper boundary condition for vad
    14 !
    15 ! Revision 1.8  2004/04/30 12:55:23  raasch
    16 ! Enlarged transposition arrays introduced
    17 !
    18 ! Revision 1.7  2003/03/16 09:49:33  raasch
    19 ! Two underscores (_) are placed in front of all define-strings
    20 !
    21 ! Revision 1.6  2001/03/30 07:54:07  raasch
    22 ! Arrays r vad, and wrk_spline changed from 3D to 2D and removed from argument
    23 ! list. Several loops over i combined to one loop.
    24 ! Application of long filter moved to this routine. All comments and identifiers
    25 ! translated into English.
    26 !
    27 ! Revision 1.5  2001/01/22 08:10:21  raasch
    28 ! Module test_variables removed
    29 !
    30 ! Revision 1.4  1999/03/25 07:34:55  raasch
    31 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,
    32 ! Ueberschwinger werden in gewissen Grenzen erlaubt
    33 !
    34 ! Revision 1.3  1999/02/26 17:55:41  schroeter
    35 ! - Gradientenkontrolle fuer den nicht-parallelen Teil
    36 ! - statistische Auswertung ueber den prozentualen Anteil des
    37 !   Upstream-Verfahrens an der Gesamtadvektion fuer nicht-
    38 !   parallelen Teil
    39 !
    40 ! Revision 1.2  1999/02/17 09:32:55  raasch
    41 ! Wertebegrenzung zur Verhinderung von Ueberschwingern
    4216!
    4317! Revision 1.1  1999/02/05 09:17:16  raasch
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: sum_up_3d_data.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.1  2006/02/23 12:55:23  raasch
    1214! Initial revision
  • palm/trunk/SOURCE/swap_timelevel.f90

    r1 r3  
    11 SUBROUTINE swap_timelevel
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: swap_timelevel.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.8  2004/01/28 15:28:18  raasch
    1214! Swaping for Runge-Kutta schemes implemented
    13 !
    14 ! Revision 1.7  2002/05/02 18:58:36  raasch
    15 ! Swaping of prognostic quantities re-changed to three timelevels
    16 !
    17 ! Revision 1.6  2001/08/21  10:01:58  10:01:58  raasch (Siegfried Raasch)
    18 ! Surface momentum fluxes are swapped also when they are prescribed
    19 !
    20 ! Revision 1.5  2001/03/30 07:54:39  raasch
    21 ! Swaping of prognostic quantities changed over from three to two timelevels,
    22 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    23 !
    24 ! Revision 1.4  2001/01/29 12:36:03  raasch
    25 ! Passive scalar is considered
    26 !
    27 ! Revision 1.3  2001/01/25 07:26:48  raasch
    28 ! Sensible and latent heat fluxes are swaped, when use_surface_fluxes = .TRUE.,
    29 ! module test_variables removed
    30 !
    31 ! Revision 1.2  2000/04/13 13:59:26  schroeter
    32 ! + swapping timelevels of moisture variables
    3315!
    3416! Revision 1.1  2000/01/10  10:08:58  10:08:58  raasch (Siegfried Raasch)
     
    3921! ------------
    4022! Swap of timelevels of variables after each timestep
    41 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    4224
    4325    USE arrays_3d
  • palm/trunk/SOURCE/temperton_fft.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: temperton_fft.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2003/04/16 12:49:25  raasch
    1214! Abort in case of illegal factors
  • palm/trunk/SOURCE/time_integration.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: time_integration.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.8  2006/08/22 14:16:05  raasch
    1214! Disturbances are imposed only for the last Runge-Kutta-substep
    13 !
    14 ! Revision 1.7  2006/08/04 15:04:39  raasch
    15 ! Output of particle timeseries
    16 !
    17 ! Revision 1.6  2006/02/23 12:58:38  raasch
    18 ! Data output (1d, 2d, 3d, sp) can be skipped for a given time interval from
    19 ! simulation start, additional argument nzb_.._inner added to call of
    20 ! disturb_field, nt_anz renamed current_timestep_number, pl.. renamed do..,
    21 ! ..anz renamed ..n, dt_average renamed dt_averaging_interval_pr,
    22 ! time_average_.. renamed time_do.._av, routines plot_.. renamed data_output_..,
    23 ! additional call of user_actions (after_integration)
    24 !
    25 ! Revision 1.5  2005/06/26 20:10:20  raasch
    26 ! +call of interaction_droplets_ptq (+module of the same name)
    27 !
    28 ! Revision 1.4  2005/04/23 09:39:01  raasch
    29 ! Revised calculation of output time counters regarding a possible decrease of
    30 ! the output time interval in case of restart runs
    31 !
    32 ! Revision 1.3  2005/03/26 21:10:26  raasch
    33 ! Pres is called by default only at the last Runge-Kutta-substep.
    34 ! Determination of the timestep moved from end to the beginning of the
    35 ! time loop in order to assure consistent output in run_control (first line
    36 ! in restart runs must be equal to the last line in the corresponding
    37 ! previous run). Exchange of ghost points for km and kh removed.
    38 ! Changes regarding non-cyclic boundary conditions:
    39 ! arguments are added to call of routines disturb_field, exchange_horiz,
    40 ! and boundary_conds, which is now called for a second time after pres in
    41 ! order to set the velocity boundary conditions at the outflow. Velocity
    42 ! perturbations are imposed near the inflow throughout the whole run.
    4315!
    4416! Revision 1.2  2004/04/30 13:03:40  raasch
     
    4719! impulse_advec renamed momentum_advec
    4820!
    49 ! Revision 1.1  2004/04/30 12:58:21  raasch
    50 ! Initial revision
    51 !
    52 ! Revision 1.44  2004/01/28 15:17:26  raasch
    53 ! Additional intermediate timesteps for the Runge-Kutta schemes, setting
    54 ! of the steering factors for prognostic equations within the new
    55 ! routine timestep_scheme_steering at the beginning of each step (former
    56 ! part of routine timestep called at the end of each step), routine
    57 ! run_control is not called for the Runge-Kutta schemes automatically in case
    58 ! of timestep changes.
    59 !
    60 ! Revision 1.43  2003/03/16 09:40:12  raasch
    61 ! Two underscores (_) are placed in front of all define-strings
    62 !
    63 ! Revision 1.42  2003/03/12 16:33:23  raasch
    64 ! New routine prognostic_equations_vec is used on NEC machines
    65 !
    66 ! Revision 1.41  2003/03/04 11:31:39  raasch
    67 ! Module DVRP added in case that dvrp graphics package is switched on
    68 !
    69 ! Revision 1.40  2002/12/19 15:45:13  raasch
    70 ! Routine check_cpu_time renamed to check_for_restart
    71 !
    72 ! Revision 1.39  2002/06/11 13:11:44  raasch
    73 ! +USE user_actions_mod, prognostic_equations_mod,
    74 ! temperature offset at cyclic boundaries in case of sloping surface moved
    75 ! from prognostic_equations to this subroutine. Particle advection moved from
    76 ! prognostic_equations to this routine. A fast version of prognostic_equations
    77 ! is called in case that the Piascek-Williams advection scheme is used.
    78 !
    79 ! Revision 1.38  2002/05/02 18:52:19  raasch
    80 ! Application of time filter and exchange of ghost points moved from
    81 ! prognostic_equations to this routine.
    82 ! Error in calculating output time intervals for various quantities removed.
    83 !
    84 ! Revision 1.37  2001/08/21  09:51:26  09:51:26  raasch (Siegfried Raasch)
    85 ! Time measurements for user_actions moved to that subroutine
    86 !
    87 ! Revision 1.36  2001/03/30 07:34:24  raasch
    88 ! Timelevel t+dt eliminated, Asselin filter now included in
    89 ! prognostic_equations, arguments removed from pres,
    90 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    91 !
    92 ! Revision 1.35  2001/01/25 07:06:43  raasch
    93 ! Time measurements with dvrp software included
    94 !
    95 ! Revision 1.34  2001/01/05 15:12:50  raasch
    96 ! Calculation and plot of spectra
    97 !
    98 ! Revision 1.33  2000/04/27 06:40:49  raasch
    99 ! Old revision remarks deleted, output of isosurfaces and particles now
    100 ! combined into new routine data_output_dvrp (includes former routines
    101 ! plot_particles and plot_isosurface)
    102 !
    103 ! Revision 1.32  2000/04/13 13:39:30  schroeter
    104 ! computation of prognostic equations moved to the new routine
    105 ! prognostic_equation, new routines calc_liquid_water_content and
    106 ! compute_vpt
    107 !
    108 ! Revision 1.31  2000/01/25  16:25:24  16:25:24  letzel (Marcus Letzel)
    109 ! All comments translated into English
    110 !
    11121! Revision 1.1  1997/08/11 06:19:04  raasch
    11222! Initial revision
  • palm/trunk/SOURCE/time_to_string.f90

    r1 r3  
    11 FUNCTION time_to_string( time )
    22
    3 !-------------------------------------------------------------------------------!
     3!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
     
    88! Former revisions:
    99! -----------------
    10 ! $Log: time_to_string.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.3  2001/01/22 08:16:04  raasch
    1214! Comments translated into English
    13 !
    14 ! Revision 1.2  1998/03/09 16:24:53  raasch
    15 ! Dreistellige Stundenausgabe ermoeglicht
    1615!
    1716! Revision 1.1  1997/08/11 06:26:08  raasch
     
    2221! ------------
    2322! Transforming the time from real to character-string hh:mm:ss
    24 !-------------------------------------------------------------------------------!
     23!------------------------------------------------------------------------------!
    2524
    2625    IMPLICIT NONE
  • palm/trunk/SOURCE/timestep.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: timestep.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.21  2006/02/23 12:59:44  raasch
    1214! nt_anz renamed current_timestep_number
    13 !
    14 ! Revision 1.20  2005/04/23 09:45:38  raasch
    15 ! fcl_factor renamed cfl_factor
    16 !
    17 ! Revision 1.19  2005/03/26 21:11:53  raasch
    18 ! In case of non-cyclic lateral boundary conditions, the increased diffusivity
    19 ! in the lateral damping layer is regarded in the minimum timestep calculation.
    20 !
    21 ! Revision 1.18  2004/01/30 10:40:02  raasch
    22 ! Timestep increment limitation in case of Runge-Kutta removed
    23 !
    24 ! Revision 1.17  2004/01/28 15:34:03  raasch
    25 ! Setting of steering factors for the prognostic equations moved to new
    26 ! routine timestep_scheme_steering now called in routine leap_frog.
    27 !
    28 ! Revision 1.16  2003/03/16 09:49:47  raasch
    29 ! Two underscores (_) are placed in front of all define-strings
    30 !
    31 ! Revision 1.15  2002/12/19 16:19:54  raasch
    32 ! Two mpi_allreduce calls for determining mean horizontal velocity reduced to
    33 ! one call. Calculation of minimum timestep for diffusion optimized by
    34 ! determining minimum grid spacing in a single k loop.
    35 !
    36 ! Revision 1.14  2001/03/30 07:55:12  raasch
    37 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    38 !
    39 ! Revision 1.13  2001/01/30 22:10:11  letzel
    40 ! All comments translated into English.
    41 !
    42 ! Revision 1.12  2001/01/22 08:17:33  raasch
    43 ! Module test_variables removed
    44 !
    45 ! Revision 1.11  2000/01/10 10:10:30  raasch
    46 ! Translationsgeschwindigkeit fuer Galilei-Transformation wird nur noch
    47 ! bedingt berechnet
    48 !
    49 ! Revision 1.10  1999/02/05 09:19:06  raasch
    50 ! Zeitschrittsteuerung beruecksichtigt jetzt auch Wahl des reinen Euler-
    51 ! Verfahrens
    52 !
    53 ! Revision 1.9  1998/09/22 17:30:41  raasch
    54 ! Ausfuehrlichere Informationen bei Abbruch durch Unterschreitung des
    55 ! minimalen Zeitschritts
    56 !
    57 ! Revision 1.8  1998/07/06 12:33:52  raasch
    58 ! Korrektur der ermittelten horizontalen Geschwindigkeitsmaxima, falls diese
    59 ! bei k=0 gefunden wurden (wegen Spiegelungsrandbedingung treten dort negative
    60 ! Werte auf, die bei Galilei-Transformation den Zeitschritt erheblich
    61 ! einschraenken). + USE test_variables
    62 !
    63 ! Revision 1.7  1998/04/21 15:57:46  raasch
    64 ! Implementierung der Galilei-Transformation
    65 !
    66 ! Revision 1.6  1998/03/25 13:57:20  raasch
    67 ! dt in dt_3d umbenannt
    68 !
    69 ! Revision 1.5  1998/03/09 16:25:10  raasch
    70 ! Schoenheitskorrekturen
    71 !
    72 ! Revision 1.4  1998/03/03 08:02:15  raasch
    73 ! Leap-Frog auch nach Zeitschrittaenderung erlaubt
    74 !
    75 ! Revision 1.3  1997/09/12 06:30:27  raasch
    76 ! Einbau des Diffusions-Kriteriums
    77 !
    78 ! Revision 1.2  1997/08/26 06:34:55  raasch
    79 ! Fehler bei w_max_ijk - Index
    8015!
    8116! Revision 1.1  1997/08/11 06:26:19  raasch
  • palm/trunk/SOURCE/timestep_scheme_steering.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: timestep_scheme_steering.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2005/03/26 21:17:06  raasch
    1214! No pressure term for Runge-Kutta-schemes (tsc(4)=0.0)
  • palm/trunk/SOURCE/transpose.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: transpose.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.2  2004/04/30 13:12:17  raasch
    1214! Switched from mpi_alltoallv to the simpler mpi_alltoall,
     
    1719! Initial revision (collection of former routines transpose_xy, transpose_xz,
    1820!                   transpose_yx, transpose_yz, transpose_zx, transpose_zy)
    19 !
    20 ! Revision 1.7  2003/03/16 09:49:53  raasch
    21 ! Two underscores (_) are placed in front of all define-strings
    22 !
    23 ! Revision 1.6  2001/03/30 07:55:50  raasch
    24 ! Array bound for work2 is explicitely calculated, all comments and identifiers
    25 ! translated into English
    26 !
    27 ! Revision 1.5  2001/01/22 08:19:14  raasch
    28 ! Module test_variables removed
    29 !
    30 ! Revision 1.4  1998/07/06 12:34:22  raasch
    31 ! + USE test_variables
    32 !
    33 ! Revision 1.3  1997/08/13 13:52:42  raasch
    34 ! Zeitmessung von MPI_ALLTOALLV
    35 !
    36 ! Revision 1.2  1997/08/11 06:27:31  raasch
    37 ! Keine Aenderungen
    3821!
    3922! Revision 1.1  1997/07/24 11:25:18  raasch
  • palm/trunk/SOURCE/user_interface.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: user_interface.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.18  2006/06/02 15:25:00  raasch
    1214! +change of grid-defining arguments in routine user_define_netcdf_grid,
    1315! new argument "found" in user_data_output_2d and user_data_output_3d
    14 !
    15 ! Revision 1.17  2006/04/11 14:59:28  raasch
    16 ! Reading and writing user-defined quantities from/to the restart file improved
    17 !
    18 ! Revision 1.16  2006/03/03 20:00:48  raasch
    19 ! Error removed from last version
    20 !
    21 ! Revision 1.15  2006/03/03 19:55:36  raasch
    22 ! Call of user_init_grid now requires the formal parameter nzb_local
    23 !
    24 ! Revision 1.14  2006/02/23 13:13:18  raasch
    25 ! Additional interfaces for the calculation and output of user-defined 2d- and
    26 ! 3d-data
    27 !
    28 ! Revision 1.13  2003/03/12 16:44:56  raasch
    29 ! Full code replaced in routine user_actions (in the call for all gridpoints)
    30 ! instead of calling the _ij version (required by NEC, because otherwise no
    31 ! vectorization)
    32 !
    33 ! Revision 1.12  2002/12/19 16:22:04  raasch
    34 ! SAVE added in module user, STOP statement replaced by call of subroutine
    35 ! local_stop
    36 !
    37 ! Revision 1.11  2002/06/11 13:22:09  raasch
    38 ! Former subroutine changed to a module which allows to be called for all grid
    39 ! points of a single vertical column with index i,j or for all grid points by
    40 ! using function overloading.
    41 ! Log point for user_actions removed from the grid point version.
    42 !
    43 ! Revision 1.10  2001/08/21 12:13:55  raasch
    44 ! Wrong module in user_actions changed
    45 !
    46 ! Revision 1.9  2001/08/21 10:04:30  raasch
    47 ! User_actions can be called from every prognostic equation,
    48 ! subroutine user_statistics, user_particle_attributes, user_dvrp_coltab added
    49 !
    50 ! Revision 1.8  2001/07/12 12:45:07  raasch
    51 ! Subroutine user_init_particles added
    52 !
    53 ! Revision 1.7  2001/03/30 07:58:52  raasch
    54 ! Translation of remaining German identifiers (variables, subroutines, etc.)
    55 !
    56 ! Revision 1.6  2000/08/07 12:43:16  raasch
    57 ! All comments translated into English
    58 !
    59 ! Revision 1.5  1998/04/06 14:42:21  raasch
    60 ! +SUBROUTINE user_last_actions
    61 !
    62 ! Revision 1.4  1998/03/30 11:38:37  raasch
    63 ! Fehlerkorrekturen und Ausgabeaenderungen
    64 !
    65 ! Revision 1.3  1998/03/25 20:11:17  raasch
    66 ! +region in userpar
    67 !
    68 ! Revision 1.2  1998/03/25 13:57:37  raasch
    69 ! userpar beinhaltet jetzt keine vordefinierten Variablen mehr
    7016!
    7117! Revision 1.1  1998/03/24 15:29:04  raasch
  • palm/trunk/SOURCE/write_3d_binary.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: write_3d_binary.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.21  2006/08/04 15:05:11  raasch
    1214! +iran, iran_part
    13 !
    14 ! Revision 1.20  2006/02/23 13:14:34  raasch
    15 ! +volume_flow_area, volume_flow_initial, -hom_sum_pr,
    16 ! output of arrays used for time-averaging
    17 !
    18 ! Revision 1.19  2005/05/18 15:55:45  raasch
    19 ! pl1d_callcount renamed dopr_time_count
    20 !
    21 ! Revision 1.18  2004/04/30 13:13:40  raasch
    22 ! Arrays of timelevel t-dt (_m) are written only in case of the leapfrog scheme
    23 !
    24 ! Revision 1.17  2003/10/29 09:14:11  raasch
    25 ! Output of random seed (random_iv, random_iy), module random_function_mod
    26 ! is used, new binary version 2.2
    27 !
    28 ! Revision 1.16  2002/12/19 16:23:08  raasch
    29 ! Writing of array hom moved to write_var_list
    30 !
    31 ! Revision 1.15  2001/03/30 07:59:18  raasch
    32 ! Translation of remaining German identifiers (variables, subroutines, etc.),
    33 ! closing of unit 14 moved to the main program in order to allow user-defined
    34 ! writing to this unit
    35 !
    36 ! Revision 1.14  2001/03/11 11:27:18  letzel
    37 ! All comments translated into English.
    38 !
    39 ! Revision 1.13  2001/01/29 12:36:56  raasch
    40 ! Passive scalar is considered
    41 !
    42 ! Revision 1.12  2001/01/22 08:26:59  raasch
    43 ! Module test_variables removed
    44 !
    45 ! Revision 1.11  2001/01/05 15:16:39  raasch
    46 ! +spectrum_x, spectrum_y, new binary version 2.1
    47 !
    48 ! Revision 1.10  2001/01/02 17:34:45  raasch
    49 ! Unit 14 closed at the end of the subroutine
    50 !
    51 ! Revision 1.9  2000/04/13 14:12:37  schroeter
    52 ! Output of water and humidity variables
    53 !
    54 ! Revision 1.8  99/03/03  09:27:43  09:27:43  raasch (Siegfried Raasch)
    55 ! Speziellen I/O-Messpunkt entfernt
    56 !
    57 ! Revision 1.7  1999/02/05 09:19:39  raasch
    58 ! program_version in binary_version umbenannt
    59 !
    60 ! Revision 1.6  1998/12/15 07:34:40  raasch
    61 ! Zusaetzliche Ermittelung der tatsaechlichen Ausgabezeiten auf T3E
    62 ! mit messpunkt_e(15) unter Verwendung von TIMEF (s. cpu_zeitmessung)
    63 !
    64 ! Revision 1.5  1998/07/06 12:36:26  raasch
    65 ! + USE test_variables
    66 !
    67 ! Revision 1.4  1998/04/15 11:24:30  raasch
    68 ! Steuerparameter-Ausgaben gehen jetzt auch auf Einheit 14 und werden von
    69 ! jedem PE ausgegeben
    70 !
    71 ! Revision 1.3  1998/04/06 14:34:21  raasch
    72 ! Ausgabe von weiteren Plotvariablen
    73 !
    74 ! Revision 1.2  1998/03/30 11:39:00  raasch
    75 ! Statistik- und Plotvariablen werden zusaetzlich ausgegeben
    7615!
    7716! Revision 1.1  1998/03/18 20:20:21  raasch
  • palm/trunk/SOURCE/write_compressed.f90

    r1 r3  
    99! Former revisions:
    1010! ---------------------
    11 ! $Log: write_compressed.f90,v $
     11! $Id$
     12! RCS Log replace by Id keyword, revision history cleaned up
     13!
    1214! Revision 1.4  2006/02/23 13:15:09  raasch
    1315! nz_plot3d renamed nz_do3d
    14 !
    15 ! Revision 1.3  2003/03/16 09:51:07  raasch
    16 ! Two underscores (_) are placed in front of all define-strings
    17 !
    18 ! Revision 1.2  2001/03/30 07:59:48  raasch
    19 ! All comments translated into English
    2016!
    2117! Revision 1.1  1999/03/02 09:25:21  raasch
  • palm/trunk/SOURCE/write_var_list.f90

    r1 r3  
    88! Former revisions:
    99! -----------------
    10 ! $Log: write_var_list.f90,v $
     10! $Id$
     11! RCS Log replace by Id keyword, revision history cleaned up
     12!
    1113! Revision 1.34  2006/08/22 14:30:52  raasch
    1214! +dz_max
    13 !
    14 ! Revision 1.33  2006/08/04 15:05:52  raasch
    15 ! +first_call_advec_particles, time_dopts, use_upstream_for_tke
    16 !
    17 ! Revision 1.32  2006/02/23 13:16:02  raasch
    18 ! +building parameters, conserve_volume_flow, dissipation_1d, e_min,
    19 ! mixing_length_1d, topography,
    20 ! -average_count_pl, time_average_pr, time_average_pl renamed time_dopr_av,
    21 ! nt_anz renamed current_timestep_number, ..pl.. renamed ..do..
    22 !
    23 ! Revision 1.31  2005/06/28 14:18:02  steinfeld
    24 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
    25 ! ug_vertical_gradient_level_ind, vg_surface, vg_vertical_gradient,
    26 ! vg_vertical_gradient_level, vg_vertical_gradient_level_ind
    27 !
    28 ! Revision 1.30  2005/06/26 20:10:54  raasch
    29 ! +cloud_droplets
    30 !
    31 ! Revision 1.29  2005/05/18 15:56:20  raasch
    32 ! +netcdf_precision
    33 !
    34 ! Revision 1.28  2005/03/26 21:18:39  raasch
    35 ! +bc_lr, bc_ns, inflow_disturbance_begin, inflow_disturbance_end, km_damp_max,
    36 ! outflow_damping_width
    37 !
    38 ! Revision 1.27  2004/04/30 13:14:31  raasch
    39 ! +grid_matching, impulse_advec renamed momentum_advec
    40 !
    41 ! Revision 1.26  2004/01/28 15:38:55  raasch
    42 ! at and bt replaced by array tsc
    43 !
    44 ! Revision 1.25  2003/10/29 09:15:12  raasch
    45 ! +last_dt_change
    46 !
    47 ! Revision 1.24  2003/03/14 13:47:55  raasch
    48 ! +random_generator
    49 !
    50 ! Revision 1.23  2002/12/19 16:24:09  raasch
    51 ! Writing of array hom moved from write_3d_binary to this routine,
    52 ! binary version incremented to 2.2, +time_restart
    53 !
    54 ! Revision 1.22  2002/04/16 08:22:02  raasch
    55 ! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,
    56 ! s_vertical_gradient, s_vertical_gradient_level
    57 !
    58 ! Revision 1.21  2001/08/21 10:05:05  raasch
    59 ! +wall_adjustment
    60 !
    61 ! Revision 1.20  2001/03/30 08:00:07  raasch
    62 ! Translation of remaining German identifiers (variables, subroutines, etc.),
    63 ! runtime and package parameters are not written on unit 14 any more
    64 !
    65 ! Revision 1.19  2001/01/29 12:38:42  raasch
    66 ! Error writing use_surface_fluxes-string corrected, +passive_scalar
    67 !
    68 ! Revision 1.18  2001/01/25 07:27:48  raasch
    69 ! +fft_method, use_surface_fluxes
    70 !
    71 ! Revision 1.17  2001/01/05 15:16:59  raasch
    72 ! +average_count_sp, time_dosp
    73 !
    74 ! Revision 1.16  2000/12/28 13:39:21  raasch
    75 ! Module dvrp_variables removed, -particle_transport
    76 !
    77 ! Revision 1.15  2000/04/27 06:58:19  raasch
    78 ! +dt_dvrp, dvrp_filecount, npex, npey, threshold, time_dvrp
    79 ! -dt_plisos, dt_plpart, time_plisos, time_plpart
    80 ! +module pegrid,
    81 ! all comments translated into English, old revision remarks deleted
    82 !
    83 ! Revision 1.14  2000/04/13 14:10:58  schroeter
    84 ! Upgrade for calculations with cloud physics
    85 !
    86 ! Revision 1.13  2000/01/10  10:11:06  10:11:06  raasch (Siegfried Raasch)
    87 ! +use_ug_for_galilei_tr
    8815!
    8916! Revision 1.1  1998/03/18 20:20:38  raasch
Note: See TracChangeset for help on using the changeset viewer.