Changeset 3 for palm/trunk/SOURCE
- Timestamp:
- Feb 13, 2007 11:30:58 AM (18 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 106 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/Makefile
r1 r3 8 8 # Former revisions: 9 9 # ----------------- 10 # $Log: Makefile,v $ 10 # $Id $ 11 # RCS Log replace by Id keyword, revision history cleaned up 12 # 11 13 # Revision 1.13 2006/08/04 14:45:31 raasch 12 14 # +data_output_ptseries 13 #14 # Revision 1.12 2006/04/26 15:10:49 raasch15 # modules.o added to dependencies of user_interface.o16 #17 # Revision 1.11 2006/02/23 13:17:11 raasch18 # +average_3d_data, data_log, sum_up_3d_data,19 # plot_.. renamed data_output_..20 #21 # Revision 1.10 2005/06/26 20:00:15 raasch22 # +interaction_droplets_ptq23 #24 # Revision 1.9 2005/05/18 14:52:20 raasch25 # +netcdf26 #27 # Revision 1.8 2004/04/30 13:15:46 raasch28 # dependencies for calc_spectra changed, +read_3d_binary,29 # leap_frog changed to time_integration, fft_for_1d_decomp to fft_xy,30 # init_pegrid and pres now depend on poisfft, all transpose routines joined31 # to one file transpose, local_getenv depends on modules32 #33 # Revision 1.7 2004/01/28 14:58:55 raasch34 # +timestep_scheme_steering35 #36 # Revision 1.6 2003/10/29 08:21:40 raasch37 # advec_particles, disturb_field, init_particles, init_3d_model, random_gauss38 # and write_3d_binary are using module random_function_mod39 #40 # Revision 1.5 2003/05/09 14:04:39 raasch41 # poisfft depends on fft_for_1d_decomp and no more on singleton,42 # init_pegrid depends on fft_for_1d_decomp, temperton_fft depends on modules43 #44 # Revision 1.4 2003/03/13 08:46:52 raasch45 # fft_for_1d_decomp has to depend on temperton_fft46 #47 # Revision 1.3 2003/03/12 15:57:37 raasch48 # modules and user_interface now depend on the f90 files, dependency of49 # singleton added, LDFLAGS moved to the end of PROG rule, +temperton_fft50 #51 # Revision 1.2 2003/01/06 08:48:19 raasch52 # +local_stop, check_cpu_time renamed to check_for_restart53 15 # 54 16 # Revision 1.1 2002/06/11 13:42:20 raasch -
palm/trunk/SOURCE/advec_particles.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_particles.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.32 2007/02/11 12:48:20 raasch 12 14 ! Allways the lower level k is used for interpolation -
palm/trunk/SOURCE/advec_s_bc.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_s_bc.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.22 2006/02/23 09:42:08 raasch 12 14 ! anz renamed ngp 13 !14 ! Revision 1.21 2003/03/16 09:22:13 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.20 2002/12/19 13:41:50 raasch18 ! STOP statement replaced by call of subroutine local_stop19 !20 ! Revision 1.19 2001/03/29 17:29:52 raasch21 ! Result is stored as tendency on tend, sk_p is now a local array, rho is22 ! replaced by d23 ! Translation of remaining German identifiers (variables, subroutines, etc.)24 !25 ! Revision 1.18 2001/01/22 04:25:52 raasch26 ! Module test_variables removed27 !28 ! Revision 1.17 2000/07/03 12:54:35 raasch29 ! Arguments rho, sk and sk_p now declared as pointers30 !31 ! Revision 1.16 2000/04/27 07:03:07 raasch32 ! Treatment at cyclic boundaries is simplified by adding one additional33 ! ghost layer to tend at each cyclic boundary. Result is that most of34 ! mpi_sendrecv-calls (for polynom-coefficients, switches and fluxes) are no35 ! more necessary so that much cpu-time is spared.36 ! Adjustment of cyclic boundaries in x-direction in case of a sloping surface37 ! (temperature only).38 ! Vertical boundary conditions for humidity added.39 !40 ! Revision 1.15 2000/01/20 09:07:01 letzel41 ! All comments translated into English42 !43 ! Revision 1.14 1999/02/17 09:01:38 raasch44 ! Verfahren erweitert fuer TKE-Advektion45 ! Dichteberechnung (Korrekturverfahren) korrigiert (Galilei-Abzug war falsch)46 !47 ! Revision 1.10 1998/04/21 15:51:34 raasch48 ! Galilei-Transformation eingebaut49 !50 ! Revision 1.4 1997/09/08 10:42:35 raasch51 ! Optimierungen im Exponentialteil, MIN/MAX-Bildung aus progn. Gleichung52 ! entfernt, sw und m1 auf t3e einfach genau definiert, ffmax aus progn.53 ! Gleichung entfernt (Skalarfeld wird glatter im gradientschwachen Bereich)54 ! Geschwindigkeitsgewinn: ca 32%55 15 ! 56 16 ! Revision 1.1 1997/08/29 08:53:46 raasch -
palm/trunk/SOURCE/advec_s_pw.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_s_pw.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.12 2006/02/23 09:42:55 raasch 12 14 ! nzb_2d replaced by nzb_s_inner 13 !14 ! Revision 1.11 2004/01/30 10:10:03 raasch15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d16 !17 ! Revision 1.10 2003/03/12 16:15:32 raasch18 ! Full code replaced in the call for all gridpoints instead of calling the19 ! _ij version (required by NEC, because otherwise no vectorization)20 !21 ! Revision 1.9 2002/06/11 12:22:14 raasch22 ! Former subroutine changed to a module which allows to be called for all grid23 ! points of a single vertical column with index i,j or for all grid points by24 ! using function overloading.25 ! Array sk now declared as a pointer.26 !27 ! Revision 1.8 2001/03/29 17:30:14 raasch28 ! Translation of remaining German identifiers (variables, subroutines, etc.)29 !30 ! Revision 1.7 2001/01/22 04:41:39 raasch31 ! Module test_variables removed32 !33 ! Revision 1.6 2000/01/20 09:13:33 letzel34 ! All comments translated into English35 !36 ! Revision 1.5 1998/07/06 12:03:09 raasch37 ! + USE test_variables38 !39 ! Revision 1.4 1998/04/21 15:51:53 raasch40 ! Galilei-Transformation eingebaut41 !42 ! Revision 1.3 1998/02/10 14:57:53 raasch43 ! Es wird jetzt grundsaetzlich die Form C3 verwendet44 !45 ! Revision 1.2 1997/09/09 08:27:10 raasch46 ! Kehrwerte der Gitterweiten implementiert47 15 ! 48 16 ! Revision 1.1 1997/08/29 08:54:20 raasch -
palm/trunk/SOURCE/advec_s_up.f90
r1 r3 8 8 ! Forner revisions: 9 9 ! ----------------- 10 ! $Log: advec_s_up.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2006/02/23 09:43:44 raasch 12 14 ! nzb_2d replaced by nzb_s_inner 13 !14 ! Revision 1.10 2004/01/30 10:10:39 raasch15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d16 !17 ! Revision 1.9 2003/03/12 16:16:25 raasch18 ! Full code replaced in the call for all gridpoints instead of calling the19 ! _ij version (required by NEC, because otherwise no vectorization)20 !21 ! Revision 1.8 2002/06/11 12:23:31 raasch22 ! Former subroutine changed to a module which allows to be called for all grid23 ! points of a single vertical column with index i,j or for all grid points by24 ! using function overloading.25 ! Array sk now declared as a pointer.26 !27 ! Revision 1.7 2001/03/29 17:30:31 raasch28 ! Translation of remaining German identifiers (variables, subroutines, etc.)29 !30 ! Revision 1.6 2001/01/22 04:43:56 raasch31 ! Module test_variables removed32 !33 ! Revision 1.5 2000/01/20 09:16:59 letzel34 ! All comments translated into English35 !36 ! Revision 1.4 1998/07/06 12:03:46 raasch37 ! + USE test_variables38 !39 ! Revision 1.3 1998/04/21 15:52:14 raasch40 ! Galilei-Transformation eingebaut41 !42 ! Revision 1.2 1997/09/09 08:27:24 raasch43 ! Kehrwerte der Gitterweiten implementiert44 15 ! 45 16 ! Revision 1.1 1997/08/29 08:54:33 raasch -
palm/trunk/SOURCE/advec_s_ups.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_s_ups.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.6 2004/04/30 08:02:43 raasch 12 14 ! Enlarged transposition arrays introduced 13 !14 ! Revision 1.5 2003/03/16 09:23:56 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.4 2001/03/29 17:34:05 raasch18 ! v_ad is now allocated locally, several temporary arrays removed from19 ! argument lists of spline_* and long_filter_*20 ! Translation of remaining German identifiers (variables, subroutines, etc.)21 !22 ! Revision 1.3 2001/01/22 04:46:10 raasch23 ! Module test_variables removed24 !25 ! Revision 1.2 2000/01/20 09:28:52 letzel26 ! All comments translated into English27 15 ! 28 16 ! Revision 1.1 1999/02/05 08:44:47 raasch -
palm/trunk/SOURCE/advec_u_pw.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_u_pw.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 09:44:21 raasch 12 14 ! nzb_2d replaced by nzb_u_inner 13 !14 ! Revision 1.14 2005/03/26 14:49:31 raasch15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions16 !17 ! Revision 1.13 2004/01/30 10:11:21 raasch18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d19 !20 ! Revision 1.12 2003/03/12 16:17:06 raasch21 ! Full code replaced in the call for all gridpoints instead of calling the22 ! _ij version (required by NEC, because otherwise no vectorization)23 !24 ! Revision 1.11 2002/06/11 12:25:00 raasch25 ! Former subroutine changed to a module which allows to be called for all grid26 ! points of a single vertical column with index i,j or for all grid points by27 ! using function overloading.28 !29 ! Revision 1.10 2001/03/29 17:34:21 raasch30 ! Translation of remaining German identifiers (variables, subroutines, etc.)31 !32 ! Revision 1.9 2001/01/22 04:55:36 raasch33 ! Module test_variables removed34 !35 ! Revision 1.8 2000/01/20 09:34:26 letzel36 ! All comments translated into English37 !38 ! Revision 1.7 1998/08/05 06:47:04 raasch39 ! Alle Gitterpunkte gleich behandelt (keine zentralen Differenzenquotienten40 ! mehr bei Vertikaladvektion bei k=1)41 15 ! 42 16 ! Revision 1.1 1997/08/11 06:09:21 raasch -
palm/trunk/SOURCE/advec_u_up.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_u_up.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.12 2006/02/23 09:45:04 raasch 12 14 ! nzb_2d replaced by nzb_u_inner 13 !14 ! Revision 1.11 2005/03/26 14:50:43 raasch15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions16 !17 ! Revision 1.10 2004/01/30 10:11:47 raasch18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d19 !20 ! Revision 1.9 2003/03/12 16:17:42 raasch21 ! Full code replaced in the call for all gridpoints instead of calling the22 ! _ij version (required by NEC, because otherwise no vectorization)23 !24 ! Revision 1.8 2002/06/11 12:25:43 raasch25 ! Former subroutine changed to a module which allows to be called for all grid26 ! points of a single vertical column with index i,j or for all grid points by27 ! using function overloading.28 !29 ! Revision 1.7 2001/03/29 17:35:08 raasch30 ! Translation of remaining German identifiers (variables, subroutines, etc.)31 !32 ! Revision 1.6 2001/01/22 04:57:23 raasch33 ! Module test_variables removed34 !35 ! Revision 1.5 2000/01/20 09:40:00 letzel36 ! All comments translated into English37 !38 ! Revision 1.4 1998/07/06 12:04:43 raasch39 ! + USE test_variables40 !41 ! Revision 1.3 1998/04/21 15:52:51 raasch42 ! Galilei-Transformation eingebaut43 !44 ! Revision 1.2 1997/09/09 08:27:49 raasch45 ! Kehrwerte der Gitterweiten implementiert46 15 ! 47 16 ! Revision 1.1 1997/08/29 08:55:25 raasch -
palm/trunk/SOURCE/advec_u_ups.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_u_ups.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.6 2004/04/30 08:03:19 raasch 12 14 ! Enlarged transposition arrays introduced 13 !14 ! Revision 1.5 2003/03/16 09:25:07 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.4 2001/03/29 17:35:32 raasch18 ! v_ad is now allocated locally, several temporary arrays removed from19 ! argument lists of spline_* and long_filter_*20 ! Translation of remaining German identifiers (variables, subroutines, etc.)21 !22 ! Revision 1.3 2001/01/22 04:59:29 raasch23 ! Module test_variables removed24 !25 ! Revision 1.2 2000/01/20 09:46:52 letzel26 ! All comments translated into English27 15 ! 28 16 ! Revision 1.1 1999/02/05 08:49:08 raasch -
palm/trunk/SOURCE/advec_v_pw.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_v_pw.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 09:46:08 raasch 12 14 ! nzb_2d replaced by nzb_v_inner 13 !14 ! Revision 1.14 2005/03/26 14:52:51 raasch15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions16 !17 ! Revision 1.13 2004/01/30 10:12:25 raasch18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d19 !20 ! Revision 1.12 2003/03/12 16:18:19 raasch21 ! Full code replaced in the call for all gridpoints instead of calling the22 ! _ij version (required by NEC, because otherwise no vectorization)23 !24 ! Revision 1.11 2002/06/11 12:26:18 raasch25 ! Former subroutine changed to a module which allows to be called for all grid26 ! points of a single vertical column with index i,j or for all grid points by27 ! using function overloading.28 !29 ! Revision 1.10 2001/03/29 17:36:06 raasch30 ! Translation of remaining German identifiers (variables, subroutines, etc.)31 !32 ! Revision 1.9 2001/01/22 05:14:48 raasch33 ! Module test_variables removed34 !35 ! Revision 1.8 2000/01/20 10:16:37 letzel36 ! All comments translated into English37 !38 ! Revision 1.7 1998/08/05 06:49:43 raasch39 ! Alle Gitterpunkte gleich behandelt (keine zentralen Differenzenquotienten mehr40 ! bei Vertikaladvektion bei k=1)41 15 ! 42 16 ! Revision 1.1 1997/08/11 06:09:57 raasch -
palm/trunk/SOURCE/advec_v_up.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_v_up.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.12 2006/02/23 09:46:37 raasch 12 14 ! nzb_2d replaced by nzb_v_inner 13 !14 ! Revision 1.11 2005/03/26 14:55:59 raasch15 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions16 !17 ! Revision 1.10 2004/01/30 10:12:45 raasch18 ! Scalar lower k index nzb replaced by 2d-array nzb_2d19 !20 ! Revision 1.9 2003/03/12 16:19:00 raasch21 ! Full code replaced in the call for all gridpoints instead of calling the22 ! _ij version (required by NEC, because otherwise no vectorization)23 !24 ! Revision 1.8 2002/06/11 12:26:48 raasch25 ! Former subroutine changed to a module which allows to be called for all grid26 ! points of a single vertical column with index i,j or for all grid points by27 ! using function overloading.28 !29 ! Revision 1.7 2001/03/29 17:36:22 raasch30 ! Translation of remaining German identifiers (variables, subroutines, etc.)31 !32 ! Revision 1.6 2001/01/22 05:16:19 raasch33 ! Module test_variables removed34 !35 ! Revision 1.5 2000/01/20 10:20:31 letzel36 ! All comments translated into English37 !38 ! Revision 1.4 1998/07/06 12:05:37 raasch39 ! + USE test_variables40 !41 ! Revision 1.3 1998/04/21 15:53:36 raasch42 ! Galilei-Transformation eingebaut43 !44 ! Revision 1.2 1997/09/09 08:28:21 raasch45 ! Kehrwerte der Gitterweiten implementiert46 15 ! 47 16 ! Revision 1.1 1997/08/29 08:56:05 raasch -
palm/trunk/SOURCE/advec_v_ups.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_v_ups.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.7 2004/04/30 08:03:52 raasch 12 14 ! Enlarged transposition arrays introduced 13 !14 ! Revision 1.6 2003/03/16 09:25:59 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.5 2001/03/29 17:36:43 raasch18 ! v_ad is now allocated locally, several temporary arrays removed from19 ! argument lists of spline_* and long_filter_*20 ! Translation of remaining German identifiers (variables, subroutines, etc.)21 !22 ! Revision 1.4 2001/01/22 08:31:10 raasch23 ! Module test_variables removed24 !25 !26 ! Revision 1.2 2000/01/20 09:52:54 letzel27 ! All comments translated into English28 15 ! 29 16 ! Revision 1.1 1999/02/05 08:50:32 raasch -
palm/trunk/SOURCE/advec_w_pw.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_w_pw.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 09:47:01 raasch 12 14 ! nzb_2d replaced by nzb_w_inner 13 !14 ! Revision 1.14 2004/01/30 10:13:07 raasch15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d16 !17 ! Revision 1.13 2003/03/12 16:19:29 raasch18 ! Full code replaced in the call for all gridpoints instead of calling the19 ! _ij version (required by NEC, because otherwise no vectorization)20 !21 ! Revision 1.12 2002/06/11 12:27:33 raasch22 ! Former subroutine changed to a module which allows to be called for all grid23 ! points of a single vertical column with index i,j or for all grid points by24 ! using function overloading.25 !26 ! Revision 1.11 2001/03/29 17:37:10 raasch27 ! Translation of remaining German identifiers (variables, subroutines, etc.)28 !29 ! Revision 1.10 2001/01/22 05:19:42 raasch30 ! Module test_variables removed31 !32 ! Revision 1.9 2000/01/20 10:18:31 letzel33 ! All comments translated into English34 !35 ! Revision 1.8 1998/08/05 06:51:10 raasch36 ! Alle Gitterpunkte werden gleich behandelt (keine zentralen Differenzenquotien-37 ! ten mehr bei Vertikaladvektion bei k=1)38 15 ! 39 16 ! Revision 1.1 1997/08/11 06:10:29 raasch -
palm/trunk/SOURCE/advec_w_up.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_w_up.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2006/02/23 09:47:23 raasch 12 14 ! *** empty log message *** 13 !14 ! Revision 1.10 2004/01/30 10:13:31 raasch15 ! Scalar lower k index nzb replaced by 2d-array nzb_2d16 !17 ! Revision 1.9 2003/03/12 16:20:07 raasch18 ! Full code replaced in the call for all gridpoints instead of calling the19 ! _ij version (required by NEC, because otherwise no vectorization)20 !21 ! Revision 1.8 2002/06/11 12:28:33 raasch22 ! Former subroutine changed to a module which allows to be called for all grid23 ! points of a single vertical column with index i,j or for all grid points by24 ! using function overloading.25 !26 ! Revision 1.7 2001/03/29 17:37:31 raasch27 ! Translation of remaining German identifiers (variables, subroutines, etc.)28 !29 ! Revision 1.6 2001/01/22 05:21:01 raasch30 ! Module test_variables removed31 !32 ! Revision 1.5 2000/01/20 10:22:08 letzel33 ! All comments translated into English34 !35 ! Revision 1.4 1998/07/06 12:06:25 raasch36 ! + USE test_variables37 !38 ! Revision 1.3 1998/04/21 15:54:20 raasch39 ! Galilei-Transformation eingebaut40 !41 ! Revision 1.2 1997/09/09 08:28:55 raasch42 ! Kehrwerte der Gitterweiten implementiert43 15 ! 44 16 ! Revision 1.1 1997/08/29 08:56:33 raasch -
palm/trunk/SOURCE/advec_w_ups.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: advec_w_ups.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.6 2004/04/30 08:05:05 raasch 12 14 ! Enlarged transposition arrays introduced 13 !14 ! Revision 1.5 2003/03/16 09:26:46 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.4 2001/03/29 17:37:53 raasch18 ! v_ad is now allocated locally, several temporary arrays removed from19 ! argument lists of spline_* and long_filter_*20 ! Translation of remaining German identifiers (variables, subroutines, etc.)21 !22 ! Revision 1.3 2001/01/22 05:23:31 raasch23 ! Module test_variables removed24 !25 ! Revision 1.2 2000/01/20 09:57:15 letzel26 ! All comments translated into English27 15 ! 28 16 ! Revision 1.1 1999/02/05 08:52:09 raasch -
palm/trunk/SOURCE/asselin_filter.f90
r1 r3 1 1 SUBROUTINE asselin_filter 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! --------------------- 10 ! $Log: asselin_filter.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.8 2004/01/30 10:14:02 raasch 12 14 ! Scalar lower k index nzb replaced by 2d-array nzb_2d 13 !14 ! Revision 1.7 2003/03/16 09:26:58 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.6 2003/03/14 13:37:38 raasch18 ! Do loop constructs only for IBM, otherwise array notation re-introduced19 !20 ! Revision 1.5 2002/12/19 13:45:57 raasch21 ! Array notation changed to do-loop constructs due to better performance,22 ! OpenMP directives added23 !24 ! Revision 1.4 2002/05/15 10:00:18 raasch25 ! Nothing changed26 !27 ! Revision 1.3 2002/05/02 19:39:07 19:39:07 raasch (Siegfried Raasch)28 ! Use of module interfaces included29 !30 ! Revision 1.2 2002/05/02 14:06:02 raasch31 ! *** empty log message ***32 15 ! 33 16 ! Revision 1.1 2002/05/02 13:43:53 13:43:53 raasch (Siegfried Raasch) … … 38 21 ! ------------- 39 22 ! Time filter needed for the leap-frog method 40 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 41 24 42 25 USE arrays_3d -
palm/trunk/SOURCE/average_3d_data.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: average_3d_data.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.1 2006/02/23 09:48:58 raasch 12 14 ! Initial revision -
palm/trunk/SOURCE/boundary_conds.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: boundary_conds.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 09:54:55 raasch 12 14 ! Surface boundary conditions in case of topography: nzb replaced by … … 14 16 ! unchanged (still using nzb) because a non-flat topography must use a 15 17 ! Prandtl-layer, which don't requires explicit setting of the surface values. 16 !17 ! Revision 1.14 2005/06/29 09:49:40 steinfeld18 ! Baroclinicity considered in the Dirichlet boundary condition for u and v at19 ! the top boundary20 !21 ! Revision 1.13 2005/04/23 08:35:30 raasch22 ! Error removed in Dirichlet bottom boundary conditions for pt and q in case23 ! of Runge-Kutta schemes. So far, new surface values have been taken from24 ! timelevel t-dt, which does not exist in case of Runge-Kutta schemes.25 !26 ! Revision 1.12 2005/03/26 14:58:24 raasch27 ! Non-cyclic boundary conditions included, argument range added28 !29 ! Revision 1.11 2001/03/30 06:54:18 raasch30 ! Timelevel t+dt replaced by timelevel t,31 ! Translation of remaining German identifiers (variables, subroutines, etc.)32 !33 ! Revision 1.10 2001/01/29 12:19:27 raasch34 ! Passive scalar is considered35 !36 ! Revision 1.9 2001/01/22 05:25:56 raasch37 ! Module test_variables removed38 !39 ! Revision 1.8 2000/07/04 14:07:59 raasch40 ! Missing diriclet boundary conditions for temperature and total water41 ! content added42 !43 ! Revision 1.7 2000/04/13 13:56:05 schroeter44 ! Boundaray conditions for total water content45 !46 ! Revision 1.6 2000/01/20 10:44:48 10:44:48 letzel (Marcus Letzel)47 ! All comments translated into English48 !49 ! Revision 1.5 2000/01/10 09:28:53 raasch50 ! Variablenuebergabe jetzt per Modul51 ! Randbedingungen fuer w im Rahmen der pointer-Einfuehrung52 !53 ! Revision 1.4 1998/07/06 12:07:05 raasch54 ! + USE test_variables55 !56 ! Revision 1.3 1998/03/10 07:19:37 raasch57 ! Beschreibung ergaenzt um Randbedingung fuer TKE58 !59 ! Revision 1.2 1997/09/19 07:39:03 raasch60 ! Randbedingungen fuer TKE61 18 ! 62 19 ! Revision 1.1 1997/09/12 06:21:34 raasch -
palm/trunk/SOURCE/buoyancy.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: buoyancy.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.19 2006/04/26 12:09:56 raasch 12 14 ! OpenMP optimization (one dimension added to sums_l) 13 !14 ! Revision 1.18 2006/02/23 09:55:33 raasch15 ! nzb_2d replaced by nzb_s_inner/outer, nanz_2dh(0) replaced by16 ! ngp_2dh_outer(:,0)17 !18 ! Revision 1.17 2005/12/06 15:29:16 raasch19 ! Mean pt profile is calculated only in case of the first intermediate timestep20 !21 ! Revision 1.16 2004/01/30 10:14:31 raasch22 ! Scalar lower k index nzb replaced by 2d-array nzb_2d23 !24 ! Revision 1.15 2003/10/29 08:41:47 raasch25 ! Horizontal mean temperature is now taken from array hom instead of array sums26 !27 ! Revision 1.14 2003/03/16 09:27:22 raasch28 ! Two underscores (_) are placed in front of all define-strings29 !30 ! Revision 1.13 2003/03/12 16:20:47 raasch31 ! Full code replaced in the call for all gridpoints instead of calling the32 ! _ij version (required by NEC, because otherwise no vectorization)33 !34 ! Revision 1.12 2002/12/19 13:47:38 raasch35 ! STOP statement replaced by call of subroutine local_stop36 !37 ! Revision 1.11 2002/06/11 12:33:07 raasch38 ! Former subroutine changed to a module which allows to be called for all grid39 ! points of a single vertical column with index i,j or for all grid points by40 ! using function overloading.41 ! Calculation of the horizontally averaged temperature profile moved to new42 ! subroutine calc_mean_pt_profile, which is part of this module.43 !44 ! Revision 1.10 2001/03/30 06:54:56 raasch45 ! Translation of remaining German identifiers (variables, subroutines, etc.)46 !47 ! Revision 1.9 2001/01/22 05:27:03 raasch48 ! Module test_variables removed49 !50 ! Revision 1.8 2000/07/03 12:55:56 raasch51 ! Argument theta declared as pointer52 !53 ! Revision 1.7 2000/04/27 07:04:19 raasch54 ! in case of a sloping surface the temperature difference is computed relative55 ! to the initial (2D) temperature field56 !57 ! Revision 1.6 2000/04/13 14:27:34 schroeter58 ! considering the influence of humidity to the buoyancy term59 ! (only for sloping_surface = false),60 !61 ! Revision 1.5 2000/01/21 16:17:23 16:17:23 letzel (Marcus Letzel)62 ! All comments translated into English63 !64 ! Revision 1.4 1998/09/22 17:16:22 raasch65 ! Auftriebsterm erweitert fuer Rechnungen mit geneigter Ebene66 !67 ! Revision 1.3 1998/07/06 12:07:27 raasch68 ! + USE test_variables69 !70 ! Revision 1.2 1998/03/30 11:33:29 raasch71 ! nanz_2dh in nanz_2dh(0) geaendert72 15 ! 73 16 ! Revision 1.1 1997/08/29 08:56:48 raasch -
palm/trunk/SOURCE/calc_liquid_water_content.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: calc_liquid_water_content.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2005/03/26 15:22:06 raasch 12 14 ! Arguments for non-cyclic boundary conditions added to argument list of 13 15 ! routine exchange_horiz, 14 16 ! ql calculated for the ghost points, exchange of ghost points removed 15 !16 ! Revision 1.4 2004/01/30 10:15:13 raasch17 ! Scalar lower k index nzb replaced by 2d-array nzb_2d18 !19 ! Revision 1.3 2001/01/22 05:30:59 raasch20 ! Module test_variables removed21 !22 ! Revision 1.2 2000/06/13 07:49:32 schroeter23 ! Horizontal boundary conditions (cyclic) of ql will be set now24 17 ! 25 18 ! Revision 1.1 2000/04/13 14:50:45 schroeter -
palm/trunk/SOURCE/calc_precipitation.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: calc_precipitation.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2004/01/30 10:15:57 raasch 12 14 ! Scalar lower k index nzb replaced by 2d-array nzb_2d 13 !14 ! Revision 1.4 2003/03/12 16:21:23 raasch15 ! Full code replaced in the call for all gridpoints instead of calling the16 ! _ij version (required by NEC, because otherwise no vectorization)17 !18 ! Revision 1.3 2002/06/11 12:38:58 raasch19 ! Former subroutine changed to a module which allows to be called for all grid20 ! points of a single vertical column with index i,j or for all grid points by21 ! using function overloading.22 !23 ! Revision 1.2 2001/01/22 05:32:25 raasch24 ! Module test_variables removed25 15 ! 26 16 ! Revision 1.1 2000/04/13 14:45:22 schroeter -
palm/trunk/SOURCE/calc_radiation.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: calc_radiation.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.6 2004/01/30 10:17:03 raasch 12 14 ! Scalar lower k index nzb replaced by 2d-array nzb_2d 13 !14 ! Revision 1.5 2003/03/12 16:21:56 raasch15 ! Full code replaced in the call for all gridpoints instead of calling the16 ! _ij version (required by NEC, because otherwise no vectorization)17 !18 ! Revision 1.4 2002/06/11 12:40:33 raasch19 ! Former subroutine changed to a module which allows to be called for all grid20 ! points of a single vertical column with index i,j or for all grid points by21 ! using function overloading.22 ! 1D-arrays lwp_groud,... are only allocated once in the first call.23 !24 ! Revision 1.3 2001/03/30 06:55:27 raasch25 ! Translation of remaining German identifiers (variables, subroutines, etc.)26 !27 ! Revision 1.2 2001/01/22 05:35:05 raasch28 ! Module test_variables removed29 15 ! 30 16 ! Revision 1.1 2000/04/13 14:42:45 schroeter -
palm/trunk/SOURCE/calc_spectra.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: calc_spectra.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.9 2006/04/11 14:56:00 raasch 12 14 ! pl_spectra renamed data_output_sp 13 !14 ! Revision 1.8 2005/05/18 15:08:18 raasch15 ! n is assigned the number of spectra (formerly, the number of spectra was16 ! given by n-1)17 !18 ! Revision 1.7 2005/04/23 08:41:39 raasch19 ! Abort in case of non-cyclic lateral boundaries20 !21 ! Revision 1.6 2004/04/30 08:23:12 raasch22 ! Usage of all available FFT-methods is possible, direct transpositions from23 ! z to x and from z to y are used, FFT routines renamed, module24 ! fft_for_1d_decomp renamed fft_xy, extensions for enlarged subdomain added25 !26 ! Revision 1.5 2003/03/16 09:27:35 raasch27 ! Two underscores (_) are placed in front of all define-strings28 !29 ! Revision 1.4 2002/06/11 12:44:15 raasch30 ! TRIM function added (needed on IBM-systems)31 !32 ! Revision 1.3 2001/03/30 06:55:58 06:55:58 raasch (Siegfried Raasch)33 ! Work array replaced by tend,34 ! Translation of remaining German identifiers (variables, subroutines, etc.)35 !36 ! Revision 1.2 2001/01/25 06:52:21 raasch37 ! Module singleton moved to modules.f9038 15 ! 39 16 ! Revision 1.1 2001/01/05 15:08:07 raasch -
palm/trunk/SOURCE/check_for_restart.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: check_for_restart.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2007/02/11 12:55:13 raasch 12 14 ! Informative output to the job protocol 13 !14 ! Revision 1.10 2003/03/16 09:27:49 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.9 2002/12/19 13:48:59 raasch18 ! Stop flag is also set in case of user forced restarts, routine name19 ! changed from check_cpu_time to check_for_restart, the run description20 ! header will be written on file CONTINUE_RUN21 !22 ! Revision 1.8 2001/08/21 08:22:23 raasch23 ! MPI_ALLREDUCE moved out of IF-clause24 !25 ! Revision 1.7 2001/03/30 06:56:28 raasch26 ! Translation of remaining German identifiers (variables, subroutines, etc.)27 !28 ! Revision 1.6 2001/01/22 05:36:33 raasch29 ! Module test_variables removed30 !31 ! Revision 1.5 2001/01/02 17:19:27 raasch32 ! Initial value of terminate_run_l is set to .FALSE. (undefined before).33 ! Flag file will only be created by PE0. Unit 90 is used instead of 91.34 !35 ! Revision 1.4 2000/01/21 16:27:59 letzel36 ! All comments translated into English37 !38 ! Revision 1.3 1998/07/06 12:07:47 raasch39 ! + USE test_variables40 !41 ! Revision 1.2 1998/03/24 15:24:54 raasch42 ! Lauf wird nur dann beendet und Fortsetzungslauf gestartet, falls auch eine43 ! Abspeicherung der entsprechenden Binaerdateien vorgesehen ist.44 ! Beendigung der Rechnung auf mehreren Prozessoren richtig implementiert.45 15 ! 46 16 ! Revision 1.1 1998/03/18 20:06:51 raasch -
palm/trunk/SOURCE/check_open.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: check_open.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.44 2006/08/22 13:48:34 raasch 12 14 ! xz and yz cross sections now up to nzt+1 13 !14 ! Revision 1.43 2006/08/04 14:05:52 raasch15 ! Opening of unit 109 (particle time series)16 !17 ! Revision 1.42 2006/03/01 16:39:47 raasch18 ! In parallel mode, all PEs are opening the same copy of file PARIN19 !20 ! Revision 1.41 2006/02/23 10:05:05 raasch21 ! New unit 20 for data logging,22 ! new units (111:113,116) for 2d/3d data averaging23 ! NetCDF files are allowed a 64 bit offset format,24 ! error number argument for handle_netcdf_error25 !26 ! Revision 1.40 2005/06/30 10:39:04 raasch27 ! No opening of units 101-108 in case of missing netcdf support28 !29 ! Revision 1.39 2005/06/26 19:51:36 raasch30 ! Data format version number for unit 85 changed31 !32 ! Revision 1.38 2005/05/18 15:09:01 raasch33 ! Extensions for NetCDF output34 !35 ! Revision 1.37 2005/04/23 08:50:18 raasch36 ! batch_scp calls is given the remote username as an additional argument37 !38 ! Revision 1.36 2004/01/28 15:03:21 raasch39 ! In case of opening unit 80, a barrier is set only for the first call of40 ! check_open (from routine init_particles), in order to avoid the possibility41 ! of hanging jobs, which may occur if unit 80 is opened within routine42 ! allocate_prt_memory.43 !44 ! Revision 1.35 2003/10/29 08:47:03 raasch45 ! Call of batch_scp on IBM corrected, output of data format version and46 ! particle group data on unit 8547 !48 ! Revision 1.34 2003/05/09 14:06:58 raasch49 ! Batch_scp paths for IBM, NEC and linux added50 !51 ! Revision 1.33 2003/03/16 09:28:00 raasch52 ! Two underscores (_) are placed in front of all define-strings53 !54 ! Revision 1.32 2002/12/19 14:03:15 raasch55 ! Unit 18 (CPU_MEASURES) will be opened on PE0 only, unit 14 is opened56 ! with special process-id-string myid_char_14 (see routine parin),57 ! STOP statement replaced by call of subroutine local_stop58 !59 ! Revision 1.31 2002/04/24 19:01:15 raasch60 ! Unit 85 is allowed to be re-opened. In case that a re-open of a file is not61 ! allowed, the subroutine will be exit.62 !63 ! Revision 1.30 2002/04/16 08:00:27 raasch64 ! ftpcopy replaced by batch_scp, retry of opening CPU_MEASURES/ in case of65 ! errors forced by adding ERR-parameter in OPEN statement.66 ! New unit 85 for particle data67 !68 ! Revision 1.29 2001/07/12 12:07:21 raasch69 ! Opening of unit 33 moved out of ELSE-block70 !71 ! Revision 1.28 2001/03/30 06:56:49 raasch72 ! Translation of remaining German identifiers (variables, subroutines, etc.)73 !74 ! Revision 1.27 2001/01/29 12:21:52 raasch75 ! In case of AVS-output: if no path for ftpcopy is given, filenames are set76 ! to "unknown"77 !78 ! Revision 1.26 2001/01/25 06:54:28 raasch79 ! Error concerning initialization of unit 23 removed (initial informations80 ! were written also when reopened)81 !82 ! Revision 1.25 2001/01/05 15:09:29 raasch83 ! New units 81-84 for spectra84 !85 ! Revision 1.24 2001/01/02 17:21:11 raasch86 ! Re-opening of files is allowed, variable openfile now is of type file_status,87 ! new unit 80 for particle informations88 ! old revison remarks deleted89 !90 ! Revision 1.23 2000/04/27 06:51:12 raasch91 ! remote_addres changed to return_addres, old revision remarks deleted92 !93 ! Revision 1.22 2000/01/21 17:07:18 letzel94 ! All comments translated into English95 15 ! 96 16 ! Revision 1.1 1997/08/11 06:10:55 raasch -
palm/trunk/SOURCE/check_parameters.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: check_parameters.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.61 2006/08/04 14:20:25 raasch 12 14 ! do2d_unit and do3d_unit now defined as 2d-arrays, check of … … 14 16 ! generation of file header moved from routines palm and header to here 15 17 ! 16 ! Revision 1.60 2006/06/02 15:09:15 raasch17 ! data_output_xy/xz/yz are set18 !19 ! Revision 1.59 2006/03/03 15:55:34 letzel20 ! Spelling inconsistency removed21 !22 ! Revision 1.58 2006/03/03 08:50:08 raasch23 ! Default setting of constant_heatflux now only depends on the value of24 ! surface_heatflux again25 !26 ! Revision 1.57 2006/02/27 10:16:44 letzel27 ! Error removed in initialization of dt_averaging_input_pr28 !29 ! Revision 1.56 2006/02/23 10:06:55 raasch30 ! Volume flow conservation is not allowed with non-cyclic lateral boundaries,31 ! calculation of nzb_diff moved to init_grid, the default setting of32 ! constant_heatflux = .F. additionally requires a surface initial temperature33 ! change /= 0, check of topography, dissipation_1d, mixing_length_1d,34 ! check and setting of new data output steering parameters and user defined35 ! output quantities,36 ! ebene renamed section, pl1d_anz, plts_anz renamed dopr_n, dots_n,37 ! plts, pl1d renamed data_output_ts/1d38 !39 ! Revision 1.55 2005/12/06 16:39:34 raasch40 ! Output of ql profile is allowed in case of using cloud droplets41 !42 ! Revision 1.54 2005/06/30 11:27:33 steinfeld43 ! Initialization of the profiles of the geostrophic wind components44 ! ug and vg;45 ! the usage of the geostrophic wind in a galilei transformation is46 ! not allowed with a simultaneous specification of a baroclinicity47 !48 ! Revision 1.53 2005/05/18 15:12:36 raasch49 ! Check of data_output_format and netcdf precision,50 ! abort if poisfft_hybrid is called in a non-parallel environment51 !52 ! Revision 1.52 2005/04/23 09:04:10 raasch53 ! fcl_factor renamed cfl_factor, default setting of outflow_damping_width54 ! corrected, additional check for non-cyclic lateral boundary conditions55 !56 ! Revision 1.51 2005/03/26 19:47:32 raasch57 ! Check of bc_lr, bc_ns, km_damp_max, inflow_disturbance_begin,58 ! inflow_disturbance_end, outflow_damping_width59 !60 ! Revision 1.50 2004/04/30 11:20:27 raasch61 ! Check of grid_matching, check of time series profiles modified,62 ! impulse_advec renamed momentum_advec63 !64 ! Revision 1.49 2004/01/28 15:04:56 raasch65 ! Additional checks and settings for timestep schemes because of the66 ! newly implemented Runge-Kutta schemes67 !68 ! Revision 1.48 2003/04/17 07:07:43 raasch69 ! Hybrid solver check in parallel environment only70 !71 ! Revision 1.47 2003/04/16 09:43:20 raasch72 ! Additional checks for hybrid solver and host, Monin Obukhov length added to73 ! time series74 !75 ! Revision 1.46 2003/03/31 15:43:19 raasch76 ! Check for horizontal temperature flux profiles cleaned up77 !78 ! Revision 1.45 2003/03/14 13:38:02 raasch79 ! Random generator check included80 !81 ! Revision 1.44 2003/03/12 16:23:10 raasch82 ! Small checks for NEC and temperton-algorithm83 !84 ! Revision 1.43 2002/12/19 14:04:26 raasch85 ! Setting default values for termination_time_needed and check the86 ! value for IBM-regatta, STOP statement replaced by call of subroutine87 ! local_stop88 !89 ! Revision 1.42 2002/09/12 13:02:25 raasch90 ! Error in calculating initial pt and q profiles removed (array bound 10 in91 ! *_ind arrays could be exceeded)92 !93 ! Revision 1.41 2002/06/11 12:49:33 raasch94 ! Hybrid solver allowed for solving the Poisson equation95 !96 ! Revision 1.40 2002/05/02 18:48:56 18:48:56 raasch (Siegfried Raasch)97 ! SELECT case of cross_normalized_x and ..y requires usage of TRIM-function98 ! on IBM99 !100 ! Revision 1.39 2002/04/16 08:03:27 raasch101 ! Parameters for scalar transport included, output of scalar profiles changed102 !103 ! Revision 1.38 2001/11/28 09:48:32 schroeter104 ! Profiles for horizontal heat fluxes and105 ! flux divergences allowed (58-72)106 !107 ! Revision 1.37 2001/09/04 11:58:55 11:58:55 raasch (Siegfried Raasch)108 ! Profiles for vertical flux of resolved/subgrid scale energy and pressure109 ! fluctuations allowed (#55-57)110 !111 ! Revision 1.36 2001/07/20 12:52:32 raasch112 ! Multigrid method allowed for solving the Poisson equation113 !114 ! Revision 1.35 2001/03/30 06:57:38 raasch115 ! Asselin filter factor is set to 0, if Euler scheme is used,116 ! Translation of remaining German identifiers (variables, subroutines, etc.)117 !118 ! Revision 1.34 2001/02/09 07:04:45 raasch119 ! Checking of boundary conditions for pt and q and their combination with120 ! prescribed surface fluxes revised121 !122 ! Revision 1.33 2001/01/29 12:23:31 raasch123 ! New parameter passive_scalar is checked124 !125 ! Revision 1.32 2001/01/25 06:55:37 raasch126 ! New control-parameters fft_method and use_surface_fluxes are checked127 !128 ! Revision 1.31 2000/04/18 11:19:43 schroeter129 ! Bott-Chlond Advektionsschema nun auch fuer Rechnungen130 ! mit Feuchte/Wolkenphysik zugelassen131 !132 ! Revision 1.30 2000/04/13 13:28:56 schroeter133 ! + initialize profiles of total water content134 ! + check boundary conditions of total water content135 ! + initialize output of cloud-physics-profiles136 !137 ! Revision 1.29 2000/01/25 14:38:37 14:38:37 letzel (Marcus Letzel)138 ! All comments translated into English139 !140 18 ! Revision 1.1 1997/08/26 06:29:23 raasch 141 19 ! Initial revision -
palm/trunk/SOURCE/close_file.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: close_file.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.10 2006/08/22 13:50:01 raasch 12 14 ! xz and yz cross sections now up to nzt+1 13 !14 ! Revision 1.9 2006/08/04 14:22:53 raasch15 ! Closing unit 109 (particle timeseries)16 !17 ! Revision 1.8 2006/02/23 10:08:23 raasch18 ! Extensions for 2d/3d data averaging19 ! pl2d_.._anz, pl3d_anz renamed do2d_.._n, do3d_avs_n, respectively20 ! error number argument for handle_netcdf_error21 !22 ! Revision 1.7 2005/05/18 15:13:12 raasch23 ! Extensions for NetCDF output24 !25 ! Revision 1.6 2005/03/26 19:54:23 raasch26 ! +datform in global-namelist (iso2d)27 !28 ! Revision 1.5 2003/03/04 11:24:46 raasch29 ! AVS data format changed from float to xdr_float (needed on linux machines30 ! due to the little/big endian problem)31 !32 ! Revision 1.4 2002/12/19 14:06:34 raasch33 ! ISO2D parameter dp set to true for ibm hosts34 !35 ! Revision 1.3 2001/03/30 06:58:36 raasch36 ! Translation of remaining German identifiers (variables, subroutines, etc.)37 !38 ! Revision 1.2 2001/01/22 05:39:06 raasch39 ! Module test_variables removed40 15 ! 41 16 ! Revision 1.1 2001/01/02 17:23:41 raasch -
palm/trunk/SOURCE/compute_vpt.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: compute_vpt.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2001/03/30 06:58:52 raasch 12 14 ! Translation of remaining German identifiers (variables, subroutines, etc.) 13 !14 ! Revision 1.4 2001/01/22 05:42:55 raasch15 ! Minor corrections within file header16 !17 ! Revision 1.2 2000/09/07 09:00:36 schroeter18 ! Now, vpt will be computed for k=nzb and nzt+119 ! additionally to keep boundary values20 15 ! 21 16 ! Revision 1.1 2000/04/13 14:40:53 schroeter -
palm/trunk/SOURCE/coriolis.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: coriolis.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.12 2006/02/23 10:08:57 raasch 12 14 ! nzb_2d replaced by nzb_u/v/w_inner 13 !14 ! Revision 1.11 2005/06/29 10:30:17 steinfeld15 ! A potential dependency of the geostrophic wind components ug and vg on height16 ! is considered17 !18 ! Revision 1.10 2005/03/26 20:02:04 raasch19 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions20 !21 ! Revision 1.9 2004/01/30 10:17:34 raasch22 ! Scalar lower k index nzb replaced by 2d-array nzb_2d23 !24 ! Revision 1.8 2003/03/12 16:23:47 raasch25 ! Full code replaced in the call for all gridpoints instead of calling the26 ! _ij version (required by NEC, because otherwise no vectorization)27 !28 ! Revision 1.7 2002/12/19 14:08:14 raasch29 ! STOP statement replaced by call of subroutine local_stop30 !31 ! Revision 1.6 2002/06/11 12:51:09 raasch32 ! Former subroutine changed to a module which allows to be called for all grid33 ! points of a single vertical column with index i,j or for all grid points by34 ! using function overloading.35 !36 ! Revision 1.5 2001/03/30 06:59:09 raasch37 ! Translation of remaining German identifiers (variables, subroutines, etc.)38 !39 ! Revision 1.4 2001/01/22 05:44:16 raasch40 ! Module test_variables removed41 !42 ! Revision 1.3 2000/12/20 09:50:27 letzel43 ! All comments translated into English.44 !45 ! Revision 1.2 1998/07/06 12:09:54 raasch46 ! + USE test_variables47 15 ! 48 16 ! Revision 1.1 1997/08/29 08:57:38 raasch -
palm/trunk/SOURCE/cpu_log.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: cpu_log.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.24 2006/06/02 15:12:17 raasch 12 14 ! cpp-directives extended for lctit 13 !14 ! Revision 1.23 2004/01/28 15:08:52 raasch15 ! Type log changed to logpoint due to name conflict with intrinsic log16 !17 ! Revision 1.22 2003/05/09 14:15:18 raasch18 ! Time measurement on linux machines included, measurements on IBM are now using19 ! function irtc, which allows correct measurements for jobs running over20 ! the 24:00 timeline21 !22 ! Revision 1.21 2003/03/16 09:30:04 raasch23 ! Two underscores (_) are placed in front of all define-strings24 !25 ! Revision 1.20 2003/03/12 16:24:18 raasch26 ! Time measurement on NEC implemented27 !28 ! Revision 1.19 2002/12/20 09:38:43 raasch29 ! variable first is defined for all hosts30 !31 ! Revision 1.18 2002/12/19 14:09:36 raasch32 ! Output of warnings in case of negative cpu-times, STOP statement replaced by33 ! call of subroutine local_stop34 !35 ! Revision 1.17 2002/05/02 18:49:37 raasch36 ! Time measurement on IBM implemented37 !38 ! Revision 1.16 2001/03/30 07:00:43 raasch39 ! Translation of remaining German identifiers (variables, subroutines, etc.)40 !41 ! Revision 1.14 2001/01/25 06:56:49 raasch42 ! +cpp-directives for dec-alpha-workstations43 !44 ! Revision 1.13 2001/01/22 08:32:46 raasch45 ! Module test_variables removed46 !47 ! Revision 1.11 2000/12/20 10:18:58 letzel48 ! All comments translated into English.49 !50 ! Revision 1.9 1999/11/25 16:19:07 raasch51 ! TIMEF jetzt als Function-Aufruf52 !53 ! Revision 1.8 1999/03/03 09:26:14 raasch54 ! Zeitmessung auf T3E grundsaetzlich auf TIMEF umgestellt55 !56 ! Revision 1.7 1998/12/15 07:32:39 raasch57 ! Messung der wahren I/O-Zeiten auf T3E mit TIMEF moeglich58 !59 ! Revision 1.4 1997/09/16 06:37:24 raasch60 ! Zeitmessungen auf t3e umgestellt auf tsecnd (schneller als mpi_wtime)61 15 ! 62 16 ! Revision 1.1 1997/07/24 11:12:29 raasch … … 67 21 ! ------------ 68 22 ! Cpu-time measurements for any program part whatever. 69 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 70 24 71 25 USE cpulog -
palm/trunk/SOURCE/cpu_statistics.f90
r1 r3 1 1 SUBROUTINE cpu_statistics 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: cpu_statistics.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.13 2006/04/26 12:10:51 raasch 12 14 ! Output of number of threads per task, max = min in case of 1 PE 13 !14 ! Revision 1.12 2003/03/16 09:30:22 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.11 2002/12/19 14:12:08 raasch18 ! PE results are collected on PE0 in order to calculate mean cpu statistics.19 ! Output format changed.20 !21 ! Revision 1.10 2001/03/30 07:01:02 raasch22 ! Translation of remaining German identifiers (variables, subroutines, etc.)23 !24 ! Revision 1.9 2001/01/22 05:49:49 raasch25 ! Module test_variables removed26 !27 ! Revision 1.8 2001/01/02 17:25:23 raasch28 ! Unit 18 closed at the end of the subroutine29 !30 ! Revision 1.7 2000/12/20 10:10:26 letzel31 ! All comments translated into English.32 15 ! 33 16 ! Revision 1.1 1997/07/24 11:11:11 raasch … … 41 24 ! statistics. The output is sorted according to the amount of cpu-time consumed 42 25 ! and output on PE0. 43 !------------------------------------------------------------------------------ -!26 !------------------------------------------------------------------------------! 44 27 45 28 USE cpulog -
palm/trunk/SOURCE/data_log.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_log.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.1 2006/02/23 10:09:29 raasch 12 14 ! Initial revision -
palm/trunk/SOURCE/data_output_2d.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_2d.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2006/08/22 13:50:29 raasch 12 14 ! xz and yz cross sections now up to nzt+1 13 !14 ! Revision 1.4 2006/06/06 07:55:18 raasch15 ! Error removed: p is assigned to to_be_resorted instead of w16 !17 ! Revision 1.3 2006/06/02 15:18:21 raasch18 ! +argument "found" in call of routine user_data_output_2d,19 ! immediate return, if no output is requested20 15 ! 21 16 ! Revision 1.2 2006/02/23 10:19:22 raasch … … 26 21 ! Former subroutine plot_2d renamed data_output_2d, pl2d.. renamed do2d.., 27 22 ! anz renamed ngp, ebene renamed section, pl2d_.._anz renamed do2d_.._n 28 !29 ! Revision 1.1 2006/02/23 10:12:31 raasch30 ! Initial revision31 !32 ! Revision 1.25 2005/10/20 15:41:41 raasch33 ! Error in output of netcdf yz-cross-sections removed34 !35 ! Revision 1.24 2005/06/26 20:07:25 raasch36 ! +ql_c, ql_v, ql_vp37 !38 ! Revision 1.23 2005/05/18 15:46:56 raasch39 ! Extensions for NetCDF output, output of vpt, qv and lwp40 !41 ! Revision 1.22 2005/03/26 20:47:45 raasch42 ! local_2d is allocated with upper bound nzt (former nzt+1 was wrong),43 ! additional barrier in case that PE0 gathers data from the other44 ! PEs, this barrier had to be set at the end of the gathering because45 ! otherwise PE0 may receive wrong data on tag 0.46 !47 ! Revision 1.21 2004/04/30 12:32:50 raasch48 ! Adjustments for enlarged grids (unequal subdomain sizes)49 !50 ! Revision 1.20 2003/03/16 09:42:27 raasch51 ! Two underscores (_) are placed in front of all define-strings52 !53 ! Revision 1.19 2002/12/19 15:56:24 raasch54 ! STOP statement replaced by call of subroutine local_stop55 !56 ! Revision 1.18 2002/06/11 13:13:00 raasch57 ! TRIM function added (needed on IBM-systems)58 !59 ! Revision 1.17 2002/04/16 08:11:16 08:11:16 raasch (Siegfried Raasch)60 ! Output of scalar field allowed with strings s_xy, s_xz, s_yz61 !62 ! Revision 1.16 2001/03/30 07:44:29 raasch63 ! Translation of remaining German identifiers (variables, subroutines, etc.)64 !65 ! Revision 1.15 2001/03/09 21:36:29 letzel66 ! All comments translated into English.67 !68 ! Revision 1.14 2001/01/29 12:32:07 raasch69 ! Passive scalar is considered, xy-plot of q and ql corrected (level_z was70 ! missing)71 !72 ! Revision 1.13 2001/01/22 07:38:20 raasch73 ! Module test_variables removed74 !75 ! Revision 1.12 2001/01/02 17:33:39 raasch76 ! Plot output file is closed at the end of the subroutine77 !78 ! Revision 1.11 2000/04/13 14:01:44 schroeter79 ! Plotausgaben fuer Gesamtwasser- und Fluessigwassergehalt80 !81 ! Revision 1.10 99/02/05 09:14:06 09:14:06 raasch (Siegfried Raasch)82 ! u=v=0 am unteren Rand bei Vertikalschnitten wird jetzt ausschliesslich83 ! bei Spiegelungsrandbedingungen verwendet84 !85 ! Revision 1.9 1998/09/22 17:25:58 raasch86 ! u*- und theta*-Ausgabe korrigiert (muessen auf nzb+1 des Hilfsfeldes local_2d87 ! umgespeichert werden.88 !89 ! Revision 1.8 1998/07/06 12:26:06 raasch90 ! Exakte Ausgabe der Lage der Schnittebenen fuer die verschiedenen Felder91 ! bei Horizontalschnitten, + USE test_variables92 !93 ! Revision 1.7 1998/03/09 16:24:19 raasch94 ! TRIM( simulated_time_chr )95 !96 ! Revision 1.6 1998/02/19 07:10:30 raasch97 ! Ausgabe der Teilfelder jetzt auch wieder wahlweise durch die einzelnen PEs98 !99 ! Revision 1.5 1998/02/10 15:08:25 raasch100 ! Spiegelungsrandbedingung fuer u und v bei Vertikalschnitten durch101 ! u=v=0 ersetzt, Ausgabe der Schnittlevel,102 ! Empfang der Teilfelder durch PE0 erfolgt in beliebiger Reihenfolge durch103 ! Verwendung von MPI_ANY_TAG, Barrieren werden gesetzt, um MPI-Fehler zu umgehen104 !105 ! Revision 1.4 1998/01/23 10:05:42 raasch106 ! Ausgabe von 2D-Feldern (u*_xy, t*_xy)107 !108 ! Revision 1.3 1997/09/08 10:47:23 raasch109 ! Prozessor 0 gibt jetzt gesamte Schnittebene des Modells aus110 !111 ! Revision 1.2 1997/08/29 09:00:36 raasch112 ! Ausgabe der potentiellen Temperatur113 23 ! 114 24 ! Revision 1.1 1997/08/11 06:24:09 raasch -
palm/trunk/SOURCE/data_output_3d.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_3d.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.3 2006/06/02 15:18:59 raasch 12 14 ! +argument "found", -argument grid in call of routine user_data_output_3d … … 16 18 ! .._anz renamed .._n, 17 19 ! output extended to (almost) all quantities, output of user-defined quantities 18 !19 ! Revision 1.1 2006/02/23 10:21:34 raasch20 ! Initial revision21 !22 ! Revision 1.13 2005/06/26 20:08:36 raasch23 ! +ql_c, ql_v, ql_vp24 !25 ! Revision 1.12 2005/05/18 15:47:34 raasch26 ! Extensions for NetCDF output27 !28 ! Revision 1.11 2003/03/16 09:42:31 raasch29 ! Two underscores (_) are placed in front of all define-strings30 !31 ! Revision 1.10 2002/06/11 13:13:27 raasch32 ! TRIM function added (needed on IBM-systems)33 !34 ! Revision 1.9 2002/04/16 08:11:50 08:11:50 raasch (Siegfried Raasch)35 ! Output of scalar field s allowed36 !37 ! Revision 1.8 2001/03/30 07:44:45 raasch38 ! Translation of remaining German identifiers (variables, subroutines, etc.)39 !40 ! Revision 1.7 2001/03/09 21:46:23 letzel41 ! All comments translated into English.42 !43 ! Revision 1.6 2001/01/22 07:40:10 raasch44 ! Module test_variables removed45 !46 ! Revision 1.5 2000/04/13 14:02:39 schroeter47 ! Plotausgabe fuer Gesamtwasser- und Fluessigwassergehalt sowie potentielle48 ! Fluessigwassertemperatur49 !50 ! Revision 1.4 99/03/02 09:24:51 09:24:51 raasch (Siegfried Raasch)51 ! Datenkompression implementiert52 !53 ! Revision 1.3 1998/07/06 12:26:35 raasch54 ! + USE test_variables55 !56 ! Revision 1.2 1998/03/09 16:24:33 raasch57 ! simulated_time_mod verlaengert58 20 ! 59 21 ! Revision 1.1 1997/09/03 06:29:36 raasch -
palm/trunk/SOURCE/data_output_dvrp.f90
r1 r3 32 32 ! Actual revisions: 33 33 ! ----------------- 34 !35 34 ! TEST: write statements 36 35 ! 37 36 ! Former revisions: 38 37 ! ----------------- 39 ! $Log: data_output_dvrp.f90,v $ 38 ! $Id$ 39 ! RCS Log replace by Id keyword, revision history cleaned up 40 ! 40 41 ! Revision 1.13 2006/02/23 10:25:12 raasch 41 42 ! Former routine plot_dvrp renamed data_output_dvrp, 42 43 ! Only a fraction of the particles may have a tail, 43 44 ! pl.. replaced by do.., %size renamed %dvrp_psize 44 !45 ! Revision 1.12 2005/03/31 07:00:46 raasch46 ! Testversion with print statements checked in47 !48 ! Revision 1.11 2003/08/01 08:44:50 raasch49 ! Routine defined as recursive, because it is called from DVRP_STEERING_UPDATE50 ! in case that the user has suspended the run by steering control.51 ! Routine color_dvrp placed at the beginning of this file. This routine is used52 ! for slicer colortable settings.53 !54 ! Revision 1.10 2003/03/16 09:42:34 raasch55 ! Two underscores (_) are placed in front of all define-strings56 !57 ! Revision 1.9 2003/03/04 11:35:07 raasch58 ! Update for new version of dvrp software (e.g. using module dvrp is now59 ! mandatory)60 !61 ! Revision 1.8 2001/11/09 13:41:21 raasch62 ! Colour can change along the particle tails63 !64 ! Revision 1.7 2001/08/21 09:55:20 raasch65 ! Superelevation for particles implemented, plotting of particle tails added,66 ! colourtable definitions moved to subroutine user_dvrp_coltab67 !68 ! Revision 1.6 2001/07/12 12:44:12 raasch69 ! Test: color definitions for particles changed70 !71 ! Revision 1.5 2001/03/30 07:45:00 raasch72 ! Translation of remaining German identifiers (variables, subroutines, etc.)73 !74 ! Revision 1.4 2001/01/29 12:32:51 raasch75 ! Passive scalar considered76 !77 ! Revision 1.3 2001/01/22 07:41:12 raasch78 ! Module test_variables removed79 !80 ! Revision 1.2 2000/12/28 13:37:19 raasch81 ! Total revision due to new features of dvrp-software,82 ! plot of slicers included,83 ! code is used only optionally (cpp-directives are added)84 45 ! 85 46 ! Revision 1.1 2000/04/27 06:27:17 raasch -
palm/trunk/SOURCE/data_output_profiles.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_profiles.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.18 2006/08/16 14:27:04 raasch 12 14 ! PRINT* statements for testing removed 13 !14 ! Revision 1.17 2006/08/08 12:45:47 raasch15 ! Immediate return in case of average_count_pr = 0.16 !17 ! Revision 1.16 2006/02/23 10:26:53 raasch18 ! Former routine plot_1d renamed data_output_profiles19 ! pl.. renamed do.., .._anz renamed .._no20 ! error number argument for handle_netcdf_error21 !22 ! Revision 1.15 2005/05/18 15:44:44 raasch23 ! Extensions for NetCDF output24 !25 ! Revision 1.14 2001/03/30 07:44:12 raasch26 ! Translation of remaining German identifiers (variables, subroutines, etc.)27 !28 ! Revision 1.13 2001/03/02 14:21:17 letzel29 ! All comments translated into English.30 !31 ! Revision 1.12 2001/01/22 07:36:55 raasch32 ! Module test_variables removed33 !34 ! Revision 1.11 1998/07/06 12:25:36 raasch35 ! + USE test_variables36 !37 ! Revision 1.10 1998/04/15 11:22:02 raasch38 ! zeitlich gemittelte Profile werden jetzt auch mit zeitlich gemittelten39 ! Groessen normiert40 !41 ! Revision 1.9 1998/04/06 14:33:53 raasch42 ! Implementierung weiterer Normierungsmethoden43 !44 ! Revision 1.8 1998/03/30 11:37:30 raasch45 ! Erweiterungen fuer Gebietsstatistiken korrigiert46 !47 ! Revision 1.7 1998/03/25 13:55:45 raasch48 ! Erweiterungen fuer Gebietsstatistiken49 !50 ! Revision 1.6 1998/03/09 16:24:04 raasch51 ! Ueberpruefung der Normierungsfaktoren auf 0.0 und evtl.52 ! Abschaltung der Normalisierung53 !54 ! Revision 1.5 1998/03/03 07:58:11 raasch55 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier56 ! Ausgabe von zeitlich gemittelten Profilen, Abspeichern von Normierungsgroessen57 !58 ! Revision 1.4 1998/02/19 07:10:04 raasch59 ! x-Wertebereichsermittelung nur bezueglich des auszugebenden Hoehenbereiches,60 ! x-Wertebereichsermittelung auch fuer Anfangsprofile61 !62 ! Revision 1.3 1998/02/04 16:09:11 raasch63 ! x-Wertebereiche fuer die einzelnen Koordinatenkreuze ermitteln64 !65 ! Revision 1.2 1997/09/16 06:38:20 raasch66 ! style- und color-Indexberechnung korrigiert67 15 ! 68 16 ! Revision 1.1 1997/09/12 06:28:48 raasch -
palm/trunk/SOURCE/data_output_ptseries.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_ptseries.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2006/08/22 13:51:13 raasch 12 14 ! Seperate output for particle groups -
palm/trunk/SOURCE/data_output_spectra.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_spectra.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.7 2006/04/11 14:56:38 raasch 12 14 ! pl_spectra renamed data_output_sp 13 15 ! 14 ! Revision 1.6 2006/02/23 10:29:38 raasch15 ! Former routine plot_spectra renamed data_output_spectra16 ! error number argument for handle_netcdf_error17 !18 ! Revision 1.5 2005/05/18 15:48:29 raasch19 ! Extensions for NetCDF output20 !21 ! Revision 1.4 2003/03/16 09:42:42 raasch22 ! Two underscores (_) are placed in front of all define-strings23 !24 ! Revision 1.3 2002/06/11 13:14:02 raasch25 ! TRIM function added (needed on IBM-systems), additional translations into26 ! English and further formatting27 !28 ! Revision 1.2 2001/03/30 07:45:18 07:45:18 raasch (Siegfried Raasch)29 ! Translation of remaining German identifiers (variables, subroutines, etc.)30 !31 16 ! Revision 1.1 2001/01/05 15:14:20 raasch 32 17 ! Initial revision -
palm/trunk/SOURCE/data_output_tseries.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: data_output_tseries.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.13 2006/03/14 12:42:51 raasch 12 14 ! Error removed: NetCDF output only if switched on 13 !14 ! Revision 1.12 2006/02/23 10:30:54 raasch15 ! Former routine plot_ts renamed data_output_tseries16 ! error number argument for handle_netcdf_error17 ! plts_anz renamed dots_n18 !19 ! Revision 1.11 2005/05/18 15:50:54 raasch20 ! Extensions for NetCDF output21 !22 ! Revision 1.10 2003/04/16 13:10:40 raasch23 ! Output of Monin Obukhov length included24 !25 ! Revision 1.9 2001/03/30 07:45:34 raasch26 ! Translation of remaining German identifiers (variables, subroutines, etc.)27 !28 ! Revision 1.8 2001/01/22 20:36:59 letzel29 ! All comments translated into English.30 !31 ! Revision 1.7 2001/01/22 07:42:50 raasch32 ! Module test_variables removed33 !34 ! Revision 1.6 1999/02/05 09:14:35 raasch35 ! Erweiterung fuer Ausgabe von Upstream-Anteilen, Erweiterung von ts_value36 ! auf 30 Elemente37 !38 ! Revision 1.5 1998/07/06 12:26:59 raasch39 ! + USE test_variables40 !41 ! Revision 1.4 1998/04/15 11:22:22 raasch42 ! Erweiterung der Zeitreihenausgabe auf oberflaechennahe Temperaturen und43 ! Waermestroeme44 !45 ! Revision 1.3 1998/03/30 11:48:03 raasch46 ! Divergenzen stehen jetzt in hom47 !48 ! Revision 1.2 1998/03/25 13:56:01 raasch49 ! Erweiterungen fuer Gebietsstatistiken,50 ! dt in dt_3d umbenannt51 15 ! 52 16 ! Revision 1.1 1998/03/03 08:00:13 raasch -
palm/trunk/SOURCE/diffusion_e.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusion_e.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.18 2006/08/04 14:29:43 raasch 12 14 ! dissipation is stored in extra array diss if needed later on for calculating 13 15 ! the sgs particle velocities 14 !15 ! Revision 1.17 2006/02/23 10:31:46 raasch16 ! nzb_2d replaced by nzb_s_inner17 !18 ! Revision 1.16 2004/01/30 10:18:18 raasch19 ! Scalar lower k index nzb replaced by 2d-array nzb_2d20 !21 ! Revision 1.15 2003/03/14 13:39:33 raasch22 ! Loop optimization for diffusion_e, l and ll are now automatic arrays23 !24 ! Revision 1.14 2003/03/12 16:25:03 raasch25 ! Full code replaced in the call for all gridpoints instead of calling the26 ! _ij version (required by NEC, because otherwise no vectorization)27 !28 ! Revision 1.13 2002/12/19 14:25:28 raasch29 ! Correction of mixing length term (l(k)/ll(k)). The condition kh=3*km in30 ! the unstable case is now also exactly met in the wall adjustment region.31 ! Factor 0.7 in wall adjustment part replaced by variable32 ! wall_adjustment_factor, which is set to ... in modules.f90.33 !34 ! Revision 1.12 2002/06/11 12:51:59 raasch35 ! Former subroutine changed to a module which allows to be called for all grid36 ! points of a single vertical column with index i,j or for all grid points by37 ! using function overloading.38 ! 1D-array l is allocated only once in the first call.39 !40 ! Revision 1.11 2001/08/21 08:24:34 raasch41 ! Wall adjustment of mixing length to 0.7 z can be switched off42 !43 ! Revision 1.10 2001/03/30 07:06:30 raasch44 ! Near surface mixing length is limited to 0.7*zu,45 ! e**1.5 replaced by e*SQRT(e) (more than 10% total increase in performance46 ! of this routine),47 ! Translation of remaining German identifiers (variables, subroutines, etc.)48 !49 ! Revision 1.9 2001/01/22 06:05:28 raasch50 ! Module test_variables removed51 !52 ! Revision 1.8 2001/01/02 17:27:00 raasch53 ! -dpt_dz_d, dpt_dz_u54 !55 ! Revision 1.7 2000/07/03 12:56:34 raasch56 ! array l changed from dummy argument to local allocatable array,57 ! dummy arguments, whose corresponding actual arguments are pointers,58 ! are now also defined as pointers59 ! all comments translated into English60 !61 ! Revision 1.6 2000/04/18 08:10:12 schroeter62 ! Revision 1.4 wieder rueckgaengig gemacht, das Stabilitaets-63 ! kriterium basiert nun wieder auf zentralen Differenzen64 !65 ! Revision 1.5 2000/04/13 14:33:08 schroeter66 ! je nach Initialisierungsmodus (trocken/feucht) fliesst in die67 ! Berechnung des Mischungsweges pt oder vpt ein, wird durch68 ! entsprechende Variablenuebergabe geregelt69 !70 ! Revision 1.4 99/02/17 09:15:52 09:15:52 raasch (Siegfried Raasch)71 ! Dissipation jetzt gemaess dem originalen Deardorff-Ansatz72 ! Kriterium fuer reduzierten Mischungsweg im stabil geschichteten Fall enger73 ! gefasst (Schichtung muss sowohl oberhalb als auch unterhalb des betrachteten74 ! Gitterpunkts stabil sein)75 !76 ! Revision 1.3 1998/07/06 12:10:56 raasch77 ! + USE test_variables78 !79 ! Revision 1.2 1998/03/11 11:48:59 raasch80 ! Anpassung des Mischungsweges an den Prandtlschen Mischungsweg moeglich81 16 ! 82 17 ! Revision 1.1 1997/09/19 07:40:24 raasch -
palm/trunk/SOURCE/diffusion_s.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusion_s.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.8 2006/02/23 10:34:17 raasch 12 14 ! nzb_2d replaced by nzb_s_outer in horizontal diffusion and by nzb_s_inner 13 15 ! or nzb_diff_s_inner, respectively, in vertical diffusion, prescribed surface 14 16 ! fluxes at vertically oriented topography 15 !16 ! Revision 1.7 2004/01/30 10:20:56 raasch17 ! Scalar lower k index nzb replaced by 2d-array nzb_2d18 !19 ! Revision 1.6 2003/03/12 16:25:32 raasch20 ! Full code replaced in the call for all gridpoints instead of calling the21 ! _ij version (required by NEC, because otherwise no vectorization)22 !23 ! Revision 1.5 2002/06/11 12:52:41 raasch24 ! Former subroutine changed to a module which allows to be called for all grid25 ! points of a single vertical column with index i,j or for all grid points by26 ! using function overloading.27 !28 ! Revision 1.4 2001/03/30 07:11:44 raasch29 ! Translation of remaining German identifiers (variables, subroutines, etc.)30 !31 ! Revision 1.3 2001/01/25 06:58:14 raasch32 ! Variable "prandtl_layer replaced by "use_surface_fluxes"33 !34 ! Revision 1.2 2000/07/03 12:57:13 raasch35 ! dummy arguments, whose corresponding actual arguments are pointers,36 ! are now also defined as pointers,37 ! all comments translated into English38 17 ! 39 18 ! Revision 1.1 2000/04/13 14:54:02 schroeter -
palm/trunk/SOURCE/diffusion_u.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusion_u.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 10:35:35 raasch 12 14 ! nzb_2d replaced by nzb_u_outer in horizontal diffusion and by nzb_u_inner … … 16 18 ! because they cause errors at the edges of topography 17 19 ! WARNING: loops containing the MAX function are still not properly vectorized! 18 !19 ! Revision 1.14 2005/03/26 20:08:21 raasch20 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,21 ! additional damping layer at the outflow in case of non-cyclic lateral22 ! boundaries, additional argument km_damp_y23 !24 ! Revision 1.13 2004/01/30 10:21:32 raasch25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d26 !27 ! Revision 1.12 2003/03/12 16:26:14 raasch28 ! Full code replaced in the call for all gridpoints instead of calling the29 ! _ij version (required by NEC, because otherwise no vectorization)30 !31 ! Revision 1.11 2002/06/11 12:53:12 raasch32 ! Former subroutine changed to a module which allows to be called for all grid33 ! points of a single vertical column with index i,j or for all grid points by34 ! using function overloading.35 !36 ! Revision 1.10 2001/08/21 08:25:55 raasch37 ! Special treatment at k=1 generally if momentum flux is prescribed (not only in38 ! case of a Prandtl layer)39 !40 ! Revision 1.9 2001/03/30 07:12:02 raasch41 ! Translation of remaining German identifiers (variables, subroutines, etc.)42 !43 ! Revision 1.8 2001/01/22 06:14:00 raasch44 ! Module test_variables removed45 !46 ! Revision 1.7 2000/07/03 12:57:46 raasch47 ! dummy arguments, whose corresponding actual arguments are pointers,48 ! are now also defined as pointers,49 ! all comments translated into English50 !51 ! Revision 1.6 2000/03/14 13:55:33 schroeter52 ! vertical impulse fluxes are computed over dz at the53 ! top of the Prandtl-layer54 !55 ! Revision 1.5 98/07/06 12:11:42 12:11:42 raasch (Siegfried Raasch)56 ! + USE test_variables57 !58 ! Revision 1.4 1998/01/23 09:56:44 raasch59 ! Einbau der Prandtl-Schicht60 !61 ! Revision 1.3 1997/09/12 07:22:50 raasch62 ! Leerzeilen mussten entfernt werden63 !64 ! Revision 1.2 1997/09/12 06:43:27 raasch65 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen66 20 ! 67 21 ! Revision 1.1 1997/09/12 06:23:51 raasch -
palm/trunk/SOURCE/diffusion_v.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusion_v.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.15 2006/02/23 10:36:00 raasch 12 14 ! nzb_2d replaced by nzb_v_outer in horizontal diffusion and by nzb_v_inner … … 16 18 ! because they cause errors at the edges of topography 17 19 ! WARNING: loops containing the MAX function are still not properly vectorized! 18 !19 ! Revision 1.14 2005/03/26 20:09:56 raasch20 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,21 ! additional damping layer at the outflow in case of non-cyclic lateral22 ! boundaries, additional argument km_damp_x23 !24 ! Revision 1.13 2004/01/30 10:21:56 raasch25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d26 !27 ! Revision 1.12 2003/03/12 16:26:31 raasch28 ! Full code replaced in the call for all gridpoints instead of calling the29 ! _ij version (required by NEC, because otherwise no vectorization)30 !31 ! Revision 1.11 2002/06/11 12:53:46 raasch32 ! Former subroutine changed to a module which allows to be called for all grid33 ! points of a single vertical column with index i,j or for all grid points by34 ! using function overloading.35 !36 ! Revision 1.10 2001/08/21 08:29:31 raasch37 ! Special treatment at k=1 generally if momentum flux is prescribed (not only in38 ! case of a Prandtl layer)39 !40 ! Revision 1.9 2001/03/30 07:12:19 raasch41 ! Translation of remaining German identifiers (variables, subroutines, etc.)42 !43 ! Revision 1.8 2001/01/22 06:18:41 raasch44 ! Module test_variables removed45 !46 ! Revision 1.7 2000/07/03 12:58:14 raasch47 ! dummy arguments, whose corresponding actual arguments are pointers,48 ! are now also defined as pointers,49 ! all comments translated into English50 !51 ! Revision 1.6 2000/03/14 13:58:18 schroeter52 ! vertical impulse fluxes are computed over dz at the53 ! top of the Prandtl-layer54 !55 ! Revision 1.5 98/07/06 12:12:05 12:12:05 raasch (Siegfried Raasch)56 ! + USE test_variables57 !58 ! Revision 1.4 1998/01/23 09:57:18 raasch59 ! Einbau der Prandtl-Schicht60 !61 ! Revision 1.3 1997/09/12 07:24:20 raasch62 ! Leerzeilen mussten entfernt werden63 !64 ! Revision 1.2 1997/09/12 06:43:46 raasch65 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen66 20 ! 67 21 ! Revision 1.1 1997/09/12 06:24:01 raasch -
palm/trunk/SOURCE/diffusion_w.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusion_w.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.12 2006/02/23 10:38:03 raasch 12 14 ! nzb_2d replaced by nzb_w_outer, wall functions added for all vertical walls, 13 15 ! +z0 in argument list 14 16 ! WARNING: loops containing the MAX function are still not properly vectorized! 15 !16 ! Revision 1.11 2005/03/26 20:11:05 raasch17 ! Additional damping layer at the outflow in case of non-cyclic lateral18 ! boundaries, additional arguments km_damp_x, km_damp_y19 !20 ! Revision 1.10 2004/01/30 10:22:32 raasch21 ! Scalar lower k index nzb replaced by 2d-array nzb_2d22 !23 ! Revision 1.9 2003/03/12 16:26:49 raasch24 ! Full code replaced in the call for all gridpoints instead of calling the25 ! _ij version (required by NEC, because otherwise no vectorization)26 !27 ! Revision 1.8 2002/06/11 12:54:16 raasch28 ! Former subroutine changed to a module which allows to be called for all grid29 ! points of a single vertical column with index i,j or for all grid points by30 ! using function overloading.31 !32 ! Revision 1.7 2001/03/30 07:12:35 raasch33 ! Translation of remaining German identifiers (variables, subroutines, etc.)34 !35 ! Revision 1.6 2001/01/22 06:20:12 raasch36 ! Module test_variables removed37 !38 ! Revision 1.5 2000/07/03 12:58:40 raasch39 ! dummy arguments, whose corresponding actual arguments are pointers,40 ! are now also defined as pointers,41 ! all comments translated into English42 !43 ! Revision 1.4 1998/07/06 12:12:24 raasch44 ! + USE test_variables45 !46 ! Revision 1.3 1997/09/12 07:24:58 raasch47 ! Leerzeilen mussten entfernt werden48 !49 ! Revision 1.2 1997/09/12 06:44:00 raasch50 ! HP-Compiler erfordert & am Beginn von Fortsetzungszeilen51 17 ! 52 18 ! Revision 1.1 1997/09/12 06:24:11 raasch -
palm/trunk/SOURCE/diffusivities.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: diffusivities.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.24 2006/04/26 12:16:26 raasch 12 14 ! OpenMP optimization (+sums_l_l_t), sqrt_e must be private 13 !14 ! Revision 1.23 2006/02/23 12:13:29 raasch15 ! nzb_2d replaced by nzb_s_inner/outer, option of a minimum tke16 ! quasi boundary conditions is now set for sums_l_l instead of sums_l17 !18 ! Revision 1.22 2005/03/26 20:16:26 raasch19 ! Eddy diffusivities are also calculated for the ghost points in order to spare20 ! the communication otherwise needed for the exchange of these points.21 ! Neumann boundary conditions set at the outflow boundaries in case of22 ! non-cyclic lateral boundaries23 !24 ! Revision 1.21 2004/01/30 10:23:12 raasch25 ! Scalar lower k index nzb replaced by 2d-array nzb_2d26 !27 ! Revision 1.20 2003/03/14 13:40:09 raasch28 ! Loop optimization (SQRT(e) solved in a seperate loop, joining of loops)29 !30 ! Revision 1.19 2003/03/12 16:28:56 raasch31 ! phi_m is initialized in declaration statement in order to avoid run time32 ! errors due to compiler problems on ibm and nec33 !34 ! Revision 1.18 2002/12/19 14:32:23 raasch35 ! Correction of mixing length term (l(k)/ll(k)). The condition kh=3*km in36 ! the unstable case is now also exactly met in the wall adjustment region.37 ! Factor 0.7 in wall adjustment part replaced by variable38 ! wall_adjustment_factor, which is set to 1.8 in modules.f90.39 ! Missing OMP DO statement added.40 !41 ! Revision 1.17 2002/06/11 12:56:49 raasch42 ! OpenMP directives added,43 ! array l changed from allocatable to automatic, because it must be on the stack44 ! in case of using OpenMP45 !46 ! Revision 1.16 2001/08/21 08:30:34 08:30:34 raasch (Siegfried Raasch)47 ! Wall adjustment of mixing length to 0.7 z can be switched off48 ! Lower boundary condition for km and kh changed to Neumann49 !50 ! Revision 1.15 2001/03/30 07:22:09 raasch51 ! Near surface mixing length is limited to 0.7*zu,52 ! if adjust_mixing_length=T, l is modified at all gridpoints (previously at53 ! k=1 only),54 ! Translation of remaining German identifiers (variables, subroutines, etc.)55 !56 ! Revision 1.14 2001/01/22 06:31:36 raasch57 ! Module test_variables removed58 !59 ! Revision 1.13 2001/01/02 17:27:41 raasch60 ! -dpt_dz_d, dpt_dz_u61 !62 ! Revision 1.12 2000/12/20 10:35:47 letzel63 ! All comments translated into English.64 !65 ! Revision 1.11 2000/04/18 08:10:42 schroeter66 ! Revision 1.4 wieder rueckgaengig gemacht, das Stabilitaets-67 ! kriterium basiert nun wieder auf zentralen Differenzen68 !69 ! Revision 1.9 2000/04/13 13:42:14 schroeter70 ! je nach Initialisierungsmodus (trocken/feucht) fliesst in71 ! Mischungswegberechnung pt oder vpt ein, wird durch entsprechende72 ! Variablenuebergabe geregelt73 !74 ! Revision 1.8 99/02/17 09:17:39 09:17:39 raasch (Siegfried Raasch)75 ! Kriterium fuer reduzierten Mischungsweg im stabil geschichteten Fall enger76 ! gefasst (Schichtung muss sowohl oberhalb als auch unterhalb des betrachteten77 ! Gitterpunkts stabil sein)78 !79 ! Revision 1.7 1998/09/22 17:18:41 raasch80 ! Testweise cm = 0.4 bei k=1, aber vorerst auskommentiert81 !82 ! Revision 1.6 1998/08/05 06:52:44 raasch83 ! Mischungsweganpassung an Prandtl nur bei k=184 !85 ! Revision 1.2 1998/03/11 11:49:26 raasch86 ! Anpassung des Mischungsweges an den Prandtlschen Mischungsweg moeglich87 15 ! 88 16 ! Revision 1.1 1997/09/19 07:41:10 raasch -
palm/trunk/SOURCE/disturb_field.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: disturb_field.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2006/08/04 14:31:59 raasch 12 14 ! izuf renamed iran 13 !14 ! Revision 1.10 2006/02/23 12:14:06 raasch15 ! Additional argument nzb_uv_inner (no perturbations below the topography)16 !17 ! Revision 1.9 2005/03/26 20:18:30 raasch18 ! Arguments for non-cyclic boundary conditions added to argument list and to19 ! argument list of routine exchange_horiz.20 ! Perturbations can be imposed for different ranges (steered by dist_range).21 !22 ! Revision 1.8 2003/10/29 08:48:51 raasch23 ! Module random_function_mod is used24 !25 ! Revision 1.7 2003/03/14 13:41:48 raasch26 ! Different random number generators available27 !28 ! Revision 1.6 2001/03/30 07:22:26 raasch29 ! Translation of remaining German identifiers (variables, subroutines, etc.)30 !31 ! Revision 1.5 2001/01/25 06:59:38 raasch32 ! Module interfaces added (missing this interface caused errors on machines33 ! other than Cray)34 !35 ! Revision 1.4 2001/01/22 06:32:49 raasch36 ! Module test_variables removed37 !38 ! Revision 1.3 2000/07/03 12:59:03 raasch39 ! All comments translated into English40 !41 ! Revision 1.2 1998/07/06 12:13:01 raasch42 ! + USE test_variables43 15 ! 44 16 ! Revision 1.1 1998/02/04 15:40:45 raasch -
palm/trunk/SOURCE/disturb_heatflux.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: disturb_heatflux.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.7 2006/08/04 14:35:07 raasch 12 14 ! Additional parameter in function random_gauss which limits the range of the 13 15 ! created random numbers, izuf renamed iran 14 !15 ! Revision 1.6 2006/02/23 12:14:38 raasch16 ! Over topography surface_heatflux is replaced by wall_heatflux(0)17 !18 ! Revision 1.5 2001/03/30 07:22:44 raasch19 ! Translation of remaining German identifiers (variables, subroutines, etc.)20 !21 ! Revision 1.4 2001/01/22 06:34:34 raasch22 ! Module test_variables removed23 !24 ! Revision 1.3 2000/12/20 11:57:38 letzel25 ! All comments translated into English.26 !27 ! Revision 1.2 1998/07/06 12:13:19 raasch28 ! + USE test_variables29 16 ! 30 17 ! Revision 1.1 1998/03/25 20:03:47 raasch -
palm/trunk/SOURCE/exchange_horiz.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: exchange_horiz.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.16 2006/02/23 12:19:08 raasch 12 14 ! anz_yz renamed ngp_yz 13 !14 ! Revision 1.15 2005/03/26 20:20:39 raasch15 ! Extensions for non-cyclic boundary conditions along x or y16 !17 ! Revision 1.14 2003/10/29 08:51:10 raasch18 ! Exchange of boundaries within the PE memory in case that multigrid method19 ! has switched to PE0 only20 !21 ! Revision 1.13 2003/03/16 09:30:38 raasch22 ! Two underscores (_) are placed in front of all define-strings23 !24 ! Revision 1.12 2002/12/20 10:40:08 raasch25 ! Integer variables will be defined in the parallel case only26 !27 ! Revision 1.11 2002/12/19 14:46:43 raasch28 ! Cyclic boundary conditions along y used instead of sendrecv in case of29 ! pdims(2)=1. SENDRECV replaced by nonblocking routines ISEND and IRECV.30 ! Array "feld" is renamed to "ar".31 !32 ! Revision 1.9 2002/04/16 08:16:50 08:16:50 raasch (Siegfried Raasch)33 ! Cyclic boundary conditions along x used instead of sendrecv in case of34 ! pdims(1)=135 !36 ! Revision 1.8 2001/07/20 13:04:59 raasch37 ! anz_yz and type_xz changed to arrays to allow for the boundary value exchange38 ! for the different grid levels used in the multigrid method,39 ! +module control_parameters40 !41 ! Revision 1.7 2001/03/30 07:23:00 raasch42 ! Translation of remaining German identifiers (variables, subroutines, etc.)43 !44 ! Revision 1.6 2001/01/22 06:36:07 raasch45 ! Module test_variables removed46 !47 ! Revision 1.5 2000/12/20 12:06:44 letzel48 ! All comments translated into English.49 !50 ! Revision 1.4 1998/07/06 12:13:36 raasch51 ! + USE test_variables52 !53 ! Revision 1.3 1998/01/23 09:57:45 raasch54 ! Umbenennung von anz_yz_p in anz_yz sowie von type_xz_p in type_xz55 !56 ! Revision 1.2 1997/08/11 06:14:55 raasch57 ! Felder werden jetzt ueber Formalparameter uebergeben.58 15 ! 59 16 ! Revision 1.1 1997/07/24 11:13:29 raasch -
palm/trunk/SOURCE/exchange_horiz_2d.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: exchange_horiz_2d.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.9 2006/05/12 19:15:52 letzel 12 14 ! MPI_REAL replaced by MPI_INTEGER in exchange_horiz_2d_int 13 !14 ! Revision 1.8 2006/02/23 12:18:32 raasch15 ! Additional subroutine exchange_horiz_2d_int for 2D integer arrays,16 ! extensions for non-cyclic boundary conditions along x or y for non-parallel17 ! case, anz_y renamed ngp_y18 !19 ! Revision 1.7 2003/03/16 09:30:43 raasch20 ! Two underscores (_) are placed in front of all define-strings21 !22 ! Revision 1.6 2002/06/11 12:59:35 raasch23 ! Cyclic boundary conditions are used instead of sendrecv in case of24 ! pdims(..)=1. Array "feld" is renamed to "ar".25 !26 ! Revision 1.5 2001/03/30 07:23:16 07:23:16 raasch (Siegfried Raasch)27 ! Translation of remaining German identifiers (variables, subroutines, etc.)28 !29 ! Revision 1.4 2001/01/22 06:43:50 raasch30 ! Module test_variables removed31 !32 ! Revision 1.3 2000/12/20 12:09:27 letzel33 ! All comments translated into English.34 !35 ! Revision 1.2 1998/07/06 12:13:53 raasch36 ! + USE test_variables37 15 ! 38 16 ! Revision 1.1 1998/01/23 09:58:21 raasch -
palm/trunk/SOURCE/fft_xy.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: fft_xy.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.4 2006/03/28 12:27:09 raasch 12 14 ! Stop when system-specific fft is selected on NEC. For unknown reasons this 13 15 ! causes a program abort during first allocation in init_grid. 14 !15 ! Revision 1.3 2006/02/23 12:22:41 raasch16 ! System specfic routines are not allowed to be used on ibmy17 16 ! 18 17 ! Revision 1.2 2004/04/30 11:44:27 raasch … … 25 24 ! decalpha), 26 25 ! non vectorized FFT for nec included 27 !28 ! Revision 1.1 2004/04/30 11:37:14 raasch29 ! Initial revision30 !31 ! Revision 1.6 2003/08/01 08:15:09 raasch32 ! No abort in fft_1dd_init on t3e-systems if system specific routines are used33 !34 ! Revision 1.5 2003/04/16 12:51:47 raasch35 ! Temperton-fft implemented in routines fft_x_1dd and fft_y_1dd36 !37 ! Revision 1.4 2003/03/16 09:37:26 raasch38 ! Two underscores (_) are placed in front of all define-strings39 !40 ! Revision 1.3 2003/03/12 16:29:39 raasch41 ! Routines fft_x_1dd_m and fft_y_1dd_m added (suitable for multiple 1d-fft on42 ! vector processors)43 !44 ! Revision 1.2 2002/12/19 14:49:18 raasch45 ! STOP statement replaced by call of subroutine local_stop46 26 ! 47 27 ! Revision 1.1 2002/06/11 13:00:49 raasch -
palm/trunk/SOURCE/flow_statistics.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: flow_statistics.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.41 2006/08/04 14:37:50 raasch 12 14 ! Error removed in non-parallel part (sums_l) 13 !14 ! Revision 1.40 2006/04/26 14:06:01 raasch15 ! Call of omp_get_thread_num written as comment line. Since lcmuk does not16 ! understand this comment, a cpp-directive is required.17 !18 ! Revision 1.39 2006/04/26 12:19:16 raasch19 ! OpenMP optimization, thread number added as third dimension of sums_l,20 ! nzb_s_inner(j,i) replaced by nzb in determination of z_i21 !22 ! Revision 1.38 2006/02/23 12:23:15 raasch23 ! nzb and nzb_diff replaced by nzb_s_outer and nzb_diff_s_outer resrectively,24 ! nanz_sums, nanz_2dh, nanz_3d renamed ngp_sums, ngp_2dh and ngp_3d,25 ! respectively,26 ! ngp_2dh and ngp_3d partially replaced by ngp_2dh_outer and ngp_3d_inner,27 ! respectively28 !29 ! Revision 1.37 2005/12/06 16:40:21 raasch30 ! Calculation of ql (54) profile in case of using cloud droplets31 !32 ! Revision 1.36 2005/04/23 09:10:27 raasch33 ! +local array sums_ll instead of temporarily using sums_l34 !35 ! Revision 1.35 2004/04/30 11:46:56 raasch36 ! leap_frog changed to time_integration37 !38 ! Revision 1.34 2003/05/09 14:25:12 raasch39 ! Parameters in MPI_ALLREDUCE call changed40 !41 ! Revision 1.33 2003/03/31 15:42:53 raasch42 ! Calculation of horizontal temperature fluxes modified and cleaned up43 !44 ! Revision 1.32 2003/03/16 09:39:07 raasch45 ! Two underscores (_) are placed in front of all define-strings46 !47 ! Revision 1.31 2003/03/12 16:30:07 raasch48 ! Calculation of perturbation and total energy of the total domain49 ! re-arranged in order to allow vectorization of the respective loop50 !51 ! Revision 1.30 2003/03/04 11:26:35 raasch52 ! Error in calculation of the vertical flux of resolved scale energy53 ! (profile 57) removed54 !55 ! Revision 1.29 2002/12/19 14:50:02 raasch56 ! STOP statement replaced by call of subroutine local_stop57 !58 ! Revision 1.28 2002/04/16 08:04:19 raasch59 ! Error in passive scalar profile (hom(..41..)) removed60 !61 ! Revision 1.27 2002/04/11 16:31:03 schroeter62 ! No changes.63 !64 ! Revision 1.26 2001/09/04 12:00:12 12:00:12 raasch (Siegfried Raasch)65 ! Calculation of the vertical flux of resolved scale energy and pressure66 ! fluctuations67 !68 ! Revision 1.25 2001/08/21 08:32:24 raasch69 ! Special calculation of surface fluxes whenever they are prescribed (not only70 ! in case of a Prandtl layer), call of user_interface71 !72 ! Revision 1.24 2001/03/30 07:23:34 raasch73 ! Translation of remaining German identifiers (variables, subroutines, etc.)74 !75 ! Revision 1.23 2001/02/07 06:48:25 raasch76 ! Error in non-parallel passive scalar code removed77 !78 ! Revision 1.22 2001/01/29 12:25:02 raasch79 ! Passive scalar is considered80 !81 ! Revision 1.21 2001/01/22 06:49:27 raasch82 ! Module test_variables removed83 !84 ! Revision 1.20 2000/04/13 13:54:15 schroeter85 ! + computation of new profiles for calculations with86 ! cloud_physics87 !88 ! Revision 1.19 2000/01/26 09:32:58 09:32:58 letzel (Marcus Letzel)89 ! All comments translated into English90 !91 ! Revision 1.16 1998/04/21 15:54:54 raasch92 ! Berechnung von w* beruecksichtigt jetzt den Fall z_i=0,93 ! Faktor 0.5 bei Energietransport hinzugefuegt94 !95 ! Revision 1.6 1998/02/10 15:06:02 raasch96 ! Berechnung des Gesamtwaermestroms, in aufgeloesten Waermestroemen wird das97 ! horizontale Temperaturmittel nicht mehr abgezogen, Berechnung der98 ! Gesamtimpulsfluesse99 15 ! 100 16 ! Revision 1.1 1997/08/11 06:15:17 raasch -
palm/trunk/SOURCE/global_min_max.f90
r1 r3 9 9 ! Former revisions: 10 10 ! ----------------- 11 ! $Log: global_min_max.f90,v $ 11 ! $Id$ 12 ! RCS Log replace by Id keyword, revision history cleaned up 13 ! 12 14 ! Revision 1.11 2003/04/16 12:56:58 raasch 13 15 ! Index values of the extrema are limited to the range 0..nx, 0..ny 14 !15 ! Revision 1.10 2003/03/16 09:39:21 raasch16 ! Two underscores (_) are placed in front of all define-strings17 !18 ! Revision 1.9 2003/03/14 13:42:50 raasch19 ! IF clause extracted from loop which determines the absolut array maximum,20 ! because otherwise it prevents vectorization of that loop21 !22 ! Revision 1.8 2002/12/19 14:51:03 raasch23 ! Speed optimization by removing MINVAL/MAXVAL calls and by handling24 ! the "abs" case in a different way than the min/max cases,25 ! translation of remaining variables26 !27 ! Revision 1.7 2002/06/11 13:08:16 raasch28 ! Last change limited to IBM since it does not work on 64-bit machines29 !30 ! Revision 1.6 2002/05/02 18:50:36 18:50:36 raasch (Siegfried Raasch)31 ! Type of fmax, fmax_l, fmin and fmin_l changed to REAL (KIND=4).32 !33 ! Revision 1.5 2001/01/22 06:51:02 raasch34 ! Module test_variables removed35 !36 ! Revision 1.4 2000/12/20 12:19:57 letzel37 ! All comments translated into English.38 !39 ! Revision 1.3 1998/07/06 12:14:37 raasch40 ! + USE test_variables41 !42 ! Revision 1.2 1997/08/11 06:16:54 raasch43 ! Indices des ermittelten Minumums/Maximums werden an PEs gesendet44 16 ! 45 17 ! Revision 1.1 1997/07/24 11:14:03 raasch -
palm/trunk/SOURCE/header.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: header.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.63 2006/08/22 13:53:13 raasch 12 14 ! Output of dz_max 13 !14 ! Revision 1.62 2006/08/04 14:38:41 raasch15 ! Generation of the run description header as well as run date and run time16 ! moved to routine check_parameters,17 ! output for additional particle quantities18 !19 ! Revision 1.61 2006/04/26 12:19:51 raasch20 ! Output of threads per task21 !22 ! Revision 1.60 2006/04/11 14:57:01 raasch23 ! pl_spectra renamed data_output_sp24 !25 ! Revision 1.59 2006/03/14 12:44:10 raasch26 ! Output of the initial geostrophic wind profile corrected27 !28 ! Revision 1.58 2006/02/23 12:24:22 raasch29 ! Output of psl, psr, pdx, etc. for each particle group, output of e_min,30 ! information about volume flow conservation and usage of NetCDF 64 bit offset31 ! format, pl.. renamed do.., skip_time control parameters for do..,32 ! ebene renamed section, .._anz renamed .._no, dt_average renamed33 ! dt_averaging_input_pr34 ! output of 2d/3d data averaging informations, buffer is forced to be written35 ! on file at the end of the routine36 !37 ! Revision 1.57 2005/10/20 14:06:44 raasch38 ! Error in output of yz-slice information removed39 !40 ! Revision 1.56 2005/07/01 07:48:14 steinfeld41 ! Information on the initial profiles of ug and vg added;42 ! dependency of ug and vg on height considered in the output of43 ! the boundary conditions for u and v44 !45 ! Revision 1.55 2005/06/26 19:53:15 raasch46 ! gas_constant renamed r_d, latent_heat renamed l_v, radius is used instead47 ! of diameter, output formats 335+337 changed48 !49 ! Revision 1.54 2005/05/19 08:32:27 raasch50 ! Error in output of particle inertia information removed51 !52 ! Revision 1.53 2005/05/18 15:31:59 raasch53 ! Informations about data output format (netcdf, iso2d, etc.) are given54 !55 ! Revision 1.52 2005/04/23 09:16:10 raasch56 ! fcl_factor renamed cfl_factor57 !58 ! Revision 1.51 2005/03/26 20:23:14 raasch59 ! Output of horizontal boundary conditions, output of additional call of60 ! pressure solver in case of Runge-Kutta schemes,61 ! calculate the number of particle groups (formerly unknown = 0 )62 !63 ! Revision 1.50 2004/04/30 11:48:02 raasch64 ! Forcing a 1d-decomposition is output on linux-cluster and decalpha,65 ! output of subdomain sizes, impulse_advec renamed momentum_advec66 !67 ! Revision 1.49 2004/01/28 15:09:15 raasch68 ! Output of timestep schemes modified69 !70 ! Revision 1.48 2003/10/29 08:52:04 raasch71 ! Additional output for multigrid method (lowest levels gathered on PE0),72 ! particle output modified for particle groups73 !74 ! Revision 1.47 2003/04/16 12:58:18 raasch75 ! Output for mixing length limitations revised76 !77 ! Revision 1.46 2003/03/16 09:39:31 raasch78 ! Two underscores (_) are placed in front of all define-strings79 !80 ! Revision 1.45 2003/03/14 13:43:21 raasch81 ! Informations about random generator82 !83 ! Revision 1.44 2003/03/12 16:31:45 raasch84 ! Small change for NEC system85 !86 ! Revision 1.43 2002/12/19 14:53:47 raasch87 ! Informations about user defined restart times added88 !89 ! Revision 1.42 2002/09/12 13:03:24 raasch90 ! Informations about particle inertia added91 !92 ! Revision 1.41 2002/06/11 13:09:11 raasch93 ! Information about hybrid solver, usage of 1d-decomposition and fast version94 ! of prognostic_equations included,95 ! no output of mixing length adjustments in case of constant diffusion96 !97 ! Revision 1.40 2002/04/16 08:06:13 08:06:13 raasch (Siegfried Raasch)98 ! Informations about particle boundary conditions and output of particle99 ! data added, direct output of informations about scalar transport (instead100 ! of using humidity)101 !102 ! Revision 1.39 2001/08/21 08:50:54 raasch103 ! Informations about mixing length adjustment, particle tails and dvrp PE usage104 ! included105 !106 ! Revision 1.38 2001/07/20 13:05:39 raasch107 ! Information about multigrid method included108 !109 ! Revision 1.37 2001/07/12 12:08:14 raasch110 ! Additional informations about dvrp-output (host and filenames) and111 ! particles (maximum age, random start positions), + module particle_attributes112 !113 ! Revision 1.36 2001/03/30 07:23:53 raasch114 ! Translation of remaining German identifiers (variables, subroutines, etc.)115 !116 ! Revision 1.35 2001/01/29 12:25:55 raasch117 ! Informations about using passive scalar118 !119 ! Revision 1.34 2001/01/25 07:00:48 raasch120 ! Information about using surface fluxes and fft-method added121 !122 ! Revision 1.33 2001/01/05 15:11:41 raasch123 ! Spectra informations added. Old revision remarks deleted.124 !125 ! Revision 1.32 2001/01/03 13:08:47 letzel126 ! File output translated into English.127 !128 ! Revision 1.31 2000/12/28 13:04:54 raasch129 ! Slicer information added to dvrp-output.130 ! Informations about optionally used packages are put in cpp-define-brackets.131 !132 ! Revision 1.30 2000/04/27 07:08:14 raasch133 ! notice if processor topology is predefined by user,134 ! general notice about dvrp-output (including isosurface and particles),135 ! all comments translated into English136 15 ! 137 16 ! Revision 1.1 1997/08/11 06:17:20 raasch -
palm/trunk/SOURCE/impact_of_latent_heat.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: impact_of_latent_heat.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2004/01/30 10:25:59 raasch 12 14 ! Scalar lower k index nzb replaced by 2d-array nzb_2d 13 !14 ! Revision 1.4 2003/03/12 16:32:14 raasch15 ! Full code replaced in the call for all gridpoints instead of calling the16 ! _ij version (required by NEC, because otherwise no vectorization)17 !18 ! Revision 1.3 2002/06/11 13:10:06 raasch19 ! Former subroutine changed to a module which allows to be called for all grid20 ! points of a single vertical column with index i,j or for all grid points by21 ! using function overloading.22 !23 ! Revision 1.2 2001/01/22 06:52:34 raasch24 ! Module test_variables removed25 15 ! 26 16 ! Revision 1.1 2000/04/13 14:48:40 schroeter -
palm/trunk/SOURCE/init_1d_model.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_1d_model.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.21 2006/06/02 15:19:57 raasch 12 14 ! cpp-directives extended for lctit 13 15 ! 14 ! Revision 1.20 2006/04/20 13:16:24 letzel15 ! Errors removed in the computation of the diabatic mixing length16 !17 ! Revision 1.19 2006/03/14 12:48:52 raasch18 ! Error in initializing l_black removed,19 ! maximum possible timestep for runge schemes reduced, output for units 15 and20 ! 17 is flushed21 !22 ! Revision 1.18 2006/02/23 12:27:38 raasch23 ! Switches for calculating the mixing length and the dissipation,24 ! remaining e1d_m**1.5 replaced by e1d_m*SQRT(e1d_m),25 ! nt_anz_1d renamed current_timestep_number_1d26 !27 ! Revision 1.17 2005/06/29 09:58:25 steinfeld28 ! A potential dependency of ug and vg on height is considered in the29 ! tendency terms of u and v; calculation of lambda is commented out30 !31 ! Revision 1.16 2005/04/23 09:22:04 raasch32 ! Initial horizontal velocities at the lowest grid levels are set to a very33 ! small value in order to avoid too small time steps caused by the diffusion34 ! criterion35 !36 ! Revision 1.15 2005/03/26 20:27:33 raasch37 ! Output format of timestep changed, e1d_m**1.5 replaced by e1d_m*SQRT(e1d_m)38 ! because of wrong results on NEC-SX639 !40 ! Revision 1.14 2004/04/30 11:55:02 raasch41 ! Runge-Kutta scheme implemented, routine leap_frog_1d renamed to42 ! time_integration_1d43 !44 ! Revision 1.13 2001/03/30 07:24:10 raasch45 ! Translation of remaining German identifiers (variables, subroutines, etc.)46 !47 ! Revision 1.12 2001/01/29 12:26:57 raasch48 ! Passive scalar is considered49 !50 ! Revision 1.11 2001/01/22 06:54:13 raasch51 ! Module test_variables removed52 !53 ! Revision 1.10 2001/01/02 17:28:35 raasch54 ! Adjustment of computations of energy dissipation, mixing length and55 ! eddy diffusivities to methods used in the 3D-model:56 ! cm changed from 0.4 to 0.1, grid dependend mixing length instead of57 ! Blackadar mixing length is used, Deardorff modification of mxixing length58 ! when stratification is stable59 !60 ! Revision 1.9 2000/12/21 14:16:15 letzel61 ! All comments translated into English.62 !63 ! Revision 1.8 2000/04/13 13:23:47 schroeter64 ! + 1D-Modell an Rechnungen mit Feuchte angepasst65 !66 ! Revision 1.7 2000/03/14 12:56:02 12:56:02 schroeter (Michael Schroeter)67 ! Asselin-filter is called before swapping timelevels68 !69 16 ! Revision 1.1 1998/03/09 16:22:10 raasch 70 17 ! Initial revision -
palm/trunk/SOURCE/init_3d_model.f90
r1 r3 11 11 ! Former revisions: 12 12 ! ----------------- 13 ! $Log: init_3d_model.f90,v $ 13 ! $Id$ 14 ! RCS Log replace by Id keyword, revision history cleaned up 15 ! 14 16 ! Revision 1.49 2006/08/22 15:59:07 raasch 15 17 ! No optimization of this file on the ibmy (Yonsei Univ.) 16 !17 ! Revision 1.48 2006/08/22 13:58:30 raasch18 ! Initialisation of all tendency arrays (t.._m) needed for the Runge-Kutta19 ! schemes20 !21 ! Revision 1.47 2006/08/04 14:39:20 raasch22 ! Allocation of diss, needed for calculating the sgs particle velocities23 !24 ! Revision 1.46 2006/05/15 08:14:56 letzel25 ! Initialize volume_flow_initial_l and volume_flow_area_l with 0.026 !27 ! Revision 1.45 2006/04/26 13:32:44 raasch28 ! Provisional correction for Piacsek & Williams advection scheme: keep u and29 ! v zero one layer below the topography (in case of set_1d-model_profiles).30 !31 ! Revision 1.44 2006/04/26 12:22:09 raasch32 ! sums_l and sums_l_l changed to three dimensions33 !34 ! Revision 1.43 2006/02/23 12:29:11 raasch35 ! Allocation and default setting of k index arrays nzb_2d and nzb_diff_2d36 ! moved to routine init_grid, nanz_sums, nanz_2dh and nanz_3d renamed ngp_sums,37 ! ngp_2dh and ngp_3d, respectively,38 ! allocation and setting of ngp_2dh and ngp_3d moved from parin and init_pegrid39 ! to here,40 ! read_fields_from_prior_run renamed read_restart_data,41 ! new grid point counters ngp_2dh_outer and ngp_3d_inner are calculated,42 ! additional argument nzb_.._inner added to call of disturb_field43 ! optional calculation of the initial volume flow through the right and north44 ! boundary45 !46 ! Revision 1.42 2005/06/26 19:55:01 raasch47 ! ql is now a pointer, init_cloud_physics is called anyway48 !49 ! Revision 1.41 2005/03/26 20:29:41 raasch50 ! Arguments for non-cyclic boundary conditions added in call of routine51 ! disturb_field, diffusivities km_damp_x and km_damp_y used within the outflow52 ! damping layer are allocated and initialized, uvmean_outflow is allocated53 !54 ! Revision 1.40 2004/04/30 11:59:05 raasch55 ! Reading of restart data is moved to new routine read_3d_binary,56 ! diffusion quantity arrays are allocated for previous timelevels only57 ! in case of leapfrog scheme, enlarged transposition array d58 !59 ! Revision 1.39 2004/01/30 10:27:20 raasch60 ! Allocation and default setting of index arrays nzb_2d and nzb_diff_2d61 !62 ! Revision 1.38 2004/01/28 15:09:56 raasch63 ! Time levels for Runge-Kutta scheme initialized, setting of bt before64 ! calling pres removed65 !66 ! Revision 1.37 2003/10/29 08:54:05 raasch67 ! Random numbers are initialized or random seed is read from file of68 ! previous run, module random_function_mod is used, new binary version 2.269 !70 ! Revision 1.36 2003/08/01 08:25:08 raasch71 ! Array sums is initialized (must be defined in the first call of pres),72 ! call of init_cloud_physics moved before init_particles73 !74 ! Revision 1.35 2003/04/09 11:50:50 raasch75 ! In case of initialisation with constant profiles, the horizontal76 ! velocities at k=0,1 are explicitly set to zero77 !78 ! Revision 1.34 2003/03/16 09:39:39 raasch79 ! Two underscores (_) are placed in front of all define-strings80 !81 ! Revision 1.33 2002/12/19 15:32:42 raasch82 ! Subroutine local_tremain_ini (initializing the remaining cpu-time)83 ! moved to main program (palm), reading of array hom is moved to routine84 ! read_var_list, STOP statement replaced by call of subroutine local_stop85 !86 ! Revision 1.32 2002/05/02 18:51:11 raasch87 ! Time level t+dt re-introduced and initialized88 !89 ! Revision 1.31 2001/11/05 09:31:38 09:31:38 weinbrec (Sonja Weinbrecht)90 ! Initialization of spectrum_x and spectrum_y for continuation runs91 !92 ! Revision 1.30 2001/03/30 07:25:25 raasch93 ! Arguments removed from pres, arrays work and work1 removed,94 ! initialization of timelevel t+dt removed,95 ! Translation of remaining German identifiers (variables, subroutines, etc.)96 !97 ! Revision 1.29 2001/01/29 12:27:49 raasch98 ! Passive scalar is considered99 !100 ! Revision 1.28 2001/01/25 07:02:19 raasch101 ! Initialization of surface fluxes changed due to new variable102 ! use_surface_fluxes103 !104 ! Revision 1.27 2001/01/05 15:12:23 raasch105 ! Allocate spectrum_x and spectrum_y, new binary version 2.1106 !107 ! Revision 1.26 2001/01/02 17:29:43 raasch108 ! Unit 13 closed after reading binary input (after call of user_init, where109 ! read from unit 13 may occur)110 !111 ! Revision 1.25 2000/12/28 13:05:53 raasch112 ! Init_particles is called unconditionally113 !114 ! Revision 1.24 2000/04/27 07:05:19 raasch115 ! Call of subroutine init_slope added to compute the initial temperature field116 ! and other constants used in case of a sloping surface,117 ! call of subroutine init_dvrp added to initialize dvrp-software,118 ! old revision remarks deleted119 !120 ! Revision 1.23 2000/04/13 13:22:09 schroeter121 ! Allokieren und Initialisierung der Feuchtefelder und Wolkenphysikparameter122 !123 ! Revision 1.22 2000/01/26 10:30:49 10:30:49 letzel (Marcus Letzel)124 ! All comments translated into English125 !126 ! Revision 1.21 2000/01/10 09:54:10 raasch127 ! Trennung der zeitabhaengigen Variablen in pointer und targets,128 ! Anfangszuordnung der pointer129 !130 ! Revision 1.17 1999/03/03 07:22:35 raasch131 ! Beim binaeren Einlesen werden alle Variablen mit exakter Laenge von Datei132 ! gelesen, da sich bei schneller I/O Records nicht einfach ueberspringen lassen133 !134 ! Revision 1.14 1998/12/15 07:33:50 raasch135 ! Zusaetzliche Ermittelung der tatsaechlichen Einlesezeiten auf T3E136 ! mit messpunkt_e(14) unter Verwendung von TIMEF (s. cpu_zeitmessung)137 18 ! 138 19 ! Revision 1.1 1998/03/09 16:22:22 raasch -
palm/trunk/SOURCE/init_advec.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_advec.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.6 2004/04/30 11:59:31 raasch 12 14 ! impulse_advec renamed momentum_advec 13 !14 ! Revision 1.5 2001/03/30 07:25:53 raasch15 ! Translation of remaining German identifiers (variables, subroutines, etc.)16 !17 ! Revision 1.4 2001/01/22 06:57:43 raasch18 ! Module test_variables removed19 !20 ! Revision 1.3 2000/12/21 16:12:52 letzel21 ! All comments translated into English.22 !23 ! Revision 1.2 1999/02/17 09:18:45 raasch24 ! Fehlerkorrektur, Verfahren muessen unabhaengig voneinander initialisierbar25 ! sein26 15 ! 27 16 ! Revision 1.1 1999/02/05 09:07:38 raasch -
palm/trunk/SOURCE/init_cloud_physics.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ------------------ 10 ! $Log: init_cloud_physics.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2005/06/26 19:55:58 raasch 12 14 ! Initialization of cloud droplet constants, gas_constant renamed r_d, 13 15 ! latent_heat renamed l_v 14 !15 ! Revision 1.4 2001/03/30 07:26:30 raasch16 ! Translation of remaining German identifiers (variables, subroutines, etc.),17 ! surface_pressure keeps unit hPa and is not converted to Pa18 !19 ! Revision 1.3 2001/01/25 07:03:10 raasch20 ! Module test_variables removed21 !22 ! Revision 1.2 2001/01/22 09:16:57 schroeter23 ! To calculate verticle pressure profile use actual surface-temperature in24 ! place of potential temperature as reference value.25 16 ! 26 17 ! Revision 1.1 2000/04/13 14:37:22 schroeter -
palm/trunk/SOURCE/init_dvrp.f90
r1 r3 4 4 ! Actual revisions: 5 5 ! ----------------- 6 !7 6 ! Test output of isosurface on camera file 8 7 ! TEST: print* statements … … 11 10 ! Former revisions: 12 11 ! ----------------- 13 ! $Log: init_dvrp.f90,v $ 12 ! $Id$ 13 ! RCS Log replace by Id keyword, revision history cleaned up 14 ! 14 15 ! Revision 1.12 2006/02/23 12:30:22 raasch 15 16 ! ebene renamed section, pl.. replaced by do.., 16 !17 ! Revision 1.11 2005/03/31 07:01:20 raasch18 ! Test version with print statements checked in19 !20 ! Revision 1.10 2003/08/01 08:36:28 raasch21 ! Module dvrp added to routine close_dvrp, set the maximum time the program22 ! can be suspended on user request, every steering variable is defined for23 ! every stream24 !25 ! Revision 1.9 2003/03/16 09:39:44 raasch26 ! Two underscores (_) are placed in front of all define-strings27 !28 ! Revision 1.8 2003/03/04 11:28:53 raasch29 ! Update for new version of dvrp software (e.g. using module dvrp is now30 ! mandatory). Initialization of dvrp steering variables included.31 !32 ! Revision 1.7 2002/12/19 15:37:57 raasch33 ! STOP statement replaced by call of subroutine local_stop34 !35 ! Revision 1.6 2001/08/21 09:45:54 raasch36 ! Superelevation implemented, additional PE can be used for handling the37 ! dvrp output (see init_dvrp_logging and close_dvrp)38 !39 ! Revision 1.5 2001/07/12 12:10:16 raasch40 ! Adjustment of camera position, output of boundingbox moved from the scene41 ! files to camera.dvr, destination for dvrp_output is now controlled via42 ! parameters given in the parameter file.43 !44 ! Revision 1.4 2001/03/30 07:31:03 raasch45 ! Translation of remaining German identifiers (variables, subroutines, etc.)46 !47 ! Revision 1.3 2001/01/25 07:04:23 raasch48 ! Subroutines init_dvrp_logging and close_dvrp attached.49 !50 ! Revision 1.2 2000/12/28 13:11:43 raasch51 ! Complete revision due to new features of dvrp-software52 17 ! 53 18 ! Revision 1.1 2000/04/27 06:24:39 raasch -
palm/trunk/SOURCE/init_grid.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_grid.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.17 2006/08/22 14:00:05 raasch 12 14 ! +dz_max to limit vertical stretching, 13 15 ! bugfix in index array initialization for line- or point-like topography 14 16 ! structures 15 !16 ! Revision 1.16 2006/03/03 14:28:31 letzel17 ! call of user_init_grid now requires the actual parameter nzb_local18 !19 ! Revision 1.15 2006/02/24 16:05:59 raasch20 ! Error removed: l_wall modification (multiplication with wall adjustment21 ! factor and ghost point exchange) must be taken out of the if-clause (see end22 ! of this routine)23 !24 ! Revision 1.14 2006/02/23 12:30:59 raasch25 ! Definition of index and wall arrays describing topography26 !27 ! Revision 1.13 2002/12/19 15:38:38 raasch28 ! Calculation of dl_grid removed, STOP statement replaced by call of subroutine29 ! local_stop30 !31 ! Revision 1.12 2001/07/20 13:06:33 raasch32 ! +dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg33 !34 ! Revision 1.11 2001/03/30 07:31:20 raasch35 ! Translation of remaining German identifiers (variables, subroutines, etc.)36 !37 ! Revision 1.10 2001/01/22 19:35:19 letzel38 ! All comments translated into English.39 !40 ! Revision 1.9 2001/01/22 07:06:16 raasch41 ! Module test_variables removed42 !43 ! Revision 1.8 1999/02/05 09:08:32 raasch44 ! Erweiterung von dzw um einen Gitterpunkt nach oben fuer Upstream-Spline-45 ! Verfahren46 !47 ! Revision 1.7 1998/07/06 12:16:43 raasch48 ! + USE test_variables49 !50 ! Revision 1.6 1998/03/18 20:13:16 raasch51 ! dz wird unveraendert gelassen (hatte bisher den Wert der groessten,52 ! gestreckten) Gitterweite53 !54 ! Revision 1.5 1997/09/19 07:43:27 raasch55 ! +dd2zu, dl_grid, l_grid56 !57 ! Revision 1.4 1997/09/12 06:26:02 raasch58 ! Kehrwerte der Gitterweitenquadrate werden berechnet59 !60 ! Revision 1.3 1997/09/09 08:29:12 raasch61 ! Kehrwerte der Gitterweiten werden berechnet62 !63 ! Revision 1.2 1997/08/26 06:31:11 raasch64 ! Berechnung der Vertikalniveaus und der Gitterweiten65 17 ! 66 18 ! Revision 1.1 1997/08/11 06:17:45 raasch -
palm/trunk/SOURCE/init_particles.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_particles.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.24 2007/02/11 13:00:17 raasch 12 14 ! Bugfix: allocation of tail_mask and new_tail_id in case of restart-runs 13 15 ! Bugfix: __ was missing in a cpp-directive 14 !15 ! Revision 1.23 2006/08/22 14:11:08 raasch16 ! Default settings of particle start positions changed17 !18 ! Revision 1.22 2006/08/04 14:44:58 raasch19 ! New particles attributes initialized (e_m, dt_sum, speed_x/y/z_sgs),20 ! total length of mpi_particle_type adjusted, more than one particle can be21 ! released per point, izuf renamed iran, set a special seed value iran_part22 ! to be used for the particle advection code, determination of particle group23 ! number removed from header to here24 !25 ! Revision 1.21 2006/03/14 12:55:24 raasch26 ! Determination of the number of particle groups removed27 !28 ! Revision 1.20 2006/02/23 12:34:00 raasch29 ! Allocation of prt_start_index and prt_count + calling subroutine30 ! sort_particles,31 ! particle source parameters psl, psr, pss, psn, psb, pst, pdx, pdy, pdz are now32 ! 1D arrays (1:max_number_of_particle_groups) and allow up to specify33 ! max_number_of_particle_groups different sources by setting34 ! number_of_particle_groups >= 1 even if density_ratio is not specified,35 ! only a fraction of the particles may have tails, therefore input from unit 9036 ! has changed,37 ! 0.5 replaced by 0.4999999999 in setting the random fluctuations,38 ! Defaults of psl, psr, pss, psn adjusted.39 ! If the number of particle groups is set by the user, density_ratio is reset to40 ! default.41 ! In case of user-specified pdx, pdy, pdz = 0.0 reset them to defaults,42 ! improve particle release at PE boundaries, nt_anz renamed43 ! current_timestep_number44 !45 ! Revision 1.19 2005/06/26 19:56:47 raasch46 ! particle_groups%radius is used instead of diameter, particle%radius is47 ! initialized, new initial informations about radius and weighting factor48 ! stored on array particles, total length of mpi_particle_type adjusted49 !50 ! Revision 1.18 2004/04/30 12:03:28 raasch51 ! MPI data type for particles modified due to modification of FORTRAN particle52 ! type53 !54 ! Revision 1.17 2003/10/29 08:55:35 raasch55 ! Module random_function_mod is used, version check of the binary file,56 ! modifications for new particle group feature57 !58 ! Revision 1.16 2003/03/16 09:39:51 raasch59 ! Two underscores (_) are placed in front of all define-strings60 !61 ! Revision 1.15 2003/03/04 11:30:48 raasch62 ! Displacement for integers in mpi_particle_type reduced from 16 to 8 on63 ! ibm, since integers seem to be only 4 byte long. Accordingly, total length64 ! reduced from 104 to 96.65 !66 ! Revision 1.14 2002/12/19 15:39:48 raasch67 ! STOP statement replaced by call of subroutine local_stop68 !69 ! Revision 1.13 2002/09/12 13:04:13 raasch70 ! Read density_ratio from restart file71 !72 ! Revision 1.12 2002/04/16 08:08:34 raasch73 ! New initial informations about speed components and start positions stored on74 ! array particles. Check boundary conditions and read them from restart file.75 !76 ! Revision 1.11 2001/11/12 16:06:06 raasch77 ! Array particle_tail_coordinates extended to contain colour informations,78 ! reading particle data from restart file can be prohibited79 !80 ! Revision 1.9 2001/08/21 09:49:05 raasch81 ! Particle tail is initialized82 !83 ! Revision 1.8 2001/07/12 12:10:58 raasch84 ! Particle sources (starting positions) are now controlled by parameters85 !86 ! Revision 1.7 2001/03/30 07:31:37 raasch87 ! Translation of remaining German identifiers (variables, subroutines, etc.)88 !89 ! Revision 1.6 2001/01/25 07:05:36 raasch90 ! Writing of particle informations is optional now91 !92 ! Revision 1.5 2001/01/22 07:16:17 raasch93 ! Module test_particles removed94 !95 ! Revision 1.4 2001/01/02 17:30:16 raasch96 ! Opening of unit 80 moved to check_open. Unit 90 instead of 81 is used97 ! for binary input.98 !99 ! Revision 1.3 2000/12/28 13:15:58 raasch100 ! Complete revision due to new features of dvrp-software,101 ! all comments translated into English,102 ! code is used only optionally (cpp-directives are added)103 !104 ! Revision 1.2 2000/04/27 06:38:48 raasch105 ! compute position of camera with new indepent routine vtk_camera and106 ! store result on a separate file107 16 ! 108 17 ! Revision 1.1 1999/11/25 16:22:38 raasch -
palm/trunk/SOURCE/init_pegrid.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_pegrid.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.28 2006/04/26 13:23:32 raasch 12 14 ! lcmuk does not understand the !$ comment so a cpp-directive is required 13 !14 ! Revision 1.27 2006/04/26 12:50:32 raasch15 ! Get the number of OpenMP threads16 !17 ! Revision 1.26 2006/02/23 12:34:57 raasch18 ! Calculation of nanz_sums, nanz_2dh and nanz_3d moved to init_3d_model,19 ! anz_y.. renamed ngp_y.., + MPI type type_x_int20 !21 ! Revision 1.25 2005/10/20 15:37:36 raasch22 ! error removed: the last PE in a row is not allowed to have more grid points23 ! than the other PES (only less)24 !25 ! Revision 1.24 2005/03/26 20:35:06 raasch26 ! Extensions for non-cyclic boundary conditions27 ! MIN function must be used instead of MAX for deciding if a switch to PE0 is28 ! useful (multigrid-method), maximuma_grid_level_l is defined in case that a29 ! pre-defined value for mg_switch_to_pe0_level is given (caused errors on NEC)30 !31 ! Revision 1.23 2004/04/30 12:27:13 raasch32 ! Setting 1d-decomposition as default on decalpha, indices for direct33 ! transpositions z --> y and y --> x added, +module poisfft_mod,34 ! data types for mpi_alltoallv removed (not neeeded any more),35 ! poisfft_hybrid is no longer enforced in case of a 1d-decomposition along x,36 ! an enlarged grid which will fit all transposition demands is calculated,37 ! impulse_advec renamed momentum_advec38 !39 ! Revision 1.22 2003/10/29 08:56:20 raasch40 ! Number of possible grid levels for the mulitgrid-method are calculated41 ! based on the total domain grid42 !43 ! Revision 1.21 2003/08/01 08:40:03 raasch44 ! Error removed: for a 1D-decomposition along x, the number PEs must also45 ! be an integral divisor of the number of gridpoints along y46 !47 ! Revision 1.20 2003/05/09 14:27:26 raasch48 ! Setting 1d-decomposition as default on linux clusters49 !50 ! Revision 1.19 2003/04/16 13:00:11 raasch51 ! Call of fft_1dd_init, +module fft_for_1d_decomp52 !53 ! Revision 1.18 2003/03/16 09:39:57 raasch54 ! Two underscores (_) are placed in front of all define-strings55 !56 ! Revision 1.17 2003/03/12 16:32:47 raasch57 ! A 1d-decomposition along x is the default on NEC hosts58 !59 ! Revision 1.16 2002/12/19 15:40:39 raasch60 ! In case of 1d-decomposition along x only a part of the integral divisor61 ! conditions is checked, on IBM hosts a 1d-decomposition along x is the62 ! default, a switch from Poisson-FFT-solver to the hybrid-solver is made63 ! in case of 1d-decomposition along x, the hybrid-solver does not force a64 ! 1d-decomposition any more, STOP replaced by call of local_stop65 !66 ! Revision 1.15 2002/06/11 13:10:58 raasch67 ! Usage of the Poisson-hybrid-solver forces a 1d-decomposition,68 ! initialization of hybrid-solver69 !70 ! Revision 1.14 2001/08/21 09:49:45 09:49:45 raasch (Siegfried Raasch)71 ! MPI_COMM_WORLD replaced by comm_palm72 !73 ! Revision 1.13 2001/07/20 13:07:26 raasch74 ! In case of using multigrid method, grid levels and additional MPI derived data75 ! types are calculated. anz_yz and type_xz are now arrays.76 !77 ! Revision 1.12 2001/03/30 07:32:01 raasch78 ! Calculations of nw and nw1 removed,79 ! Translation of remaining German identifiers (variables, subroutines, etc.)80 !81 ! Revision 1.11 2001/01/22 08:33:48 raasch82 ! Module test_variables removed83 !84 ! Revision 1.9 2000/12/28 13:16:37 raasch85 ! definition of mpi_particle_type moved to subroutine init_particles86 !87 ! Revision 1.8 2000/04/27 06:54:03 raasch88 ! Processor topology can now be prescribed by the user using npex and npey,89 ! all comments translated into English90 !91 ! Revision 1.7 1999/11/25 16:23:36 raasch92 ! MPI-Datentyp fuer Partikeltransfer vereinbart93 !94 ! Revision 1.6 1998/07/16 06:50:23 raasch95 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert96 !97 ! Revision 1.5 1998/07/06 12:17:06 raasch98 ! + USE test_variables99 !100 ! Revision 1.4 1998/03/18 20:13:42 raasch101 ! myid_char auf 5 Zeichen erweitert102 !103 ! Revision 1.3 1998/01/23 10:02:08 raasch104 ! Gitterpunktzahlen und Datentyp fuer 2D-Randaustausch105 ! anz_yz_p, type_xz_p in anz_yz und type_xz umbenannt106 !107 ! Revision 1.2 1997/08/11 06:18:28 raasch108 ! Gitterpunktanzahl von globalen Feldern werden berechnet109 15 ! 110 16 ! Revision 1.1 1997/07/24 11:15:09 raasch -
palm/trunk/SOURCE/init_pt_anomaly.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_pt_anomaly.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.7 2005/03/26 20:36:55 raasch 12 14 ! Arguments for non-cyclic boundary conditions added to argument list of 13 15 ! routine exchange_horiz 14 !15 ! Revision 1.6 2001/03/30 07:32:28 raasch16 ! Translation of remaining German identifiers (variables, subroutines, etc.)17 !18 ! Revision 1.5 2001/01/22 19:40:09 letzel19 ! All comments translated into English.20 !21 ! Revision 1.4 2001/01/22 07:19:14 raasch22 ! Module test_variables removed23 !24 ! Revision 1.3 1998/07/06 12:17:25 raasch25 ! + USE test_variables26 !27 ! Revision 1.2 1997/09/12 06:26:16 raasch28 ! Stoerungsradius jetzt proportional zur Gitterweite29 16 ! 30 17 ! Revision 1.1 1997/08/29 08:58:56 raasch -
palm/trunk/SOURCE/init_rankine.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_rankine.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2005/03/26 20:38:49 raasch 12 14 ! Arguments for non-cyclic boundary conditions added to argument list of 13 15 ! routine exchange_horiz 14 !15 ! Revision 1.10 2004/01/28 15:11:27 raasch16 ! Setting of bt before calling pres removed17 !18 ! Revision 1.9 2001/03/30 07:32:49 raasch19 ! Arguments removed from pres,20 ! Translation of remaining German identifiers (variables, subroutines, etc.)21 !22 ! Revision 1.8 2001/01/22 19:44:52 letzel23 ! All comments translated into English.24 !25 ! Revision 1.7 2001/01/22 07:20:30 raasch26 ! Module test_variables removed27 !28 ! Revision 1.6 1998/07/06 12:17:44 raasch29 ! + USE test_variables30 !31 ! Revision 1.5 1998/04/15 11:20:38 raasch32 ! pt, usws und vsws werden an pres uebergeben33 !34 ! Revision 1.4 1997/09/12 06:26:30 raasch35 ! Wirbelradius jetzt proportional zur Gitterweite36 !37 ! Revision 1.3 1997/09/09 08:29:24 raasch38 ! Kehrwerte der Gitterweiten implementiert39 !40 ! Revision 1.2 1997/08/26 06:31:33 raasch41 ! +tend in pres-Aufruf42 16 ! 43 17 ! Revision 1.1 1997/08/11 06:18:43 raasch -
palm/trunk/SOURCE/init_slope.f90
r1 r3 1 1 SUBROUTINE init_slope 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: init_slope.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2006/02/23 12:35:34 raasch 12 14 ! nanz_2dh renamed ngp_2dh 13 !14 ! Revision 1.4 2003/03/16 09:40:02 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.3 2001/03/30 07:33:59 raasch18 ! Translation of remaining German identifiers (variables, subroutines, etc.)19 !20 ! Revision 1.2 2001/01/22 07:21:35 raasch21 ! Module test_variables removed22 15 ! 23 16 ! Revision 1.1 2000/04/27 07:06:24 raasch … … 31 24 ! Remember: when a sloping surface is used, only one constant temperature 32 25 ! gradient is allowed! 33 !------------------------------------------------------------------------------ -!26 !------------------------------------------------------------------------------! 34 27 35 28 USE arrays_3d -
palm/trunk/SOURCE/interaction_droplets_ptq.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: interaction_droplets_ptq.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.1 2005/06/26 19:57:47 raasch 12 14 ! Initial revision -
palm/trunk/SOURCE/local_getenv.f90
r1 r3 1 1 SUBROUTINE local_getenv( var, ivar, value, ivalue ) 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: local_getenv.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.5 2003/05/09 14:37:07 raasch 12 14 ! On the MUK cluster, only PE0 is able to read environment variables. 13 15 ! Therefore, they have to be communicated via broadcast to the other PEs. 14 !15 ! Revision 1.4 2003/03/16 09:40:22 raasch16 ! Two underscores (_) are placed in front of all define-strings17 !18 ! Revision 1.3 2001/01/22 07:23:55 raasch19 ! Comments translated into English20 !21 ! Revision 1.2 1998/07/16 06:50:50 raasch22 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert23 16 ! 24 17 ! Revision 1.1 1997/08/11 06:21:01 raasch … … 29 22 ! ------------ 30 23 ! Getting the values of environment-variabls (for different operating-systems) 31 !------------------------------------------------------------------------------ -!24 !------------------------------------------------------------------------------! 32 25 33 26 #if defined( __lcmuk ) -
palm/trunk/SOURCE/local_stop.f90
r1 r3 1 1 SUBROUTINE local_stop 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: local_stop.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2003/03/16 09:40:28 raasch 12 14 ! Two underscores (_) are placed in front of all define-strings … … 19 21 ! ------------ 20 22 ! Stop program execution 21 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 22 24 23 25 USE pegrid -
palm/trunk/SOURCE/local_system.f90
r1 r3 1 1 SUBROUTINE local_system( command ) 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: local_system.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.4 2003/03/16 09:40:33 raasch 12 14 ! Two underscores (_) are placed in front of all define-strings 13 !14 ! Revision 1.3 2001/01/22 07:25:35 raasch15 ! Comments translated into English16 !17 ! Revision 1.2 1998/07/16 06:51:19 raasch18 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert19 15 ! 20 16 ! Revision 1.1 1997/09/03 06:27:27 raasch … … 25 21 ! ------------ 26 22 ! System calls for different operating systems 27 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 28 24 29 25 CHARACTER (LEN=*) :: command -
palm/trunk/SOURCE/local_tremain.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: local_tremain.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.14 2006/06/02 15:20:33 raasch 12 14 ! Extended to TIT Sun Fire X4600 System (lctit) 13 !14 ! Revision 1.13 2003/05/09 14:39:55 raasch15 ! Time measurement on IBM is now done using function irtc, which allows correct16 ! measurements for jobs running over the 24:00 timeline17 !18 ! Revision 1.12 2003/03/16 09:41:49 raasch19 ! Two underscores (_) are placed in front of all define-strings20 !21 ! Revision 1.11 2003/03/14 13:44:20 raasch22 ! Print statement for testing removed23 !24 ! Revision 1.10 2003/03/12 16:33:55 raasch25 ! Extended to NEC-SX6 systems26 !27 ! Revision 1.9 2002/12/19 15:49:17 raasch28 ! Extended to IBM-Regatta systems29 !30 ! Revision 1.8 2001/03/30 07:35:00 raasch31 ! Translation of remaining German identifiers (variables, subroutines, etc.)32 !33 ! Revision 1.7 2001/01/22 08:34:53 raasch34 ! Module test_variables removed35 !36 ! Revision 1.5 1998/07/16 06:51:44 raasch37 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert38 !39 ! Revision 1.4 1998/07/06 12:19:02 raasch40 ! + USE test_variables41 !42 ! Revision 1.3 1998/03/24 15:27:23 raasch43 ! Zeitbestimmung auf t3eh aufgenommen44 !45 ! Revision 1.2 1998/03/23 08:39:53 raasch46 ! Bestimmung der verbleibenden Zeit auf vpp ergaenzt47 15 ! 48 16 ! Revision 1.1 1998/03/18 20:14:47 raasch -
palm/trunk/SOURCE/local_tremain_ini.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: local_tremain_ini.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.13 2007/02/11 13:07:03 raasch 12 14 ! Allowed cpu limit is now read from file instead of reading the value from 13 15 ! environment variable (see routine parin) 14 !15 ! Revision 1.12 2006/06/02 15:21:09 raasch16 ! Extended to TIT Sun Fire X4600 System (lctit)17 !18 ! Revision 1.11 2003/05/09 14:40:20 raasch19 ! Time measurement on IBM is now done using function irtc, which allows correct20 ! measurements for jobs running over the 24:00 timeline21 !22 ! Revision 1.10 2003/03/16 09:42:01 raasch23 ! Two underscores (_) are placed in front of all define-strings24 !25 ! Revision 1.9 2002/12/19 15:50:31 raasch26 ! Measurements extended for IBM27 !28 ! Revision 1.8 2001/03/30 07:35:16 raasch29 ! Translation of remaining German identifiers (variables, subroutines, etc.)30 !31 ! Revision 1.7 2001/01/22 08:35:38 raasch32 ! Module test_variables removed33 !34 ! Revision 1.5 1998/07/16 06:52:11 raasch35 ! cpp-Direktiven fuer t3ej2 und t3ej5 erweitert36 !37 ! Revision 1.4 1998/07/06 12:19:21 raasch38 ! + USE test_variables39 !40 ! Revision 1.3 1998/03/24 15:27:38 raasch41 ! Initialisierung auf t3eh ergaenzt42 !43 ! Revision 1.2 1998/03/23 08:40:23 raasch44 ! Initialisierung auf vpp ergaenzt45 16 ! 46 17 ! Revision 1.1 1998/03/18 20:15:05 raasch -
palm/trunk/SOURCE/modules.f90
r1 r3 9 9 ! Former revisions: 10 10 ! ----------------- 11 ! $Log: modules.f90,v $ 11 ! $Id$ 12 ! RCS Log replace by Id keyword, revision history cleaned up 13 ! 12 14 ! Revision 1.95 2007/02/11 13:18:30 raasch 13 15 ! version 3.1b (last under RCS control) 14 !15 ! Revision 1.94 2006/08/22 14:12:18 raasch16 ! +dz_max,17 ! id_var_dopts extended for second dimension18 !19 ! Revision 1.93 2006/08/04 14:50:58 raasch20 ! +c_0, diss, dt_min_part, end_time_prel, particles_per_point, sgs_wfu/v/w_part,21 ! use_sgs_for_particles, use_upstream_for_tke needed for the particle package,22 ! +dt_sum, e_m, speed_x/y/z/_sgs in particle_type,23 ! +dt_dopts, time_dopts, netcdf ids for output of particle time series,24 ! do2d_unit and do3d_unit now defined as 2d-arrays, izuf renamed iran25 !26 ! Revision 1.92 2006/06/02 15:21:45 raasch27 ! +NetCDF ids for the shifted u- and v-grid, +data_output_xy/xz/yz28 !29 ! Revision 1.91 2006/04/26 12:22:52 raasch30 ! +threads_per_task,31 ! sums_l, sums_l_l arrays changed to three dimensions32 !33 ! Revision 1.90 2006/04/11 14:58:27 raasch34 ! pl_spectra renamed data_output_sp, length of run_identifier increased to 8035 !36 ! Revision 1.89 2006/03/14 12:55:52 raasch37 ! Default for number_of_particle_groups changed to 138 !39 ! Revision 1.88 2006/03/03 14:26:53 letzel40 ! control parameter topography now has length 4041 !42 ! Revision 1.87 2006/03/03 08:50:25 raasch43 ! Default value for surface_heatflux changed to 9999999.944 !45 ! Revision 1.86 2006/02/23 12:39:14 raasch46 ! Version 3.0 with topography47 ! index arrays nzb_2d replaced by several new variable dependend arrays48 ! (e.g. nzb_u_inner), nzb_diff_2d removed, +new grid point49 ! counters ngp_2dh_outer and ngp_3d_inner,50 ! +skip time control parameters skip_time_..., +log_message,51 ! particle parameters psl, psr, pss, psn, psb, pst, pdx, pdy, pdz52 ! are now 1D arrays (1:max_number_of_particle_groups), default of53 ! number_of_particle_groups changed to 999999, idum in particle_type renamed54 ! tail_id,55 ! +data_output, data_output_user, do2d, do3d, doav for new steering of data56 ! output,57 ! +type_x_int for ghost exchange of integer arrays,58 ! +netcdf_64bit for supporting NetCDF 64-bit offset format,59 ! +mixing_length_1d, dissipation_1d for steering the 1d-model,60 ! +e_min (control parameter for minimum tke), +parameters for assuring a61 ! constant volume flow, +prt_count, prt_start_index, particle_advection,62 ! ebene renamed section, nt_anz renamed current_timestep_number, anz_y..63 ! renamed ngp_y.., nanz_sums, nanz_2dh, nanz_3d renamed ngp_sums, ngp_2dh,64 ! ngp_3d, respectively, pl1d_anz, pl2d_.._anz, plts_anz, pl3d_anz renamed65 ! dopr_n, pl2d_.._no, dots_n, do3d_n, respectively,66 ! ql_v and ql_vp now have the target attribute,67 ! average_period_pr1d removed68 !69 ! Revision 1.85 2005/12/06 16:40:56 raasch70 ! Version 2.11c71 !72 ! Revision 1.84 2005/10/20 15:38:23 raasch73 ! +id_var_rnop_prt, some values of prt_var_names and prt_var_units changed74 !75 ! Revision 1.83 2005/06/29 10:40:30 steinfeld76 ! Scalars ug and vg changed into one-dimensional arrays;77 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,78 ! ug_vertical_gradient_level_ind,79 ! vg_surface, vg_vertical_gradient, vg_vertical_gradient_level,80 ! vg_vertical_gradient_level_ind81 !82 ! Revision 1.82 2005/06/26 20:03:33 raasch83 ! +ql_1, ql_2, ql_m, ql is now a pointer, +variables for cloud droplets (arrays84 ! ql_v,ql_vp, pointer ql_c), gas_constant renamed r_d, latent_heat renamed l_v,85 ! radius instead of diameter, radius + weighting factor included in86 ! particle_type, pl2d and pl3d string length increased from 6 to 10,87 ! default of call_psolver_at_all_substeps is TRUE88 !89 ! Revision 1.81 2005/05/18 15:36:46 raasch90 ! New module netcdf_control, +data_output_format, pl1d_callcount renamed91 ! dopr_time_count,92 ! default values of pl_spectra and spectra_direction changed93 !94 ! Revision 1.80 2005/04/23 09:31:14 raasch95 ! Implicit counters i renamed i9 due to declaration conflicts with i in other96 ! subroutines (reported by Portland compiler),97 ! crmax increased to 100, fcl_factor renamed cfl_factor98 !99 ! Revision 1.79 2005/03/26 20:40:48 raasch100 ! +variables for non-cyclic boundary conditions, call_psolver_at_all_substeps,101 ! default value for residual limit increased from 1E-6 to 1E-4102 !103 ! Revision 1.78 2004/04/30 12:30:01 raasch104 ! idum added in particle_type for better alignment (decalpha),105 ! DATA statements replaced by (/ .../) because of internal compiler errors106 ! on decalpha, + indices for direct transpositions z --> y,107 ! dt_max_1d increased to 300.0 (from 60.0), +grid_matching,108 ! arrays cross_ts_numbers and cross_ts_number_count enlarged,109 ! impulse_advec renamed momentum_advec110 !111 ! Revision 1.77 2004/01/30 10:30:26 raasch112 ! +index arrays nzb_2d, nzb_diff_2d113 !114 ! Revision 1.76 2004/01/28 15:25:15 raasch115 ! +array tsc for timestep steering, +tu_m, tv_m etc. pointer for additional116 ! tendencies needed for Runge-Kutta scheme, +intermediate_timestep_count(_max)117 ! type log changed to logpoint due to name conflict with intrinsic log118 !119 ! Revision 1.75 2003/10/29 08:58:54 raasch120 ! +mg_cycles, mg_switch_to_pe0, mg_switch_to_pe0_level in module121 ! control_parameters,122 ! +mg_loc_ind in module indices, variables related to particle groups added123 ! to module particle_attributes124 !125 ! Revision 1.74 2003/08/01 08:40:59 raasch126 ! +isclice_dvrp, lock_steering_update, slicer_position_dvrp,127 ! slicer_range_limits_dvrp in module dvrp_variables128 !129 ! Revision 1.73 2003/05/09 14:41:11 raasch130 ! New version number 2.7131 !132 ! Revision 1.72 2003/04/25 18:54:39 raasch133 ! Module mpi added in module pegrid in case of linux cluster134 !135 ! Revision 1.71 2003/04/16 13:27:33 raasch136 ! Arrays cross_ts_uy... enlarged137 !138 ! Revision 1.70 2003/04/16 13:07:43 raasch139 ! Array cross_ts_profiles enlarged140 !141 ! Revision 1.69 2003/03/16 09:42:07 raasch142 ! Two underscores (_) are placed in front of all define-strings143 !144 ! Revision 1.68 2003/03/14 13:44:52 raasch145 ! +random_generator in control_parameters146 !147 ! Revision 1.67 2003/03/12 16:34:19 raasch148 ! +tasks_per_node in pegrid149 !150 ! Revision 1.66 2003/03/04 11:32:50 raasch151 ! +type steering, +steering_dvrp in dvrp_variables152 !153 ! Revision 1.65 2002/12/19 15:52:37 raasch154 ! +initial_wallclock_time in cpulog, +myid_char_14 in pegrid,155 ! default value for termination_time_needed now set in check_parameters,156 ! array dl_grid removed, +wall_adjustment_factor, dt_restart, restart_time,157 ! time_restart158 !159 ! Revision 1.64 2002/09/12 13:05:44 raasch160 ! Default values of *_vertical_gradient_level changed from 100000.0 to -1.0,161 ! +molecular_viscosity, density_ratio162 !163 ! Revision 1.63 2002/06/11 13:12:29 raasch164 ! Entrys buoyancy, diffusion_e, diffusion_s, diffusion_u, diffusion_v and165 ! diffusion_w removed from module pointer_interfaces,166 ! length of character variables host + psolver increased to 64 and 16167 !168 ! Revision 1.62 2002/05/02 18:53:18 raasch169 ! Module singleton moved to a separate file. Arrays e_3, pt_3, q_3, u_3, v_3,170 ! w_3 re-introduced. Pointer e_p, pt_p, q_p, u_p, v_p, w_p re-introduced.171 !172 ! Revision 1.61 2002/04/16 08:09:32 raasch173 ! Information about speed components and start positions included in174 ! particle_type, +dt_write_particle_data, particle_advection_start,175 ! time_write_particle_data and parameters for particle boundary conditions,176 ! +vertical_particle_advection,177 ! + variables for scalar boundary conditions and initialization178 !179 ! Revision 1.60 2001/11/12 16:05:01 raasch180 ! +read_particles_from_restartfile181 !182 ! Revision 1.58 2001/09/04 12:00:52 raasch183 ! New version number 2.3b184 !185 ! Revision 1.57 2001/08/21 09:52:08 raasch186 ! +comm_palm, maximum_number_of_tailpoints, maximum_tailpoint_age,187 ! minimum_tailpoint_distance, particle_tail_coordinates, superelevation_x,188 ! superelevation_y, wall_adjustment189 !190 ! Revision 1.56 2001/07/20 13:09:07 raasch191 ! +grid_level, maximum_grid_level, cycle_mg, gamma_mg, ngsrb,192 ! residual_limit in control_parameters, +dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg193 ! in arrays_3d, +superelevation in dvrp_variables, +nxl_mg .. nzt_mg in194 ! indices, +ddx2_mg, ddy2_mg in grid_variables195 !196 ! Revision 1.55 2001/07/12 12:13:24 raasch197 ! default of bc_e_b changed from (u*)**2+neumann to neumann,198 ! +particle_maximum_age, psl-pst, pdx-pdz,199 ! default values defined for dvrp_directory, dvrp_file, dvrp_host, dvrp_password200 ! SAVE-statements added to all modules201 !202 ! Revision 1.54 2001/03/30 07:35:57 raasch203 ! Interface for advec_s_bc changed, interface for production_e eliminated,204 ! arrays work, work1, e_3, pt_3, q_3, u_3, v_3, w_3 removed,205 ! pointer e_p, pt_p, q_p, u_p, v_p, w_p removed,206 ! -nw, nw1,207 ! Translation of remaining German identifiers (variables, subroutines, etc.)208 !209 ! Revision 1.53 2001/01/30 01:57:16 raasch210 ! +passive_scalar211 !212 ! Revision 1.51 2001/01/26 00:31:49 raasch213 ! +fft_method, use_surface_fluxes in steuer, +module singleton,214 ! +write_particle_informations in particle_attributes215 !216 ! Revision 1.49 2001/01/05 15:13:16 raasch217 ! new module spectrum for computing spectra, +average_count_sp,218 ! averaging_interval_sp, dt_dosp, time_average_sp, time_dosp in module steuer,219 ! +spectrum_x, spectrum_y in module statistik220 !221 ! Revision 1.48 2001/01/02 17:31:36 raasch222 ! openfile has the new type file_status, version updated to 2.1a223 !224 ! Revision 1.47 2000/12/30 18:48:53 raasch225 ! Default value of mode_dvrp changed226 !227 ! Revision 1.46 2000/12/28 13:21:35 raasch228 ! version updated to 2.1229 !230 ! Revision 1.45 2000/12/28 13:19:35 raasch231 ! +mode_dvrp, +interval_..._dvrp, REAL*4 changed to REAL in particle_attributes,232 ! +dt_prel, time_prel, -particle_transport233 !234 ! Revision 1.44 2000/07/03 13:02:32 raasch235 ! version updated to 2.0b236 !237 ! Revision 1.43 2000/07/03 12:59:53 raasch238 ! +dvrp_directory, dvrp_host, dvrp_password, dvrp_username,239 ! new module pointer_interfaces which includes interfaces for subroutines with240 ! pointer arguments called in prognostic_equations241 !242 ! Revision 1.42 2000/04/27 07:35:56 raasch243 ! npex and npey must be defined in general (not only for parallel case)244 !245 ! Revision 1.41 2000/04/27 07:07:14 raasch246 ! +pt_slope_offset, pt_slope_ref, npex, npey, dt_dvrp, time_dvrp,247 ! new module dvrp_variables248 ! -dt_plisos, dt_plpart, isosurface_filecount, time_plisos, time_plpart249 ! Version 2.0a, all comments translated into English, old revision remarks250 ! deleted251 !252 ! Revision 1.40 2000/04/13 13:17:25 schroeter253 ! Felder fuer Wolkenphysik zu Modul arrays3d hinzugefuegt, neues Modul254 ! (cloud_parameters) fuer Wolkenphysikparameter eingefuegt255 !256 ! Revision 1.39 2000/01/10 10:05:11 10:05:11 raasch (Siegfried Raasch)257 ! +use_ug_for_galilei_tr, Aufteilung der zeitabhaengigen Groessen in pointer258 ! und targets, + timestep_count259 !260 ! Revision 1.38 1999/12/28 10:06:52 raasch261 ! +Variablen fuer Plot von Isooberflaechen262 !263 ! Revision 1.37 1999/11/25 16:26:38 raasch264 ! Neues Modul particle_attributes,265 ! +Variablen fuer Partikelverlagerung: particle_transport, dt_plpart,266 ! Version 1.1e267 16 ! 268 17 ! Revision 1.1 1997/07/24 11:21:26 raasch -
palm/trunk/SOURCE/netcdf.f90
r1 r3 11 11 ! Former revisions: 12 12 ! ----------------- 13 ! $Log: netcdf.f90,v $ 13 ! $Id$ 14 ! RCS Log replace by Id keyword, revision history cleaned up 15 ! 14 16 ! Revision 1.12 2006/09/26 19:35:16 raasch 15 17 ! Bugfix yv coordinates for yz cross sections 16 !17 ! Revision 1.11 2006/08/22 14:13:11 raasch18 ! Extension for timeseries information of particle groups,19 ! xz and yz cross sections now up to nzt+120 !21 ! Revision 1.10 2006/08/04 14:51:42 raasch22 ! Extensions for output of particle time series23 ! do2d_unit and do3d_unit now defined as 2d-arrays24 !25 ! Revision 1.9 2006/06/13 13:01:08 raasch26 ! Some errors of revision 1.8 removed27 !28 ! Revision 1.8 2006/06/02 15:23:13 raasch29 ! x- and y-coordinates of 2d- and 3d-variables are now precisely defined30 ! depending on which grid (u-, v- or scalar) the respective variable is defined31 !32 ! Revision 1.7 2006/04/21 14:23:47 letzel33 ! Error removed in extend mode when checking whether the selected cross sections34 ! match those in the already existing NetCDF file.35 !36 ! Revision 1.6 2006/04/11 14:58:47 raasch37 ! pl_spectra renamed data_output_sp38 !39 ! Revision 1.5 2006/02/23 12:42:07 raasch40 ! Extensions for output of 2d/3d time averaged data, definition of grids for41 ! user-defined output quantities,42 ! error removed getting the variable ids for _ext for profiles and time series,43 ! ebene renamed section, pl1d_anz renamed dopr_n,44 ! error number argument for handle_netcdf_error,45 ! no optimization on ibmy due to compiler error46 !47 ! Revision 1.4 2005/12/06 16:42:57 raasch48 ! Error in netcdf variable declaration for statistic_regions /= 0 removed49 !50 ! Revision 1.3 2005/10/20 15:40:17 raasch51 ! real number of particles is additionally contained as a variable on the52 ! particle data file,53 ! pr_xz added to vertical cross sections,54 ! error removed in calculating y axis data for yz cross sections55 !56 ! Revision 1.2 2005/06/26 20:04:31 raasch57 ! Particle output extended, 2d- and 3d-files are not opened if there is no58 ! output for these files (pl2d=pl3d=' ')59 18 ! 60 19 ! Revision 1.1 2005/05/18 15:37:16 raasch -
palm/trunk/SOURCE/package_parin.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: package_parin.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.18 2006/08/04 14:52:23 raasch 12 14 ! +dt_dopts, dt_min_part, end_time_prel, particles_per_point, 13 15 ! use_sgs_for_particles in particles_par 14 !15 ! Revision 1.17 2006/04/11 14:59:05 raasch16 ! pl_spectra renamed data_output_sp17 !18 ! Revision 1.16 2006/02/23 12:43:24 raasch19 ! +number_of_particle_groups, skip_particles_for_tail, use_particle_tails in20 ! particles_par,21 ! +skip_time_dosp in spectra_par,22 ! dt_plsp renamed dt_dosp, default setting of dt_dops,23 ! set flag particle_advection24 !25 ! Revision 1.15 2005/06/26 20:05:42 raasch26 ! diameter renamed radius27 !28 ! Revision 1.14 2003/10/29 09:01:29 raasch29 ! +diameter in particles_par, uniform_psize renamed dvrp_psize30 !31 ! Revision 1.13 2003/08/01 08:43:33 raasch32 ! +slicer_range_limits_dvrp in dvrp_graphics_par33 !34 ! Revision 1.12 2003/03/16 09:42:11 raasch35 ! Two underscores (_) are placed in front of all define-strings36 !37 ! Revision 1.11 2002/09/12 13:06:50 raasch38 ! +density_ratio in particles_par39 !40 ! Revision 1.10 2002/04/16 08:10:08 raasch41 ! +variables for particle boundary conditions, +vertical_particle_advection42 ! in particles_par43 !44 ! Revision 1.9 2001/11/12 16:07:04 raasch45 ! +read_particles_from_restartfile46 !47 ! Revision 1.7 2001/08/21 09:52:55 raasch48 ! +superelevation, superelevation_x, superelevation_y in dvrp_graphics_par,49 ! +maximum_number_of_particles, maximum_number_of_tailpoints,50 ! maximum_tailpoint_age, minimum_tailpoint_distance in particles_par51 !52 ! Revision 1.6 2001/07/12 12:42:49 raasch53 ! +several new parameters in dvrp_graphics_par and particles_par54 !55 ! Revision 1.5 2001/03/30 07:37:56 raasch56 ! Namelist names changed,57 ! Translation of remaining German identifiers (variables, subroutines, etc.)58 !59 ! Revision 1.4 2001/01/25 07:10:48 raasch60 ! +write_particle_informations61 !62 ! Revision 1.2 2001/01/05 15:13:34 raasch63 ! Spectra package added64 16 ! 65 17 ! Revision 1.1 2000/12/28 13:21:57 raasch -
palm/trunk/SOURCE/palm.f90
r1 r3 12 12 ! Former revisions: 13 13 ! ----------------- 14 ! $Log: palm.f90,v $ 14 ! $Id$ 15 ! RCS Log replace by Id keyword, revision history cleaned up 16 ! 15 17 ! Revision 1.10 2006/08/04 14:53:12 raasch 16 18 ! Distibution of run description header removed, call of header moved behind 17 19 ! init_3d_model 18 !19 ! Revision 1.9 2006/02/23 12:44:43 raasch20 ! plot_.. renamed data_output_..21 !22 ! Revision 1.8 2004/04/30 12:31:07 raasch23 ! Missing cpp-instructions added, leap_frog changed to time_integration24 !25 ! Revision 1.7 2003/05/09 14:46:49 raasch26 ! String converted to integer before send via broadcast due to Linux MPICH27 ! problem28 !29 ! Revision 1.6 2003/03/16 09:42:17 raasch30 ! Two underscores (_) are placed in front of all define-strings31 !32 ! Revision 1.5 2002/12/19 15:53:59 raasch33 ! Routine local_tremain_ini moved from init_3d_model to here34 !35 ! Revision 1.4 2001/08/21 09:53:35 raasch36 ! comm_palm is now used as the global communicator since one PE maybe split from37 ! MPI_COMM_WORLD for usage in dvrp-graphics38 !39 ! Revision 1.3 2001/03/30 07:39:41 raasch40 ! Translation of remaining German identifiers (variables, subroutines, etc.),41 ! closing of unit 14 moved from routine write_3d_binary to this routine42 20 ! 43 21 ! Revision 1.2 2001/01/25 07:15:06 raasch … … 45 23 ! Initialization of dvrp logging as well as exit of dvrp moved to new 46 24 ! subroutines init_dvrp_logging and close_dvrp (file init_dvrp.f90) 47 !48 ! Revision 1.23 2001/01/05 15:13:51 raasch49 ! +module spectrum50 !51 ! Revision 1.22 2001/01/02 17:32:55 raasch52 ! Subroutine close_files renamed to close_file, now called with an argument53 !54 ! Revision 1.21 2000/12/28 13:32:16 raasch55 ! Call of MPI_FINALIZE moved to the end of the program. At the beginning56 ! comm2d is set to MPI_COMM_WORLD (it was undefined before).57 ! Call of write_particles at the end. Open and close of DVRP-Software at the58 ! beginning and at the end (optionally, by cpp-directives)59 !60 ! Revision 1.20 2000/04/27 06:42:43 raasch61 ! closing call to dvrp-software (if used) at the end of the run,62 ! old revision remarks deleted63 !64 ! Revision 1.19 2000/01/25 15:29:23 letzel65 ! All comments translated into English66 !67 ! Revision 1.18 1999/11/25 16:27:49 raasch68 ! Bei Benutzung von VTK-Grafik wird Modell als UP in C++ - Programm eingebunden69 25 ! 70 26 ! Revision 1.1 1997/07/24 11:23:35 raasch -
palm/trunk/SOURCE/parin.f90
r2 r3 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +netcdf_64bit_3d in d3par, revision history cleaned up6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! +netcdf_64bit_3d in d3par, 12 ! RCS Log replace by Id keyword, revision history cleaned up 11 13 ! 12 14 ! Revision 1.57 2007/02/11 13:11:22 raasch -
palm/trunk/SOURCE/poisfft.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: poisfft.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.24 2006/08/04 15:00:24 raasch 12 14 ! Default setting of the thread number tn in case of not using OpenMP … … 15 17 ! Additional compiler directive in routine tridia_1dd for preventing loop 16 18 ! exchange on NEC-SX6 17 !18 ! Revision 1.22 2005/05/18 15:53:00 raasch19 ! Missing argument in ffty20 !21 ! Revision 1.21 2005/03/26 20:51:59 raasch22 ! Use of module indices removed from routine split_1dd since it caused errors23 ! in case of nx /= ny24 19 ! 25 20 ! Revision 1.20 2004/04/30 12:38:09 raasch … … 30 25 ! is now a 1d-array, variables passed by modules instead of using parameter 31 26 ! lists, enlarged transposition arrays introduced 32 !33 ! Revision 1.19 2003/04/16 13:11:16 raasch34 ! Internal fft method calls substituted by calls of fft_x_1dd and fft_y_1dd35 !36 ! Revision 1.18 2003/03/16 09:44:45 raasch37 ! Two underscores (_) are placed in front of all define-strings38 !39 ! Revision 1.17 2003/03/12 16:36:34 raasch40 ! System specific version for NEC included (equivalent to CRAY scilib41 ! calls, but different subroutine names; the CRAY names are replaced with42 ! NEC names by cpp), size of array work1 is expanded because NEC routines43 ! need more memory44 !45 ! Revision 1.15 2002/06/11 13:14:31 raasch46 ! System specific version for IBM included (parallel case only)47 !48 ! Revision 1.13 2001/01/25 22:19:18 letzel49 ! All comments translated into English.50 !51 ! Revision 1.12 2001/01/23 01:53:04 raasch52 ! Total revision of parallel and non-paralle fft-routines. Additional53 ! implementation of the singleton algorithm.54 !55 ! Revision 1.9 1998/07/27 08:12:56 schroeter56 ! Berechnung von l in maketri geschieht nun ausserhalb57 ! der k-Schleife, da unabhaengig von k, dadurch58 ! Rechenzeitersparnis59 !60 ! Revision 1.8 1998/07/24 07:52:09 schroeter61 ! Anpassung des Druckloesealgorithmusses an Vorgehensweise62 ! von Schmidt, Schumann und Volkert63 ! (DFVLR-Mitteilung 84-15, S. 20-26)64 27 ! 65 28 ! Revision 1.1 1997/07/24 11:24:14 raasch -
palm/trunk/SOURCE/poisfft_hybrid.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: poisfft_hybrid.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2004/04/30 12:43:14 raasch 12 14 ! Renaming of fft routines, additional argument in calls of fft_y_m 13 !14 ! Revision 1.10 2003/08/01 08:49:30 raasch15 ! Error removed: reqa=req after calling routine MPI_IRECV (cascade)16 !17 ! Revision 1.9 2003/05/09 14:50:47 raasch18 ! Parameters in calls of MPI_ALLTOALL, MPI_WAITALL changed in order to fulfill19 ! f95 type and rank requirements20 !21 ! Revision 1.8 2003/03/16 11:17:04 raasch22 ! +module pegrid in poisfft_hybrid_ini23 !24 ! Revision 1.7 2003/03/16 09:45:52 raasch25 ! Two underscores (_) are placed in front of all define-strings26 !27 ! Revision 1.6 2003/03/16 09:11:42 raasch28 ! Fortran error removed29 !30 ! Revision 1.5 2003/03/14 13:45:53 raasch31 ! Loop optimization in maketri_hybrid32 !33 ! Revision 1.4 2003/03/12 16:37:22 raasch34 ! Vectorized version (poisfft_hybrid_omp_vec) added, arguments ar and tri35 ! passed to routine substi_hybrid (otherwise no vectorization and very36 ! bad performance on NEC), reading of environment variable tasks_per_node37 ! moved to routine parin38 !39 ! Revision 1.3 2002/12/20 15:18:02 raasch40 ! Variable statement restricted for parallel environment in subroutine41 ! cascade42 15 ! 43 16 ! Revision 1.2 2002/12/19 16:08:31 raasch -
palm/trunk/SOURCE/poismg.f90
r1 r3 12 12 ! Former revisions: 13 13 ! ----------------- 14 ! $Log: poismg.f90,v $ 14 ! $Id$ 15 ! RCS Log replace by Id keyword, revision history cleaned up 16 ! 15 17 ! Revision 1.6 2005/03/26 20:55:54 raasch 16 18 ! Implementation of non-cyclic (Neumann) horizontal boundary conditions, 17 19 ! routine prolong simplified (one call of exchange_horiz spared) 18 !19 ! Revision 1.5 2003/10/29 09:04:34 raasch20 ! On a defined level, data are gathered on PE0 and further calculations are21 ! carried out only on this PE.22 !23 ! Revision 1.4 2003/03/16 09:45:57 raasch24 ! Two underscores (_) are placed in front of all define-strings25 !26 ! Revision 1.3 2001/07/23 12:00:31 raasch27 ! Neumann boundary conditions for the upper boundary and Diriclet conditions28 ! for the lower boundary implemented29 !30 ! Revision 1.2 2001/07/20 15:21:55 raasch31 ! Praeprocessor directives included for allreduce32 20 ! 33 21 ! Revision 1.1 2001/07/20 13:10:51 raasch -
palm/trunk/SOURCE/prandtl_fluxes.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: prandtl_fluxes.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.19 2006/04/26 12:24:35 raasch 12 14 ! +OpenMP directives and optimization (array assignments replaced by DO loops) 13 !14 ! Revision 1.18 2006/02/23 12:49:32 raasch15 ! shf, ts, rif, us, usws, vsws, qs and qsws are now defined at the actual16 ! surface of the model domain, i.e. either at the bottom or at the height17 ! of the topography. Thus, zu(1) is replaced18 ! by zp = ( zu(nzb_[variable]_inner(j,i)+1 - zw(nzb_[variable]_inner(j,i) ),19 ! [variable](1,j,i) by [variable](nzb_[variable]_inner(j,i)+1,j,i) and20 ! [variable](0,j,i) by [variable](nzb_[variable]_inner(j,i) ,j,i).21 !22 ! Revision 1.17 2004/04/30 12:43:50 raasch23 ! rif_m replaced by rif (they are the same when used here)24 !25 ! Revision 1.16 2003/11/20 15:13:26 raasch26 ! Variable name (B) changed from capital to small27 !28 ! Revision 1.15 2001/03/30 07:46:47 raasch29 ! Translation of remaining German identifiers (variables, subroutines, etc.)30 !31 ! Revision 1.14 2001/01/29 12:33:38 raasch32 ! Passive scalar is considered33 !34 ! Revision 1.13 2001/01/25 07:23:22 raasch35 ! Range of ts is limited, since otherwise overflow may occur on REAL*436 ! machines37 !38 ! Revision 1.11 2001/01/22 07:50:34 raasch39 ! Module test_variables removed40 !41 ! Revision 1.10 2000/04/13 13:44:41 schroeter42 ! + implentation of the Prandtl layer for the total water content43 !44 ! Revision 1.9 2000/01/26 14:48:48 14:48:48 letzel (Marcus Letzel)45 ! correct wrong time level of rif in computation of theta*,46 ! all comments translated into English47 !48 ! Revision 1.8 1998/09/22 17:27:15 raasch49 ! Testweise Randbedingung fuer TKE mit cm = 0.4, aber vorerst auskommentiert50 !51 ! Revision 1.7 1998/08/05 06:54:23 raasch52 ! Begrenzung von rif ist jetzt variabel (rif_max, rif_min)53 !54 ! Revision 1.6 1998/07/06 12:30:01 raasch55 ! + USE test_variables56 !57 ! Revision 1.5 1998/04/15 11:22:58 raasch58 ! Berechnung von theta* ueber Oberflaechentemperatur moeglich59 !60 ! Revision 1.4 1998/03/11 11:53:22 raasch61 ! Zusaetzliche untere Randbedingung fuer TKE ( (u*/0.1)**2 )62 !63 ! Revision 1.3 1998/02/10 15:08:53 raasch64 ! Grenzfall bei labiler Schichtung ( a=1, b=1 ) wird stabil gerechnet65 ! Begrenzung von rif aktiviert66 !67 ! Revision 1.2 1998/02/04 16:09:29 raasch68 ! Berechnung der Waermefluesse69 15 ! 70 16 ! Revision 1.1 1998/01/23 10:06:06 raasch -
palm/trunk/SOURCE/pres.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: pres.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.25 2006/04/26 13:26:12 raasch 12 14 ! OpenMP optimization (+localsum, threadsum) 13 !14 ! Revision 1.24 2006/02/23 12:50:08 raasch15 ! nzb replaced by nzb_.._inner, optional volume flow control, calculation of16 ! divergence sum reordered due to optimization17 !18 ! Revision 1.23 2005/05/18 15:53:44 raasch19 ! call conditions for poisfft_hybrid modified20 !21 ! Revision 1.22 2005/03/26 20:57:24 raasch22 ! Adjustments for non-cyclic boundary conditions: arguments added to routine23 ! exchange_horiz, mean horizontal wind parallel to the outflow is calculated.24 !25 ! Revision 1.21 2004/04/30 12:45:00 raasch26 ! Most of the arguments removed from poisfft call, +module poisfft_mod,27 ! poisfft is also used for 1d-decompositions, more memory has to be28 ! allocated for tend, when an enlarged domain is used for transposition29 !30 ! Revision 1.20 2003/03/16 09:46:02 raasch31 ! Two underscores (_) are placed in front of all define-strings32 !33 ! Revision 1.19 2003/03/14 13:46:32 raasch34 ! Optimization of loops (IBM-version and vectorizable version)35 !36 ! Revision 1.18 2003/03/12 16:38:20 raasch37 ! Simulated_time has been replaced by sums(nzb+1,4) in first if clause38 ! in order to avoid run time errors due to compiler bugs on ibm and nec39 ! machines, due to optimization, dividing d by dt_3d is now done in the40 ! first loop (not in a seperated loop),41 ! error removed: tend=p added after calling sor method42 !43 ! Revision 1.17 2002/06/11 13:17:03 raasch44 ! Hybrid solver included, array operations on d removed due to bad performance45 ! on IBM, loops including the array d are combined, OpenMP directives added.46 ! Total perturbation pressure is no more obtained iteratively (p=p+tend) from47 ! the current perturbation pressure.48 !49 ! Revision 1.16 2001/09/04 12:05:52 12:05:52 raasch (Siegfried Raasch)50 ! Value of perturbation pressure is no more obtained iteratively (p = p + tend,51 ! FFT-solver only)52 !53 ! Revision 1.15 2001/07/20 13:12:54 raasch54 ! Multigrid method included55 !56 ! Revision 1.14 2001/03/30 07:47:29 raasch57 ! All arguments removed, dp replaced by tend, array work removed from58 ! sor argument list,59 ! Translation of remaining German identifiers (variables, subroutines, etc.)60 !61 ! Revision 1.13 2001/01/22 07:53:53 raasch62 ! Module test_variables removed63 !64 ! Revision 1.12 2000/01/26 13:23:32 letzel65 ! All comments translated into English66 !67 ! Revision 1.11 1998/09/22 17:28:14 raasch68 ! dp war bei SOR-Aufrufen bisher nicht initialisiert, entsprechend geaendert69 !70 ! Revision 1.10 1998/07/06 12:30:22 raasch71 ! + USE test_variables72 !73 ! Revision 1.9 1998/04/15 11:23:20 raasch74 ! Testweise Zusatzbedingung fuer den Druck am unteren Rand bei inhomogener75 ! Oberflaechentemperatur (pt wird uebergeben)76 !77 ! Revision 1.8 1998/03/30 11:37:59 raasch78 ! Divergenzsummenberechnung (PEs) nach flow_statistics verlagert79 !80 ! Revision 1.7 1998/03/25 13:56:16 raasch81 ! dt in dt_3d umbenannt82 !83 ! Revision 1.6 1997/09/16 06:38:37 raasch84 ! Ueberfluessige Initialisierung von d entfernt85 !86 ! Revision 1.5 1997/09/12 06:29:56 raasch87 ! Randbedingungen umbenannt88 !89 ! Revision 1.4 1997/09/09 08:30:27 raasch90 ! Kehrwerte der Gitterweiten implementiert91 !92 ! Revision 1.3 1997/08/26 06:34:16 raasch93 ! Gesamtstoerdruck p wird in Bew.gleichungen mitgerechnet und ergibt sich94 ! durch Summation des aktuellen Stoerdrucks dp, mit dem hier korrigiert wird95 !96 ! Revision 1.2 1997/08/11 06:25:24 raasch97 ! Felder werden ueber Formalparameterliste uebergeben.98 15 ! 99 16 ! Revision 1.1 1997/07/24 11:24:44 raasch -
palm/trunk/SOURCE/print_1d.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: print_1d.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.11 2006/02/23 12:50:43 raasch 12 14 ! Preliminary no output of time-averaged data 13 !14 ! Revision 1.10 2001/03/30 07:48:24 raasch15 ! Translation of remaining German identifiers (variables, subroutines, etc.)16 !17 ! Revision 1.9 2001/01/22 20:51:40 letzel18 ! All comments translated into English.19 !20 ! Revision 1.8 2001/01/22 07:55:14 raasch21 ! Module test_variables removed22 !23 ! Revision 1.7 1998/07/06 12:31:17 raasch24 ! Mittelungszeitraum wird direkt durch average_period_pr1d angegeben,25 ! + USE test_variables26 !27 ! Revision 1.6 1998/03/30 11:38:16 raasch28 ! Ausgaben fuer die (Teil-) Gebiete29 !30 ! Revision 1.5 1998/03/25 13:56:31 raasch31 ! Erweiterungen fuer Gebietsstatistiken32 !33 ! Revision 1.4 1998/03/18 20:18:33 raasch34 ! Zusaetzliche Header werden ausgegeben35 !36 ! Revision 1.3 1998/03/03 08:01:28 raasch37 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier38 ! Ausgabe von zeitlich gemittelten Profilen39 !40 ! Revision 1.2 1998/02/10 15:09:16 raasch41 ! Ausdruck von Werten im w-Niveau42 15 ! 43 16 ! Revision 1.1 1997/09/19 07:45:22 raasch -
palm/trunk/SOURCE/production_e.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: production_e.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.21 2006/04/26 12:45:35 raasch 12 14 ! OpenMP parallelization of production_e_init 13 !14 ! Revision 1.20 2006/02/23 12:51:25 raasch15 ! nzb_2d and nzb_diff_2d replaced by nzb_s_outer and nzb_diff_s_outer16 ! respectively, momentum flux calculation at vertical walls using profile17 ! functions, loop rearrangement for better vectorization18 !19 ! Revision 1.19 2005/06/29 08:15:33 steinfeld20 ! Error in calculating u_0 and v_0 removed (multiplication instead of21 ! division by 2dz)22 !23 ! Revision 1.18 2004/04/30 12:45:54 raasch24 ! dptdz eliminated25 !26 ! Revision 1.17 2004/01/30 10:36:00 raasch27 ! Velocity gradients at the surface limited (see u_0, v_0),28 ! scalar lower k index nzb_diff replaced by 2d-array nzb_diff_2d29 !30 ! Revision 1.16 2003/03/12 16:40:22 raasch31 ! Full code replaced in the call for all gridpoints instead of calling the32 ! _ij version (required by NEC, because otherwise no vectorization)33 !34 ! Revision 1.15 2002/12/19 16:16:34 raasch35 ! Calculation of deformation tensor re-designed (most of the finite36 ! differences are now formed over two grid spacings, errors in derivatives37 ! along y removed)38 !39 ! Revision 1.14 2002/09/12 13:09:44 raasch40 ! Error in calculating v_0 removed41 !42 ! Revision 1.13 2002/06/11 13:19:04 raasch43 ! Former subroutine changed to a module which allows to be called for all grid44 ! points of a single vertical column with index i,j or for all grid points by45 ! using function overloading.46 ! Calculation of u_0 and v_0 moved to the new subroutine production_e_init47 ! due to the global communication necessary. These arrays are allocated only48 ! once during the first call.49 !50 ! Revision 1.12 2001/08/21 09:59:47 raasch51 ! Special treatment at k=1 generally if surface fluxes are prescribed (not only52 ! in case of a Prandtl layer). In these cases, the boundary conditions are53 ! also included in the shear production term.54 !55 ! Revision 1.11 2001/03/30 07:48:54 raasch56 ! Calculation of shear production and buoyancy production is split into57 ! two loops, arguments k1 and k2 are eliminated and now used as scalars,58 ! Translation of remaining German identifiers (variables, subroutines, etc.)59 !60 ! Revision 1.10 2001/01/22 07:56:19 raasch61 ! Module test_variables removed62 !63 ! Revision 1.9 2001/01/02 17:34:12 raasch64 ! -dpt_dz_d, dpt_dz_u65 !66 ! Revision 1.8 2000/07/03 13:00:24 raasch67 ! Arguments k1 and k2 declared as pointers,68 ! all comments tranlated into English69 !70 ! Revision 1.7 2000/04/18 08:13:00 schroeter71 ! Revision 1.5 rueckgaengig gemacht, Temperaturgrandient im72 ! Auftriebtrem wird wieder durch zentrale Differenzen gebildet73 !74 ! Revision 1.6 2000/04/13 14:24:02 schroeter75 ! condsidering the influence of humidity to TKE-production76 !77 ! Revision 1.5 99/02/17 09:29:18 09:29:18 raasch (Siegfried Raasch)78 ! Vertikalscherung des Windes exakter formuliert79 ! Temperaturgradient im Auftriebsterm enger gefasst, im stabilen Fall wird80 ! jetzt das Minimum vom oberen und unteren Differenzenquotienten verwendet81 !82 ! Revision 1.4 1998/07/06 12:31:52 raasch83 ! + USE test_variables84 !85 ! Revision 1.3 1998/04/15 11:23:49 raasch86 ! Im Fall einer Prandtl-Schicht wird die Energieproduktion bei nzb+187 ! jetzt immer mittels shf und sonst ueber dpt/dz berechnet (bisher nur bei88 ! vorgegebenem Waermestrom)89 !90 ! Revision 1.2 1998/02/19 07:10:51 raasch91 ! vorgegebener Waermestrom wird gegebenenfalls im Auftriebsterm bei92 ! k=nzb+1 direkt angegeben93 15 ! 94 16 ! Revision 1.1 1997/09/19 07:45:35 raasch -
palm/trunk/SOURCE/prognostic_equations.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: prognostic_equations.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.21 2006/08/04 15:01:07 raasch 12 14 ! upstream scheme can be forced to be used for tke (use_upstream_for_tke) 13 15 ! regardless of the timestep scheme used for the other quantities, 14 16 ! new argument diss in call of diffusion_e 15 !16 ! Revision 1.20 2006/02/23 12:52:08 raasch17 ! nzb_2d replaced by nzb_u/v/w/s_inner, +z0 in argument list of diffusion_u/v/w18 !19 ! Revision 1.19 2005/06/29 10:33:41 steinfeld20 ! Scalars ug and vg are replaced by homonymous arrays in order to allow for21 ! the consideration of a potential dependency of the geostrophic wind22 ! components on height23 !24 ! Revision 1.18 2005/03/26 20:58:45 raasch25 ! Extension of horizontal loop upper bounds for non-cyclic boundary conditions,26 ! additional arguments km_damp_x/y for diffusion_u/v/w27 !28 ! Revision 1.17 2004/04/30 12:46:40 raasch29 ! impulse_advec renamed momentum_advec30 !31 ! Revision 1.16 2004/01/30 10:36:57 raasch32 ! Scalar lower k index nzb replaced by 2d-array nzb_2d33 !34 ! Revision 1.15 2004/01/28 15:24:35 raasch35 ! Runge-Kutta schemes available, steering variables at and bt in equations36 ! replaced by array sct37 !38 ! Revision 1.14 2003/03/12 16:40:48 raasch39 ! New routine prognostic_equations_vec, optimized for vector processors40 !41 ! Revision 1.13 2002/06/11 13:20:21 raasch42 ! Single node optimization: loops i and j extracted from all those43 ! tendency-subroutines which do not contain communication. Some communication44 ! parts moved before the i,j loops. Call of particle advection moved to45 ! subroutine leap_frog.46 ! Former subroutine changed to a module which contains two versions: one47 ! with loop optimization for the single equations and one version with one48 ! big optimized loop containing all prognostic equations.49 !50 ! Revision 1.12 2002/05/02 18:54:12 raasch51 ! Timelevel t+dt re-introduced, Asselin filter and exchange of ghost points52 ! moved to routine leap_frog53 !54 ! Revision 1.11 2002/04/16 08:12:40 08:12:40 raasch (Siegfried Raasch)55 ! Particle advection will be performed only after the simulated time has56 ! exceeded a user defined limit57 !58 ! Revision 1.10 2001/09/04 12:10:28 raasch59 ! Exchange of ghost points added for the time filtered arrays60 !61 ! Revision 1.9 2001/08/21 09:59:03 raasch62 ! Calls to user-interface for tendency terms added63 !64 ! Revision 1.8 2001/03/30 07:50:26 raasch65 ! Timelevel t+dt eliminated, Asselin filter included in prognostic equations,66 ! calling arguments eliminated or cut down from advec_s_bc, production_e,67 ! advec_*_ups,68 ! Translation of remaining German identifiers (variables, subroutines, etc.)69 !70 ! Revision 1.7 2001/01/29 12:34:20 raasch71 ! Passive scalar is considered72 !73 ! Revision 1.6 2001/01/22 07:57:20 raasch74 ! Module test_variables removed75 !76 ! Revision 1.5 2000/12/28 13:38:00 raasch77 ! Advec_particles is called unconditionally78 !79 ! Revision 1.4 2000/07/03 13:01:02 raasch80 ! module pointer_interfaces added to reduce cpu-time and memory demands81 ! of the compilation process, actual argument "work" eliminated from82 ! parameter list of diffusion_e83 !84 ! Revision 1.3 2000/04/27 07:09:57 raasch85 ! Temperature offset is added at cyclic boundaries (x-direction) when using86 ! a sloping surface87 !88 ! Revision 1.2 2000/04/13 15:04:28 schroeter89 ! prognostic equation for the total water content added, new routines:90 ! calc_precipitation, impact_of_latent_heat, calc_radiation,91 ! diffusion_pt renamed to diffusion_s, changes in the parameter list of92 ! subroutines diffusion_e, buyoancy, production_e93 17 ! 94 18 ! Revision 1.1 2000/04/13 14:56:27 schroeter -
palm/trunk/SOURCE/random_function.f90
r1 r3 1 1 MODULE random_function_mod 2 !------------------------------------------------------------------------------ -!2 !------------------------------------------------------------------------------! 3 3 ! Actual revisions: 4 4 ! ----------------- … … 7 7 ! Former revisions: 8 8 ! ----------------- 9 ! $Log: random_function.f90,v $ 9 ! $Id$ 10 ! RCS Log replace by Id keyword, revision history cleaned up 11 ! 10 12 ! Revision 1.3 2003/10/29 09:06:57 raasch 11 13 ! Former function changed to a module. 12 !13 ! Revision 1.2 2001/01/22 08:00:23 raasch14 ! Comments translated into English15 14 ! 16 15 ! Revision 1.1 1998/02/04 16:09:45 raasch … … 22 21 ! Random number generator, produces numbers equally distributed in interval [0,1] 23 22 ! This routine is taken from the "numerical recipies" 24 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 25 24 26 25 IMPLICIT NONE -
palm/trunk/SOURCE/random_gauss.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: random_gauss.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.4 2006/08/04 15:01:48 raasch 12 14 ! Range of random number is limited by an upper limit (new second parameter) 13 !14 ! Revision 1.3 2003/10/29 09:08:04 raasch15 ! Module random_function_mod is used16 !17 ! Revision 1.2 2001/01/22 08:03:31 raasch18 ! Comments translated into English19 15 ! 20 16 ! Revision 1.1 1998/03/25 20:09:47 raasch -
palm/trunk/SOURCE/read_3d_binary.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: read_3d_binary.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.4 2006/08/04 15:02:32 raasch 12 14 ! +iran, iran_part 13 !14 ! Revision 1.3 2006/02/23 12:52:45 raasch15 ! +volume_flow_area, volume_flow_initial, -hom_sum_pr,16 ! reading arrays used for time-averaging17 !18 ! Revision 1.2 2005/05/18 15:54:33 raasch19 ! pl1d_callcount renamed dopr_time_count20 15 ! 21 16 ! Revision 1.1 2004/04/30 12:47:27 raasch -
palm/trunk/SOURCE/read_var_list.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: read_var_list.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.34 2006/08/22 14:14:27 raasch 12 14 ! +dz_max 13 15 ! 14 ! Revision 1.33 2006/08/04 15:04:08 raasch15 ! +first_call_advec_particles, time_dopts, use_upstream_for_tke16 !17 ! Revision 1.32 2006/02/23 12:53:31 raasch18 ! +building parameters, conserve_volume_flow, dissipation_1d, e_min,19 ! mixing_length_1d, topography,20 ! -average_count_pl, time_average_pr,21 ! nt_anz renamed current_timestep_number, ..pl.. renamed ..do..22 !23 ! Revision 1.31 2005/06/29 09:12:49 steinfeld24 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,25 ! ug_vertical_gradient_level_ind, vg_surface, vg_vertical_gradient,26 ! vg_vertical_gradient_level, vg_vertical_gradient_level_ind;27 ! allocation of the new arrays ug and vg28 !29 ! Revision 1.30 2005/06/26 20:09:42 raasch30 ! +cloud_droplets31 !32 ! Revision 1.29 2005/05/18 15:55:12 raasch33 ! +netcdf_precision34 !35 ! Revision 1.28 2005/03/26 21:01:00 raasch36 ! +bc_lr, bc_ns, inflow_disturbance_begin, inflow_disturbance_end, km_damp_max,37 ! outflow_damping_width38 !39 ! Revision 1.27 2004/04/30 12:49:56 raasch40 ! +grid_matching, impulse_advec renamed momentum_advec41 !42 ! Revision 1.26 2004/01/28 15:27:33 raasch43 ! at and bt replaced by array tsc44 !45 ! Revision 1.25 2003/10/29 09:10:41 raasch46 ! +last_dt_change47 !48 ! Revision 1.24 2003/03/14 13:47:15 raasch49 ! +random_generator50 !51 ! Revision 1.23 2002/12/19 16:18:27 raasch52 ! Array hom, which was formerly read in init_3d_model, will now be allocated53 ! and read here, binary version incremented to 2.2, +time_restart,54 ! STOP statement replaced by call of subroutine local_stop55 !56 ! Revision 1.22 2002/04/16 08:21:06 raasch57 ! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,58 ! s_vertical_gradient, s_vertical_gradient_level59 !60 ! Revision 1.21 2001/08/21 10:01:23 raasch61 ! +wall_adjustment62 !63 ! Revision 1.20 2001/03/30 07:51:25 raasch64 ! Translation of remaining German identifiers (variables, subroutines, etc.),65 ! runtime and package parameters are not read from unit 13 any more66 !67 ! Revision 1.19 2001/01/29 12:35:03 raasch68 ! +passive_scalar69 !70 ! Revision 1.18 2001/01/25 07:24:28 raasch71 ! +fft_method, use_surface_fluxes72 !73 ! Revision 1.17 2001/01/05 15:16:18 raasch74 ! +average_count_sp, time_dosp75 !76 ! Revision 1.16 2000/12/28 13:38:46 raasch77 ! Module dvrp_variables removed, -particle_transport78 !79 ! Revision 1.15 2000/04/27 06:57:04 raasch80 ! +dt_dvrp, dvrp_filecount, npex, npey, threshold, time_dvrp81 ! -dt_plisos, dt_plpart, time_plisos, time_plpart,82 ! all comments translated into English, old revision remarks deleted83 !84 ! Revision 1.14 2000/04/13 13:06:16 schroeter85 ! Einlesen der Initialisierungsparameter fuer Rechnungen86 ! mit Feuchte/Wolkenphysik87 !88 ! Revision 1.13 2000/01/10 10:06:32 10:06:32 raasch (Siegfried Raasch)89 ! +use_ug_for_galilei_tr90 !91 16 ! Revision 1.1 1998/03/18 20:18:48 raasch 92 17 ! Initial revision -
palm/trunk/SOURCE/run_control.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: run_control.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.20 2006/06/02 15:23:47 raasch 12 14 ! cpp-directives extended for lctit 13 !14 ! Revision 1.19 2006/02/23 12:54:21 raasch15 ! Output buffer is flushed also on NEC, nt_anz renamed current_timestep_number16 !17 ! Revision 1.18 2004/01/30 10:38:26 raasch18 ! Euler-step information removed in case of Runge-Kutta19 !20 ! Revision 1.17 2003/10/29 09:13:18 raasch21 ! Data are forced to be immediately written on disc on IBM and Linux clusters22 !23 ! Revision 1.16 2003/04/16 13:12:39 raasch24 ! Output format of iteration count enlarged25 !26 ! Revision 1.15 2001/07/20 13:13:45 raasch27 ! Additional output of the number of multigrid cycles used28 !29 ! Revision 1.14 2001/03/30 07:52:24 raasch30 ! Translation of remaining German identifiers (variables, subroutines, etc.)31 !32 ! Revision 1.13 2001/01/22 08:04:54 raasch33 ! Module test_variables removed34 !35 ! Revision 1.12 2001/01/03 13:46:12 letzel36 ! File output and all comments translated into English.37 !38 ! Revision 1.11 1998/09/22 17:29:49 raasch39 ! Ausgabeformate t.w. um einige Kommastellen erweitert40 !41 ! Revision 1.10 1998/07/06 12:32:56 raasch42 ! + USE test_variables43 !44 ! Revision 1.9 1998/04/21 16:23:48 raasch45 ! Ausgabe von advected_distance46 !47 ! Revision 1.8 1998/03/30 11:50:16 raasch48 ! Divergenzen stehen jetzt in hom49 !50 ! Revision 1.7 1998/03/25 13:57:06 raasch51 ! Erweiterungen fuer Gebietsstatistiken,52 ! dt in dt_3d umbenannt53 !54 ! Revision 1.6 1998/03/18 20:19:50 raasch55 ! Pruefen, ob Einheit 15 bereits geoeffnet ist56 !57 ! Revision 1.5 1998/03/03 08:01:49 raasch58 ! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier,59 ! Ausgabe von w* und z_i60 !61 ! Revision 1.4 1998/02/04 16:10:46 raasch62 ! e_total und e_pertu stehen jetzt im Feld hom63 !64 ! Revision 1.3 1998/01/23 10:08:01 raasch65 ! Ausgabe von u* und theta*66 !67 ! Revision 1.2 1997/08/26 06:34:36 raasch68 ! Schoenheitskorrekturen69 15 ! 70 16 ! Revision 1.1 1997/08/11 06:25:38 raasch -
palm/trunk/SOURCE/singleton.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: singleton.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2004/04/30 12:52:09 raasch 12 14 ! Shape of arrays is explicitly stored in ishape and handled to the -
palm/trunk/SOURCE/sor.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: sor.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.9 2005/03/26 21:02:23 raasch 12 14 ! Implementation of non-cyclic (Neumann) horizontal boundary conditions, 13 15 ! dx2,dy2 replaced by ddx2,ddy2 14 !15 ! Revision 1.8 2001/03/30 07:52:44 raasch16 ! Arrays f1-f3 changed from dummy arguments to allocatable arrays,17 ! Translation of remaining German identifiers (variables, subroutines, etc.)18 !19 ! Revision 1.7 2001/01/30 21:41:25 letzel20 ! All comments translated into English.21 !22 ! Revision 1.6 2001/01/22 08:06:14 raasch23 ! Module test_variables removed24 !25 ! Revision 1.5 1998/07/06 12:33:13 raasch26 ! + USE test_variables27 !28 ! Revision 1.4 1997/09/12 06:30:12 raasch29 ! Randbedingungen umbenannt30 !31 ! Revision 1.3 1997/09/09 08:30:44 raasch32 ! Kehrwerte der Gitterweiten implementiert33 !34 ! Revision 1.2 1997/08/29 09:00:58 raasch35 ! omega --> omega_sor36 16 ! 37 17 ! Revision 1.1 1997/08/11 06:25:56 raasch … … 42 22 ! ------------ 43 23 ! Solve the Poisson-equation with the SOR-Red/Black-scheme. 44 !------------------------------------------------------------------------------ -!24 !------------------------------------------------------------------------------! 45 25 46 26 USE grid_variables -
palm/trunk/SOURCE/spline_x.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: spline_x.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.8 2004/04/30 12:54:20 raasch 12 14 ! Names of transpose indices changed, enlarged transposition arrays introduced 13 !14 ! Revision 1.7 2003/03/16 09:47:17 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.6 2001/03/30 07:53:19 raasch18 ! Arrays r, vad and wrk_spline changed from 3D to 2D and removed from argument19 ! list. Several loops over k/j (parallel/non-parallel) combined to one loop.20 ! Application of long filter moved to this routine. All comments and identifiers21 ! translated into English.22 !23 ! Revision 1.5 2001/01/22 08:07:41 raasch24 ! Module test_variables removed25 !26 ! Revision 1.4 1999/03/25 07:32:31 raasch27 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,28 ! Ueberschwinger werden in gewissen Grenzen erlaubt29 !30 ! Revision 1.3 1999/02/26 17:53:26 schroeter31 ! - Gradientenkontrolle fuer den nicht-parallelen Teil32 ! - statistische Auswertung ueber den prozentualen Anteil des33 ! Upstream-Verfahrens an der Gesamtadvektion fuer nicht-34 ! parallelen Teil35 !36 ! Revision 1.2 1999/02/17 09:30:19 raasch37 ! Wertebegrenzung zur Verhinderung von Ueberschwingern38 15 ! 39 16 ! Revision 1.1 1999/02/05 09:15:59 raasch -
palm/trunk/SOURCE/spline_y.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: spline_y.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.9 2004/04/30 12:54:37 raasch 12 14 ! Names of transpose indices changed, enlarged transposition arrays introduced 13 !14 ! Revision 1.8 2003/03/16 09:48:41 raasch15 ! Two underscores (_) are placed in front of all define-strings16 !17 ! Revision 1.7 2001/03/30 07:53:40 raasch18 ! Arrays r and wrk_spline changed from 3D to 2D and removed from argument19 ! list. Several loops over k/i (parallel/non-parallel) combined to one loop.20 ! Application of long filter moved to this routine. All comments and identifiers21 ! translated into English.22 !23 ! Revision 1.6 2001/01/22 08:09:04 raasch24 ! Module test_variables removed25 !26 ! Revision 1.5 1999/11/25 16:30:42 raasch27 ! Laufindex-Fehler korrigiert28 !29 ! Revision 1.4 1999/03/25 07:33:16 raasch30 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,31 ! Ueberschwinger werden in gewissen Grenzen erlaubt32 !33 ! Revision 1.3 1999/02/26 17:54:28 schroeter34 ! - Gradientenkontrolle fuer den nicht-parallelen Teil35 ! - statistische Auswertung ueber den prozentualen Anteil des36 ! Upstream-Verfahrens an der Gesamtadvektion fuer nicht-37 ! parallelen Teil38 !39 ! Revision 1.2 1999/02/17 09:31:34 raasch40 ! Test einer Wertebegrenzung zur Verhinderung von Ueberschwingern41 15 ! 42 16 ! Revision 1.1 1999/02/05 09:16:31 raasch -
palm/trunk/SOURCE/spline_z.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: spline_z.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.9 2005/06/29 08:22:56 steinfeld 12 14 ! Dependency of ug and vg on height considered in the determination of the 13 15 ! upper boundary condition for vad 14 !15 ! Revision 1.8 2004/04/30 12:55:23 raasch16 ! Enlarged transposition arrays introduced17 !18 ! Revision 1.7 2003/03/16 09:49:33 raasch19 ! Two underscores (_) are placed in front of all define-strings20 !21 ! Revision 1.6 2001/03/30 07:54:07 raasch22 ! Arrays r vad, and wrk_spline changed from 3D to 2D and removed from argument23 ! list. Several loops over i combined to one loop.24 ! Application of long filter moved to this routine. All comments and identifiers25 ! translated into English.26 !27 ! Revision 1.5 2001/01/22 08:10:21 raasch28 ! Module test_variables removed29 !30 ! Revision 1.4 1999/03/25 07:34:55 raasch31 ! Filterung der Ueberschwinger geschieht optional, ups_limit_e eingefuehrt,32 ! Ueberschwinger werden in gewissen Grenzen erlaubt33 !34 ! Revision 1.3 1999/02/26 17:55:41 schroeter35 ! - Gradientenkontrolle fuer den nicht-parallelen Teil36 ! - statistische Auswertung ueber den prozentualen Anteil des37 ! Upstream-Verfahrens an der Gesamtadvektion fuer nicht-38 ! parallelen Teil39 !40 ! Revision 1.2 1999/02/17 09:32:55 raasch41 ! Wertebegrenzung zur Verhinderung von Ueberschwingern42 16 ! 43 17 ! Revision 1.1 1999/02/05 09:17:16 raasch -
palm/trunk/SOURCE/sum_up_3d_data.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: sum_up_3d_data.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.1 2006/02/23 12:55:23 raasch 12 14 ! Initial revision -
palm/trunk/SOURCE/swap_timelevel.f90
r1 r3 1 1 SUBROUTINE swap_timelevel 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: swap_timelevel.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.8 2004/01/28 15:28:18 raasch 12 14 ! Swaping for Runge-Kutta schemes implemented 13 !14 ! Revision 1.7 2002/05/02 18:58:36 raasch15 ! Swaping of prognostic quantities re-changed to three timelevels16 !17 ! Revision 1.6 2001/08/21 10:01:58 10:01:58 raasch (Siegfried Raasch)18 ! Surface momentum fluxes are swapped also when they are prescribed19 !20 ! Revision 1.5 2001/03/30 07:54:39 raasch21 ! Swaping of prognostic quantities changed over from three to two timelevels,22 ! Translation of remaining German identifiers (variables, subroutines, etc.)23 !24 ! Revision 1.4 2001/01/29 12:36:03 raasch25 ! Passive scalar is considered26 !27 ! Revision 1.3 2001/01/25 07:26:48 raasch28 ! Sensible and latent heat fluxes are swaped, when use_surface_fluxes = .TRUE.,29 ! module test_variables removed30 !31 ! Revision 1.2 2000/04/13 13:59:26 schroeter32 ! + swapping timelevels of moisture variables33 15 ! 34 16 ! Revision 1.1 2000/01/10 10:08:58 10:08:58 raasch (Siegfried Raasch) … … 39 21 ! ------------ 40 22 ! Swap of timelevels of variables after each timestep 41 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 42 24 43 25 USE arrays_3d -
palm/trunk/SOURCE/temperton_fft.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: temperton_fft.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2003/04/16 12:49:25 raasch 12 14 ! Abort in case of illegal factors -
palm/trunk/SOURCE/time_integration.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: time_integration.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.8 2006/08/22 14:16:05 raasch 12 14 ! Disturbances are imposed only for the last Runge-Kutta-substep 13 !14 ! Revision 1.7 2006/08/04 15:04:39 raasch15 ! Output of particle timeseries16 !17 ! Revision 1.6 2006/02/23 12:58:38 raasch18 ! Data output (1d, 2d, 3d, sp) can be skipped for a given time interval from19 ! simulation start, additional argument nzb_.._inner added to call of20 ! disturb_field, nt_anz renamed current_timestep_number, pl.. renamed do..,21 ! ..anz renamed ..n, dt_average renamed dt_averaging_interval_pr,22 ! time_average_.. renamed time_do.._av, routines plot_.. renamed data_output_..,23 ! additional call of user_actions (after_integration)24 !25 ! Revision 1.5 2005/06/26 20:10:20 raasch26 ! +call of interaction_droplets_ptq (+module of the same name)27 !28 ! Revision 1.4 2005/04/23 09:39:01 raasch29 ! Revised calculation of output time counters regarding a possible decrease of30 ! the output time interval in case of restart runs31 !32 ! Revision 1.3 2005/03/26 21:10:26 raasch33 ! Pres is called by default only at the last Runge-Kutta-substep.34 ! Determination of the timestep moved from end to the beginning of the35 ! time loop in order to assure consistent output in run_control (first line36 ! in restart runs must be equal to the last line in the corresponding37 ! previous run). Exchange of ghost points for km and kh removed.38 ! Changes regarding non-cyclic boundary conditions:39 ! arguments are added to call of routines disturb_field, exchange_horiz,40 ! and boundary_conds, which is now called for a second time after pres in41 ! order to set the velocity boundary conditions at the outflow. Velocity42 ! perturbations are imposed near the inflow throughout the whole run.43 15 ! 44 16 ! Revision 1.2 2004/04/30 13:03:40 raasch … … 47 19 ! impulse_advec renamed momentum_advec 48 20 ! 49 ! Revision 1.1 2004/04/30 12:58:21 raasch50 ! Initial revision51 !52 ! Revision 1.44 2004/01/28 15:17:26 raasch53 ! Additional intermediate timesteps for the Runge-Kutta schemes, setting54 ! of the steering factors for prognostic equations within the new55 ! routine timestep_scheme_steering at the beginning of each step (former56 ! part of routine timestep called at the end of each step), routine57 ! run_control is not called for the Runge-Kutta schemes automatically in case58 ! of timestep changes.59 !60 ! Revision 1.43 2003/03/16 09:40:12 raasch61 ! Two underscores (_) are placed in front of all define-strings62 !63 ! Revision 1.42 2003/03/12 16:33:23 raasch64 ! New routine prognostic_equations_vec is used on NEC machines65 !66 ! Revision 1.41 2003/03/04 11:31:39 raasch67 ! Module DVRP added in case that dvrp graphics package is switched on68 !69 ! Revision 1.40 2002/12/19 15:45:13 raasch70 ! Routine check_cpu_time renamed to check_for_restart71 !72 ! Revision 1.39 2002/06/11 13:11:44 raasch73 ! +USE user_actions_mod, prognostic_equations_mod,74 ! temperature offset at cyclic boundaries in case of sloping surface moved75 ! from prognostic_equations to this subroutine. Particle advection moved from76 ! prognostic_equations to this routine. A fast version of prognostic_equations77 ! is called in case that the Piascek-Williams advection scheme is used.78 !79 ! Revision 1.38 2002/05/02 18:52:19 raasch80 ! Application of time filter and exchange of ghost points moved from81 ! prognostic_equations to this routine.82 ! Error in calculating output time intervals for various quantities removed.83 !84 ! Revision 1.37 2001/08/21 09:51:26 09:51:26 raasch (Siegfried Raasch)85 ! Time measurements for user_actions moved to that subroutine86 !87 ! Revision 1.36 2001/03/30 07:34:24 raasch88 ! Timelevel t+dt eliminated, Asselin filter now included in89 ! prognostic_equations, arguments removed from pres,90 ! Translation of remaining German identifiers (variables, subroutines, etc.)91 !92 ! Revision 1.35 2001/01/25 07:06:43 raasch93 ! Time measurements with dvrp software included94 !95 ! Revision 1.34 2001/01/05 15:12:50 raasch96 ! Calculation and plot of spectra97 !98 ! Revision 1.33 2000/04/27 06:40:49 raasch99 ! Old revision remarks deleted, output of isosurfaces and particles now100 ! combined into new routine data_output_dvrp (includes former routines101 ! plot_particles and plot_isosurface)102 !103 ! Revision 1.32 2000/04/13 13:39:30 schroeter104 ! computation of prognostic equations moved to the new routine105 ! prognostic_equation, new routines calc_liquid_water_content and106 ! compute_vpt107 !108 ! Revision 1.31 2000/01/25 16:25:24 16:25:24 letzel (Marcus Letzel)109 ! All comments translated into English110 !111 21 ! Revision 1.1 1997/08/11 06:19:04 raasch 112 22 ! Initial revision -
palm/trunk/SOURCE/time_to_string.f90
r1 r3 1 1 FUNCTION time_to_string( time ) 2 2 3 !------------------------------------------------------------------------------ -!3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- … … 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: time_to_string.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.3 2001/01/22 08:16:04 raasch 12 14 ! Comments translated into English 13 !14 ! Revision 1.2 1998/03/09 16:24:53 raasch15 ! Dreistellige Stundenausgabe ermoeglicht16 15 ! 17 16 ! Revision 1.1 1997/08/11 06:26:08 raasch … … 22 21 ! ------------ 23 22 ! Transforming the time from real to character-string hh:mm:ss 24 !------------------------------------------------------------------------------ -!23 !------------------------------------------------------------------------------! 25 24 26 25 IMPLICIT NONE -
palm/trunk/SOURCE/timestep.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: timestep.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.21 2006/02/23 12:59:44 raasch 12 14 ! nt_anz renamed current_timestep_number 13 !14 ! Revision 1.20 2005/04/23 09:45:38 raasch15 ! fcl_factor renamed cfl_factor16 !17 ! Revision 1.19 2005/03/26 21:11:53 raasch18 ! In case of non-cyclic lateral boundary conditions, the increased diffusivity19 ! in the lateral damping layer is regarded in the minimum timestep calculation.20 !21 ! Revision 1.18 2004/01/30 10:40:02 raasch22 ! Timestep increment limitation in case of Runge-Kutta removed23 !24 ! Revision 1.17 2004/01/28 15:34:03 raasch25 ! Setting of steering factors for the prognostic equations moved to new26 ! routine timestep_scheme_steering now called in routine leap_frog.27 !28 ! Revision 1.16 2003/03/16 09:49:47 raasch29 ! Two underscores (_) are placed in front of all define-strings30 !31 ! Revision 1.15 2002/12/19 16:19:54 raasch32 ! Two mpi_allreduce calls for determining mean horizontal velocity reduced to33 ! one call. Calculation of minimum timestep for diffusion optimized by34 ! determining minimum grid spacing in a single k loop.35 !36 ! Revision 1.14 2001/03/30 07:55:12 raasch37 ! Translation of remaining German identifiers (variables, subroutines, etc.)38 !39 ! Revision 1.13 2001/01/30 22:10:11 letzel40 ! All comments translated into English.41 !42 ! Revision 1.12 2001/01/22 08:17:33 raasch43 ! Module test_variables removed44 !45 ! Revision 1.11 2000/01/10 10:10:30 raasch46 ! Translationsgeschwindigkeit fuer Galilei-Transformation wird nur noch47 ! bedingt berechnet48 !49 ! Revision 1.10 1999/02/05 09:19:06 raasch50 ! Zeitschrittsteuerung beruecksichtigt jetzt auch Wahl des reinen Euler-51 ! Verfahrens52 !53 ! Revision 1.9 1998/09/22 17:30:41 raasch54 ! Ausfuehrlichere Informationen bei Abbruch durch Unterschreitung des55 ! minimalen Zeitschritts56 !57 ! Revision 1.8 1998/07/06 12:33:52 raasch58 ! Korrektur der ermittelten horizontalen Geschwindigkeitsmaxima, falls diese59 ! bei k=0 gefunden wurden (wegen Spiegelungsrandbedingung treten dort negative60 ! Werte auf, die bei Galilei-Transformation den Zeitschritt erheblich61 ! einschraenken). + USE test_variables62 !63 ! Revision 1.7 1998/04/21 15:57:46 raasch64 ! Implementierung der Galilei-Transformation65 !66 ! Revision 1.6 1998/03/25 13:57:20 raasch67 ! dt in dt_3d umbenannt68 !69 ! Revision 1.5 1998/03/09 16:25:10 raasch70 ! Schoenheitskorrekturen71 !72 ! Revision 1.4 1998/03/03 08:02:15 raasch73 ! Leap-Frog auch nach Zeitschrittaenderung erlaubt74 !75 ! Revision 1.3 1997/09/12 06:30:27 raasch76 ! Einbau des Diffusions-Kriteriums77 !78 ! Revision 1.2 1997/08/26 06:34:55 raasch79 ! Fehler bei w_max_ijk - Index80 15 ! 81 16 ! Revision 1.1 1997/08/11 06:26:19 raasch -
palm/trunk/SOURCE/timestep_scheme_steering.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: timestep_scheme_steering.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2005/03/26 21:17:06 raasch 12 14 ! No pressure term for Runge-Kutta-schemes (tsc(4)=0.0) -
palm/trunk/SOURCE/transpose.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: transpose.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.2 2004/04/30 13:12:17 raasch 12 14 ! Switched from mpi_alltoallv to the simpler mpi_alltoall, … … 17 19 ! Initial revision (collection of former routines transpose_xy, transpose_xz, 18 20 ! transpose_yx, transpose_yz, transpose_zx, transpose_zy) 19 !20 ! Revision 1.7 2003/03/16 09:49:53 raasch21 ! Two underscores (_) are placed in front of all define-strings22 !23 ! Revision 1.6 2001/03/30 07:55:50 raasch24 ! Array bound for work2 is explicitely calculated, all comments and identifiers25 ! translated into English26 !27 ! Revision 1.5 2001/01/22 08:19:14 raasch28 ! Module test_variables removed29 !30 ! Revision 1.4 1998/07/06 12:34:22 raasch31 ! + USE test_variables32 !33 ! Revision 1.3 1997/08/13 13:52:42 raasch34 ! Zeitmessung von MPI_ALLTOALLV35 !36 ! Revision 1.2 1997/08/11 06:27:31 raasch37 ! Keine Aenderungen38 21 ! 39 22 ! Revision 1.1 1997/07/24 11:25:18 raasch -
palm/trunk/SOURCE/user_interface.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: user_interface.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.18 2006/06/02 15:25:00 raasch 12 14 ! +change of grid-defining arguments in routine user_define_netcdf_grid, 13 15 ! new argument "found" in user_data_output_2d and user_data_output_3d 14 !15 ! Revision 1.17 2006/04/11 14:59:28 raasch16 ! Reading and writing user-defined quantities from/to the restart file improved17 !18 ! Revision 1.16 2006/03/03 20:00:48 raasch19 ! Error removed from last version20 !21 ! Revision 1.15 2006/03/03 19:55:36 raasch22 ! Call of user_init_grid now requires the formal parameter nzb_local23 !24 ! Revision 1.14 2006/02/23 13:13:18 raasch25 ! Additional interfaces for the calculation and output of user-defined 2d- and26 ! 3d-data27 !28 ! Revision 1.13 2003/03/12 16:44:56 raasch29 ! Full code replaced in routine user_actions (in the call for all gridpoints)30 ! instead of calling the _ij version (required by NEC, because otherwise no31 ! vectorization)32 !33 ! Revision 1.12 2002/12/19 16:22:04 raasch34 ! SAVE added in module user, STOP statement replaced by call of subroutine35 ! local_stop36 !37 ! Revision 1.11 2002/06/11 13:22:09 raasch38 ! Former subroutine changed to a module which allows to be called for all grid39 ! points of a single vertical column with index i,j or for all grid points by40 ! using function overloading.41 ! Log point for user_actions removed from the grid point version.42 !43 ! Revision 1.10 2001/08/21 12:13:55 raasch44 ! Wrong module in user_actions changed45 !46 ! Revision 1.9 2001/08/21 10:04:30 raasch47 ! User_actions can be called from every prognostic equation,48 ! subroutine user_statistics, user_particle_attributes, user_dvrp_coltab added49 !50 ! Revision 1.8 2001/07/12 12:45:07 raasch51 ! Subroutine user_init_particles added52 !53 ! Revision 1.7 2001/03/30 07:58:52 raasch54 ! Translation of remaining German identifiers (variables, subroutines, etc.)55 !56 ! Revision 1.6 2000/08/07 12:43:16 raasch57 ! All comments translated into English58 !59 ! Revision 1.5 1998/04/06 14:42:21 raasch60 ! +SUBROUTINE user_last_actions61 !62 ! Revision 1.4 1998/03/30 11:38:37 raasch63 ! Fehlerkorrekturen und Ausgabeaenderungen64 !65 ! Revision 1.3 1998/03/25 20:11:17 raasch66 ! +region in userpar67 !68 ! Revision 1.2 1998/03/25 13:57:37 raasch69 ! userpar beinhaltet jetzt keine vordefinierten Variablen mehr70 16 ! 71 17 ! Revision 1.1 1998/03/24 15:29:04 raasch -
palm/trunk/SOURCE/write_3d_binary.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: write_3d_binary.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.21 2006/08/04 15:05:11 raasch 12 14 ! +iran, iran_part 13 !14 ! Revision 1.20 2006/02/23 13:14:34 raasch15 ! +volume_flow_area, volume_flow_initial, -hom_sum_pr,16 ! output of arrays used for time-averaging17 !18 ! Revision 1.19 2005/05/18 15:55:45 raasch19 ! pl1d_callcount renamed dopr_time_count20 !21 ! Revision 1.18 2004/04/30 13:13:40 raasch22 ! Arrays of timelevel t-dt (_m) are written only in case of the leapfrog scheme23 !24 ! Revision 1.17 2003/10/29 09:14:11 raasch25 ! Output of random seed (random_iv, random_iy), module random_function_mod26 ! is used, new binary version 2.227 !28 ! Revision 1.16 2002/12/19 16:23:08 raasch29 ! Writing of array hom moved to write_var_list30 !31 ! Revision 1.15 2001/03/30 07:59:18 raasch32 ! Translation of remaining German identifiers (variables, subroutines, etc.),33 ! closing of unit 14 moved to the main program in order to allow user-defined34 ! writing to this unit35 !36 ! Revision 1.14 2001/03/11 11:27:18 letzel37 ! All comments translated into English.38 !39 ! Revision 1.13 2001/01/29 12:36:56 raasch40 ! Passive scalar is considered41 !42 ! Revision 1.12 2001/01/22 08:26:59 raasch43 ! Module test_variables removed44 !45 ! Revision 1.11 2001/01/05 15:16:39 raasch46 ! +spectrum_x, spectrum_y, new binary version 2.147 !48 ! Revision 1.10 2001/01/02 17:34:45 raasch49 ! Unit 14 closed at the end of the subroutine50 !51 ! Revision 1.9 2000/04/13 14:12:37 schroeter52 ! Output of water and humidity variables53 !54 ! Revision 1.8 99/03/03 09:27:43 09:27:43 raasch (Siegfried Raasch)55 ! Speziellen I/O-Messpunkt entfernt56 !57 ! Revision 1.7 1999/02/05 09:19:39 raasch58 ! program_version in binary_version umbenannt59 !60 ! Revision 1.6 1998/12/15 07:34:40 raasch61 ! Zusaetzliche Ermittelung der tatsaechlichen Ausgabezeiten auf T3E62 ! mit messpunkt_e(15) unter Verwendung von TIMEF (s. cpu_zeitmessung)63 !64 ! Revision 1.5 1998/07/06 12:36:26 raasch65 ! + USE test_variables66 !67 ! Revision 1.4 1998/04/15 11:24:30 raasch68 ! Steuerparameter-Ausgaben gehen jetzt auch auf Einheit 14 und werden von69 ! jedem PE ausgegeben70 !71 ! Revision 1.3 1998/04/06 14:34:21 raasch72 ! Ausgabe von weiteren Plotvariablen73 !74 ! Revision 1.2 1998/03/30 11:39:00 raasch75 ! Statistik- und Plotvariablen werden zusaetzlich ausgegeben76 15 ! 77 16 ! Revision 1.1 1998/03/18 20:20:21 raasch -
palm/trunk/SOURCE/write_compressed.f90
r1 r3 9 9 ! Former revisions: 10 10 ! --------------------- 11 ! $Log: write_compressed.f90,v $ 11 ! $Id$ 12 ! RCS Log replace by Id keyword, revision history cleaned up 13 ! 12 14 ! Revision 1.4 2006/02/23 13:15:09 raasch 13 15 ! nz_plot3d renamed nz_do3d 14 !15 ! Revision 1.3 2003/03/16 09:51:07 raasch16 ! Two underscores (_) are placed in front of all define-strings17 !18 ! Revision 1.2 2001/03/30 07:59:48 raasch19 ! All comments translated into English20 16 ! 21 17 ! Revision 1.1 1999/03/02 09:25:21 raasch -
palm/trunk/SOURCE/write_var_list.f90
r1 r3 8 8 ! Former revisions: 9 9 ! ----------------- 10 ! $Log: write_var_list.f90,v $ 10 ! $Id$ 11 ! RCS Log replace by Id keyword, revision history cleaned up 12 ! 11 13 ! Revision 1.34 2006/08/22 14:30:52 raasch 12 14 ! +dz_max 13 !14 ! Revision 1.33 2006/08/04 15:05:52 raasch15 ! +first_call_advec_particles, time_dopts, use_upstream_for_tke16 !17 ! Revision 1.32 2006/02/23 13:16:02 raasch18 ! +building parameters, conserve_volume_flow, dissipation_1d, e_min,19 ! mixing_length_1d, topography,20 ! -average_count_pl, time_average_pr, time_average_pl renamed time_dopr_av,21 ! nt_anz renamed current_timestep_number, ..pl.. renamed ..do..22 !23 ! Revision 1.31 2005/06/28 14:18:02 steinfeld24 ! +ug_surface, ug_vertical_gradient, ug_vertical_gradient_level,25 ! ug_vertical_gradient_level_ind, vg_surface, vg_vertical_gradient,26 ! vg_vertical_gradient_level, vg_vertical_gradient_level_ind27 !28 ! Revision 1.30 2005/06/26 20:10:54 raasch29 ! +cloud_droplets30 !31 ! Revision 1.29 2005/05/18 15:56:20 raasch32 ! +netcdf_precision33 !34 ! Revision 1.28 2005/03/26 21:18:39 raasch35 ! +bc_lr, bc_ns, inflow_disturbance_begin, inflow_disturbance_end, km_damp_max,36 ! outflow_damping_width37 !38 ! Revision 1.27 2004/04/30 13:14:31 raasch39 ! +grid_matching, impulse_advec renamed momentum_advec40 !41 ! Revision 1.26 2004/01/28 15:38:55 raasch42 ! at and bt replaced by array tsc43 !44 ! Revision 1.25 2003/10/29 09:15:12 raasch45 ! +last_dt_change46 !47 ! Revision 1.24 2003/03/14 13:47:55 raasch48 ! +random_generator49 !50 ! Revision 1.23 2002/12/19 16:24:09 raasch51 ! Writing of array hom moved from write_3d_binary to this routine,52 ! binary version incremented to 2.2, +time_restart53 !54 ! Revision 1.22 2002/04/16 08:22:02 raasch55 ! +bc_s_b, bc_s_t, surface_scalarflux, s_surface, s_surface_initial_change,56 ! s_vertical_gradient, s_vertical_gradient_level57 !58 ! Revision 1.21 2001/08/21 10:05:05 raasch59 ! +wall_adjustment60 !61 ! Revision 1.20 2001/03/30 08:00:07 raasch62 ! Translation of remaining German identifiers (variables, subroutines, etc.),63 ! runtime and package parameters are not written on unit 14 any more64 !65 ! Revision 1.19 2001/01/29 12:38:42 raasch66 ! Error writing use_surface_fluxes-string corrected, +passive_scalar67 !68 ! Revision 1.18 2001/01/25 07:27:48 raasch69 ! +fft_method, use_surface_fluxes70 !71 ! Revision 1.17 2001/01/05 15:16:59 raasch72 ! +average_count_sp, time_dosp73 !74 ! Revision 1.16 2000/12/28 13:39:21 raasch75 ! Module dvrp_variables removed, -particle_transport76 !77 ! Revision 1.15 2000/04/27 06:58:19 raasch78 ! +dt_dvrp, dvrp_filecount, npex, npey, threshold, time_dvrp79 ! -dt_plisos, dt_plpart, time_plisos, time_plpart80 ! +module pegrid,81 ! all comments translated into English, old revision remarks deleted82 !83 ! Revision 1.14 2000/04/13 14:10:58 schroeter84 ! Upgrade for calculations with cloud physics85 !86 ! Revision 1.13 2000/01/10 10:11:06 10:11:06 raasch (Siegfried Raasch)87 ! +use_ug_for_galilei_tr88 15 ! 89 16 ! Revision 1.1 1998/03/18 20:20:38 raasch
Note: See TracChangeset
for help on using the changeset viewer.