Changeset 2270 for palm/trunk/SOURCE/flow_statistics.f90
- Timestamp:
- Jun 9, 2017 12:18:47 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/flow_statistics.f90
r2252 r2270 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Revised numbering (removed 2 timeseries) 28 ! 29 ! 2252 2017-06-07 09:35:37Z knoop 27 30 ! perturbation pressure now depending on flux_output_mode 28 31 ! … … 399 402 IF ( humidity ) sums_l(:,49,i) = sums_wsqs_ws_l(:,i) & 400 403 * waterflux_output_conversion ! w*q* 401 IF ( passive_scalar ) sums_l(:,11 6,i) = sums_wsss_ws_l(:,i) ! w*s*404 IF ( passive_scalar ) sums_l(:,114,i) = sums_wsss_ws_l(:,i) ! w*s* 402 405 ENDDO 403 406 … … 484 487 DO j = nys, nyn 485 488 DO k = nzb, nzt+1 486 sums_l(k,11 7,tn) = sums_l(k,117,tn) + s(k,j,i) &489 sums_l(k,115,tn) = sums_l(k,115,tn) + s(k,j,i) & 487 490 * rmask(j,i,sr) & 488 491 * MERGE( 1.0_wp, 0.0_wp, & … … 512 515 ENDIF 513 516 IF ( passive_scalar ) THEN 514 sums_l(:,11 7,0) = sums_l(:,117,0) + sums_l(:,117,i)517 sums_l(:,115,0) = sums_l(:,115,0) + sums_l(:,115,i) 515 518 ENDIF 516 519 ENDDO … … 553 556 IF ( passive_scalar ) THEN 554 557 IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr ) 555 CALL MPI_ALLREDUCE( sums_l(nzb,11 7,0), sums(nzb,117), nzt+2-nzb, &558 CALL MPI_ALLREDUCE( sums_l(nzb,115,0), sums(nzb,115), nzt+2-nzb, & 556 559 MPI_REAL, MPI_SUM, comm2d, ierr ) 557 560 ENDIF … … 569 572 ENDIF 570 573 ENDIF 571 IF ( passive_scalar ) sums(:,11 7) = sums_l(:,117,0)574 IF ( passive_scalar ) sums(:,115) = sums_l(:,115,0) 572 575 #endif 573 576 … … 607 610 ! 608 611 !-- Passive scalar 609 IF ( passive_scalar ) hom(:,1,11 7,sr) = sums(:,117) / &612 IF ( passive_scalar ) hom(:,1,115,sr) = sums(:,115) / & 610 613 ngp_2dh_s_inner(:,sr) ! s 611 614 … … 654 657 ENDIF 655 658 IF ( passive_scalar ) THEN 656 sums_l(k,11 8,tn) = sums_l(k,118,tn) + &657 ( s(k,j,i)-hom(k,1,11 7,sr) )**2 * rmask(j,i,sr)&659 sums_l(k,116,tn) = sums_l(k,116,tn) + & 660 ( s(k,j,i)-hom(k,1,115,sr) )**2 * rmask(j,i,sr)& 658 661 * flag 659 662 ENDIF … … 889 892 !-- Passive scalar flux 890 893 IF ( passive_scalar ) THEN 891 sums_l(k,11 9,tn) = sums_l(k,119,tn) &894 sums_l(k,117,tn) = sums_l(k,117,tn) & 892 895 - 0.5_wp * ( kh(k,j,i) + kh(k+1,j,i) )& 893 896 * ( s(k+1,j,i) - s(k,j,i) ) & … … 949 952 ENDIF 950 953 IF ( passive_scalar ) THEN 951 sums_l(k+ki,11 9,tn) = sums_l(k+ki,119,tn) + &954 sums_l(k+ki,117,tn) = sums_l(k+ki,117,tn) + & 952 955 surf_def_h(l)%ssws(m) * rmask(j,i,sr) ! w"s" 953 956 ENDIF … … 1001 1004 ENDIF 1002 1005 IF ( passive_scalar ) THEN 1003 sums_l(nzb,11 9,tn) = sums_l(nzb,119,tn) + &1006 sums_l(nzb,117,tn) = sums_l(nzb,117,tn) + & 1004 1007 surf_lsm_h%ssws(m) * rmask(j,i,sr) ! w"s" 1005 1008 ENDIF … … 1051 1054 ENDIF 1052 1055 IF ( passive_scalar ) THEN 1053 sums_l(nzb,11 9,tn) = sums_l(nzb,119,tn) + &1056 sums_l(nzb,117,tn) = sums_l(nzb,117,tn) + & 1054 1057 surf_usm_h%ssws(m) * rmask(j,i,sr) ! w"s" 1055 1058 ENDIF … … 1063 1066 IF ( surf_def_h(0)%ns >= 1 ) THEN 1064 1067 m = surf_def_h(0)%start_index(j,i) 1065 sums_l(nzb,11 4,tn) = sums_l(nzb,114,tn) + &1068 sums_l(nzb,112,tn) = sums_l(nzb,112,tn) + & 1066 1069 surf_def_h(0)%ol(m) * rmask(j,i,sr) ! L 1067 1070 ENDIF 1068 1071 IF ( surf_lsm_h%ns >= 1 ) THEN 1069 1072 m = surf_lsm_h%start_index(j,i) 1070 sums_l(nzb,11 4,tn) = sums_l(nzb,114,tn) + &1073 sums_l(nzb,112,tn) = sums_l(nzb,112,tn) + & 1071 1074 surf_lsm_h%ol(m) * rmask(j,i,sr) ! L 1072 1075 ENDIF 1073 1076 IF ( surf_usm_h%ns >= 1 ) THEN 1074 1077 m = surf_usm_h%start_index(j,i) 1075 sums_l(nzb,11 4,tn) = sums_l(nzb,114,tn) + &1078 sums_l(nzb,112,tn) = sums_l(nzb,112,tn) + & 1076 1079 surf_usm_h%ol(m) * rmask(j,i,sr) ! L 1077 1080 ENDIF … … 1079 1082 1080 1083 IF ( radiation .AND. radiation_scheme /= 'constant' ) THEN 1081 sums_l(nzb, 101,tn) = sums_l(nzb,101,tn) + rad_net(j,i)1082 sums_l(nzb,10 2,tn) = sums_l(nzb,102,tn) + rad_lw_in(nzb,j,i)1083 sums_l(nzb,10 3,tn) = sums_l(nzb,103,tn) + rad_lw_out(nzb,j,i)1084 sums_l(nzb,10 4,tn) = sums_l(nzb,104,tn) + rad_sw_in(nzb,j,i)1085 sums_l(nzb,10 5,tn) = sums_l(nzb,105,tn) + rad_sw_out(nzb,j,i)1084 sums_l(nzb,99,tn) = sums_l(nzb,99,tn) + rad_net(j,i) 1085 sums_l(nzb,100,tn) = sums_l(nzb,100,tn) + rad_lw_in(nzb,j,i) 1086 sums_l(nzb,101,tn) = sums_l(nzb,101,tn) + rad_lw_out(nzb,j,i) 1087 sums_l(nzb,102,tn) = sums_l(nzb,102,tn) + rad_sw_in(nzb,j,i) 1088 sums_l(nzb,103,tn) = sums_l(nzb,103,tn) + rad_sw_out(nzb,j,i) 1086 1089 1087 1090 #if defined ( __rrtmg ) 1088 1091 IF ( radiation_scheme == 'rrtmg' ) THEN 1089 sums_l(nzb,1 10,tn) = sums_l(nzb,110,tn) + rrtm_aldif(0,j,i)1090 sums_l(nzb,1 11,tn) = sums_l(nzb,111,tn) + rrtm_aldir(0,j,i)1091 sums_l(nzb,11 2,tn) = sums_l(nzb,112,tn) + rrtm_asdif(0,j,i)1092 sums_l(nzb,11 3,tn) = sums_l(nzb,113,tn) + rrtm_asdir(0,j,i)1092 sums_l(nzb,108,tn) = sums_l(nzb,108,tn) + rrtm_aldif(0,j,i) 1093 sums_l(nzb,109,tn) = sums_l(nzb,109,tn) + rrtm_aldir(0,j,i) 1094 sums_l(nzb,110,tn) = sums_l(nzb,110,tn) + rrtm_asdif(0,j,i) 1095 sums_l(nzb,111,tn) = sums_l(nzb,111,tn) + rrtm_asdir(0,j,i) 1093 1096 ENDIF 1094 1097 #endif … … 1144 1147 ENDIF 1145 1148 IF ( passive_scalar ) THEN 1146 sums_l(nzt,11 9,tn) = sums_l(nzt,119,tn) + &1149 sums_l(nzt,117,tn) = sums_l(nzt,117,tn) + & 1147 1150 surf_def_h(2)%ssws(m) * rmask(j,i,sr) ! w"s" 1148 1151 ENDIF … … 1247 1250 IF ( passive_scalar .AND. ( .NOT. ws_scheme_sca & 1248 1251 .OR. sr /= 0 ) ) THEN 1249 pts = 0.5_wp * ( s(k,j,i) - hom(k,1,11 7,sr) + &1250 s(k+1,j,i) - hom(k+1,1,11 7,sr) )1251 sums_l(k,11 6,tn) = sums_l(k,116,tn) + pts * w(k,j,i) * &1252 pts = 0.5_wp * ( s(k,j,i) - hom(k,1,115,sr) + & 1253 s(k+1,j,i) - hom(k+1,1,115,sr) ) 1254 sums_l(k,114,tn) = sums_l(k,114,tn) + pts * w(k,j,i) * & 1252 1255 rmask(j,i,sr) * flag 1253 1256 ENDIF … … 1277 1280 IF ( i >= nxl .AND. i <= nxr .AND. & 1278 1281 j >= nys .AND. j <= nyn ) THEN 1279 sums_l(nzb,93,tn) = sums_l(nzb,93,tn) + surf_lsm_h%ghf_eb(m) 1280 sums_l(nzb,94,tn) = sums_l(nzb,94,tn) + surf_lsm_h%shf_eb(m) 1281 sums_l(nzb,95,tn) = sums_l(nzb,95,tn) + surf_lsm_h%qsws_eb(m) 1282 sums_l(nzb,96,tn) = sums_l(nzb,96,tn) + surf_lsm_h%qsws_liq_eb(m) 1283 sums_l(nzb,97,tn) = sums_l(nzb,97,tn) + surf_lsm_h%qsws_soil_eb(m) 1284 sums_l(nzb,98,tn) = sums_l(nzb,98,tn) + surf_lsm_h%qsws_veg_eb(m) 1285 sums_l(nzb,99,tn) = sums_l(nzb,99,tn) + surf_lsm_h%r_a(m) 1286 sums_l(nzb,100,tn) = sums_l(nzb,100,tn)+ surf_lsm_h%r_s(m) 1282 sums_l(nzb,93,tn) = sums_l(nzb,93,tn) + surf_lsm_h%ghf(m) 1283 sums_l(nzb,94,tn) = sums_l(nzb,94,tn) + surf_lsm_h%qsws_liq(m) 1284 sums_l(nzb,95,tn) = sums_l(nzb,95,tn) + surf_lsm_h%qsws_soil(m) 1285 sums_l(nzb,96,tn) = sums_l(nzb,96,tn) + surf_lsm_h%qsws_veg(m) 1286 sums_l(nzb,97,tn) = sums_l(nzb,97,tn) + surf_lsm_h%r_a(m) 1287 sums_l(nzb,98,tn) = sums_l(nzb,98,tn)+ surf_lsm_h%r_s(m) 1287 1288 ENDIF 1288 1289 ENDDO … … 1379 1380 ENDIF 1380 1381 IF ( passive_scalar ) THEN 1381 pts = 0.5_wp * ( s(k,j,i) - hom(k,1,11 7,sr) + &1382 s(k+1,j,i) - hom(k+1,1,11 7,sr) )1383 sums_l(k,11 6,tn) = sums_l(k,116,tn) + pts * w(k,j,i) * &1382 pts = 0.5_wp * ( s(k,j,i) - hom(k,1,115,sr) + & 1383 s(k+1,j,i) - hom(k+1,1,115,sr) ) 1384 sums_l(k,114,tn) = sums_l(k,114,tn) + pts * w(k,j,i) * & 1384 1385 rmask(j,i,sr) * flag 1385 1386 ENDIF … … 1581 1582 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 1582 1583 1583 sums_l(k,10 2,tn) = sums_l(k,102,tn) + rad_lw_in(k,j,i) &1584 sums_l(k,100,tn) = sums_l(k,100,tn) + rad_lw_in(k,j,i) & 1584 1585 * rmask(j,i,sr) * flag 1585 sums_l(k,10 3,tn) = sums_l(k,103,tn) + rad_lw_out(k,j,i) &1586 sums_l(k,101,tn) = sums_l(k,101,tn) + rad_lw_out(k,j,i) & 1586 1587 * rmask(j,i,sr) * flag 1587 sums_l(k,10 4,tn) = sums_l(k,104,tn) + rad_sw_in(k,j,i) &1588 sums_l(k,102,tn) = sums_l(k,102,tn) + rad_sw_in(k,j,i) & 1588 1589 * rmask(j,i,sr) * flag 1589 sums_l(k,10 5,tn) = sums_l(k,105,tn) + rad_sw_out(k,j,i) &1590 sums_l(k,103,tn) = sums_l(k,103,tn) + rad_sw_out(k,j,i) & 1590 1591 * rmask(j,i,sr) * flag 1591 sums_l(k,10 6,tn) = sums_l(k,106,tn) + rad_lw_cs_hr(k,j,i) &1592 sums_l(k,104,tn) = sums_l(k,104,tn) + rad_lw_cs_hr(k,j,i) & 1592 1593 * rmask(j,i,sr) * flag 1593 sums_l(k,10 7,tn) = sums_l(k,107,tn) + rad_lw_hr(k,j,i) &1594 sums_l(k,105,tn) = sums_l(k,105,tn) + rad_lw_hr(k,j,i) & 1594 1595 * rmask(j,i,sr) * flag 1595 sums_l(k,10 8,tn) = sums_l(k,108,tn) + rad_sw_cs_hr(k,j,i) &1596 sums_l(k,106,tn) = sums_l(k,106,tn) + rad_sw_cs_hr(k,j,i) & 1596 1597 * rmask(j,i,sr) * flag 1597 sums_l(k,10 9,tn) = sums_l(k,109,tn) + rad_sw_hr(k,j,i) &1598 sums_l(k,107,tn) = sums_l(k,107,tn) + rad_sw_hr(k,j,i) & 1598 1599 * rmask(j,i,sr) * flag 1599 1600 ENDDO … … 1655 1656 sums(k,55:63) = sums(k,55:63) / ngp_2dh(sr) 1656 1657 sums(k,81:88) = sums(k,81:88) / ngp_2dh(sr) 1657 sums(k,89:11 4) = sums(k,89:114) / ngp_2dh(sr)1658 sums(k,11 6) = sums(k,116) / ngp_2dh(sr)1659 sums(k,11 9) = sums(k,119) / ngp_2dh(sr)1658 sums(k,89:112) = sums(k,89:112) / ngp_2dh(sr) 1659 sums(k,114) = sums(k,114) / ngp_2dh(sr) 1660 sums(k,117) = sums(k,117) / ngp_2dh(sr) 1660 1661 IF ( ngp_2dh_s_inner(k,sr) /= 0 ) THEN 1661 1662 sums(k,8:11) = sums(k,8:11) / ngp_2dh_s_inner(k,sr) … … 1666 1667 sums(k,64) = sums(k,64) / ngp_2dh_s_inner(k,sr) 1667 1668 sums(k,70:80) = sums(k,70:80) / ngp_2dh_s_inner(k,sr) 1668 sums(k,11 8) = sums(k,118) / ngp_2dh_s_inner(k,sr)1669 sums(k,1 20:pr_palm-2) = sums(k,120:pr_palm-2) / ngp_2dh_s_inner(k,sr)1669 sums(k,116) = sums(k,116) / ngp_2dh_s_inner(k,sr) 1670 sums(k,118:pr_palm-2) = sums(k,118:pr_palm-2) / ngp_2dh_s_inner(k,sr) 1670 1671 ENDIF 1671 1672 ENDDO … … 1788 1789 hom(:,1,91,sr) = sums(:,91) ! m_soil 1789 1790 ! 92 is initial m_soil profile 1790 hom(:,1,93,sr) = sums(:,93) ! ghf_eb 1791 hom(:,1,94,sr) = sums(:,94) ! shf_eb 1792 hom(:,1,95,sr) = sums(:,95) ! qsws_eb 1793 hom(:,1,96,sr) = sums(:,96) ! qsws_liq_eb 1794 hom(:,1,97,sr) = sums(:,97) ! qsws_soil_eb 1795 hom(:,1,98,sr) = sums(:,98) ! qsws_veg_eb 1796 hom(:,1,99,sr) = sums(:,99) ! r_a 1797 hom(:,1,100,sr) = sums(:,100) ! r_s 1791 hom(:,1,93,sr) = sums(:,93) ! ghf 1792 hom(:,1,94,sr) = sums(:,94) ! qsws_liq 1793 hom(:,1,95,sr) = sums(:,95) ! qsws_soil 1794 hom(:,1,96,sr) = sums(:,96) ! qsws_veg 1795 hom(:,1,97,sr) = sums(:,97) ! r_a 1796 hom(:,1,98,sr) = sums(:,98) ! r_s 1798 1797 1799 1798 ENDIF 1800 1799 1801 1800 IF ( radiation ) THEN 1802 hom(:,1, 101,sr) = sums(:,101) ! rad_net1803 hom(:,1,10 2,sr) = sums(:,102) ! rad_lw_in1804 hom(:,1,10 3,sr) = sums(:,103) ! rad_lw_out1805 hom(:,1,10 4,sr) = sums(:,104) ! rad_sw_in1806 hom(:,1,10 5,sr) = sums(:,105) ! rad_sw_out1801 hom(:,1,99,sr) = sums(:,99) ! rad_net 1802 hom(:,1,100,sr) = sums(:,100) ! rad_lw_in 1803 hom(:,1,101,sr) = sums(:,101) ! rad_lw_out 1804 hom(:,1,102,sr) = sums(:,102) ! rad_sw_in 1805 hom(:,1,103,sr) = sums(:,103) ! rad_sw_out 1807 1806 1808 1807 IF ( radiation_scheme == 'rrtmg' ) THEN 1809 hom(:,1,10 6,sr) = sums(:,106) ! rad_lw_cs_hr1810 hom(:,1,10 7,sr) = sums(:,107) ! rad_lw_hr1811 hom(:,1,10 8,sr) = sums(:,108) ! rad_sw_cs_hr1812 hom(:,1,10 9,sr) = sums(:,109) ! rad_sw_hr1813 1814 hom(:,1,1 10,sr) = sums(:,110) ! rrtm_aldif1815 hom(:,1,1 11,sr) = sums(:,111) ! rrtm_aldir1816 hom(:,1,11 2,sr) = sums(:,112) ! rrtm_asdif1817 hom(:,1,11 3,sr) = sums(:,113) ! rrtm_asdir1808 hom(:,1,104,sr) = sums(:,104) ! rad_lw_cs_hr 1809 hom(:,1,105,sr) = sums(:,105) ! rad_lw_hr 1810 hom(:,1,106,sr) = sums(:,106) ! rad_sw_cs_hr 1811 hom(:,1,107,sr) = sums(:,107) ! rad_sw_hr 1812 1813 hom(:,1,108,sr) = sums(:,108) ! rrtm_aldif 1814 hom(:,1,109,sr) = sums(:,109) ! rrtm_aldir 1815 hom(:,1,110,sr) = sums(:,110) ! rrtm_asdif 1816 hom(:,1,111,sr) = sums(:,111) ! rrtm_asdir 1818 1817 ENDIF 1819 1818 ENDIF 1820 1819 1821 hom(:,1,11 4,sr) = sums(:,114) !: L1820 hom(:,1,112,sr) = sums(:,112) !: L 1822 1821 1823 1822 IF ( passive_scalar ) THEN 1824 hom(:,1,11 9,sr) = sums(:,119) ! w"s"1825 hom(:,1,11 6,sr) = sums(:,116) ! w*s*1826 hom(:,1,1 20,sr) = sums(:,119) + sums(:,116) ! ws1827 hom(:,1,11 8,sr) = sums(:,118) ! s*21828 ENDIF 1829 1830 hom(:,1,1 21,sr) = rho_air ! rho_air in Kg/m^31831 hom(:,1,12 2,sr) = rho_air_zw ! rho_air_zw in Kg/m^31823 hom(:,1,117,sr) = sums(:,117) ! w"s" 1824 hom(:,1,114,sr) = sums(:,114) ! w*s* 1825 hom(:,1,118,sr) = sums(:,117) + sums(:,114) ! ws 1826 hom(:,1,116,sr) = sums(:,116) ! s*2 1827 ENDIF 1828 1829 hom(:,1,119,sr) = rho_air ! rho_air in Kg/m^3 1830 hom(:,1,120,sr) = rho_air_zw ! rho_air_zw in Kg/m^3 1832 1831 1833 1832 hom(:,1,pr_palm,sr) = sums(:,pr_palm) … … 1972 1971 1973 1972 IF ( .NOT. neutral ) THEN 1974 ts_value(22,sr) = hom(nzb,1,11 4,sr) ! L1973 ts_value(22,sr) = hom(nzb,1,112,sr) ! L 1975 1974 ELSE 1976 1975 ts_value(22,sr) = 1.0E10_wp … … 1980 1979 1981 1980 IF ( passive_scalar ) THEN 1982 ts_value(24,sr) = hom(nzb+13,1,11 9,sr) ! w"s" ( to do ! )1981 ts_value(24,sr) = hom(nzb+13,1,117,sr) ! w"s" ( to do ! ) 1983 1982 ts_value(25,sr) = hom(nzb+13,1,pr_palm,sr) ! s* 1984 1983 ENDIF … … 1987 1986 !-- Collect land surface model timeseries 1988 1987 IF ( land_surface ) THEN 1989 ts_value(dots_soil ,sr) = hom(nzb,1,93,sr) ! ghf_eb 1990 ts_value(dots_soil+1,sr) = hom(nzb,1,94,sr) ! shf_eb 1991 ts_value(dots_soil+2,sr) = hom(nzb,1,95,sr) ! qsws_eb 1992 ts_value(dots_soil+3,sr) = hom(nzb,1,96,sr) ! qsws_liq_eb 1993 ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr) ! qsws_soil_eb 1994 ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr) ! qsws_veg_eb 1995 ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr) ! r_a 1996 ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr) ! r_s 1988 ts_value(dots_soil ,sr) = hom(nzb,1,93,sr) ! ghf 1989 ts_value(dots_soil+1,sr) = hom(nzb,1,94,sr) ! qsws_liq 1990 ts_value(dots_soil+2,sr) = hom(nzb,1,95,sr) ! qsws_soil 1991 ts_value(dots_soil+3,sr) = hom(nzb,1,96,sr) ! qsws_veg 1992 ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr) ! r_a 1993 ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr) ! r_s 1997 1994 ENDIF 1998 1995 ! 1999 1996 !-- Collect radiation model timeseries 2000 1997 IF ( radiation ) THEN 2001 ts_value(dots_rad,sr) = hom(nzb,1, 101,sr)! rad_net2002 ts_value(dots_rad+1,sr) = hom(nzb,1,10 2,sr) ! rad_lw_in2003 ts_value(dots_rad+2,sr) = hom(nzb,1,10 3,sr) ! rad_lw_out2004 ts_value(dots_rad+3,sr) = hom(nzb,1,10 4,sr) ! rad_sw_in2005 ts_value(dots_rad+4,sr) = hom(nzb,1,10 5,sr) ! rad_sw_out1998 ts_value(dots_rad,sr) = hom(nzb,1,99,sr) ! rad_net 1999 ts_value(dots_rad+1,sr) = hom(nzb,1,100,sr) ! rad_lw_in 2000 ts_value(dots_rad+2,sr) = hom(nzb,1,101,sr) ! rad_lw_out 2001 ts_value(dots_rad+3,sr) = hom(nzb,1,102,sr) ! rad_sw_in 2002 ts_value(dots_rad+4,sr) = hom(nzb,1,103,sr) ! rad_sw_out 2006 2003 2007 2004 IF ( radiation_scheme == 'rrtmg' ) THEN 2008 ts_value(dots_rad+5,sr) = hom(nzb,1,1 10,sr) ! rrtm_aldif2009 ts_value(dots_rad+6,sr) = hom(nzb,1,1 11,sr) ! rrtm_aldir2010 ts_value(dots_rad+7,sr) = hom(nzb,1,11 2,sr) ! rrtm_asdif2011 ts_value(dots_rad+8,sr) = hom(nzb,1,11 3,sr) ! rrtm_asdir2005 ts_value(dots_rad+5,sr) = hom(nzb,1,108,sr) ! rrtm_aldif 2006 ts_value(dots_rad+6,sr) = hom(nzb,1,109,sr) ! rrtm_aldir 2007 ts_value(dots_rad+7,sr) = hom(nzb,1,110,sr) ! rrtm_asdif 2008 ts_value(dots_rad+8,sr) = hom(nzb,1,111,sr) ! rrtm_asdir 2012 2009 ENDIF 2013 2010
Note: See TracChangeset
for help on using the changeset viewer.