Ignore:
Timestamp:
Jan 16, 2018 2:16:49 PM (3 years ago)
Author:
suehring
Message:

Tile approach for spectral albedo implemented.

File:
1 edited

Legend:

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

    r2718 r2753  
    2525! -----------------
    2626! $Id$
     27! Tile approach for spectral albedo implemented.
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    11651168                      m = surf_def_h(0)%start_index(j,i)
    11661169                      sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
    1167                                            surf_def_h(0)%rrtm_aldif(m) * rmask(j,i,sr)
     1170                                   surf_def_h(0)%rrtm_aldif(0,m) * rmask(j,i,sr)
    11681171                      sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
    1169                                            surf_def_h(0)%rrtm_aldir(m) * rmask(j,i,sr)
     1172                                   surf_def_h(0)%rrtm_aldir(0,m) * rmask(j,i,sr)
    11701173                      sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
    1171                                            surf_def_h(0)%rrtm_asdif(m) * rmask(j,i,sr)
     1174                                   surf_def_h(0)%rrtm_asdif(0,m) * rmask(j,i,sr)
    11721175                      sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
    1173                                            surf_def_h(0)%rrtm_asdir(m) * rmask(j,i,sr)
     1176                                   surf_def_h(0)%rrtm_asdir(0,m) * rmask(j,i,sr)
    11741177                   ENDIF
    11751178                   IF ( surf_lsm_h%end_index(j,i) >=                           &
     
    11771180                      m = surf_lsm_h%start_index(j,i)
    11781181                      sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
    1179                                            surf_lsm_h%rrtm_aldif(m) * rmask(j,i,sr)
     1182                               SUM( surf_lsm_h%frac(:,m) *                     &
     1183                                    surf_lsm_h%rrtm_aldif(:,m) ) * rmask(j,i,sr)
    11801184                      sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
    1181                                            surf_lsm_h%rrtm_aldir(m) * rmask(j,i,sr)
     1185                               SUM( surf_lsm_h%frac(:,m) *                     &
     1186                                    surf_lsm_h%rrtm_aldir(:,m) ) * rmask(j,i,sr)
    11821187                      sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
    1183                                            surf_lsm_h%rrtm_asdif(m) * rmask(j,i,sr)
     1188                               SUM( surf_lsm_h%frac(:,m) *                     &
     1189                                    surf_lsm_h%rrtm_asdif(:,m) ) * rmask(j,i,sr)
    11841190                      sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
    1185                                            surf_lsm_h%rrtm_asdir(m) * rmask(j,i,sr)
     1191                               SUM( surf_lsm_h%frac(:,m) *                     &
     1192                                    surf_lsm_h%rrtm_asdir(:,m) ) * rmask(j,i,sr)
    11861193                   ENDIF
    11871194                   IF ( surf_usm_h%end_index(j,i) >=                           &
     
    11891196                      m = surf_usm_h%start_index(j,i)
    11901197                      sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
    1191                                            surf_usm_h%rrtm_aldif(m) * rmask(j,i,sr)
     1198                               SUM( surf_usm_h%frac(:,m) *                     &
     1199                                    surf_usm_h%rrtm_aldif(:,m) ) * rmask(j,i,sr)
    11921200                      sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
    1193                                            surf_usm_h%rrtm_aldir(m) * rmask(j,i,sr)
     1201                               SUM( surf_usm_h%frac(:,m) *                     &
     1202                                    surf_usm_h%rrtm_aldir(:,m) ) * rmask(j,i,sr)
    11941203                      sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
    1195                                            surf_usm_h%rrtm_asdif(m) * rmask(j,i,sr)
     1204                               SUM( surf_usm_h%frac(:,m) *                     &
     1205                                    surf_usm_h%rrtm_asdif(:,m) ) * rmask(j,i,sr)
    11961206                      sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
    1197                                            surf_usm_h%rrtm_asdir(m) * rmask(j,i,sr)
     1207                               SUM( surf_usm_h%frac(:,m) *                     &
     1208                                    surf_usm_h%rrtm_asdir(:,m) ) * rmask(j,i,sr)
    11981209                   ENDIF
    11991210
Note: See TracChangeset for help on using the changeset viewer.