Changeset 1337 for palm/trunk/SOURCE/prognostic_equations.f90
- Timestamp:
- Mar 25, 2014 3:11:48 PM (11 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/prognostic_equations.f90 (modified) (68 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/prognostic_equations.f90
r1333 r1337 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! Bugfix: REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 299 299 IF ( .NOT. outflow_l .OR. i > nxl ) THEN 300 300 301 tend(:,j,i) = 0.0 301 tend(:,j,i) = 0.0_wp 302 302 IF ( timestep_scheme(1:5) == 'runge' ) THEN 303 303 IF ( ws_scheme_mom ) THEN … … 354 354 intermediate_timestep_count_max ) THEN 355 355 DO k = nzb_u_inner(j,i)+1, nzt 356 tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tu_m(k,j,i)356 tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i) 357 357 ENDDO 358 358 ENDIF … … 365 365 IF ( .NOT. outflow_s .OR. j > nys ) THEN 366 366 367 tend(:,j,i) = 0.0 367 tend(:,j,i) = 0.0_wp 368 368 IF ( timestep_scheme(1:5) == 'runge' ) THEN 369 369 IF ( ws_scheme_mom ) THEN … … 413 413 intermediate_timestep_count_max ) THEN 414 414 DO k = nzb_v_inner(j,i)+1, nzt 415 tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tv_m(k,j,i)415 tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i) 416 416 ENDDO 417 417 ENDIF … … 422 422 ! 423 423 !-- Tendency terms for w-velocity component 424 tend(:,j,i) = 0.0 424 tend(:,j,i) = 0.0_wp 425 425 IF ( timestep_scheme(1:5) == 'runge' ) THEN 426 426 IF ( ws_scheme_mom ) THEN … … 471 471 intermediate_timestep_count_max ) THEN 472 472 DO k = nzb_w_inner(j,i)+1, nzt-1 473 tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tw_m(k,j,i)473 tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i) 474 474 ENDDO 475 475 ENDIF … … 484 484 ! 485 485 !-- Tendency terms for potential temperature 486 tend(:,j,i) = 0.0 486 tend(:,j,i) = 0.0_wp 487 487 IF ( timestep_scheme(1:5) == 'runge' ) THEN 488 488 IF ( ws_scheme_sca ) THEN … … 516 516 ! 517 517 !-- Consideration of heat sources within the plant canopy 518 IF ( plant_canopy .AND. cthf /= 0.0 ) THEN518 IF ( plant_canopy .AND. cthf /= 0.0_wp ) THEN 519 519 CALL plant_canopy_model( i, j, 4 ) 520 520 ENDIF … … 551 551 intermediate_timestep_count_max ) THEN 552 552 DO k = nzb_s_inner(j,i)+1, nzt 553 tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + &554 5.3125 * tpt_m(k,j,i)553 tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 554 5.3125_wp * tpt_m(k,j,i) 555 555 ENDDO 556 556 ENDIF … … 565 565 ! 566 566 !-- Tendency-terms for salinity 567 tend(:,j,i) = 0.0 567 tend(:,j,i) = 0.0_wp 568 568 IF ( timestep_scheme(1:5) == 'runge' ) & 569 569 THEN … … 588 588 - tsc(5) * rdf_sc(k) * & 589 589 ( sa(k,j,i) - sa_init(k) ) 590 IF ( sa_p(k,j,i) < 0.0 ) sa_p(k,j,i) = 0.1* sa(k,j,i)590 IF ( sa_p(k,j,i) < 0.0_wp ) sa_p(k,j,i) = 0.1_wp * sa(k,j,i) 591 591 ENDDO 592 592 … … 601 601 intermediate_timestep_count_max ) THEN 602 602 DO k = nzb_s_inner(j,i)+1, nzt 603 tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + &604 5.3125 * tsa_m(k,j,i)603 tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 604 5.3125_wp * tsa_m(k,j,i) 605 605 ENDDO 606 606 ENDIF … … 620 620 ! 621 621 !-- Tendency-terms for total water content / scalar 622 tend(:,j,i) = 0.0 622 tend(:,j,i) = 0.0_wp 623 623 IF ( timestep_scheme(1:5) == 'runge' ) & 624 624 THEN … … 666 666 - tsc(5) * rdf_sc(k) * & 667 667 ( q(k,j,i) - q_init(k) ) 668 IF ( q_p(k,j,i) < 0.0 ) q_p(k,j,i) = 0.1* q(k,j,i)668 IF ( q_p(k,j,i) < 0.0_wp ) q_p(k,j,i) = 0.1_wp * q(k,j,i) 669 669 ENDDO 670 670 … … 679 679 intermediate_timestep_count_max ) THEN 680 680 DO k = nzb_s_inner(j,i)+1, nzt 681 tq_m(k,j,i) = -9.5625 * tend(k,j,i) + &682 5.3125 * tq_m(k,j,i)681 tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 682 5.3125_wp * tq_m(k,j,i) 683 683 ENDDO 684 684 ENDIF … … 692 692 ! 693 693 !-- Calculate prognostic equation for rain water content 694 tend(:,j,i) = 0.0 694 tend(:,j,i) = 0.0_wp 695 695 IF ( timestep_scheme(1:5) == 'runge' ) & 696 696 THEN … … 718 718 tsc(3) * tqr_m(k,j,i) ) & 719 719 - tsc(5) * rdf_sc(k) * qr(k,j,i) 720 IF ( qr_p(k,j,i) < 0.0 ) qr_p(k,j,i) = 0.0720 IF ( qr_p(k,j,i) < 0.0_wp ) qr_p(k,j,i) = 0.0_wp 721 721 ENDDO 722 722 ! … … 730 730 intermediate_timestep_count_max ) THEN 731 731 DO k = nzb_s_inner(j,i)+1, nzt 732 tqr_m(k,j,i) = -9.5625 * tend(k,j,i) + &733 5.3125 * tqr_m(k,j,i)732 tqr_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 733 5.3125_wp * tqr_m(k,j,i) 734 734 ENDDO 735 735 ENDIF … … 738 738 ! 739 739 !-- Calculate prognostic equation for rain drop concentration. 740 tend(:,j,i) = 0.0 740 tend(:,j,i) = 0.0_wp 741 741 IF ( timestep_scheme(1:5) == 'runge' ) THEN 742 742 IF ( ws_scheme_sca ) THEN … … 764 764 tsc(3) * tnr_m(k,j,i) ) & 765 765 - tsc(5) * rdf_sc(k) * nr(k,j,i) 766 IF ( nr_p(k,j,i) < 0.0 ) nr_p(k,j,i) = 0.0766 IF ( nr_p(k,j,i) < 0.0_wp ) nr_p(k,j,i) = 0.0_wp 767 767 ENDDO 768 768 ! … … 776 776 intermediate_timestep_count_max ) THEN 777 777 DO k = nzb_s_inner(j,i)+1, nzt 778 tnr_m(k,j,i) = -9.5625 * tend(k,j,i) + &779 5.3125 * tnr_m(k,j,i)778 tnr_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 779 5.3125_wp * tnr_m(k,j,i) 780 780 ENDDO 781 781 ENDIF … … 793 793 ! 794 794 !-- Tendency-terms for TKE 795 tend(:,j,i) = 0.0 795 tend(:,j,i) = 0.0_wp 796 796 IF ( timestep_scheme(1:5) == 'runge' & 797 797 .AND. .NOT. use_upstream_for_tke ) THEN … … 830 830 e_p(k,j,i) = e(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) + & 831 831 tsc(3) * te_m(k,j,i) ) 832 IF ( e_p(k,j,i) < 0.0 ) e_p(k,j,i) = 0.1* e(k,j,i)832 IF ( e_p(k,j,i) < 0.0_wp ) e_p(k,j,i) = 0.1_wp * e(k,j,i) 833 833 ENDDO 834 834 … … 843 843 intermediate_timestep_count_max ) THEN 844 844 DO k = nzb_s_inner(j,i)+1, nzt 845 te_m(k,j,i) = -9.5625 * tend(k,j,i) + &846 5.3125 * te_m(k,j,i)845 te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 846 5.3125_wp * te_m(k,j,i) 847 847 ENDDO 848 848 ENDIF … … 880 880 CALL cpu_log( log_point(5), 'u-equation', 'start' ) 881 881 882 tend = 0.0 882 tend = 0.0_wp 883 883 IF ( timestep_scheme(1:5) == 'runge' ) THEN 884 884 IF ( ws_scheme_mom ) THEN … … 946 946 DO j = nys, nyn 947 947 DO k = nzb_u_inner(j,i)+1, nzt 948 tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tu_m(k,j,i)948 tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i) 949 949 ENDDO 950 950 ENDDO … … 959 959 CALL cpu_log( log_point(6), 'v-equation', 'start' ) 960 960 961 tend = 0.0 961 tend = 0.0_wp 962 962 IF ( timestep_scheme(1:5) == 'runge' ) THEN 963 963 IF ( ws_scheme_mom ) THEN … … 1022 1022 DO j = nysv, nyn 1023 1023 DO k = nzb_v_inner(j,i)+1, nzt 1024 tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tv_m(k,j,i)1024 tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i) 1025 1025 ENDDO 1026 1026 ENDDO … … 1094 1094 DO j = nys, nyn 1095 1095 DO k = nzb_w_inner(j,i)+1, nzt-1 1096 tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tw_m(k,j,i)1096 tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i) 1097 1097 ENDDO 1098 1098 ENDDO … … 1120 1120 sbt = 1.0 1121 1121 ENDIF 1122 tend = 0.0 1122 tend = 0.0_wp 1123 1123 CALL advec_s_bc( pt, 'pt' ) 1124 1124 … … 1128 1128 !-- pt-tendency terms with no communication 1129 1129 IF ( scalar_advec /= 'bc-scheme' ) THEN 1130 tend = 0.0 1130 tend = 0.0_wp 1131 1131 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1132 1132 IF ( ws_scheme_sca ) THEN … … 1156 1156 ! 1157 1157 !-- Consideration of heat sources within the plant canopy 1158 IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN1158 IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN 1159 1159 CALL plant_canopy_model( 4 ) 1160 1160 ENDIF … … 1201 1201 DO j = nys, nyn 1202 1202 DO k = nzb_s_inner(j,i)+1, nzt 1203 tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + &1204 5.3125 * tpt_m(k,j,i)1203 tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 1204 5.3125_wp * tpt_m(k,j,i) 1205 1205 ENDDO 1206 1206 ENDDO … … 1229 1229 sbt = 1.0 1230 1230 ENDIF 1231 tend = 0.0 1231 tend = 0.0_wp 1232 1232 CALL advec_s_bc( sa, 'sa' ) 1233 1233 … … 1262 1262 - tsc(5) * rdf_sc(k) * & 1263 1263 ( sa(k,j,i) - sa_init(k) ) 1264 IF ( sa_p(k,j,i) < 0.0 ) sa_p(k,j,i) = 0.1* sa(k,j,i)1264 IF ( sa_p(k,j,i) < 0.0_wp ) sa_p(k,j,i) = 0.1_wp * sa(k,j,i) 1265 1265 ENDDO 1266 1266 ENDDO … … 1283 1283 DO j = nys, nyn 1284 1284 DO k = nzb_s_inner(j,i)+1, nzt 1285 tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + &1286 5.3125 * tsa_m(k,j,i)1285 tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + & 1286 5.3125_wp * tsa_m(k,j,i) 1287 1287 ENDDO 1288 1288 ENDDO … … 1315 1315 ! 1316 1316 !-- Bott-Chlond scheme always uses Euler time step. Thus: 1317 sbt = 1.0 1318 ENDIF 1319 tend = 0.0 1317 sbt = 1.0_wp 1318 ENDIF 1319 tend = 0.0_wp 1320 1320 CALL advec_s_bc( q, 'q' ) 1321 1321 … … 1325 1325 !-- Scalar/q-tendency terms with no communication 1326 1326 IF ( scalar_advec /= 'bc-scheme' ) THEN 1327 tend = 0.0 1327 tend = 0.0_wp 1328 1328 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1329 1329 IF ( ws_scheme_sca ) THEN … … 1371 1371 - tsc(5) * rdf_sc(k) * & 1372 1372 ( q(k,j,i) - q_init(k) ) 1373 IF ( q_p(k,j,i) < 0.0 ) q_p(k,j,i) = 0.1* q(k,j,i)1373 IF ( q_p(k,j,i) < 0.0_wp ) q_p(k,j,i) = 0.1_wp * q(k,j,i) 1374 1374 ENDDO 1375 1375 ENDDO … … 1392 1392 DO j = nys, nyn 1393 1393 DO k = nzb_s_inner(j,i)+1, nzt 1394 tq_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tq_m(k,j,i)1394 tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tq_m(k,j,i) 1395 1395 ENDDO 1396 1396 ENDDO … … 1417 1417 ! 1418 1418 !-- Bott-Chlond scheme always uses Euler time step. Thus: 1419 sbt = 1.0 1420 ENDIF 1421 tend = 0.0 1419 sbt = 1.0_wp 1420 ENDIF 1421 tend = 0.0_wp 1422 1422 CALL advec_s_bc( e, 'e' ) 1423 1423 … … 1429 1429 IF ( scalar_advec /= 'bc-scheme' .OR. use_upstream_for_tke ) THEN 1430 1430 IF ( use_upstream_for_tke ) THEN 1431 tend = 0.0 1431 tend = 0.0_wp 1432 1432 CALL advec_s_up( e ) 1433 1433 ELSE 1434 tend = 0.0 1434 tend = 0.0_wp 1435 1435 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1436 1436 IF ( ws_scheme_sca ) THEN … … 1472 1472 e_p(k,j,i) = e(k,j,i) + dt_3d * ( sbt * tend(k,j,i) + & 1473 1473 tsc(3) * te_m(k,j,i) ) 1474 IF ( e_p(k,j,i) < 0.0 ) e_p(k,j,i) = 0.1* e(k,j,i)1474 IF ( e_p(k,j,i) < 0.0_wp ) e_p(k,j,i) = 0.1_wp * e(k,j,i) 1475 1475 ENDDO 1476 1476 ENDDO … … 1493 1493 DO j = nys, nyn 1494 1494 DO k = nzb_s_inner(j,i)+1, nzt 1495 te_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* te_m(k,j,i)1495 te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * te_m(k,j,i) 1496 1496 ENDDO 1497 1497 ENDDO … … 1545 1545 CALL advec_u_ws_acc 1546 1546 ELSE 1547 tend = 0.0 1547 tend = 0.0_wp ! to be removed later?? 1548 1548 CALL advec_u_pw 1549 1549 ENDIF … … 1597 1597 tu_m(k,j,i) = tend(k,j,i) 1598 1598 ELSEIF ( runge_step == 2 ) THEN 1599 tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tu_m(k,j,i)1599 tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i) 1600 1600 ENDIF 1601 1601 ENDIF … … 1615 1615 CALL advec_v_ws_acc 1616 1616 ELSE 1617 tend = 0.0 ! to be removed later??1617 tend = 0.0_wp ! to be removed later?? 1618 1618 CALL advec_v_pw 1619 1619 END IF … … 1664 1664 tv_m(k,j,i) = tend(k,j,i) 1665 1665 ELSEIF ( runge_step == 2 ) THEN 1666 tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tv_m(k,j,i)1666 tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i) 1667 1667 ENDIF 1668 1668 ENDIF … … 1682 1682 CALL advec_w_ws_acc 1683 1683 ELSE 1684 tend = 0.0 ! to be removed later??1684 tend = 0.0_wp ! to be removed later?? 1685 1685 CALL advec_w_pw 1686 1686 ENDIF … … 1727 1727 tw_m(k,j,i) = tend(k,j,i) 1728 1728 ELSEIF ( runge_step == 2 ) THEN 1729 tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tw_m(k,j,i)1729 tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i) 1730 1730 ENDIF 1731 1731 ENDIF … … 1754 1754 sbt = 1.0 1755 1755 ENDIF 1756 tend = 0.0 1756 tend = 0.0_wp 1757 1757 CALL advec_s_bc( pt, 'pt' ) 1758 1758 … … 1762 1762 !-- pt-tendency terms with no communication 1763 1763 IF ( scalar_advec /= 'bc-scheme' ) THEN 1764 tend = 0.0 1764 tend = 0.0_wp 1765 1765 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1766 1766 IF ( ws_scheme_sca ) THEN 1767 1767 CALL advec_s_ws_acc( pt, 'pt' ) 1768 1768 ELSE 1769 tend = 0.0 ! to be removed later??1769 tend = 0.0_wp ! to be removed later?? 1770 1770 CALL advec_s_pw( pt ) 1771 1771 ENDIF … … 1791 1791 ! 1792 1792 !-- Consideration of heat sources within the plant canopy 1793 IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN1793 IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN 1794 1794 CALL plant_canopy_model( 4 ) 1795 1795 ENDIF … … 1827 1827 tpt_m(k,j,i) = tend(k,j,i) 1828 1828 ELSEIF ( runge_step == 2 ) THEN 1829 tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tpt_m(k,j,i)1829 tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tpt_m(k,j,i) 1830 1830 ENDIF 1831 1831 ENDIF … … 1853 1853 ! 1854 1854 !-- Bott-Chlond scheme always uses Euler time step. Thus: 1855 sbt = 1.0 1856 ENDIF 1857 tend = 0.0 1855 sbt = 1.0_wp 1856 ENDIF 1857 tend = 0.0_wp 1858 1858 CALL advec_s_bc( sa, 'sa' ) 1859 1859 … … 1863 1863 !-- sa-tendency terms with no communication 1864 1864 IF ( scalar_advec /= 'bc-scheme' ) THEN 1865 tend = 0.0 1865 tend = 0.0_wp 1866 1866 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1867 1867 IF ( ws_scheme_sca ) THEN … … 1888 1888 - tsc(5) * rdf_sc(k) * & 1889 1889 ( sa(k,j,i) - sa_init(k) ) 1890 IF ( sa_p(k,j,i) < 0.0 ) sa_p(k,j,i) = 0.1* sa(k,j,i)1890 IF ( sa_p(k,j,i) < 0.0_wp ) sa_p(k,j,i) = 0.1_wp * sa(k,j,i) 1891 1891 ! 1892 1892 !-- Tendencies for the next Runge-Kutta step … … 1894 1894 tsa_m(k,j,i) = tend(k,j,i) 1895 1895 ELSEIF ( runge_step == 2 ) THEN 1896 tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tsa_m(k,j,i)1896 tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tsa_m(k,j,i) 1897 1897 ENDIF 1898 1898 ENDDO … … 1924 1924 ! 1925 1925 !-- Bott-Chlond scheme always uses Euler time step. Thus: 1926 sbt = 1.0 1927 ENDIF 1928 tend = 0.0 1926 sbt = 1.0_wp 1927 ENDIF 1928 tend = 0.0_wp 1929 1929 CALL advec_s_bc( q, 'q' ) 1930 1930 … … 1934 1934 !-- Scalar/q-tendency terms with no communication 1935 1935 IF ( scalar_advec /= 'bc-scheme' ) THEN 1936 tend = 0.0 1936 tend = 0.0_wp 1937 1937 IF ( timestep_scheme(1:5) == 'runge' ) THEN 1938 1938 IF ( ws_scheme_sca ) THEN … … 1980 1980 - tsc(5) * rdf_sc(k) * & 1981 1981 ( q(k,j,i) - q_init(k) ) 1982 IF ( q_p(k,j,i) < 0.0 ) q_p(k,j,i) = 0.1* q(k,j,i)1982 IF ( q_p(k,j,i) < 0.0_wp ) q_p(k,j,i) = 0.1_wp * q(k,j,i) 1983 1983 ! 1984 1984 !-- Tendencies for the next Runge-Kutta step … … 1986 1986 tq_m(k,j,i) = tend(k,j,i) 1987 1987 ELSEIF ( runge_step == 2 ) THEN 1988 tq_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* tq_m(k,j,i)1988 tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tq_m(k,j,i) 1989 1989 ENDIF 1990 1990 ENDDO … … 2010 2010 ! 2011 2011 !-- Bott-Chlond scheme always uses Euler time step. Thus: 2012 sbt = 1.0 2013 ENDIF 2014 tend = 0.0 2012 sbt = 1.0_wp 2013 ENDIF 2014 tend = 0.0_wp 2015 2015 CALL advec_s_bc( e, 'e' ) 2016 2016 … … 2022 2022 IF ( scalar_advec /= 'bc-scheme' .OR. use_upstream_for_tke ) THEN 2023 2023 IF ( use_upstream_for_tke ) THEN 2024 tend = 0.0 2024 tend = 0.0_wp 2025 2025 CALL advec_s_up( e ) 2026 2026 ELSE … … 2029 2029 CALL advec_s_ws_acc( e, 'e' ) 2030 2030 ELSE 2031 tend = 0.0 ! to be removed later??2031 tend = 0.0_wp ! to be removed later?? 2032 2032 CALL advec_s_pw( e ) 2033 2033 ENDIF 2034 2034 ELSE 2035 tend = 0.0 ! to be removed later??2035 tend = 0.0_wp ! to be removed later?? 2036 2036 CALL advec_s_up( e ) 2037 2037 ENDIF … … 2071 2071 e_p(k,j,i) = e(k,j,i) + dt_3d * ( sbt * tend(k,j,i) + & 2072 2072 tsc(3) * te_m(k,j,i) ) 2073 IF ( e_p(k,j,i) < 0.0 ) e_p(k,j,i) = 0.1* e(k,j,i)2073 IF ( e_p(k,j,i) < 0.0_wp ) e_p(k,j,i) = 0.1_wp * e(k,j,i) 2074 2074 ! 2075 2075 !-- Tendencies for the next Runge-Kutta step … … 2077 2077 te_m(k,j,i) = tend(k,j,i) 2078 2078 ELSEIF ( runge_step == 2 ) THEN 2079 te_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125* te_m(k,j,i)2079 te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * te_m(k,j,i) 2080 2080 ENDIF 2081 2081 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
