Changeset 1757 for palm/trunk/SOURCE/surface_layer_fluxes.f90

Timestamp:

Feb 22, 2016 3:49:32 PM (8 years ago)

Author:

maronga

Message:

some changes in land surface model, radiation model, nudging and some minor updates

File:

• palm/trunk/SOURCE/surface_layer_fluxes.f90 (modified) (4 diffs)

palm/trunk/SOURCE/surface_layer_fluxes.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Minor fixes.
 !
 ! Former revisions:
@@ -358 +358 @@
 
 !
-!--       Invert array and switch sign, then calculate Obukhov length and bulk
-!--       Richardson number
-          z0h_min = MINVAL(z0h)
-          z0_min  = MINVAL(z0)
+!--       Save roughness lengths to temporary variables
+          z0h_min = z0h(nys,nxl)
+          z0_min  = z0(nys,nxl)
           
 !
@@ -641 +640 @@
              DO  j = nys, nyn
 
-                k   = nzb_s_inner(j,i)
-
 !
 !--             If the bulk Richardson number is outside the range of the lookup
@@ -1055 +1052 @@
              !$acc loop independent
              DO  j = nysg, nyng
-                e(nzb_s_inner(j,i)+1,j,i) = ( us(j,i) / 0.1_wp )**2
+                k = nzb_s_inner(j,i)
+                e(k+1,j,i) = ( us(j,i) / 0.1_wp )**2
 !
 !--             As a test: cm = 0.4
-!               e(nzb_s_inner(j,i)+1,j,i) = ( us(j,i) / 0.4_wp )**2
-                e(nzb_s_inner(j,i),j,i)   = e(nzb_s_inner(j,i)+1,j,i)
+!               e(k+1,j,i) = ( us(j,i) / 0.4_wp )**2
+                e(k,j,i)   = e(k+1,j,i)
              ENDDO
           ENDDO