Changeset 1658 for palm/trunk/SOURCE/flow_statistics.f90
 Timestamp:
 Sep 18, 2015 10:52:53 AM (6 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/flow_statistics.f90
r1655 r1658 21 21 ! Current revisions: 22 22 !  23 ! 23 ! bugfix: reduction variables in the openacc branch are now initialzed to zero 24 24 ! 25 25 ! Former revisions: … … 1781 1781 !$acc parallel loop gang present( pt, rflags_invers, rmask, sums_l, u, v ) create( s1, s2, s3 ) 1782 1782 DO k = nzb, nzt+1 1783 s1 = 0 1784 s2 = 0 1785 s3 = 0 1783 1786 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3 ) 1784 1787 DO i = nxl, nxr … … 1803 1806 !$acc parallel loop gang present( rflags_invers, rmask, sums_l, sa ) create( s1 ) 1804 1807 DO k = nzb, nzt+1 1808 s1 = 0 1805 1809 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 1806 1810 DO i = nxl, nxr … … 1823 1827 !$acc parallel loop gang present( q, rflags_invers, rmask, sums_l, vpt ) create( s1, s2 ) 1824 1828 DO k = nzb, nzt+1 1829 s1 = 0 1830 s2 = 0 1825 1831 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 1826 1832 DO i = nxl, nxr … … 1839 1845 !$acc parallel loop gang present( pt, q, ql, rflags_invers, rmask, sums_l ) create( s1, s2 ) 1840 1846 DO k = nzb, nzt+1 1847 s1 = 0 1848 s2 = 0 1841 1849 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 1842 1850 DO i = nxl, nxr … … 1861 1869 !$acc parallel loop gang present( q, rflags_invers, rmask, sums_l ) create( s1 ) 1862 1870 DO k = nzb, nzt+1 1871 s1 = 0 1863 1872 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 1864 1873 DO i = nxl, nxr … … 2048 2057 !$acc parallel loop gang present( e, hom, kh, km, p, pt, w, rflags_invers, rmask, sums_l ) create( s1, s2, s3, s4, s5, s6, s7 ) 2049 2058 DO k = nzb, nzt+1 2059 s1 = 0 2060 s2 = 0 2061 s3 = 0 2062 s4 = 0 2063 s5 = 0 2064 s6 = 0 2065 s7 = 0 2050 2066 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3, s4, s5, s6, s7 ) 2051 2067 DO i = nxl, nxr … … 2080 2096 !$acc parallel loop gang present( hom, q, rflags_invers, rmask, sums_l ) create( s1 ) 2081 2097 DO k = nzb, nzt+1 2098 s1 = 0 2082 2099 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2083 2100 DO i = nxl, nxr … … 2095 2112 ! Total and perturbation energy for the total domain (being 2096 2113 ! collected in the last column of sums_l). 2114 s1 = 0 2097 2115 !$OMP DO 2098 2116 !$acc parallel loop collapse(3) present( rflags_invers, rmask, u, v, w ) reduction(+:s1) … … 2113 2131 !$OMP DO 2114 2132 !$acc parallel present( rmask, sums_l, us, usws, vsws, ts ) create( s1, s2, s3, s4 ) 2133 s1 = 0 2134 s2 = 0 2135 s3 = 0 2136 s4 = 0 2115 2137 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3, s4 ) 2116 2138 DO i = nxl, nxr … … 2132 2154 IF ( humidity ) THEN 2133 2155 !$acc parallel present( qs, rmask, sums_l ) create( s1 ) 2156 s1 = 0 2134 2157 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2135 2158 DO i = nxl, nxr … … 2150 2173 !$acc parallel loop gang present( u, v, w, rflags_invers, rmask, sums_l ) create( s1, s2, s3, s4, ust2, vst2, w2 ) 2151 2174 DO k = nzb, nzt+1 2175 s1 = 0 2176 s2 = 0 2177 s3 = 0 2178 s4 = 0 2152 2179 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3, s4 ) 2153 2180 DO i = nxl, nxr … … 2176 2203 !$OMP DO 2177 2204 !$acc parallel present( sums_l ) create( s1 ) 2205 s1 = 0 2178 2206 !$acc loop reduction( +: s1 ) 2179 2207 DO k = nzb, nzt+1 … … 2203 2231 !$acc parallel loop gang present( ddzu, kh, km, pt, u, v, w, rflags_invers, rmask, sums_l ) create( s1, s2, s3 ) 2204 2232 DO k = nzb, nzt_diff 2233 s1 = 0 2234 s2 = 0 2235 s3 = 0 2205 2236 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3 ) 2206 2237 DO i = nxl, nxr … … 2244 2275 !$acc parallel loop gang present( ddzu, kh, sa, rflags_invers, rmask, sums_l ) create( s1 ) 2245 2276 DO k = nzb, nzt_diff 2277 s1 = 0 2246 2278 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2247 2279 DO i = nxl, nxr … … 2264 2296 !$acc parallel loop gang present( ddzu, kh, q, vpt, rflags_invers, rmask, sums_l ) create( s1, s2 ) 2265 2297 DO k = nzb, nzt_diff 2298 s1 = 0 2299 s2 = 0 2266 2300 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2267 2301 DO i = nxl, nxr … … 2286 2320 !$acc parallel loop gang present( ddzu, kh, q, ql, rflags_invers, rmask, sums_l ) create( s1 ) 2287 2321 DO k = nzb, nzt_diff 2322 s1 = 0 2288 2323 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2289 2324 DO i = nxl, nxr … … 2309 2344 !$acc parallel loop gang present( ddzu, kh, q, rflags_invers, rmask, sums_l ) create( s1 ) 2310 2345 DO k = nzb, nzt_diff 2346 s1 = 0 2311 2347 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2312 2348 DO i = nxl, nxr … … 2328 2364 !$OMP DO 2329 2365 !$acc parallel present( rmask, shf, sums_l, usws, vsws ) create( s1, s2, s3, s4, s5 ) 2366 s1 = 0 2367 s2 = 0 2368 s3 = 0 2369 s4 = 0 2370 s5 = 0 2330 2371 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3, s4, s5 ) 2331 2372 DO i = nxl, nxr … … 2351 2392 !$OMP DO 2352 2393 !$acc parallel present( rmask, saswsb, sums_l ) create( s1 ) 2394 s1 = 0 2353 2395 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2354 2396 DO i = nxl, nxr … … 2366 2408 !$OMP DO 2367 2409 !$acc parallel present( pt, q, qsws, rmask, shf, sums_l ) create( s1, s2 ) 2410 s1 = 0 2411 s2 = 0 2368 2412 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2369 2413 DO i = nxl, nxr … … 2382 2426 !$OMP DO 2383 2427 !$acc parallel present( pt, q, ql, qsws, rmask, shf, sums_l ) create( s1 ) 2428 s1 = 0 2384 2429 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2385 2430 DO i = nxl, nxr … … 2399 2444 !$OMP DO 2400 2445 !$acc parallel present( qsws, rmask, sums_l ) create( s1 ) 2446 s1 = 0 2401 2447 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2402 2448 DO i = nxl, nxr … … 2418 2464 !$OMP DO 2419 2465 !$acc parallel present( qsws, rmask, sums_l ) create( s1 ) 2466 s1 = 0 2420 2467 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2421 2468 DO i = nxl, nxr … … 2437 2484 !$OMP DO 2438 2485 !$acc parallel present( rmask, sums_l, tswst, uswst, vswst ) create( s1, s2, s3, s4, s5 ) 2486 s1 = 0 2487 s2 = 0 2488 s3 = 0 2489 s4 = 0 2490 s5 = 0 2439 2491 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3, s4, s5 ) 2440 2492 DO i = nxl, nxr … … 2458 2510 !$OMP DO 2459 2511 !$acc parallel present( rmask, saswst, sums_l ) create( s1 ) 2512 s1 = 0 2460 2513 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2461 2514 DO i = nxl, nxr … … 2473 2526 !$OMP DO 2474 2527 !$acc parallel present( pt, q, qswst, rmask, tswst, sums_l ) create( s1, s2 ) 2528 s1 = 0 2529 s2 = 0 2475 2530 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2476 2531 DO i = nxl, nxr … … 2489 2544 !$OMP DO 2490 2545 !$acc parallel present( pt, q, ql, qswst, rmask, tswst, sums_l ) create( s1 ) 2546 s1 = 0 2491 2547 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2492 2548 DO i = nxl, nxr … … 2507 2563 !$OMP DO 2508 2564 !$acc parallel present( qswst, rmask, sums_l ) create( s1 ) 2565 s1 = 0 2509 2566 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2510 2567 DO i = nxl, nxr … … 2526 2583 !$OMP DO 2527 2584 !$acc parallel present( qswst, rmask, sums_l ) create( s1 ) 2585 s1 = 0 2528 2586 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2529 2587 DO i = nxl, nxr … … 2546 2604 !$acc parallel loop gang present( hom, pt, rflags_invers, rmask, sums_l, u, v, w ) create( s1, s2, s3 ) 2547 2605 DO k = nzb, nzt_diff 2606 s1 = 0 2607 s2 = 0 2608 s3 = 0 2548 2609 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3 ) 2549 2610 DO i = nxl, nxr … … 2580 2641 !$acc parallel loop gang present( hom, rflags_invers, rmask, sa, sums_l, w ) create( s1 ) 2581 2642 DO k = nzb, nzt_diff 2643 s1 = 0 2582 2644 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2583 2645 DO i = nxl, nxr … … 2597 2659 !$acc parallel loop gang present( rflags_invers, rho, prho, rmask, sums_l ) create( s1, s2 ) 2598 2660 DO k = nzb, nzt_diff 2661 s1 = 0 2662 s2 = 0 2599 2663 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2600 2664 DO i = nxl, nxr … … 2620 2684 !$acc parallel loop gang present( hom, rflags_invers, rmask, sums_l, vpt, w ) create( s1 ) 2621 2685 DO k = nzb, nzt_diff 2686 s1 = 0 2622 2687 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2623 2688 DO i = nxl, nxr … … 2636 2701 !$acc parallel loop gang present( hom, q, ql, rflags_invers, rmask, sums_l, w ) create( s1 ) 2637 2702 DO k = nzb, nzt_diff 2703 s1 = 0 2638 2704 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2639 2705 DO i = nxl, nxr … … 2652 2718 !$acc parallel loop gang present( qc, ql, rflags_invers, rmask, sums_l ) create( s1, s2 ) 2653 2719 DO k = nzb, nzt_diff 2720 s1 = 0 2654 2721 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2655 2722 DO i = nxl, nxr … … 2668 2735 !$acc parallel loop gang present( nr, qr, prr, rflags_invers, rmask, sums_l ) create( s1, s2, s3 ) 2669 2736 DO k = nzb, nzt_diff 2737 s1 = 0 2738 s2 = 0 2739 s3 = 0 2670 2740 !$acc loop vector collapse( 2 ) reduction( +: s1, s2, s3 ) 2671 2741 DO i = nxl, nxr … … 2688 2758 !$acc parallel loop gang present( ql, rflags_invers, rmask, sums_l ) create( s1 ) 2689 2759 DO k = nzb, nzt_diff 2760 s1 = 0 2690 2761 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2691 2762 DO i = nxl, nxr … … 2704 2775 !$acc parallel loop gang present( ql, rflags_invers, rmask, sums_l ) create( s1 ) 2705 2776 DO k = nzb, nzt_diff 2777 s1 = 0 2706 2778 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2707 2779 DO i = nxl, nxr … … 2722 2794 !$acc parallel loop gang present( hom, rflags_invers, rmask, sums_l, vpt, w ) create( s1 ) 2723 2795 DO k = nzb, nzt_diff 2796 s1 = 0 2724 2797 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2725 2798 DO i = nxl, nxr … … 2754 2827 !$acc parallel loop gang present( hom, q, rflags_invers, rmask, sums_l, w ) create( s1 ) 2755 2828 DO k = nzb, nzt_diff 2829 s1 = 0 2756 2830 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2757 2831 DO i = nxl, nxr … … 2778 2852 !$acc parallel loop gang present( hom, rflags_invers, rmask, sums_l, u, v, w ) create( s1, s2 ) 2779 2853 DO k = nzb, nzt_diff 2854 s1 = 0 2855 s2 = 0 2780 2856 !$acc loop vector collapse( 2 ) reduction( +: s1, s2 ) 2781 2857 DO i = nxl, nxr … … 2809 2885 !$acc parallel loop gang present( hom, pt, rflags_invers, rmask, sums_l, w ) create( s1 ) 2810 2886 DO k = nzb, nzt_diff 2887 s1 = 0 2811 2888 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2812 2889 DO i = nxl, nxr … … 2828 2905 !$acc parallel loop gang present( hom, q, rflags_invers, rmask, sums_l, w ) create( s1 ) 2829 2906 DO k = nzb, nzt_diff 2907 s1 = 0 2830 2908 !$acc loop vector collapse( 2 ) reduction( +: s1 ) 2831 2909 DO i = nxl, nxr
Note: See TracChangeset
for help on using the changeset viewer.