Changeset 3744 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- Feb 15, 2019 6:38:58 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r3685 r3744 21 21 ! Current revisions: 22 22 ! ------------------ 23 ! 23 ! Bugfix, missing calculation of 10cm temperature at vertical building walls, 24 ! required for indoor model 24 25 ! 25 26 ! Former revisions: … … 443 444 CALL calc_pt_near_surface ( '2m' ) 444 445 ENDIF 445 IF ( indoor_model ) THEN 446 CALL calc_pt_near_surface ( '10cm' ) 447 ENDIF446 ! 447 !-- Calculate 10cm temperature, required in indoor model 448 IF ( indoor_model ) CALL calc_pt_near_surface ( '10cm' ) 448 449 ENDIF 449 450 … … 548 549 !-- Compute respective surface fluxes for momentum and TKE 549 550 CALL calc_surface_fluxes 551 ! 552 !-- Calculate 10cm temperature, required in indoor model 553 IF ( indoor_model ) CALL calc_pt_near_surface ( '10cm' ) 550 554 ENDIF 551 555 ENDDO … … 566 570 !-- Compute respective surface fluxes for momentum and TKE 567 571 CALL calc_surface_fluxes 572 568 573 ENDIF 569 574 ENDDO … … 625 630 !-- Compute respective surface fluxes for momentum and TKE 626 631 CALL calc_surface_fluxes 632 ! 633 !-- Calculate 10cm temperature, required in indoor model 634 IF ( indoor_model ) CALL calc_pt_near_surface ( '10cm' ) 627 635 ENDIF 628 636 ENDDO … … 2020 2028 k = surf%k(m) 2021 2029 2022 surf%pt_10cm(m) = surf%pt_surface(m) + surf%ts(m) / kappa &2023 * ( log( 0.1_wp / surf%z0h(m) )&2024 - psi_h( 0.1_wp / surf%ol(m) ) 2030 surf%pt_10cm(m) = surf%pt_surface(m) + surf%ts(m) / kappa & 2031 * ( LOG( 0.1_wp / surf%z0h(m) ) & 2032 - psi_h( 0.1_wp / surf%ol(m) ) & 2025 2033 + psi_h( surf%z0h(m) / surf%ol(m) ) ) 2026 2034 2027 2035 ENDDO 2028 2036 … … 2036 2044 k = surf%k(m) 2037 2045 2038 surf%pt_2m(m) = surf%pt_surface(m) + surf%ts(m) / kappa &2039 * ( log( 2.0_wp / surf%z0h(m) )&2040 - psi_h( 2.0_wp / surf%ol(m) ) 2046 surf%pt_2m(m) = surf%pt_surface(m) + surf%ts(m) / kappa & 2047 * ( LOG( 2.0_wp / surf%z0h(m) ) & 2048 - psi_h( 2.0_wp / surf%ol(m) ) & 2041 2049 + psi_h( surf%z0h(m) / surf%ol(m) ) ) 2042 2050
Note: See TracChangeset
for help on using the changeset viewer.