Diff [3:1] for / – PALM

/palm/trunk/SOURCE/Makefile

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

/palm/trunk/SOURCE/advec_particles.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: advec_particles.f90,v $
 ! Revision 1.32  2007/02/11 12:48:20  raasch
 ! Allways the lower level k is used for interpolation

/palm/trunk/SOURCE/advec_s_bc.f90

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

/palm/trunk/SOURCE/advec_s_pw.f90

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

/palm/trunk/SOURCE/advec_s_up.f90

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

/palm/trunk/SOURCE/advec_s_ups.f90

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

/palm/trunk/SOURCE/advec_u_pw.f90

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

/palm/trunk/SOURCE/advec_u_up.f90

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

/palm/trunk/SOURCE/advec_u_ups.f90

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

/palm/trunk/SOURCE/advec_v_pw.f90

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

/palm/trunk/SOURCE/advec_v_up.f90

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

/palm/trunk/SOURCE/advec_v_ups.f90

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

/palm/trunk/SOURCE/advec_w_pw.f90

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

/palm/trunk/SOURCE/advec_w_up.f90

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

/palm/trunk/SOURCE/advec_w_ups.f90

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

/palm/trunk/SOURCE/asselin_filter.f90

-                      r3
+                      r1
  SUBROUTINE asselin_filter
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! ---------------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: asselin_filter.f90,v $
 ! Revision 1.8  2004/01/30 10:14:02  raasch
 ! Scalar lower k index nzb replaced by 2d-array nzb_2d
+!
+! Revision 1.7  2003/03/16 09:26:58  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.6  2003/03/14 13:37:38  raasch
+! Do loop constructs only for IBM, otherwise array notation re-introduced
+!
+! Revision 1.5  2002/12/19 13:45:57  raasch
+! Array notation changed to do-loop constructs due to better performance,
+! OpenMP directives added
+!
+! Revision 1.4  2002/05/15 10:00:18  raasch
+! Nothing changed
+!
+! Revision 1.3  2002/05/02  19:39:07  19:39:07  raasch (Siegfried Raasch)
+! Use of module interfaces included
+!
+! Revision 1.2  2002/05/02 14:06:02  raasch
+! *** empty log message ***
+!
 ! Revision 1.1  2002/05/02  13:43:53  13:43:53  raasch (Siegfried Raasch)
 …
 ! -------------
 ! Time filter needed for the leap-frog method
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     USE arrays_3d

/palm/trunk/SOURCE/average_3d_data.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: average_3d_data.f90,v $
 ! Revision 1.1  2006/02/23 09:48:58  raasch
 ! Initial revision

/palm/trunk/SOURCE/boundary_conds.f90

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

/palm/trunk/SOURCE/buoyancy.f90

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

/palm/trunk/SOURCE/calc_liquid_water_content.f90

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

/palm/trunk/SOURCE/calc_precipitation.f90

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

/palm/trunk/SOURCE/calc_radiation.f90

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

/palm/trunk/SOURCE/calc_spectra.f90

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

/palm/trunk/SOURCE/check_for_restart.f90

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

/palm/trunk/SOURCE/check_open.f90

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

/palm/trunk/SOURCE/check_parameters.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: check_parameters.f90,v $
 ! Revision 1.61  2006/08/04 14:20:25  raasch
 ! do2d_unit and do3d_unit now defined as 2d-arrays, check of
 …
 ! generation of file header moved from routines palm and header to here
