Changeset 4716
- Timestamp:
- Sep 30, 2020 10:06:37 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model_mod.f90
r4713 r4716 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Revert change at water surfaces (previous) 28 ! 29 ! 4713 2020-09-29 12:02:05Z pavelkrc 27 30 ! - Optimize calculation in case of calc_soil_moisture=false 28 31 ! - Improve OMP parallelization … … 2045 2048 !-- Implicit solution when the surface layer has no heat capacity, 2046 2049 !-- otherwise use RK3 scheme. 2047 IF ( surf%water_surface(m) ) THEN 2048 ! 2049 !-- water temperature is fixed (avoid numerical inaccuracies) 2050 surf_t_surface_p%var_1d(m) = surf_t_surface%var_1d(m) 2051 ELSE 2052 surf_t_surface_p%var_1d(m) = ( coef_1 * dt_3d * tsc(2) + c_surface_tmp * & 2053 surf_t_surface%var_1d(m) ) / ( c_surface_tmp + coef_2 & 2050 surf_t_surface_p%var_1d(m) = ( coef_1 * dt_3d * tsc(2) + c_surface_tmp *& 2051 surf_t_surface%var_1d(m) ) / ( c_surface_tmp + coef_2& 2054 2052 * dt_3d * tsc(2) ) 2055 ENDIF 2053 2056 2054 ! 2057 2055 !-- Add RK3 term
Note: See TracChangeset
for help on using the changeset viewer.