Ignore:
Timestamp:
Oct 26, 2016 11:15:40 AM (5 years ago)
Author:
knoop
Message:

Anelastic approximation implemented

File:
1 edited

Legend:

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

    r2001 r2037  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Anelastic approximation implemented
    2323!
    2424! Former revisions:
     
    6565
    6666    USE arrays_3d,                                                             &
    67         ONLY:  shf
     67        ONLY:  shf, heatflux_input_conversion
    6868       
    6969    USE control_parameters,                                                    &
     
    7676   
    7777    USE indices,                                                               &
    78         ONLY:  nx, nxl, nxr, ny, nyn, nys, nzb_s_inner
     78        ONLY:  nx, nxl, nxr, ny, nyn, nys, nzb, nzb_s_inner
    7979
    8080    IMPLICIT NONE
     
    9797          THEN
    9898             IF ( nzb_s_inner(j,i) == 0 )  THEN
    99                 shf(j,i) = randomnumber * surface_heatflux
     99                shf(j,i) = randomnumber * surface_heatflux                     &
     100                           * heatflux_input_conversion(nzb)
    100101             ELSE
    101102!
    102103!--             Over topography surface_heatflux is replaced by wall_heatflux(0)
    103                 shf(j,i) = randomnumber * wall_heatflux(0)
     104                shf(j,i) = randomnumber * wall_heatflux(0)                     &
     105                           * heatflux_input_conversion(nzb_s_inner(j,i))
    104106             ENDIF
    105107          ENDIF
Note: See TracChangeset for help on using the changeset viewer.