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

Timestamp:

Dec 16, 2019 10:43:49 AM (5 years ago)

Author:

motisi

Message:

plant_canopy_model: unification of variable names: pc_ variabels now pcm_ variables
removal of bowenratio output
canopy mode 'block' now 'homogeneous'
'read_from_file_3d' now 'read_from_file'
removal of confusing comment lines
replacement of k_wall by topo_top_ind
removal of else-statement in tendency-calculation

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Renamed pc_heating_rate, pc_transpiration_rate, pc_transpiration_rate to
+! pcm_heating_rate, pcm_latent_rate, pcm_transpiration_rate
+! 
+! 4340 2019-12-16 08:17:03Z Giersch
 ! Albedo indices for building_surface_pars are now declared as parameters to 
 ! prevent an error if the gfortran compiler with -Werror=unused-value is used
@@ -302 +306 @@
 
     USE plant_canopy_model_mod,                                                &
-        ONLY:  lad_s, pc_heating_rate, pc_transpiration_rate, pc_latent_rate,  &
-               plant_canopy_transpiration, pcm_calc_transpiration_rate
+        ONLY:  lad_s,                                                          &
+               pcm_heating_rate,                                               &
+               pcm_transpiration_rate,                                         &
+               pcm_latent_rate,                                                &
+               plant_canopy_transpiration,                                     &
+               pcm_calc_transpiration_rate
 
     USE pegrid
@@ -6150 +6158 @@
 !--  push heat flux absorbed by plant canopy to respective 3D arrays
      IF ( npcbl > 0 )  THEN
-         pc_heating_rate(:,:,:) = 0.0_wp
+         pcm_heating_rate(:,:,:) = 0.0_wp
          DO ipcgb = 1, npcbl
              j = pcbl(iy, ipcgb)
@@ -6158 +6166 @@
 !--          Following expression equals former kk = k - nzb_s_inner(j,i)
              kk = k - topo_top_ind(j,i,0)  !- lad arrays are defined flat
-             pc_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
+             pcm_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
                  * pchf_prep(k) * pt(k, j, i) !-- = dT/dt
          ENDDO
@@ -6164 +6172 @@
          IF ( humidity .AND. plant_canopy_transpiration ) THEN
 !--          Calculation of plant canopy transpiration rate and correspondidng latent heat rate
-             pc_transpiration_rate(:,:,:) = 0.0_wp
-             pc_latent_rate(:,:,:) = 0.0_wp
+             pcm_transpiration_rate(:,:,:) = 0.0_wp
+             pcm_latent_rate(:,:,:) = 0.0_wp
              DO ipcgb = 1, npcbl
                  i = pcbl(ix, ipcgb)
@@ -6172 +6180 @@
                  kk = k - topo_top_ind(j,i,0)  !- lad arrays are defined flat
                  CALL pcm_calc_transpiration_rate( i, j, k, kk, pcbinsw(ipcgb), pcbinlw(ipcgb), &
-                                                   pc_transpiration_rate(kk,j,i), pc_latent_rate(kk,j,i) )
+                                                   pcm_transpiration_rate(kk,j,i), pcm_latent_rate(kk,j,i) )
               ENDDO
          ENDIF