Changeset 2213 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Apr 24, 2017 3:10:35 PM (7 years ago)

Author:

kanani

Message:

bugfix in PCM output, minor formatting and clean-up

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -21 +21 @@
 ! Current revisions:
 ! ------------------
-! 
+! Removal of output quantities usm_lad and usm_canopy_hr
 ! 
 ! Former revisions:
@@ -2005 +2005 @@
                   var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis')  THEN
             unit = '1'
-        ELSE IF ( plant_canopy  .AND.  var(1:7) == 'usm_lad' )  THEN
-            unit = 'm2/m3'
-        ELSE IF ( plant_canopy  .AND.  var(1:13) == 'usm_canopy_hr' )  THEN
-            unit = 'K/s'
         ELSE
             unit = 'illegal'
@@ -2401 +2397 @@
                  ENDIF
               ENDDO
-
-          CASE ( 'usm_lad' )
-!--           leaf area density
-              DO i = nxl, nxr
-                 DO j = nys, nyn
-                     DO k = nzb_s_inner(j,i), nzut
-                         temp_pf(k,j,i) = lad_s(k-nzb_s_inner(j,i),j,i)
-                     ENDDO
-                 ENDDO
-              ENDDO
-              
-          CASE ( 'usm_canopy_hr' )
-!--           canopy heating rate
-              DO i = nxl, nxr
-                 DO j = nys, nyn
-                     DO k = nzb_s_inner(j,i), nzut
-                         temp_pf(k,j,i) = pc_heating_rate(k-nzb_s_inner(j,i),j,i)
-                     ENDDO
-                 ENDDO
-              ENDDO
              
           CASE DEFAULT
@@ -2471 +2447 @@
              var(1:9) == 'usm_surfz'  .OR.  var(1:7) == 'usm_svf'  .OR.                     & 
              var(1:7) == 'usm_dif'  .OR.  var(1:11) == 'usm_surfcat'  .OR.                  &
-             var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis'  .OR.            &
-             var(1:7) == 'usm_lad'  .OR.  var(1:13) == 'usm_canopy_hr' )  THEN
+             var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis' )  THEN
 
             found = .TRUE.