Changeset 4180 for palm/trunk/SOURCE/advec_ws.f90
 Timestamp:
 Aug 21, 2019 2:37:54 PM (2 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/advec_ws.f90
r4110 r4180 77 77 ! OpenACC port for SPEC 78 78 ! 79 ! 3589 20181130 15:09:51Z suehring80 ! Move the control parameter "salsa" from salsa_mod to control_parameters81 ! (M. Kurppa)82 !83 ! 3582 20181129 19:16:36Z suehring84 !  Computation of vertical fluxes separated from computation of horizontal85 ! fluxes. Loops are splitted in order to avoid indirect indexing and allow86 ! for better vectorization.87 !  Accelerate code by remove typeconversions of ibits88 !  replace pointer definition in scalar routine by simple explicit definition89 ! of the passed array90 !91 ! 3547 20181121 13:21:24Z suehring92 ! variables documented93 !94 ! 3467 20181030 19:05:21Z suehring95 ! Implementation of a new aerosol module salsa.96 !97 ! 3302 20181003 02:39:40Z raasch98 ! Stokes drift velocity added in scalar advection99 !100 ! 3298 20181002 12:21:11Z kanani101 ! Add formatted note from ketelsen about statistics for chemistry102 !103 ! 3294 20181001 02:37:10Z raasch104 ! ocean renamed ocean_mode105 !106 ! 3274 20180924 15:42:55Z knoop107 ! Modularization of all bulk cloud physics code components108 !109 ! 3241 20180912 15:02:00Z raasch110 ! unused variables removed111 !112 ! 3183 20180727 14:25:55Z suehring113 ! Rename variables for boundarycondition flags and for offline nesting mode114 !115 ! 3182 20180727 13:36:03Z suehring116 ! Bugfix in calculation of leftsided fluxes for ucomponent in OpenMP case.117 !118 ! 2731 20180109 17:44:02Z suehring119 ! Enable loop vectorization by splitting the kloop120 !121 ! 2718 20180102 08:49:38Z maronga122 ! Corrected "Former revisions" section123 !124 ! 2696 20171214 17:12:51Z kanani125 ! Change in file header (GPL part)126 ! Implement advection for TKEdissipation in case of RANSTKEe closure (TG)127 ! Allocate advc_flags_1/2 within ws_init_flags instead of init_grid128 ! Change argument list for exchange_horiz_2d_int (MS)129 !130 ! 2329 20170803 14:24:56Z knoop131 ! Bugfix concerning density in divergence correction close to buildings132 !133 ! 2292 20170620 09:51:42Z schwenkel134 ! Implementation of new microphysic scheme: cloud_scheme = 'morrison'135 ! includes two more prognostic equations for cloud drop concentration (nc)136 ! and cloud water content (qc).137 !138 ! 2233 20170530 18:08:54Z suehring139 !140 ! 2232 20170530 17:47:52Z suehring141 ! Rename wall_flags_0 and wall_flags_00 into advc_flags_m and advc_flags_1,142 ! respectively.143 ! Set advc_flags_1/2 on basis of wall_flags_0/00 instead of nzb_s/u/v/w_inner.144 ! Setting advc_flags_1/2 also for downwardfacing walls145 !146 ! 2200 20170411 11:37:51Z suehring147 ! monotonic_adjustment removed148 !149 ! 2118 20170117 16:38:49Z raasch150 ! OpenACC version of subroutines removed151 !152 ! 2037 20161026 11:15:40Z knoop153 ! Anelastic approximation implemented154 !155 ! 2000 20160820 18:09:15Z knoop156 ! Forced header and separation lines into 80 columns157 !158 ! 1996 20160818 11:42:29Z suehring159 ! Bugfix concerning calculation of turbulent of turbulent fluxes160 !161 ! 1960 20160712 16:34:24Z suehring162 ! Separate humidity and passive scalar163 !164 ! 1942 20160614 12:18:18Z suehring165 ! Initialization of flags for wsscheme moved from init_grid.166 !167 ! 1873 20160418 14:50:06Z maronga168 ! Module renamed (removed _mod)169 !170 !171 ! 1850 20160408 13:29:27Z maronga172 ! Module renamed173 !174 !175 ! 1822 20160407 07:49:42Z hoffmann176 ! icloud_scheme removed, microphysics_seifert added177 !178 ! 1682 20151007 23:56:08Z knoop179 ! Code annotations made doxygen readable180 !181 ! 1630 20150826 16:57:23Z suehring182 !183 !184 ! 1629 20150826 16:56:11Z suehring185 ! Bugfix concerning wall_flags at left and south PE boundaries186 !187 ! 1581 20150410 13:45:59Z suehring188 !189 !190 ! 1580 20150410 13:43:49Z suehring191 ! Bugfix: statistical evaluation of scalar fluxes in case of monotonic limiter192 !193 ! 1567 20150310 17:57:55Z suehring194 ! Bugfixes in monotonic limiter.195 !196 ! 20150309 13:10:37Z heinze197 ! Bugfix: REAL constants provided with KINDattribute in call of198 ! intrinsic functions like MAX and MIN199 !200 ! 1557 20150305 16:43:04Z suehring201 ! Enable monotone advection for scalars using monotonic limiter202 !203 ! 1374 20140425 12:55:07Z raasch204 ! missing variables added to ONLY list205 !206 ! 1361 20140416 15:17:48Z hoffmann207 ! accelerator and vector version for qr and nr added208 !209 ! 1353 20140408 15:21:23Z heinze210 ! REAL constants provided with KINDattribute,211 ! module kinds added212 ! some formatting adjustments213 !214 ! 1322 20140320 16:38:49Z raasch215 ! REAL constants defined as wpkind216 !217 ! 1320 20140320 08:40:49Z raasch218 ! ONLYattribute added to USEstatements,219 ! kindparameters added to all INTEGER and REAL declaration statements,220 ! kinds are defined in new module kinds,221 ! old module precision_kind is removed,222 ! revision history before 2012 removed,223 ! comment fields (!:) to be used for variable explanations added to224 ! all variable declaration statements225 !226 ! 1257 20131108 15:18:40Z raasch227 ! accelerator loop directives removed228 !229 ! 1221 20130910 08:59:13Z raasch230 ! wall_flags_00 introduced, which holds bits 32...231 !232 ! 1128 20130412 06:19:32Z raasch233 ! loop index bounds in accelerator version replaced by i_left, i_right, j_south,234 ! j_north235 !236 ! 1115 20130326 18:16:16Z hoffmann237 ! calculation of qr and nr is restricted to precipitation238 !239 ! 1053 20121113 17:11:03Z hoffmann240 ! necessary expansions according to the two new prognostic equations (nr, qr)241 ! of the twomoment cloud physics scheme:242 ! +flux_l_*, flux_s_*, diss_l_*, diss_s_*, sums_ws*s_ws_l243 !244 ! 1036 20121022 13:43:42Z raasch245 ! code put under GPL (PALM 3.9)246 !247 ! 1027 20121015 17:18:39Z suehring248 ! Bugfix in calculation indices k_mm, k_pp in accelerator version249 !250 ! 1019 20120928 06:46:45Z raasch251 ! small change in comment lines252 !253 ! 1015 20120927 09:23:24Z raasch254 ! accelerator versions (*_acc) added255 !256 ! 1010 20120920 07:59:54Z raasch257 ! cpp switch __nopointer added for pointer free version258 !259 ! 888 20120420 15:03:46Z suehring260 ! Number of REAL( IBITS() calls with identical arguments is reduced.261 !262 ! 862 20120326 14:21:38Z suehring263 ! wsscheme also work with topography in combination with vector version.264 ! wsscheme also work with outflow boundaries in combination with265 ! vector version.266 ! Degradation of the applied order of scheme is now steered by multiplying with267 ! Integer advc_flags_1. 2nd order scheme, WS3 and WS5 are calculated on each268 ! grid point and mulitplied with the appropriate flag.269 ! 2nd order numerical dissipation term changed. Now the appropriate 2nd order270 ! term derived according to the 4th and 6th order terms is applied. It turns271 ! out that diss_2nd does not provide sufficient dissipation near walls.272 ! Therefore, the function diss_2nd is removed.273 ! Near walls a divergence correction is necessary to overcome numerical274 ! instabilities due to too less divergence reduction of the velocity field.275 ! boundary_flags and logicals steering the degradation are removed.276 ! Empty SUBROUTINE local_diss removed.277 ! Further formatting adjustments.278 !279 ! 801 20120110 17:30:36Z suehring280 ! Bugfix concerning OpenMP parallelization. Summation of sums_wsus_ws_l,281 ! sums_wsvs_ws_l, sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wspts_ws_l,282 ! sums_wsqs_ws_l, sums_wssas_ws_l is now threadsafe by adding an additional283 ! dimension.284 !285 ! Initial revision286 !287 ! 411 20091211 12:31:43 Z suehring288 !289 !290 !291 ! @author Matthias Suehring292 !293 79 ! 294 80 ! Description:
Note: See TracChangeset
for help on using the changeset viewer.