Changeset 3156 for palm/trunk


Ignore:
Timestamp:
Jul 19, 2018 4:30:54 PM (6 years ago)
Author:
knoop
Message:

Bugfix: replaced usage of the pt array with the surf%pt_surface array in radiation model

File:
1 edited

Legend:

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

    r3155 r3156  
    30993099!
    31003100!--             Prepare profiles of temperature and H2O volume mixing ratio
    3101                 rrtm_tlev(0,nzb+1) = pt(nzb,j,i) * ( surface_pressure          &
    3102                                                      / 1000.0_wp )**0.286_wp
     3101                DO  m = surf_lsm_h%start_index(j,i), surf_lsm_h%end_index(j,i)
     3102                   rrtm_tlev(0,nzb+1) = surf_lsm_h%pt_surface(m) * &
     3103                                    ( surface_pressure / 1000.0_wp )**0.286_wp
     3104                ENDDO
     3105                DO  m = surf_usm_h%start_index(j,i), surf_usm_h%end_index(j,i)
     3106                   rrtm_tlev(0,nzb+1) = surf_usm_h%pt_surface(m) * &
     3107                                    ( surface_pressure / 1000.0_wp )**0.286_wp
     3108                ENDDO
    31033109
    31043110
     
    32403246                               surf_lsm_h%frac(ind_wat_win,m)   *              &
    32413247                               surf_lsm_h%emissivity(ind_wat_win,m)
    3242                    rrtm_tsfc = pt(surf_lsm_h%k(m)+surf_lsm_h%koff,j,i) *       &
     3248                   rrtm_tsfc = surf_lsm_h%pt_surface(m) *                      &
    32433249                                       (surface_pressure / 1000.0_wp )**0.286_wp
    32443250                ENDDO             
     
    32503256                               surf_usm_h%frac(ind_wat_win,m)   *              &
    32513257                               surf_usm_h%emissivity(ind_wat_win,m)
    3252                    rrtm_tsfc = pt(surf_usm_h%k(m)+surf_usm_h%koff,j,i) *       &
     3258                   rrtm_tsfc = surf_usm_h%pt_surface(m) *                      &
    32533259                                       (surface_pressure / 1000.0_wp )**0.286_wp
    32543260                ENDDO
Note: See TracChangeset for help on using the changeset viewer.