Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r3685 r3744  
    2121! Current revisions:
    2222! ------------------
    23 !
     23! Bugfix, missing calculation of 10cm temperature at vertical building walls,
     24! required for indoor model
    2425!
    2526! Former revisions:
     
    443444             CALL calc_pt_near_surface ( '2m' )
    444445          ENDIF
    445           IF ( indoor_model )  THEN
    446              CALL calc_pt_near_surface ( '10cm' )
    447           ENDIF
     446!
     447!--       Calculate 10cm temperature, required in indoor model
     448          IF ( indoor_model )  CALL calc_pt_near_surface ( '10cm' )
    448449       ENDIF
    449450
     
    548549!--          Compute respective surface fluxes for momentum and TKE
    549550             CALL calc_surface_fluxes
     551!
     552!--          Calculate 10cm temperature, required in indoor model
     553             IF ( indoor_model )  CALL calc_pt_near_surface ( '10cm' )
    550554          ENDIF
    551555       ENDDO
     
    566570!--          Compute respective surface fluxes for momentum and TKE
    567571             CALL calc_surface_fluxes
     572             
    568573          ENDIF
    569574       ENDDO
     
    625630!--          Compute respective surface fluxes for momentum and TKE
    626631             CALL calc_surface_fluxes
     632!
     633!--          Calculate 10cm temperature, required in indoor model             
     634             IF ( indoor_model )  CALL calc_pt_near_surface ( '10cm' )
    627635          ENDIF
    628636       ENDDO
     
    20202028                k = surf%k(m)
    20212029
    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) )               &
    20252033                                     + psi_h( surf%z0h(m) / surf%ol(m) ) )
    2026 
     2034                                   
    20272035             ENDDO
    20282036
     
    20362044                k = surf%k(m)
    20372045
    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) )            &
    20412049                                     + psi_h( surf%z0h(m) / surf%ol(m) ) )
    20422050
Note: See TracChangeset for help on using the changeset viewer.