Index: palm/trunk/SOURCE/urban_surface_mod.f90
===================================================================
--- palm/trunk/SOURCE/urban_surface_mod.f90	(revision 3769)
+++ palm/trunk/SOURCE/urban_surface_mod.f90	(revision 3802)
@@ -28,4 +28,7 @@
 ! -----------------
 ! $Id$
+! unused subroutine commented out
+! 
+! 3769 2019-02-28 10:16:49Z moh.hefny
 ! removed unused variables
 !
@@ -8629,99 +8632,99 @@
  
  
-     CONTAINS
- !------------------------------------------------------------------------------!
- ! Description:
- ! ------------
- !> Calculation of specific humidity of the skin layer (surface). It is assumend
- !> that the skin is always saturated.
- !------------------------------------------------------------------------------!
-        SUBROUTINE calc_q_surface_usm
- 
-           IMPLICIT NONE
- 
-           REAL(wp) :: resistance    !< aerodynamic and soil resistance term
- 
-           DO  m = 1, surf_usm_h%ns
- 
-              i   = surf_usm_h%i(m)            
-              j   = surf_usm_h%j(m)
-              k   = surf_usm_h%k(m)
- 
-!
-!--          Calculate water vapour pressure at saturation
-              e_s = 0.01_wp * 610.78_wp * EXP( 17.269_wp *                  &
-                                     ( t_surf_green_h_p(m) - 273.16_wp ) /  &
-                                     ( t_surf_green_h_p(m) - 35.86_wp  )    &
-                                          )
- 
-!
-!--          Calculate specific humidity at saturation
-              q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
- 
-!              surf_usm_h%r_a_green(m) = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exner(k) ) /  &
-!                    ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-10_wp )
-!                 
-! !--          make sure that the resistance does not drop to zero
-!              IF ( ABS(surf_usm_h%r_a_green(m)) < 1.0E-10_wp )  surf_usm_h%r_a_green(m) = 1.0E-10_wp
- 
-              resistance = surf_usm_h%r_a_green(m) / ( surf_usm_h%r_a_green(m) + surf_usm_h%r_s(m) + 1E-5_wp )
- 
-!
-!--          Calculate specific humidity at surface
-              IF ( bulk_cloud_model )  THEN
-                 q(k,j,i) = resistance * q_s +                   &
-                                            ( 1.0_wp - resistance ) *              &
-                                            ( q(k,j,i) - ql(k,j,i) )
-              ELSE
-                 q(k,j,i) = resistance * q_s +                   &
-                                            ( 1.0_wp - resistance ) *              &
-                                              q(k,j,i)
-              ENDIF
- 
-!
-!--          Update virtual potential temperature
-              vpt(k,j,i) = pt(k,j,i) *         &
-                         ( 1.0_wp + 0.61_wp * q(k,j,i) )
- 
-           ENDDO
-
-!
-!--       Now, treat vertical surface elements
-           DO  l = 0, 3
-              DO  m = 1, surf_usm_v(l)%ns
-!
-!--             Get indices of respective grid point
-                 i = surf_usm_v(l)%i(m)
-                 j = surf_usm_v(l)%j(m)
-                 k = surf_usm_v(l)%k(m)
- 
-!
-!--             Calculate water vapour pressure at saturation
-                 e_s = 0.01_wp * 610.78_wp * EXP( 17.269_wp *                       &
-                                        ( t_surf_green_v_p(l)%t(m) - 273.16_wp ) /  &
-                                        ( t_surf_green_v_p(l)%t(m) - 35.86_wp  )    &
-                                             )
- 
-!
-!--             Calculate specific humidity at saturation
-                 q_s = 0.622_wp * e_s / ( surface_pressure -e_s )
- 
-!
-!--             Calculate specific humidity at surface
-                 IF ( bulk_cloud_model )  THEN
-                    q(k,j,i) = ( q(k,j,i) - ql(k,j,i) )
-                 ELSE
-                    q(k,j,i) = q(k,j,i)
-                 ENDIF
-!
-!--             Update virtual potential temperature
-                 vpt(k,j,i) = pt(k,j,i) *         &
-                            ( 1.0_wp + 0.61_wp * q(k,j,i) )
- 
-              ENDDO
- 
-           ENDDO
- 
-        END SUBROUTINE calc_q_surface_usm
+!     CONTAINS
+! !------------------------------------------------------------------------------!
+! ! Description:
+! ! ------------
+! !> Calculation of specific humidity of the skin layer (surface). It is assumend
+! !> that the skin is always saturated.
+! !------------------------------------------------------------------------------!
+!        SUBROUTINE calc_q_surface_usm
+! 
+!           IMPLICIT NONE
+! 
+!           REAL(wp) :: resistance    !< aerodynamic and soil resistance term
+! 
+!           DO  m = 1, surf_usm_h%ns
+! 
+!              i   = surf_usm_h%i(m)            
+!              j   = surf_usm_h%j(m)
+!              k   = surf_usm_h%k(m)
+! 
+!!
+!!--          Calculate water vapour pressure at saturation
+!              e_s = 0.01_wp * 610.78_wp * EXP( 17.269_wp *                  &
+!                                     ( t_surf_green_h_p(m) - 273.16_wp ) /  &
+!                                     ( t_surf_green_h_p(m) - 35.86_wp  )    &
+!                                          )
+! 
+!!
+!!--          Calculate specific humidity at saturation
+!              q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
+! 
+!!              surf_usm_h%r_a_green(m) = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exner(k) ) /  &
+!!                    ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-10_wp )
+!!                 
+!! !--          make sure that the resistance does not drop to zero
+!!              IF ( ABS(surf_usm_h%r_a_green(m)) < 1.0E-10_wp )  surf_usm_h%r_a_green(m) = 1.0E-10_wp
+! 
+!              resistance = surf_usm_h%r_a_green(m) / ( surf_usm_h%r_a_green(m) + surf_usm_h%r_s(m) + 1E-5_wp )
+! 
+!!
+!!--          Calculate specific humidity at surface
+!              IF ( bulk_cloud_model )  THEN
+!                 q(k,j,i) = resistance * q_s +                   &
+!                                            ( 1.0_wp - resistance ) *              &
+!                                            ( q(k,j,i) - ql(k,j,i) )
+!              ELSE
+!                 q(k,j,i) = resistance * q_s +                   &
+!                                            ( 1.0_wp - resistance ) *              &
+!                                              q(k,j,i)
+!              ENDIF
+! 
+!!
+!!--          Update virtual potential temperature
+!              vpt(k,j,i) = pt(k,j,i) *         &
+!                         ( 1.0_wp + 0.61_wp * q(k,j,i) )
+! 
+!           ENDDO
+!
+!!
+!!--       Now, treat vertical surface elements
+!           DO  l = 0, 3
+!              DO  m = 1, surf_usm_v(l)%ns
+!!
+!!--             Get indices of respective grid point
+!                 i = surf_usm_v(l)%i(m)
+!                 j = surf_usm_v(l)%j(m)
+!                 k = surf_usm_v(l)%k(m)
+! 
+!!
+!!--             Calculate water vapour pressure at saturation
+!                 e_s = 0.01_wp * 610.78_wp * EXP( 17.269_wp *                       &
+!                                        ( t_surf_green_v_p(l)%t(m) - 273.16_wp ) /  &
+!                                        ( t_surf_green_v_p(l)%t(m) - 35.86_wp  )    &
+!                                             )
+! 
+!!
+!!--             Calculate specific humidity at saturation
+!                 q_s = 0.622_wp * e_s / ( surface_pressure -e_s )
+! 
+!!
+!!--             Calculate specific humidity at surface
+!                 IF ( bulk_cloud_model )  THEN
+!                    q(k,j,i) = ( q(k,j,i) - ql(k,j,i) )
+!                 ELSE
+!                    q(k,j,i) = q(k,j,i)
+!                 ENDIF
+!!
+!!--             Update virtual potential temperature
+!                 vpt(k,j,i) = pt(k,j,i) *         &
+!                            ( 1.0_wp + 0.61_wp * q(k,j,i) )
+! 
+!              ENDDO
+! 
+!           ENDDO
+! 
+!        END SUBROUTINE calc_q_surface_usm
         
      END SUBROUTINE usm_surface_energy_balance
