Changeset 1845 for palm/trunk/SOURCE/calc_radiation.f90

Timestamp:

Apr 8, 2016 8:29:13 AM (8 years ago)

Author:

raasch

Message:

nzb_2d replaced by nzb_..._inner, Kessler precipitation stored on surface grid point

File:

• palm/trunk/SOURCE/calc_radiation.f90 (modified) (6 diffs)

palm/trunk/SOURCE/calc_radiation.f90

@@ -93 +93 @@
 
        USE indices,                                                            &
-           ONLY:  nxl, nxr, nyn, nys, nzb, nzb_2d, nzt
+           ONLY:  nxl, nxr, nyn, nys, nzb, nzb_s_inner, nzt
 
        USE kinds
@@ -141 +141 @@
              blackbody_emission(nzb) = sigma * temperature**4
 
-             DO  k = nzb_2d(j,i)+1, nzt
+             DO  k = nzb_s_inner(j,i)+1, nzt
 
                 k_help = ( nzt+nzb+1 ) - k
@@ -167 +167 @@
              impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
 
-             DO  k = nzb_2d(j,i)+1, nzt
+             DO  k = nzb_s_inner(j,i)+1, nzt
 !
 !--             Save some computational time, but this may cause load
@@ -238 +238 @@
 
        USE indices,                                                            &
-           ONLY:  nzb, nzb_2d, nzt
+           ONLY:  nzb, nzb_s_inner, nzt
 
        USE kinds
@@ -283 +283 @@
        blackbody_emission(nzb) = sigma * temperature**4
 
-       DO  k = nzb_2d(j,i)+1, nzt
+       DO  k = nzb_s_inner(j,i)+1, nzt
           k_help = ( nzt+nzb+1 ) - k
           lwp_ground(k)   = lwp_ground(k-1) + rho_surface * ql(k,j,i) * dzw(k)
@@ -306 +306 @@
        impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
 
-       DO  k = nzb_2d(j,i)+1, nzt
+       DO  k = nzb_s_inner(j,i)+1, nzt
 !
 !--       Store some computational time,