- Timestamp:
- Feb 19, 2019 7:29:39 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/biometeorology_mod.f90
r3749 r3750 3020 3020 reg_a(i) = coeff(i,1) + coeff(i,3) * ta + coeff(i,4) * & 3021 3021 sqrt_ws + coeff(i,5)*dtmrt 3022 delta_cold(i -1) = reg_a(i) + coeff(i,2) * pmvc3022 delta_cold(i) = reg_a(i) + coeff(i,2) * pmvc 3023 3023 ENDDO 3024 3024 ! … … 3027 3027 r_nenner = coeff(i,2) - coeff(i+1,2) 3028 3028 IF ( ABS ( r_nenner ) > 0.00001_wp ) THEN 3029 pmv_cross(i -1) = ( reg_a(i+1) - reg_a(i) ) / r_nenner3029 pmv_cross(i) = ( reg_a(i+1) - reg_a(i) ) / r_nenner 3030 3030 ELSE 3031 3031 nerr = 1_iwp … … 3036 3036 i_bin = 3 3037 3037 DO i = 1, 2 3038 IF ( pmva > pmv_cross(i -1) ) THEN3038 IF ( pmva > pmv_cross(i) ) THEN 3039 3039 i_bin = i 3040 3040 EXIT … … 3044 3044 !-- Adjust to operative temperature scaled according 3045 3045 !-- to classical PMV (Fanger) 3046 dpmv_cold_res = delta_cold(i_bin -1) - dpmv_adj(pmva)3046 dpmv_cold_res = delta_cold(i_bin) - dpmv_adj(pmva) 3047 3047 3048 3048 END SUBROUTINE dpmv_cold
Note: See TracChangeset
for help on using the changeset viewer.