Changeset 4555


Ignore:
Timestamp:
Jun 5, 2020 9:52:00 PM (8 months ago)
Author:
moh.hefny
Message:

Bugfix in averaging PC and MRT related quantities

File:
1 edited

Legend:

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

    r4552 r4555  
    2828! -----------------
    2929! $Id$
     30! Bugfix in averaging PC and MRT related quantities
     31!
     32! 4552 2020-06-02 20:33:29Z moh.hefny
    3033! Bugfix in IF statement in the emissivity coupling parameter for radiation-RTM
    3134!
     
    98639866
    98649867    INTEGER(iwp) ::  i !<
     9868    INTEGER(iwp) ::  imrt !< index of MRT
    98659869    INTEGER(iwp) ::  j !<
    98669870    INTEGER(iwp) ::  k !<
    98679871    INTEGER(iwp) ::  l, m !< index of current surface element
    9868 
     9872   
    98699873    INTEGER(iwp)                                       :: ids, idsint_u, idsint_l, isurf
    98709874    CHARACTER(LEN=varnamelength)                       :: var
     
    1071910723          CASE ( 'rtm_rad_pc_inlw' )
    1072010724              DO l = 1, npcbl
    10721                  pcbinlw_av(:) = pcbinlw_av(:) / REAL( average_count_3d, kind=wp )
     10725                 pcbinlw_av(l) = pcbinlw_av(l) / REAL( average_count_3d, kind=wp )
    1072210726              ENDDO
    1072310727
    1072410728          CASE ( 'rtm_rad_pc_insw' )
    1072510729              DO l = 1, npcbl
    10726                  pcbinsw_av(:) = pcbinsw_av(:) / REAL( average_count_3d, kind=wp )
     10730                 pcbinsw_av(l) = pcbinsw_av(l) / REAL( average_count_3d, kind=wp )
    1072710731              ENDDO
    1072810732
    1072910733          CASE ( 'rtm_rad_pc_inswdir' )
    1073010734              DO l = 1, npcbl
    10731                  pcbinswdir_av(:) = pcbinswdir_av(:) / REAL( average_count_3d, kind=wp )
     10735                 pcbinswdir_av(l) = pcbinswdir_av(l) / REAL( average_count_3d, kind=wp )
    1073210736              ENDDO
    1073310737
    1073410738          CASE ( 'rtm_rad_pc_inswdif' )
    1073510739              DO l = 1, npcbl
    10736                  pcbinswdif_av(:) = pcbinswdif_av(:) / REAL( average_count_3d, kind=wp )
     10740                 pcbinswdif_av(l) = pcbinswdif_av(l) / REAL( average_count_3d, kind=wp )
    1073710741              ENDDO
    1073810742
    1073910743          CASE ( 'rtm_rad_pc_inswref' )
    1074010744              DO l = 1, npcbl
    10741                  pcbinswref_av(:) = pcbinswref_av(:) / REAL( average_count_3d, kind=wp )
     10745                 pcbinswref_av(l) = pcbinswref_av(l) / REAL( average_count_3d, kind=wp )
    1074210746              ENDDO
    1074310747
    1074410748          CASE ( 'rad_mrt_lw' )
    1074510749             IF ( ALLOCATED( mrtinlw_av ) )  THEN
    10746                 mrtinlw_av(:) = mrtinlw_av(:) / REAL( average_count_3d, KIND=wp )
     10750                DO imrt = 1, nmrtbl
     10751                   mrtinlw_av(imrt) = mrtinlw_av(imrt) / REAL( average_count_3d, KIND=wp )
     10752                ENDDO
    1074710753             ENDIF
    1074810754
    1074910755          CASE ( 'rad_mrt' )
    1075010756             IF ( ALLOCATED( mrt_av ) )  THEN
    10751                 mrt_av(:) = mrt_av(:) / REAL( average_count_3d, KIND=wp )
     10757                DO imrt = 1, nmrtbl
     10758                   mrt_av(imrt) = mrt_av(imrt) / REAL( average_count_3d, KIND=wp )
     10759                ENDDO
    1075210760             ENDIF
    1075310761
Note: See TracChangeset for help on using the changeset viewer.