+!
+! Revision 1.60  2006/06/02 15:09:15  raasch
+! data_output_xy/xz/yz are set
+!
+! Revision 1.59  2006/03/03 15:55:34  letzel
+! Spelling inconsistency removed
+!
+! Revision 1.58  2006/03/03 08:50:08  raasch
+! Default setting of constant_heatflux now only depends on the value of
+! surface_heatflux again
+!
+! Revision 1.57  2006/02/27 10:16:44  letzel
+! Error removed in initialization of dt_averaging_input_pr
+!
+! Revision 1.56  2006/02/23 10:06:55  raasch
+! Volume flow conservation is not allowed with non-cyclic lateral boundaries,
+! calculation of nzb_diff moved to init_grid, the default setting of
+! constant_heatflux = .F. additionally requires a surface initial temperature
+! change /= 0, check of topography, dissipation_1d, mixing_length_1d,
+! check and setting of new data output steering parameters and user defined
+! output quantities,
+! ebene renamed section, pl1d_anz, plts_anz renamed dopr_n, dots_n,
+! plts, pl1d renamed data_output_ts/1d
+!
+! Revision 1.55  2005/12/06 16:39:34  raasch
+! Output of ql profile is allowed in case of using cloud droplets
+!
+! Revision 1.54  2005/06/30 11:27:33  steinfeld
+! Initialization of the profiles of the geostrophic wind components
+! ug and vg;
+! the usage of the geostrophic wind in a galilei transformation is
+! not allowed with a simultaneous specification of a baroclinicity
+!
+! Revision 1.53  2005/05/18 15:12:36  raasch
+! Check of data_output_format and netcdf precision,
+! abort if poisfft_hybrid is called in a non-parallel environment
+!
+! Revision 1.52  2005/04/23 09:04:10  raasch
+! fcl_factor renamed cfl_factor, default setting of outflow_damping_width
+! corrected, additional check for non-cyclic lateral boundary conditions
+!
+! Revision 1.51  2005/03/26 19:47:32  raasch
+! Check of bc_lr, bc_ns, km_damp_max, inflow_disturbance_begin,
+! inflow_disturbance_end, outflow_damping_width
+!
+! Revision 1.50  2004/04/30 11:20:27  raasch
+! Check of grid_matching, check of time series profiles modified,
+! impulse_advec renamed momentum_advec
+!
+! Revision 1.49  2004/01/28 15:04:56  raasch
+! Additional checks and settings for timestep schemes because of the
+! newly implemented Runge-Kutta schemes
+!
+! Revision 1.48  2003/04/17 07:07:43  raasch
+! Hybrid solver check in parallel environment only
+!
+! Revision 1.47  2003/04/16 09:43:20  raasch
+! Additional checks for hybrid solver and host, Monin Obukhov length added to
+! time series
+!
+! Revision 1.46  2003/03/31 15:43:19  raasch
+! Check for horizontal temperature flux profiles cleaned up
+!
+! Revision 1.45  2003/03/14 13:38:02  raasch
+! Random generator check included
+!
+! Revision 1.44  2003/03/12 16:23:10  raasch
+! Small checks for NEC and temperton-algorithm
+!
+! Revision 1.43  2002/12/19 14:04:26  raasch
+! Setting default values for termination_time_needed and check the
+! value for IBM-regatta, STOP statement replaced by call of subroutine
+! local_stop
+!
+! Revision 1.42  2002/09/12 13:02:25  raasch
+! Error in calculating initial pt and q profiles removed (array bound 10 in
+! *_ind arrays could be exceeded)
+!
+! Revision 1.41  2002/06/11 12:49:33  raasch
+! Hybrid solver allowed for solving the Poisson equation
+!
+! Revision 1.40  2002/05/02  18:48:56  18:48:56  raasch (Siegfried Raasch)
+! SELECT case of cross_normalized_x and ..y requires usage of TRIM-function
+! on IBM
+!
+! Revision 1.39  2002/04/16 08:03:27  raasch
+! Parameters for scalar transport included, output of scalar profiles changed
+!
+! Revision 1.38  2001/11/28 09:48:32  schroeter
+! Profiles for horizontal heat fluxes and
+! flux divergences allowed (58-72)
+!
+! Revision 1.37  2001/09/04  11:58:55  11:58:55  raasch (Siegfried Raasch)
+! Profiles for vertical flux of resolved/subgrid scale energy and pressure
+! fluctuations allowed (#55-57)
+!
+! Revision 1.36  2001/07/20 12:52:32  raasch
+! Multigrid method allowed for solving the Poisson equation
+!
+! Revision 1.35  2001/03/30 06:57:38  raasch
+! Asselin filter factor is set to 0, if Euler scheme is used,
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.34  2001/02/09 07:04:45  raasch
+! Checking of boundary conditions for pt and q and their combination with
+! prescribed surface fluxes revised
+!
+! Revision 1.33  2001/01/29 12:23:31  raasch
+! New parameter passive_scalar is checked
+!
+! Revision 1.32  2001/01/25 06:55:37  raasch
+! New control-parameters fft_method and use_surface_fluxes are checked
+!
+! Revision 1.31  2000/04/18 11:19:43  schroeter
+! Bott-Chlond Advektionsschema nun auch fuer Rechnungen
+! mit Feuchte/Wolkenphysik zugelassen
+!
+! Revision 1.30  2000/04/13 13:28:56  schroeter
+! + initialize profiles of total water content
+! + check boundary conditions of total water content
+! + initialize output of cloud-physics-profiles
+!
+! Revision 1.29  2000/01/25  14:38:37  14:38:37  letzel (Marcus Letzel)
+! All comments translated into English
+!
 ! Revision 1.1  1997/08/26 06:29:23  raasch
 ! Initial revision

/palm/trunk/SOURCE/close_file.f90

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

/palm/trunk/SOURCE/compute_vpt.f90

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

/palm/trunk/SOURCE/coriolis.f90

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

/palm/trunk/SOURCE/cpu_log.f90

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

/palm/trunk/SOURCE/cpu_statistics.f90

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

/palm/trunk/SOURCE/data_log.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: data_log.f90,v $
 ! Revision 1.1  2006/02/23 10:09:29  raasch
 ! Initial revision

/palm/trunk/SOURCE/data_output_2d.f90

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

/palm/trunk/SOURCE/data_output_3d.f90

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

/palm/trunk/SOURCE/data_output_dvrp.f90

-                      r3
+                      r1
 ! Actual revisions:
 ! -----------------
+!
 ! TEST: write statements
+!
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: data_output_dvrp.f90,v $
 ! Revision 1.13  2006/02/23 10:25:12  raasch
 ! Former routine plot_dvrp renamed data_output_dvrp,
 ! Only a fraction of the particles may have a tail,
 ! pl.. replaced by do.., %size renamed %dvrp_psize
+!
+! Revision 1.12  2005/03/31 07:00:46  raasch
+! Testversion with print statements checked in
+!
+! Revision 1.11  2003/08/01 08:44:50  raasch
+! Routine defined as recursive, because it is called from DVRP_STEERING_UPDATE
+! in case that the user has suspended the run by steering control.
+! Routine color_dvrp placed at the beginning of this file. This routine is used
+! for slicer colortable settings.
+!
+! Revision 1.10  2003/03/16 09:42:34  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.9  2003/03/04 11:35:07  raasch
+! Update for new version of dvrp software (e.g. using module dvrp is now
+! mandatory)
+!
+! Revision 1.8  2001/11/09 13:41:21  raasch
+! Colour can change along the particle tails
+!
+! Revision 1.7  2001/08/21 09:55:20  raasch
+! Superelevation for particles implemented, plotting of particle tails added,
+! colourtable definitions moved to subroutine user_dvrp_coltab
+!
+! Revision 1.6  2001/07/12 12:44:12  raasch
+! Test: color definitions for particles changed
+!
+! Revision 1.5  2001/03/30 07:45:00  raasch
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.4  2001/01/29 12:32:51  raasch
+! Passive scalar considered
+!
+! Revision 1.3  2001/01/22 07:41:12  raasch
+! Module test_variables removed
+!
+! Revision 1.2  2000/12/28 13:37:19  raasch
+! Total revision due to new features of dvrp-software,
+! plot of slicers included,
+! code is used only optionally (cpp-directives are added)
+!
 ! Revision 1.1  2000/04/27 06:27:17  raasch

/palm/trunk/SOURCE/data_output_profiles.f90

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

/palm/trunk/SOURCE/data_output_ptseries.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: data_output_ptseries.f90,v $
 ! Revision 1.2  2006/08/22 13:51:13  raasch
 ! Seperate output for particle groups

/palm/trunk/SOURCE/data_output_spectra.f90

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

/palm/trunk/SOURCE/data_output_tseries.f90

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

/palm/trunk/SOURCE/diffusion_e.f90

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

/palm/trunk/SOURCE/diffusion_s.f90

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

/palm/trunk/SOURCE/diffusion_u.f90

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

/palm/trunk/SOURCE/diffusion_v.f90

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

/palm/trunk/SOURCE/diffusion_w.f90

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

/palm/trunk/SOURCE/diffusivities.f90

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

/palm/trunk/SOURCE/disturb_field.f90

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

/palm/trunk/SOURCE/disturb_heatflux.f90

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

/palm/trunk/SOURCE/exchange_horiz.f90

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

/palm/trunk/SOURCE/exchange_horiz_2d.f90

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

/palm/trunk/SOURCE/fft_xy.f90

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

/palm/trunk/SOURCE/flow_statistics.f90

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

/palm/trunk/SOURCE/global_min_max.f90

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

/palm/trunk/SOURCE/header.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: header.f90,v $
 ! Revision 1.63  2006/08/22 13:53:13  raasch
 ! Output of dz_max
+!
+! Revision 1.62  2006/08/04 14:38:41  raasch
+! Generation of the run description header as well as run date and run time
+! moved to routine check_parameters,
+! output for additional particle quantities
+!
+! Revision 1.61  2006/04/26 12:19:51  raasch
+! Output of threads per task
+!
+! Revision 1.60  2006/04/11 14:57:01  raasch
+! pl_spectra renamed data_output_sp
+!
+! Revision 1.59  2006/03/14 12:44:10  raasch
+! Output of the initial geostrophic wind profile corrected
+!
+! Revision 1.58  2006/02/23 12:24:22  raasch
+! Output of psl, psr, pdx, etc. for each particle group, output of e_min,
+! information about volume flow conservation and usage of NetCDF 64 bit offset
+! format, pl.. renamed do.., skip_time control parameters for do..,
+! ebene renamed section, .._anz renamed .._no, dt_average renamed
+! dt_averaging_input_pr
+! output of 2d/3d data averaging informations, buffer is forced to be written
+! on file at the end of the routine
+!
+! Revision 1.57  2005/10/20 14:06:44  raasch
+! Error in output of yz-slice information removed
+!
+! Revision 1.56  2005/07/01 07:48:14  steinfeld
+! Information on the initial profiles of ug and vg added;
+! dependency of ug and vg on height considered in the output of
+! the boundary conditions for u and v
+!
+! Revision 1.55  2005/06/26 19:53:15  raasch
+! gas_constant renamed r_d, latent_heat renamed l_v, radius is used instead
+! of diameter, output formats 335+337 changed
+!
+! Revision 1.54  2005/05/19 08:32:27  raasch
+! Error in output of particle inertia information removed
+!
+! Revision 1.53  2005/05/18 15:31:59  raasch
+! Informations about data output format (netcdf, iso2d, etc.) are given
+!
+! Revision 1.52  2005/04/23 09:16:10  raasch
+! fcl_factor renamed cfl_factor
+!
+! Revision 1.51  2005/03/26 20:23:14  raasch
+! Output of horizontal boundary conditions, output of additional call of
+! pressure solver in case of Runge-Kutta schemes,
+! calculate the number of particle groups (formerly unknown = 0 )
+!
+! Revision 1.50  2004/04/30 11:48:02  raasch
+! Forcing a 1d-decomposition is output on linux-cluster and decalpha,
+! output of subdomain sizes, impulse_advec renamed momentum_advec
+!
+! Revision 1.49  2004/01/28 15:09:15  raasch
+! Output of timestep schemes modified
+!
+! Revision 1.48  2003/10/29 08:52:04  raasch
+! Additional output for multigrid method (lowest levels gathered on PE0),
+! particle output modified for particle groups
+!
+! Revision 1.47  2003/04/16 12:58:18  raasch
+! Output for mixing length limitations revised
+!
+! Revision 1.46  2003/03/16 09:39:31  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.45  2003/03/14 13:43:21  raasch
+! Informations about random generator
+!
+! Revision 1.44  2003/03/12 16:31:45  raasch
+! Small change for NEC system
+!
+! Revision 1.43  2002/12/19 14:53:47  raasch
+! Informations about user defined restart times added
+!
+! Revision 1.42  2002/09/12 13:03:24  raasch
+! Informations about particle inertia added
+!
+! Revision 1.41  2002/06/11 13:09:11  raasch
+! Information about hybrid solver, usage of 1d-decomposition and fast version
+! of prognostic_equations included,
+! no output of mixing length adjustments in case of constant diffusion
+!
+! Revision 1.40  2002/04/16  08:06:13  08:06:13  raasch (Siegfried Raasch)
+! Informations about particle boundary conditions and output of particle
+! data added, direct output of informations about scalar transport (instead
+! of using humidity)
+!
+! Revision 1.39  2001/08/21 08:50:54  raasch
+! Informations about mixing length adjustment, particle tails and dvrp PE usage
+! included
+!
+! Revision 1.38  2001/07/20 13:05:39  raasch
+! Information about multigrid method included
+!
+! Revision 1.37  2001/07/12 12:08:14  raasch
+! Additional informations about dvrp-output (host and filenames) and
+! particles (maximum age, random start positions), + module particle_attributes
+!
+! Revision 1.36  2001/03/30 07:23:53  raasch
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.35  2001/01/29 12:25:55  raasch
+! Informations about using passive scalar
+!
+! Revision 1.34  2001/01/25 07:00:48  raasch
+! Information about using surface fluxes and fft-method added
+!
+! Revision 1.33  2001/01/05 15:11:41  raasch
+! Spectra informations added. Old revision remarks deleted.
+!
+! Revision 1.32  2001/01/03 13:08:47  letzel
+! File output translated into English.
+!
+! Revision 1.31  2000/12/28 13:04:54  raasch
+! Slicer information added to dvrp-output.
+! Informations about optionally used packages are put in cpp-define-brackets.
+!
+! Revision 1.30  2000/04/27 07:08:14  raasch
+! notice if processor topology is predefined by user,
+! general notice about dvrp-output (including isosurface and particles),
+! all comments translated into English
+!
 ! Revision 1.1  1997/08/11 06:17:20  raasch

/palm/trunk/SOURCE/impact_of_latent_heat.f90

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

/palm/trunk/SOURCE/init_1d_model.f90

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

/palm/trunk/SOURCE/init_3d_model.f90

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

/palm/trunk/SOURCE/init_advec.f90

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

/palm/trunk/SOURCE/init_cloud_physics.f90

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

/palm/trunk/SOURCE/init_dvrp.f90

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

/palm/trunk/SOURCE/init_grid.f90

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

/palm/trunk/SOURCE/init_particles.f90

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

/palm/trunk/SOURCE/init_pegrid.f90

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

/palm/trunk/SOURCE/init_pt_anomaly.f90

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

/palm/trunk/SOURCE/init_rankine.f90

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

/palm/trunk/SOURCE/init_slope.f90

-                      r3
+                      r1
  SUBROUTINE init_slope
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: init_slope.f90,v $
 ! Revision 1.5  2006/02/23 12:35:34  raasch
 ! nanz_2dh renamed ngp_2dh
+!
+! Revision 1.4  2003/03/16 09:40:02  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.3  2001/03/30 07:33:59  raasch
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.2  2001/01/22 07:21:35  raasch
+! Module test_variables removed
+!
 ! Revision 1.1  2000/04/27 07:06:24  raasch
 …
 ! Remember: when a sloping surface is used, only one constant temperature
 !           gradient is allowed!
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     USE arrays_3d

/palm/trunk/SOURCE/interaction_droplets_ptq.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: interaction_droplets_ptq.f90,v $
 ! Revision 1.1  2005/06/26 19:57:47  raasch
 ! Initial revision

/palm/trunk/SOURCE/local_getenv.f90

-                      r3
+                      r1
  SUBROUTINE local_getenv( var, ivar, value, ivalue )
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: local_getenv.f90,v $
 ! Revision 1.5  2003/05/09 14:37:07  raasch
 ! On the MUK cluster, only PE0 is able to read environment variables.
 ! Therefore, they have to be communicated via broadcast to the other PEs.
+!
+! Revision 1.4  2003/03/16 09:40:22  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.3  2001/01/22 07:23:55  raasch
+! Comments translated into English
+!
+! Revision 1.2  1998/07/16 06:50:50  raasch
+! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
+!
 ! Revision 1.1  1997/08/11 06:21:01  raasch
 …
 ! ------------
 ! Getting the values of environment-variabls (for different operating-systems)
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 #if defined( __lcmuk )

/palm/trunk/SOURCE/local_stop.f90

-                      r3
+                      r1
  SUBROUTINE local_stop
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: local_stop.f90,v $
 ! Revision 1.2  2003/03/16 09:40:28  raasch
 ! Two underscores (_) are placed in front of all define-strings
 …
 ! ------------
 ! Stop program execution
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     USE pegrid

/palm/trunk/SOURCE/local_system.f90

-                      r3
+                      r1
  SUBROUTINE local_system( command )
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: local_system.f90,v $
 ! Revision 1.4  2003/03/16 09:40:33  raasch
 ! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.3  2001/01/22 07:25:35  raasch
+! Comments translated into English
+!
+! Revision 1.2  1998/07/16 06:51:19  raasch
+! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert
+!
 ! Revision 1.1  1997/09/03 06:27:27  raasch
 …
 ! ------------
 ! System calls for different operating systems
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     CHARACTER (LEN=*) ::  command

/palm/trunk/SOURCE/local_tremain.f90

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

/palm/trunk/SOURCE/local_tremain_ini.f90

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

/palm/trunk/SOURCE/modules.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: modules.f90,v $
 ! Revision 1.95  2007/02/11 13:18:30  raasch
 ! version 3.1b (last under RCS control)
+!
+! Revision 1.94  2006/08/22 14:12:18  raasch
+! +dz_max,
+! id_var_dopts extended for second dimension
+!
+! Revision 1.93  2006/08/04 14:50:58  raasch
+! +c_0, diss, dt_min_part, end_time_prel, particles_per_point, sgs_wfu/v/w_part,
+! use_sgs_for_particles, use_upstream_for_tke needed for the particle package,
+! +dt_sum, e_m, speed_x/y/z/_sgs in particle_type,
+! +dt_dopts, time_dopts, netcdf ids for output of particle time series,
+! do2d_unit and do3d_unit now defined as 2d-arrays, izuf renamed iran
+!
+! Revision 1.92  2006/06/02 15:21:45  raasch
+! +NetCDF ids for the shifted u- and v-grid, +data_output_xy/xz/yz
+!
+! Revision 1.91  2006/04/26 12:22:52  raasch
+! +threads_per_task,
+! sums_l, sums_l_l arrays changed to three dimensions
+!
+! Revision 1.90  2006/04/11 14:58:27  raasch
+! pl_spectra renamed data_output_sp, length of run_identifier increased to 80
+!
+! Revision 1.89  2006/03/14 12:55:52  raasch
+! Default for number_of_particle_groups changed to 1
+!
+! Revision 1.88  2006/03/03 14:26:53  letzel
+! control parameter topography now has length 40
+!
+! Revision 1.87  2006/03/03 08:50:25  raasch
+! Default value for surface_heatflux changed to 9999999.9
+!
+! Revision 1.86  2006/02/23 12:39:14  raasch
+! Version 3.0 with topography
+! index arrays nzb_2d replaced by several new variable dependend arrays
+! (e.g. nzb_u_inner), nzb_diff_2d removed, +new grid point
+! counters ngp_2dh_outer and ngp_3d_inner,
+! +skip time control parameters skip_time_..., +log_message,
+! particle parameters psl, psr, pss, psn, psb, pst, pdx, pdy, pdz
+! are now 1D arrays (1:max_number_of_particle_groups), default of
+! number_of_particle_groups changed to 999999, idum in particle_type renamed
+! tail_id,
+! +data_output, data_output_user, do2d, do3d, doav for new steering of data
+! output,
+! +type_x_int for ghost exchange of integer arrays,
+! +netcdf_64bit for supporting NetCDF 64-bit offset format,
+! +mixing_length_1d, dissipation_1d for steering the 1d-model,
+! +e_min (control parameter for minimum tke), +parameters for assuring a
+! constant volume flow, +prt_count, prt_start_index, particle_advection,
+! ebene renamed section, nt_anz renamed current_timestep_number, anz_y..
+! renamed ngp_y.., nanz_sums, nanz_2dh, nanz_3d renamed ngp_sums, ngp_2dh,
+! ngp_3d, respectively, pl1d_anz, pl2d_.._anz, plts_anz, pl3d_anz renamed
+! dopr_n, pl2d_.._no, dots_n, do3d_n, respectively,
+! ql_v and ql_vp now have the target attribute,
+! average_period_pr1d removed
+!
+! Revision 1.85  2005/12/06 16:40:56  raasch
+! Version 2.11c
+!
+! Revision 1.84  2005/10/20 15:38:23  raasch
+! +id_var_rnop_prt, some values of prt_var_names and prt_var_units changed
+!
+! Revision 1.83  2005/06/29 10:40:30  steinfeld
+! Scalars ug and vg changed into one-dimensional arrays;
+! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
+!  ug_vertical_gradient_level_ind,
+!  vg_surface, vg_vertical_gradient, vg_vertical_gradient_level,
+!  vg_vertical_gradient_level_ind
+!
+! Revision 1.82  2005/06/26 20:03:33  raasch
+! +ql_1, ql_2, ql_m, ql is now a pointer, +variables for cloud droplets (arrays
+! ql_v,ql_vp, pointer ql_c), gas_constant renamed r_d, latent_heat renamed l_v,
+! radius instead of diameter, radius + weighting factor included in
+! particle_type, pl2d and pl3d string length increased from 6 to 10,
+! default of call_psolver_at_all_substeps is TRUE
+!
+! Revision 1.81  2005/05/18 15:36:46  raasch
+! New module netcdf_control, +data_output_format, pl1d_callcount renamed
+! dopr_time_count,
+! default values of pl_spectra and spectra_direction changed
+!
+! Revision 1.80  2005/04/23 09:31:14  raasch
+! Implicit counters i renamed i9 due to declaration conflicts with i in other
+! subroutines (reported by Portland compiler),
+! crmax increased to 100, fcl_factor renamed cfl_factor
+!
+! Revision 1.79  2005/03/26 20:40:48  raasch
+! +variables for non-cyclic boundary conditions, call_psolver_at_all_substeps,
+! default value for residual limit increased from 1E-6 to 1E-4
+!
+! Revision 1.78  2004/04/30 12:30:01  raasch
+! idum added in particle_type for better alignment (decalpha),
+! DATA statements replaced by (/ .../) because of internal compiler errors
+! on decalpha, + indices for direct transpositions z --> y,
+! dt_max_1d increased to 300.0 (from 60.0), +grid_matching,
+! arrays cross_ts_numbers and cross_ts_number_count enlarged,
+! impulse_advec renamed momentum_advec
+!
+! Revision 1.77  2004/01/30 10:30:26  raasch
+! +index arrays nzb_2d, nzb_diff_2d
+!
+! Revision 1.76  2004/01/28 15:25:15  raasch
+! +array tsc for timestep steering, +tu_m, tv_m etc. pointer for additional
+! tendencies needed for Runge-Kutta scheme, +intermediate_timestep_count(_max)
+! type log changed to logpoint due to name conflict with intrinsic log
+!
+! Revision 1.75  2003/10/29 08:58:54  raasch
+! +mg_cycles, mg_switch_to_pe0, mg_switch_to_pe0_level in module
+! control_parameters,
+! +mg_loc_ind in module indices, variables related to particle groups added
+! to module particle_attributes
+!
+! Revision 1.74  2003/08/01 08:40:59  raasch
+! +isclice_dvrp, lock_steering_update, slicer_position_dvrp,
+! slicer_range_limits_dvrp in module dvrp_variables
+!
+! Revision 1.73  2003/05/09 14:41:11  raasch
+! New version number 2.7
+!
+! Revision 1.72  2003/04/25 18:54:39  raasch
+! Module mpi added in module pegrid in case of linux cluster
+!
+! Revision 1.71  2003/04/16 13:27:33  raasch
+! Arrays cross_ts_uy... enlarged
+!
+! Revision 1.70  2003/04/16 13:07:43  raasch
+! Array cross_ts_profiles enlarged
+!
+! Revision 1.69  2003/03/16 09:42:07  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.68  2003/03/14 13:44:52  raasch
+! +random_generator in control_parameters
+!
+! Revision 1.67  2003/03/12 16:34:19  raasch
+! +tasks_per_node in pegrid
+!
+! Revision 1.66  2003/03/04 11:32:50  raasch
+! +type steering, +steering_dvrp in dvrp_variables
+!
+! Revision 1.65  2002/12/19 15:52:37  raasch
+! +initial_wallclock_time in cpulog, +myid_char_14 in pegrid,
+! default value for termination_time_needed now set in check_parameters,
+! array dl_grid removed, +wall_adjustment_factor, dt_restart, restart_time,
+! time_restart
+!
+! Revision 1.64  2002/09/12 13:05:44  raasch
+! Default values of *_vertical_gradient_level changed from 100000.0 to -1.0,
+! +molecular_viscosity, density_ratio
+!
+! Revision 1.63  2002/06/11 13:12:29  raasch
+! Entrys buoyancy, diffusion_e, diffusion_s, diffusion_u, diffusion_v and
+! diffusion_w removed from module pointer_interfaces,
+! length of character variables host + psolver increased to 64 and 16
+!
+! Revision 1.62  2002/05/02 18:53:18  raasch
+! Module singleton moved to a separate file. Arrays e_3, pt_3, q_3, u_3, v_3,
+! w_3 re-introduced. Pointer e_p, pt_p, q_p, u_p, v_p, w_p re-introduced.
+!
+! Revision 1.61  2002/04/16 08:09:32  raasch
+! Information about speed components and start positions included in
+! particle_type, +dt_write_particle_data, particle_advection_start,
+! time_write_particle_data and parameters for particle boundary conditions,
+! +vertical_particle_advection,
+! + variables for scalar boundary conditions and initialization
+!
+! Revision 1.60  2001/11/12 16:05:01  raasch
+! +read_particles_from_restartfile
+!
+! Revision 1.58  2001/09/04 12:00:52  raasch
+! New version number 2.3b
+!
+! Revision 1.57  2001/08/21 09:52:08  raasch
+! +comm_palm, maximum_number_of_tailpoints, maximum_tailpoint_age,
+! minimum_tailpoint_distance, particle_tail_coordinates, superelevation_x,
+! superelevation_y, wall_adjustment
+!
+! Revision 1.56  2001/07/20 13:09:07  raasch
+! +grid_level, maximum_grid_level, cycle_mg, gamma_mg, ngsrb,
+! residual_limit in control_parameters, +dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg
+! in arrays_3d, +superelevation in dvrp_variables, +nxl_mg .. nzt_mg in
+! indices, +ddx2_mg, ddy2_mg in grid_variables
+!
+! Revision 1.55  2001/07/12 12:13:24  raasch
+! default of bc_e_b changed from (u*)**2+neumann to neumann,
+! +particle_maximum_age, psl-pst, pdx-pdz,
+! default values defined for dvrp_directory, dvrp_file, dvrp_host, dvrp_password
+! SAVE-statements added to all modules
+!
+! Revision 1.54  2001/03/30 07:35:57  raasch
+! Interface for advec_s_bc changed, interface for production_e eliminated,
+! arrays work, work1, e_3, pt_3, q_3, u_3, v_3, w_3 removed,
+! pointer e_p, pt_p, q_p, u_p, v_p, w_p removed,
+! -nw, nw1,
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.53  2001/01/30 01:57:16  raasch
+! +passive_scalar
+!
+! Revision 1.51  2001/01/26 00:31:49  raasch
+! +fft_method, use_surface_fluxes in steuer, +module singleton,
+! +write_particle_informations in particle_attributes
+!
+! Revision 1.49  2001/01/05 15:13:16  raasch
+! new module spectrum for computing spectra, +average_count_sp,
+! averaging_interval_sp, dt_dosp, time_average_sp, time_dosp in module steuer,
+! +spectrum_x, spectrum_y in module statistik
+!
+! Revision 1.48  2001/01/02 17:31:36  raasch
+! openfile has the new type file_status, version updated to 2.1a
+!
+! Revision 1.47  2000/12/30 18:48:53  raasch
+! Default value of mode_dvrp changed
+!
+! Revision 1.46  2000/12/28 13:21:35  raasch
+! version updated to 2.1
+!
+! Revision 1.45  2000/12/28 13:19:35  raasch
+! +mode_dvrp, +interval_..._dvrp, REAL*4 changed to REAL in particle_attributes,
+! +dt_prel, time_prel, -particle_transport
+!
+! Revision 1.44  2000/07/03 13:02:32  raasch
+! version updated to 2.0b
+!
+! Revision 1.43  2000/07/03 12:59:53  raasch
+! +dvrp_directory, dvrp_host, dvrp_password, dvrp_username,
+! new module pointer_interfaces which includes interfaces for subroutines with
+! pointer arguments called in prognostic_equations
+!
+! Revision 1.42  2000/04/27 07:35:56  raasch
+! npex and npey must be defined in general (not only for parallel case)
+!
+! Revision 1.41  2000/04/27 07:07:14  raasch
+! +pt_slope_offset, pt_slope_ref, npex, npey, dt_dvrp, time_dvrp,
+! new module dvrp_variables
+! -dt_plisos, dt_plpart, isosurface_filecount, time_plisos, time_plpart
+! Version 2.0a, all comments translated into English, old revision remarks
+! deleted
+!
+! Revision 1.40  2000/04/13 13:17:25  schroeter
+! Felder fuer Wolkenphysik zu Modul arrays3d hinzugefuegt, neues Modul
+! (cloud_parameters) fuer Wolkenphysikparameter eingefuegt
+!
+! Revision 1.39  2000/01/10  10:05:11  10:05:11  raasch (Siegfried Raasch)
+! +use_ug_for_galilei_tr, Aufteilung der zeitabhaengigen Groessen in pointer
+! und targets, + timestep_count
+!
+! Revision 1.38  1999/12/28 10:06:52  raasch
+! +Variablen fuer Plot von Isooberflaechen
+!
+! Revision 1.37  1999/11/25 16:26:38  raasch
+! Neues Modul particle_attributes,
+! +Variablen fuer Partikelverlagerung: particle_transport, dt_plpart,
+! Version 1.1e
+!
 ! Revision 1.1  1997/07/24 11:21:26  raasch

/palm/trunk/SOURCE/netcdf.f90

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

/palm/trunk/SOURCE/package_parin.f90

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

/palm/trunk/SOURCE/palm.f90

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

/palm/trunk/SOURCE/parin.f90

Property svn:keywords deleted

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! +netcdf_64bit_3d in d3par,
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: parin.f90,v $
 ! Revision 1.57  2007/02/11 13:11:22  raasch
 ! Values of environment variables are now read from file ENVPAR instead of
 ! reading them with a system call, + NAMELIST envpar
+!
+! Revision 1.56  2006/08/22 14:13:56  raasch
+! +dz_max in inipar
+!
+! Revision 1.55  2006/08/04 14:53:43  raasch
+! +use_upstream_for_tke
+!
+! Revision 1.54  2006/03/17 08:45:44  raasch
+! +skip_time_data_output_av
+!
+! Revision 1.53  2006/02/23 12:46:19  raasch
+! parameters for 2d/3d averaging,
+! +conserve_volume_flow, dissipation_1d, e_min, mixing_length_1d, topography, &
+! wall_heatflux in inipar,
+! +data_output, netcdf_64bit, skip_time_dosp, skip_time_dopr, etc. in d3par,
+! allocation of nanz_2dh and nanz_3d moved to init_3d_model, ebene renamed
+! section, average_period_pr1d removed, average_period_pl1d renamed
+! averaging_interval_pr, dt_average renamed dt_averaging_input_pr,
+! +parameters for single building topography settings
+!
+! Revision 1.52  2005/06/29 10:21:15  steinfeld
+! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
+!  vg_surface, vg_vertical_gradient, vg_vertical_gradient_level in inipar;
+! arrays ug and vg are allocated during initialization runs
+!
+! Revision 1.51  2005/06/26 20:06:26  raasch
+! +cloud_droplets in inipar
+!
+! Revision 1.50  2005/05/18 15:40:49  raasch
+! +netcdf_precision in inipar, +data_output_format in d3par
+!
+! Revision 1.49  2005/04/23 09:35:23  raasch
+! fcl_factor renamed cfl_factor
+!
+! Revision 1.48  2005/03/26 20:42:17  raasch
+! +bc_lr, bc_ns, call_psolver_at_all_substeps, inflow_disturbance_begin, &
+! inflow_disturbance_end, km_damp_max, outflow_damping_width in inipar
+!
+! Revision 1.47  2004/04/30 12:32:06  raasch
+! +grid_matching in inipar, impulse_advec renamed momentum_advec
+!
+! Revision 1.46  2003/10/29 09:03:06  raasch
+! +mg_cycles, mg_switch_to_pe0_level in d3par
+!
+! Revision 1.45  2003/03/16 09:42:23  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.44  2003/03/14 13:45:24  raasch
+! +random_generator in inipar
+!
+! Revision 1.43  2003/03/12 16:35:55  raasch
+! Reading of environment variable tasks_per_node moved from poisfft_hybrid
+! to here
+!
+! Revision 1.42  2002/12/19 15:55:31  raasch
+! Saving preliminary process-id-string for opening unit 14 (binary output for
+! restart). Array hom is allocated during first run only, allocation during
+! restart runs will now be done in read_var_list.
+! +dt_restart, restart_time in d3par, STOP statement replaced by call of
+! subroutine local_stop
+!
+! Revision 1.41  2002/04/16 08:10:50  raasch
+! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,
+! s_vertical_gradient, s_vertical_gradient_level in inipar
+!
+! Revision 1.40  2001/08/21 09:54:27  raasch
+! MPI_COMM_WORLD replaced by comm_palm, +wall_adjustment in inipar
+!
+! Revision 1.39  2001/07/20 13:10:07  raasch
+! +cycle_mg, ngsrb, residual_limit in d3par
+!
+! Revision 1.38  2001/03/30 07:43:55  raasch
+! Translation of remaining German identifiers (variables, subroutines, etc.),
+! namelist d3par can also be omitted in the initial run
+!
+! Revision 1.37  2001/01/29 12:30:26  raasch
+! +passive_scalar in inipar
+!
+! Revision 1.36  2001/01/25 07:17:39  raasch
+! +fft_method, use_surface_fluxes in inipar
+!
+! Revision 1.35  2001/01/02 17:32:13  raasch
+! Unit 11 closed at the end of the subroutine
+!
+! Revision 1.34  2000/12/28 13:24:51  raasch
+! Call of package_parin. dt_dvrp and threshold moved
+! from d3par to other namelists in subroutine package_parin
+!
+! Revision 1.33  2000/04/27 06:55:20  raasch
+! +npex, npey in initial parameter list, comments have to be translated yet,
+! -dt_plisos, dt_plpart, +dt_dvrp, threshold, in run-time parameter list,
+! all comments translated into English, old revision remarks deleted
+!
+! Revision 1.32  2000/04/13 13:03:06  schroeter
+! Erweitern der Parameterlist fuer Wolkenphysik
+!
+! Revision 1.31  2000/01/10  10:05:55  10:05:55  raasch (Siegfried Raasch)
+! +use_ug_for_galilei_tr
+!
 ! Revision 1.1  1997/07/24 11:22:50  raasch
 ! Initial revision
 …
                        dt_do2d_yz, dt_do3d, dt_restart, dt_run_control, &
                        end_time, force_print_header, mg_cycles, &
                        mg_switch_to_pe0_level, netcdf_64bit, netcdf_64bit_3d, &
                        ngsrb, normalizing_region, nsor, nz_do3d, omega_sor, &
+                       mg_switch_to_pe0_level, netcdf_64bit, ngsrb, &
+                       normalizing_region, nsor, nz_do3d, omega_sor, &
                        prandtl_number, profile_columns, profile_rows, psolver, &
                        rayleigh_damping_factor, rayleigh_damping_height, &

/palm/trunk/SOURCE/poisfft.f90

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

/palm/trunk/SOURCE/poisfft_hybrid.f90

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

/palm/trunk/SOURCE/poismg.f90

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

/palm/trunk/SOURCE/prandtl_fluxes.f90

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

/palm/trunk/SOURCE/pres.f90

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

/palm/trunk/SOURCE/print_1d.f90

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

/palm/trunk/SOURCE/production_e.f90

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

/palm/trunk/SOURCE/prognostic_equations.f90

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

/palm/trunk/SOURCE/random_function.f90

-                      r3
+                      r1
  MODULE random_function_mod
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: random_function.f90,v $
 ! Revision 1.3  2003/10/29 09:06:57  raasch
 ! Former function changed to a module.
+!
+! Revision 1.2  2001/01/22 08:00:23  raasch
+! Comments translated into English
+!
 ! Revision 1.1  1998/02/04 16:09:45  raasch
 …
 ! Random number generator, produces numbers equally distributed in interval [0,1]
 ! This routine is taken from the "numerical recipies"
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     IMPLICIT NONE

/palm/trunk/SOURCE/random_gauss.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: random_gauss.f90,v $
 ! Revision 1.4  2006/08/04 15:01:48  raasch
 ! Range of random number is limited by an upper limit (new second parameter)
+!
+! Revision 1.3  2003/10/29 09:08:04  raasch
+! Module random_function_mod is used
+!
+! Revision 1.2  2001/01/22 08:03:31  raasch
+! Comments translated into English
+!
 ! Revision 1.1  1998/03/25 20:09:47  raasch

/palm/trunk/SOURCE/read_3d_binary.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: read_3d_binary.f90,v $
 ! Revision 1.4  2006/08/04 15:02:32  raasch
 ! +iran, iran_part
+!
+! Revision 1.3  2006/02/23 12:52:45  raasch
+! +volume_flow_area, volume_flow_initial, -hom_sum_pr,
+! reading arrays used for time-averaging
+!
+! Revision 1.2  2005/05/18 15:54:33  raasch
+! pl1d_callcount renamed dopr_time_count
+!
 ! Revision 1.1  2004/04/30 12:47:27  raasch

/palm/trunk/SOURCE/read_var_list.f90

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

/palm/trunk/SOURCE/run_control.f90

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

/palm/trunk/SOURCE/singleton.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: singleton.f90,v $
 ! Revision 1.2  2004/04/30 12:52:09  raasch
 ! Shape of arrays is explicitly stored in ishape and handled to the

/palm/trunk/SOURCE/sor.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: sor.f90,v $
 ! Revision 1.9  2005/03/26 21:02:23  raasch
 ! Implementation of non-cyclic (Neumann) horizontal boundary conditions,
 ! dx2,dy2 replaced by ddx2,ddy2
+!
+! Revision 1.8  2001/03/30 07:52:44  raasch
+! Arrays f1-f3 changed from dummy arguments to allocatable arrays,
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.7  2001/01/30 21:41:25  letzel
+! All comments translated into English.
+!
+! Revision 1.6  2001/01/22 08:06:14  raasch
+! Module test_variables removed
+!
+! Revision 1.5  1998/07/06 12:33:13  raasch
+! + USE test_variables
+!
+! Revision 1.4  1997/09/12 06:30:12  raasch
+! Randbedingungen umbenannt
+!
+! Revision 1.3  1997/09/09 08:30:44  raasch
+! Kehrwerte der Gitterweiten implementiert
+!
+! Revision 1.2  1997/08/29 09:00:58  raasch
+! omega --> omega_sor
+!
 ! Revision 1.1  1997/08/11 06:25:56  raasch
 ! Initial revision
 …
 ! ------------
 ! Solve the Poisson-equation with the SOR-Red/Black-scheme.
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     USE grid_variables

/palm/trunk/SOURCE/spline_x.f90

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

/palm/trunk/SOURCE/spline_y.f90

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

/palm/trunk/SOURCE/spline_z.f90

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

/palm/trunk/SOURCE/sum_up_3d_data.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: sum_up_3d_data.f90,v $
 ! Revision 1.1  2006/02/23 12:55:23  raasch
 ! Initial revision

/palm/trunk/SOURCE/swap_timelevel.f90

-                      r3
+                      r1
  SUBROUTINE swap_timelevel
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: swap_timelevel.f90,v $
 ! Revision 1.8  2004/01/28 15:28:18  raasch
 ! Swaping for Runge-Kutta schemes implemented
+!
+! Revision 1.7  2002/05/02 18:58:36  raasch
+! Swaping of prognostic quantities re-changed to three timelevels
+!
+! Revision 1.6  2001/08/21  10:01:58  10:01:58  raasch (Siegfried Raasch)
+! Surface momentum fluxes are swapped also when they are prescribed
+!
+! Revision 1.5  2001/03/30 07:54:39  raasch
+! Swaping of prognostic quantities changed over from three to two timelevels,
+! Translation of remaining German identifiers (variables, subroutines, etc.)
+!
+! Revision 1.4  2001/01/29 12:36:03  raasch
+! Passive scalar is considered
+!
+! Revision 1.3  2001/01/25 07:26:48  raasch
+! Sensible and latent heat fluxes are swaped, when use_surface_fluxes = .TRUE.,
+! module test_variables removed
+!
+! Revision 1.2  2000/04/13 13:59:26  schroeter
+! + swapping timelevels of moisture variables
+!
 ! Revision 1.1  2000/01/10  10:08:58  10:08:58  raasch (Siegfried Raasch)
 …
 ! ------------
 ! Swap of timelevels of variables after each timestep
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     USE arrays_3d

/palm/trunk/SOURCE/temperton_fft.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: temperton_fft.f90,v $
 ! Revision 1.2  2003/04/16 12:49:25  raasch
 ! Abort in case of illegal factors

/palm/trunk/SOURCE/time_integration.f90

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

/palm/trunk/SOURCE/time_to_string.f90

-                      r3
+                      r1
  FUNCTION time_to_string( time )
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
 ! Actual revisions:
 ! -----------------
 …
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: time_to_string.f90,v $
 ! Revision 1.3  2001/01/22 08:16:04  raasch
 ! Comments translated into English
+!
+! Revision 1.2  1998/03/09 16:24:53  raasch
+! Dreistellige Stundenausgabe ermoeglicht
+!
 ! Revision 1.1  1997/08/11 06:26:08  raasch
 …
 ! ------------
 ! Transforming the time from real to character-string hh:mm:ss
 !------------------------------------------------------------------------------!
+!-------------------------------------------------------------------------------!
     IMPLICIT NONE

/palm/trunk/SOURCE/timestep.f90

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

/palm/trunk/SOURCE/timestep_scheme_steering.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: timestep_scheme_steering.f90,v $
 ! Revision 1.2  2005/03/26 21:17:06  raasch
 ! No pressure term for Runge-Kutta-schemes (tsc(4)=0.0)

/palm/trunk/SOURCE/transpose.f90

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

/palm/trunk/SOURCE/user_interface.f90

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

/palm/trunk/SOURCE/write_3d_binary.f90

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

/palm/trunk/SOURCE/write_compressed.f90

-                      r3
+                      r1
 ! Former revisions:
 ! ---------------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: write_compressed.f90,v $
 ! Revision 1.4  2006/02/23 13:15:09  raasch
 ! nz_plot3d renamed nz_do3d
+!
+! Revision 1.3  2003/03/16 09:51:07  raasch
+! Two underscores (_) are placed in front of all define-strings
+!
+! Revision 1.2  2001/03/30 07:59:48  raasch
+! All comments translated into English
+!
 ! Revision 1.1  1999/03/02 09:25:21  raasch

/palm/trunk/SOURCE/write_var_list.f90

-                      r3
+                      r1
 ! Former revisions:
 ! -----------------
+! $Id$
+! RCS Log replace by Id keyword, revision history cleaned up
+!
+! $Log: write_var_list.f90,v $
 ! Revision 1.34  2006/08/22 14:30:52  raasch
 ! +dz_max
+!
+! Revision 1.33  2006/08/04 15:05:52  raasch
+! +first_call_advec_particles, time_dopts, use_upstream_for_tke
+!
+! Revision 1.32  2006/02/23 13:16:02  raasch
+! +building parameters, conserve_volume_flow, dissipation_1d, e_min,
+! mixing_length_1d, topography,
+! -average_count_pl, time_average_pr, time_average_pl renamed time_dopr_av,
+! nt_anz renamed current_timestep_number, ..pl.. renamed ..do..
+!
+! Revision 1.31  2005/06/28 14:18:02  steinfeld
+! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,
+! ug_vertical_gradient_level_ind, vg_surface, vg_vertical_gradient,
+! vg_vertical_gradient_level, vg_vertical_gradient_level_ind
+!
+! Revision 1.30  2005/06/26 20:10:54  raasch
+! +cloud_droplets
+!
+! Revision 1.29  2005/05/18 15:56:20  raasch
+! +netcdf_precision
+!
+! Revision 1.28  2005/03/26 21:18:39  raasch
+! +bc_lr, bc_ns, inflow_disturbance_begin, inflow_disturbance_end, km_damp_max,
+! outflow_damping_width
+!
+! Revision 1.27  2004/04/30 13:14:31  raasch
+! +grid_matching, impulse_advec renamed momentum_advec
+!
+! Revision 1.26  2004/01/28 15:38:55  raasch
+! at and bt replaced by array tsc
+!
+! Revision 1.25  2003/10/29 09:15:12  raasch
+! +last_dt_change
+!
+! Revision 1.24  2003/03/14 13:47:55  raasch
+! +random_generator
+!
+! Revision 1.23  2002/12/19 16:24:09  raasch
+! Writing of array hom moved from write_3d_binary to this routine,
+! binary version incremented to 2.2, +time_restart
+!
+! Revision 1.22  2002/04/16 08:22:02  raasch
+! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,
+! s_vertical_gradient, s_vertical_gradient_level
+!
+! Revision 1.21  2001/08/21 10:05:05  raasch
+! +wall_adjustment
+!
+! Revision 1.20  2001/03/30 08:00:07  raasch
+! Translation of remaining German identifiers (variables, subroutines, etc.),
+! runtime and package parameters are not written on unit 14 any more
+!
+! Revision 1.19  2001/01/29 12:38:42  raasch
+! Error writing use_surface_fluxes-string corrected, +passive_scalar
+!
+! Revision 1.18  2001/01/25 07:27:48  raasch
+! +fft_method, use_surface_fluxes
+!
+! Revision 1.17  2001/01/05 15:16:59  raasch
+! +average_count_sp, time_dosp
+!
+! Revision 1.16  2000/12/28 13:39:21  raasch
+! Module dvrp_variables removed, -particle_transport
+!
+! Revision 1.15  2000/04/27 06:58:19  raasch
+! +dt_dvrp, dvrp_filecount, npex, npey, threshold, time_dvrp
+! -dt_plisos, dt_plpart, time_plisos, time_plpart
+! +module pegrid,
+! all comments translated into English, old revision remarks deleted
+!
+! Revision 1.14  2000/04/13 14:10:58  schroeter
+! Upgrade for calculations with cloud physics
+!
+! Revision 1.13  2000/01/10  10:11:06  10:11:06  raasch (Siegfried Raasch)
+! +use_ug_for_galilei_tr
+!
 ! Revision 1.1  1998/03/18 20:20:38  raasch

Context Navigation

Changes in / [3:1]

Legend:

Download in other formats: