Changeset 1513
- Timestamp:
- Dec 19, 2014 9:14:10 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model.f90
r1501 r1513 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! Bugfix: REAL constants provided with KIND-attribute in call of 23 ! intrinsic function MAX and MIN 23 24 ! 24 25 ! Former revisions: … … 808 809 809 810 f2 = ( m_soil(0,j,i) - m_min ) / ( m_fc(j,i) - m_min ) 810 f2 = MAX(f2,1.0E-20 )811 f2 = MAX(f2,1.0E-20_wp) 811 812 812 813 r_soil(j,i) = r_soil_min(j,i) / f2 … … 815 816 !-- Calculate fraction of liquid water reservoir 816 817 m_liq_max = m_max_depth * LAI(j,i) 817 c_liq(j,i) = MIN(1.0 , m_liq(j,i)/m_liq_max)818 c_liq(j,i) = MIN(1.0_wp, m_liq(j,i)/m_liq_max) 818 819 819 820 q_s = 0.622_wp * e_s / surface_pressure … … 1068 1069 !-- Calculate soil heat conductivity at the center of the soil 1069 1070 !-- layers 1070 Ke = 1.0 + LOG10(MAX(0.1 ,m_soil(k,j,i) / m_sat(j,i)))1071 Ke = 1.0 + LOG10(MAX(0.1_wp,m_soil(k,j,i) / m_sat(j,i))) 1071 1072 lambda_temp(k) = Ke * (lambda_h_sat(j,i) + lambda_h_dry) + & 1072 1073 lambda_h_dry
Note: See TracChangeset
for help on using the changeset viewer.