Ignore:
Timestamp:
Nov 2, 2015 7:43:04 AM (6 years ago)
Author:
maronga
Message:

minor bugfixes for radiation model. bugfix in subjob

File:
1 edited

Legend:

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

    r1692 r1701  
    11961196                   sums_l(k,108,tn)  = sums_l(k,108,tn)  + rad_sw_cs_hr(k,j,i) &
    11971197                                       * rmask(j,i,sr)
    1198                    sums_l(k,109,tn)  = sums_l(k,108,tn)  + rad_sw_hr(k,j,i)    &
     1198                   sums_l(k,109,tn)  = sums_l(k,109,tn)  + rad_sw_hr(k,j,i)    &
    11991199                                       * rmask(j,i,sr)
    12001200                ENDDO
     
    14011401
    14021402          IF ( radiation_scheme == 'rrtmg' )  THEN
     1403#if defined ( __rrtmg )
    14031404             hom(:,1,106,sr) = sums(:,106)            ! rad_lw_cs_hr
    14041405             hom(:,1,107,sr) = sums(:,107)            ! rad_lw_hr
     
    14061407             hom(:,1,109,sr) = sums(:,109)            ! rad_sw_hr
    14071408
    1408 #if defined ( __rrtmg )
    14091409             hom(:,1,110,sr) = sums(:,110)            ! rrtm_aldif
    14101410             hom(:,1,111,sr) = sums(:,111)            ! rrtm_aldir
     
    15881588          ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
    15891589          ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
    1590           ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
    1591           ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
     1590          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
     1591          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
    15921592
    15931593#if defined ( __rrtmg )
     
    35293529          ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
    35303530          ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
    3531           ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
    3532           ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
     3531          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
     3532          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
    35333533
    35343534#if defined ( __rrtmg )
Note: See TracChangeset for help on using the changeset viewer.