Changeset 3 for palm/trunk/SOURCE/production_e.f90
- Timestamp:
- Feb 13, 2007 11:30:58 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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
Note: See TracChangeset
for help on using the changeset viewer.