Changeset 2317 for palm/trunk/SOURCE/radiation_model_mod.f90

Timestamp:

Jul 20, 2017 5:27:19 PM (7 years ago)

Author:

suehring

Message:

get topograpyh top index via function call

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Get topography top index via Function call 
 ! 
 ! Former revisions:
@@ -169 +169 @@
 
     USE indices,                                                               &
-        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,           &
-               wall_flags_0
+        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
 
     USE kinds
@@ -200 +199 @@
         ONLY:  rrtmg_sw
 #endif
-
-
+    USE surface_mod,                                                           &
+        ONLY:  get_topography_top_index
 
     IMPLICIT NONE
@@ -1186 +1185 @@
 !
 !--          Obtain vertical index of topography top 
-             k = MAXLOC(                                                       &
-                        MERGE( 1, 0,                                           &
-                               BTEST( wall_flags_0(:,j,i), 12 )                &
-                             ), DIM = 1                                        &
-                       ) - 1
+             k = get_topography_top_index( j, i, 's' )
 
              exn1 = (hyp(k+1) / 100000.0_wp )**0.286_wp
@@ -1243 +1238 @@
 !--          Obtain vertical index of topography top. So far it is identical to
 !--          nzb.
-             k = MAXLOC(                                                       &
-                        MERGE( 1, 0,                                           &
-                               BTEST( wall_flags_0(:,j,i), 12 )                &
-                             ), DIM = 1                                        &
-                       ) - 1
+             k = get_topography_top_index( j, i, 's' )
 
              rad_net(j,i)      = net_radiation