Changeset 3744 for palm/trunk/SOURCE/plant_canopy_model_mod.f90

Timestamp:

Feb 15, 2019 6:38:58 PM (5 years ago)

Author:

suehring

Message:

Coupling of indoor model to atmosphere; output of indoor temperatures and waste heat; enable restarts with indoor model; bugfix plant transpiration; bugfix - missing calculation of 10cm temperature

File:

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

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -22 +22 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix in transpiration, floating invalid when temperature 
+! becomes > 40 degrees
 ! 
 ! Former revisions:
@@ -236 +237 @@
 !> 2) Calculation of sinks and sources of momentum, heat and scalar concentration 
 !> due to canopy elements (subroutine pcm_tendency).
+!
+! @todo - precalculate constant terms in pcm_calc_transpiration_rate
 !------------------------------------------------------------------------------!
  MODULE plant_canopy_model_mod
@@ -422 +425 @@
      sat_press_d = sat_press * 17.27_wp * 237.3_wp / (temp + 237.3_wp)**2
 !--  Wind speed
-     wind_speed = SQRT( ((u(k,j,i) + u(k,j,i+1))/2.0_wp )**2 +           &
-                        ((v(k,j,i) + v(k,j,i+1))/2.0_wp )**2 +           &
-                        ((w(k,j,i) + w(k,j,i+1))/2.0_wp )**2 )
+     wind_speed = SQRT( ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) )**2 +            &
+                        ( 0.5_wp * ( v(k,j,i) + v(k,j,i+1) ) )**2 +            &
+                        ( 0.5_wp * ( w(k,j,i) + w(k,j,i+1) ) )**2 )
 !--  Aerodynamic conductivity (Daudet et al. (1999) eq. 14
      g_b = 0.01_wp * wind_speed + 0.0071_wp
@@ -434 +437 @@
 !--  Second function for calculation of stomatal conductivity (temperature dependency)
 !--  Stewart (1988; Agric. and Forest. Meteorol. 43) eq. 21
-     f2 = MAX(t2_min, (temp-t2_min) * (t2_max-temp)**((t2_max-16.9_wp)/(16.9_wp-t2_min)) / &
+     f2 = MAX(t2_min, (temp-t2_min) * MAX(0.0_wp,t2_max-temp)**((t2_max-16.9_wp)/(16.9_wp-t2_min)) / &
               ((16.9_wp-t2_min) * (t2_max-16.9_wp)**((t2_max-16.9_wp)/(16.9_wp-t2_min))) )
 !--  Water pressure deficit