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

Timestamp:

Oct 21, 2016 3:11:58 PM (8 years ago)

Author:

knoop

Message:

Renamed variable rho to rho_ocean, rho_init to rho_ocean_init and rho_av to rho_ocean_av

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -21 +21 @@
 ! Current revisions:
 ! ------------------
-! 
+! renamed variable rho to rho_ocean
 ! 
 ! Former revisions:
@@ -419 +419 @@
     INTEGER(iwp), PARAMETER                        :: icsurf   = 5         !< Surface skin layer heat capacity (J m−2 K−1 )
     INTEGER(iwp), PARAMETER                        :: ithick   = 6         !< thickness of the surface (wall, roof, land)  ( m )
-    INTEGER(iwp), PARAMETER                        :: irhoC    = 7         !< volumetric heat capacity rho*C of the material ( J m−3 K−1 )
+    INTEGER(iwp), PARAMETER                        :: irhoC    = 7         !< volumetric heat capacity rho_ocean*C of the material ( J m−3 K−1 )
     INTEGER(iwp), PARAMETER                        :: ilambdah = 8         !< thermal conductivity λH of the wall (W m−1 K−1 )
     CHARACTER(12), DIMENSION(:), ALLOCATABLE       :: surface_type_names   !< names of wall types (used only for reports)
@@ -2886 +2886 @@
         IF ( plant_canopy )  THEN
             pchf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp &
-                        / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho * c_p * Vbox * T)
+                        / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho_ocean * c_p * Vbox * T)
         ENDIF
 
@@ -3863 +3863 @@
             ENDIF
             
-!--         volumetric heat capacity rho*C of the wall ( J m−3 K−1 )
+!--         volumetric heat capacity rho_ocean*C of the wall ( J m−3 K−1 )
             rho_c_wall(:,l) = surface_params(irhoC, ip)
             
@@ -3930 +3930 @@
             pt1  = pt(k,j,i)
 
-!--         calculate rho * cp coefficient at surface layer
+!--         calculate rho_ocean * cp coefficient at surface layer
             rho_cp  = cp * hyp(k) / ( r_d * pt1 * exn(k) )