Changeset 2334


Ignore:
Timestamp:
Aug 4, 2017 11:57:04 AM (4 years ago)
Author:
gronemeier
Message:

set c_m to 0.4 within 1D model

File:
1 edited

Legend:

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

    r2299 r2334  
    2525! -----------------
    2626! $Id$
     27! set c_m = 0.4 according to Detering and Etling (1985)
     28!
     29! 2299 2017-06-29 10:14:38Z maronga
    2730! Removed german text
    2831!
     
    273276    REAL(wp) ::  a            !<
    274277    REAL(wp) ::  b            !<
     278    REAL(wp) ::  c_m = 0.4_wp !< model constant, 0.4 according to Detering and Etling (1985)
    275279    REAL(wp) ::  dissipation  !<
    276280    REAL(wp) ::  dpt_dz       !<
     
    473477!--       Boundary conditions for the prognostic variables.
    474478!--       At the top boundary (nzt+1) u,v and e keep their initial values
    475 !--       (ug(nzt+1), vg(nzt+1), 0), at the bottom boundary the mirror
    476 !--       boundary condition applies to u and v.
    477 !--       The boundary condition for e is set further below ( (u*/cm)**2 ).
    478          ! u1d_p(nzb) = -u1d_p(nzb+1)
    479          ! v1d_p(nzb) = -v1d_p(nzb+1)
    480 
     479!--       (ug(nzt+1), vg(nzt+1), 0).
     480!--       At the bottom boundary, Dirichlet condition is used for u and v (0)
     481!--       and Neumann condition for e (e(nzb)=e(nzb+1)).
    481482          u1d_p(nzb) = 0.0_wp
    482483          v1d_p(nzb) = 0.0_wp
     
    499500!
    500501!--             Compute theta* using Rif numbers of the previous time step
    501                 IF ( rif1d(1) >= 0.0_wp )  THEN
     502                IF ( rif1d(nzb+1) >= 0.0_wp )  THEN
    502503!
    503504!--                Stable stratification
     
    628629!--             compatibility with the 3D model.
    629630                IF ( ibc_e_b == 2 )  THEN
    630                    e1d(nzb+1) = ( us1d / 0.1_wp )**2
    631 !                  e1d(nzb+1) = ( us1d / 0.4_wp )**2  !not used so far, see also
    632                                                       !prandtl_fluxes
     631                   e1d(nzb+1) = ( us1d / c_m )**2
    633632                ENDIF
    634633                e1d(nzb) = e1d(nzb+1)
     
    637636!
    638637!--                Compute q*
    639                    IF ( rif1d(1) >= 0.0_wp )  THEN
     638                   IF ( rif1d(nzb+1) >= 0.0_wp )  THEN
    640639!
    641640!--                   Stable stratification
     
    714713             ENDIF
    715714             DO  k = nzb_diff, nzt
    716 !                km1d(k) = 0.4 * SQRT( e1d(k) ) !changed: adjustment to 3D-model
    717                 km1d(k) = 0.1_wp * SQRT( e1d(k) )
     715                km1d(k) = c_m * SQRT( e1d(k) )
     716
    718717                IF ( rif1d(k) >= 0.0_wp )  THEN
    719718                   km1d(k) = km1d(k) * l1d(k)
Note: See TracChangeset for help on using the changeset viewer.