Ignore:
Timestamp:
Jun 9, 2017 12:18:47 PM (4 years ago)
Author:
maronga
Message:

major revisions in land surface model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r2252 r2270  
    2525! -----------------
    2626! $Id$
     27! Revised numbering (removed 2 timeseries)
     28!
     29! 2252 2017-06-07 09:35:37Z knoop
    2730! perturbation pressure now depending on flux_output_mode
    2831!
     
    399402             IF ( humidity       ) sums_l(:,49,i)  = sums_wsqs_ws_l(:,i)       &
    400403                                           * waterflux_output_conversion  ! w*q*
    401              IF ( passive_scalar ) sums_l(:,116,i) = sums_wsss_ws_l(:,i)  ! w*s*
     404             IF ( passive_scalar ) sums_l(:,114,i) = sums_wsss_ws_l(:,i)  ! w*s*
    402405          ENDDO
    403406
     
    484487             DO  j =  nys, nyn
    485488                DO  k = nzb, nzt+1
    486                    sums_l(k,117,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)             &
    487490                                    * rmask(j,i,sr)                            &
    488491                                    * MERGE( 1.0_wp, 0.0_wp,                   &
     
    512515             ENDIF
    513516             IF ( passive_scalar )  THEN
    514                 sums_l(:,117,0) = sums_l(:,117,0) + sums_l(:,117,i)
     517                sums_l(:,115,0) = sums_l(:,115,0) + sums_l(:,115,i)
    515518             ENDIF
    516519          ENDDO
     
    553556       IF ( passive_scalar )  THEN
    554557          IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    555           CALL MPI_ALLREDUCE( sums_l(nzb,117,0), sums(nzb,117), nzt+2-nzb,       &
     558          CALL MPI_ALLREDUCE( sums_l(nzb,115,0), sums(nzb,115), nzt+2-nzb,       &
    556559                              MPI_REAL, MPI_SUM, comm2d, ierr )
    557560       ENDIF
     
    569572          ENDIF
    570573       ENDIF
    571        IF ( passive_scalar )  sums(:,117) = sums_l(:,117,0)
     574       IF ( passive_scalar )  sums(:,115) = sums_l(:,115,0)
    572575#endif
    573576
     
    607610!
    608611!--    Passive scalar
    609        IF ( passive_scalar )  hom(:,1,117,sr) = sums(:,117) /                  &
     612       IF ( passive_scalar )  hom(:,1,115,sr) = sums(:,115) /                  &
    610613            ngp_2dh_s_inner(:,sr)                    ! s
    611614
     
    654657                ENDIF
    655658                IF ( passive_scalar )  THEN
    656                    sums_l(k,118,tn) = sums_l(k,118,tn) + &
    657                                   ( s(k,j,i)-hom(k,1,117,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)&
    658661                                                                  * flag
    659662                ENDIF
     
    889892!--             Passive scalar flux
    890893                IF ( passive_scalar )  THEN
    891                    sums_l(k,119,tn) = sums_l(k,119,tn)                         &
     894                   sums_l(k,117,tn) = sums_l(k,117,tn)                         &
    892895                                         - 0.5_wp * ( kh(k,j,i) + kh(k+1,j,i) )&
    893896                                                  * ( s(k+1,j,i) - s(k,j,i) )  &
     
    949952                         ENDIF
    950953                         IF ( passive_scalar )  THEN
    951                             sums_l(k+ki,119,tn) = sums_l(k+ki,119,tn) +                     &
     954                            sums_l(k+ki,117,tn) = sums_l(k+ki,117,tn) +                     &
    952955                                        surf_def_h(l)%ssws(m) * rmask(j,i,sr) ! w"s"
    953956                         ENDIF
     
    10011004                   ENDIF
    10021005                   IF ( passive_scalar )  THEN
    1003                       sums_l(nzb,119,tn) = sums_l(nzb,119,tn) +                     &
     1006                      sums_l(nzb,117,tn) = sums_l(nzb,117,tn) +                     &
    10041007                                        surf_lsm_h%ssws(m) * rmask(j,i,sr) ! w"s"
    10051008                   ENDIF
     
    10511054                   ENDIF
    10521055                   IF ( passive_scalar )  THEN
    1053                       sums_l(nzb,119,tn) = sums_l(nzb,119,tn) +                     &
     1056                      sums_l(nzb,117,tn) = sums_l(nzb,117,tn) +                     &
    10541057                                        surf_usm_h%ssws(m) * rmask(j,i,sr) ! w"s"
    10551058                   ENDIF
     
    10631066                IF ( surf_def_h(0)%ns >= 1 )  THEN
    10641067                   m = surf_def_h(0)%start_index(j,i)
    1065                    sums_l(nzb,114,tn) = sums_l(nzb,114,tn) +                      &
     1068                   sums_l(nzb,112,tn) = sums_l(nzb,112,tn) +                      &
    10661069                                        surf_def_h(0)%ol(m)  * rmask(j,i,sr) ! L
    10671070                ENDIF
    10681071                IF ( surf_lsm_h%ns >= 1 )  THEN
    10691072                   m = surf_lsm_h%start_index(j,i)
    1070                    sums_l(nzb,114,tn) = sums_l(nzb,114,tn) +                      &
     1073                   sums_l(nzb,112,tn) = sums_l(nzb,112,tn) +                      &
    10711074                                        surf_lsm_h%ol(m)  * rmask(j,i,sr) ! L
    10721075                ENDIF
    10731076                IF ( surf_usm_h%ns >= 1 )  THEN
    10741077                   m = surf_usm_h%start_index(j,i)
    1075                    sums_l(nzb,114,tn) = sums_l(nzb,114,tn) +                      &
     1078                   sums_l(nzb,112,tn) = sums_l(nzb,112,tn) +                      &
    10761079                                        surf_usm_h%ol(m)  * rmask(j,i,sr) ! L
    10771080                ENDIF
     
    10791082
    10801083             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,102,tn)  = sums_l(nzb,102,tn)  + rad_lw_in(nzb,j,i)
    1083                 sums_l(nzb,103,tn)  = sums_l(nzb,103,tn)  + rad_lw_out(nzb,j,i)
    1084                 sums_l(nzb,104,tn)  = sums_l(nzb,104,tn)  + rad_sw_in(nzb,j,i)
    1085                 sums_l(nzb,105,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)
    10861089
    10871090#if defined ( __rrtmg )
    10881091                IF ( radiation_scheme == 'rrtmg' )  THEN
    1089                    sums_l(nzb,110,tn)  = sums_l(nzb,110,tn)  + rrtm_aldif(0,j,i)
    1090                    sums_l(nzb,111,tn)  = sums_l(nzb,111,tn)  + rrtm_aldir(0,j,i)
    1091                    sums_l(nzb,112,tn)  = sums_l(nzb,112,tn)  + rrtm_asdif(0,j,i)
    1092                    sums_l(nzb,113,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)
    10931096                ENDIF
    10941097#endif
     
    11441147                ENDIF
    11451148                IF ( passive_scalar )  THEN
    1146                    sums_l(nzt,119,tn) = sums_l(nzt,119,tn) + &
     1149                   sums_l(nzt,117,tn) = sums_l(nzt,117,tn) + &
    11471150                                        surf_def_h(2)%ssws(m) * rmask(j,i,sr) ! w"s"
    11481151                ENDIF
     
    12471250                IF ( passive_scalar .AND. ( .NOT. ws_scheme_sca                &
    12481251                     .OR. sr /= 0 ) )  THEN
    1249                    pts = 0.5_wp * ( s(k,j,i)   - hom(k,1,117,sr) +             &
    1250                                     s(k+1,j,i) - hom(k+1,1,117,sr) )
    1251                    sums_l(k,116,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) *      &
    12521255                                                       rmask(j,i,sr) * flag
    12531256                ENDIF
     
    12771280             IF ( i >= nxl  .AND.  i <= nxr  .AND.                             &
    12781281                  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)
    12871288             ENDIF
    12881289          ENDDO
     
    13791380                   ENDIF
    13801381                   IF ( passive_scalar )  THEN
    1381                       pts = 0.5_wp * ( s(k,j,i)   - hom(k,1,117,sr) +           &
    1382                                       s(k+1,j,i) - hom(k+1,1,117,sr) )
    1383                       sums_l(k,116,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) *    &
    13841385                                        rmask(j,i,sr) * flag
    13851386                   ENDIF
     
    15811582                   flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    15821583
    1583                    sums_l(k,102,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)    &
    15841585                                       * rmask(j,i,sr) * flag
    1585                    sums_l(k,103,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)   &
    15861587                                       * rmask(j,i,sr) * flag
    1587                    sums_l(k,104,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)    &
    15881589                                       * rmask(j,i,sr) * flag
    1589                    sums_l(k,105,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)   &
    15901591                                       * rmask(j,i,sr) * flag
    1591                    sums_l(k,106,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) &
    15921593                                       * rmask(j,i,sr) * flag
    1593                    sums_l(k,107,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)    &
    15941595                                       * rmask(j,i,sr) * flag
    1595                    sums_l(k,108,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) &
    15961597                                       * rmask(j,i,sr) * flag
    1597                    sums_l(k,109,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)    &
    15981599                                       * rmask(j,i,sr) * flag
    15991600                ENDDO
     
    16551656          sums(k,55:63)         = sums(k,55:63)         / ngp_2dh(sr)
    16561657          sums(k,81:88)         = sums(k,81:88)         / ngp_2dh(sr)
    1657           sums(k,89:114)        = sums(k,89:114)        / ngp_2dh(sr)
    1658           sums(k,116)           = sums(k,116)           / ngp_2dh(sr)
    1659           sums(k,119)           = 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)
    16601661          IF ( ngp_2dh_s_inner(k,sr) /= 0 )  THEN
    16611662             sums(k,8:11)          = sums(k,8:11)          / ngp_2dh_s_inner(k,sr)
     
    16661667             sums(k,64)            = sums(k,64)            / ngp_2dh_s_inner(k,sr)
    16671668             sums(k,70:80)         = sums(k,70:80)         / ngp_2dh_s_inner(k,sr)
    1668              sums(k,118)           = sums(k,118)           / ngp_2dh_s_inner(k,sr)
    1669              sums(k,120: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)
    16701671          ENDIF
    16711672       ENDDO
     
    17881789          hom(:,1,91,sr) = sums(:,91)              ! m_soil
    17891790                                                   ! 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
    17981797
    17991798       ENDIF
    18001799
    18011800       IF ( radiation )  THEN
    1802           hom(:,1,101,sr) = sums(:,101)            ! rad_net
    1803           hom(:,1,102,sr) = sums(:,102)            ! rad_lw_in
    1804           hom(:,1,103,sr) = sums(:,103)            ! rad_lw_out
    1805           hom(:,1,104,sr) = sums(:,104)            ! rad_sw_in
    1806           hom(:,1,105,sr) = sums(:,105)            ! rad_sw_out
     1801          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
    18071806
    18081807          IF ( radiation_scheme == 'rrtmg' )  THEN
    1809              hom(:,1,106,sr) = sums(:,106)            ! rad_lw_cs_hr
    1810              hom(:,1,107,sr) = sums(:,107)            ! rad_lw_hr
    1811              hom(:,1,108,sr) = sums(:,108)            ! rad_sw_cs_hr
    1812              hom(:,1,109,sr) = sums(:,109)            ! rad_sw_hr
    1813 
    1814              hom(:,1,110,sr) = sums(:,110)            ! rrtm_aldif
    1815              hom(:,1,111,sr) = sums(:,111)            ! rrtm_aldir
    1816              hom(:,1,112,sr) = sums(:,112)            ! rrtm_asdif
    1817              hom(:,1,113,sr) = sums(:,113)            ! rrtm_asdir
     1808             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
    18181817          ENDIF
    18191818       ENDIF
    18201819
    1821        hom(:,1,114,sr) = sums(:,114)            !: L
     1820       hom(:,1,112,sr) = sums(:,112)            !: L
    18221821
    18231822       IF ( passive_scalar )  THEN
    1824           hom(:,1,119,sr) = sums(:,119)     ! w"s"
    1825           hom(:,1,116,sr) = sums(:,116)     ! w*s*
    1826           hom(:,1,120,sr) = sums(:,119) + sums(:,116)    ! ws
    1827           hom(:,1,118,sr) = sums(:,118)     ! s*2
    1828        ENDIF
    1829 
    1830        hom(:,1,121,sr) = rho_air       ! rho_air in Kg/m^3
    1831        hom(:,1,122,sr) = rho_air_zw    ! rho_air_zw in Kg/m^3
     1823          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
    18321831
    18331832       hom(:,1,pr_palm,sr) =   sums(:,pr_palm)
     
    19721971
    19731972       IF ( .NOT. neutral )  THEN
    1974           ts_value(22,sr) = hom(nzb,1,114,sr)          ! L
     1973          ts_value(22,sr) = hom(nzb,1,112,sr)          ! L
    19751974       ELSE
    19761975          ts_value(22,sr) = 1.0E10_wp
     
    19801979
    19811980       IF ( passive_scalar )  THEN
    1982           ts_value(24,sr) = hom(nzb+13,1,119,sr)       ! w"s" ( to do ! )
     1981          ts_value(24,sr) = hom(nzb+13,1,117,sr)       ! w"s" ( to do ! )
    19831982          ts_value(25,sr) = hom(nzb+13,1,pr_palm,sr)   ! s*
    19841983       ENDIF
     
    19871986!--    Collect land surface model timeseries
    19881987       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
    19971994       ENDIF
    19981995!
    19991996!--    Collect radiation model timeseries
    20001997       IF ( radiation )  THEN
    2001           ts_value(dots_rad,sr)   = hom(nzb,1,101,sr)          ! rad_net
    2002           ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
    2003           ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
    2004           ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
    2005           ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
     1998          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
    20062003
    20072004          IF ( radiation_scheme == 'rrtmg' )  THEN
    2008              ts_value(dots_rad+5,sr) = hom(nzb,1,110,sr)          ! rrtm_aldif
    2009              ts_value(dots_rad+6,sr) = hom(nzb,1,111,sr)          ! rrtm_aldir
    2010              ts_value(dots_rad+7,sr) = hom(nzb,1,112,sr)          ! rrtm_asdif
    2011              ts_value(dots_rad+8,sr) = hom(nzb,1,113,sr)          ! rrtm_asdir
     2005             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
    20122009          ENDIF
    20132010
Note: See TracChangeset for help on using the changeset viewer.