Changeset 3156 for palm

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:

• palm/trunk/SOURCE/radiation_model_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -3099 +3099 @@
 !
 !--             Prepare profiles of temperature and H2O volume mixing ratio
-                rrtm_tlev(0,nzb+1) = pt(nzb,j,i) * ( surface_pressure          &
-                                                     / 1000.0_wp )**0.286_wp
+                DO  m = surf_lsm_h%start_index(j,i), surf_lsm_h%end_index(j,i)
+                   rrtm_tlev(0,nzb+1) = surf_lsm_h%pt_surface(m) * &
+                                    ( surface_pressure / 1000.0_wp )**0.286_wp
+                ENDDO
+                DO  m = surf_usm_h%start_index(j,i), surf_usm_h%end_index(j,i)
+                   rrtm_tlev(0,nzb+1) = surf_usm_h%pt_surface(m) * &
+                                    ( surface_pressure / 1000.0_wp )**0.286_wp
+                ENDDO
 
 
@@ -3240 +3246 @@
                                surf_lsm_h%frac(ind_wat_win,m)   *              &
                                surf_lsm_h%emissivity(ind_wat_win,m)
-                   rrtm_tsfc = pt(surf_lsm_h%k(m)+surf_lsm_h%koff,j,i) *       &
+                   rrtm_tsfc = surf_lsm_h%pt_surface(m) *                      &
                                        (surface_pressure / 1000.0_wp )**0.286_wp
                 ENDDO             
@@ -3250 +3256 @@
                                surf_usm_h%frac(ind_wat_win,m)   *              &
                                surf_usm_h%emissivity(ind_wat_win,m)
-                   rrtm_tsfc = pt(surf_usm_h%k(m)+surf_usm_h%koff,j,i) *       &
+                   rrtm_tsfc = surf_usm_h%pt_surface(m) *                      &
                                        (surface_pressure / 1000.0_wp )**0.286_wp
                 ENDDO