Changeset 107 for palm/trunk/SOURCE


Ignore:
Timestamp:
Aug 17, 2007 1:54:45 PM (17 years ago)
Author:
raasch
Message:

further bugfix for non-cyclic BCs

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r106 r107  
    1717Quantities w*p* and w"e can be output as vertical profiles.
    1818
    19 boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_3d_model, modules, palm, parin, pres, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
     19Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)
     20
     21boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_rankine, init_3d_model, modules, palm, parin, pres, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
    2022
    2123New:
  • palm/trunk/SOURCE/boundary_conds.f90

    r106 r107  
    285285                                       ( u(k,-1,i) - u(k,0,i) ) * ddy
    286286
    287              v_p(k,-1,i) = v(k,-1,i) - dt_3d * tsc(2) * c_v(k,i) * &
     287             v_p(k,0,i)  = v(k,0,i) - dt_3d * tsc(2) * c_v(k,i) * &
    288288                                       ( v(k,0,i) - v(k,1,i) ) * ddy
    289289
     
    433433                IF ( denom /= 0.0 )  THEN
    434434                   c_u(k,j) = -c_max * ( u(k,j,1) - u_m_l(k,j,1) ) / denom
    435                    IF ( c_u(k,j) > 0.0 )  THEN
     435                   IF ( c_u(k,j) < 0.0 )  THEN
    436436                      c_u(k,j) = 0.0
    437                    ELSEIF ( c_u(k,j) < -c_max )  THEN
    438                       c_u(k,j) = -c_max
    439                    ENDIF
    440                 ELSE
    441                    c_u(k,j) = -c_max
     437                   ELSEIF ( c_u(k,j) > c_max )  THEN
     438                      c_u(k,j) = c_max
     439                   ENDIF
     440                ELSE
     441                   c_u(k,j) = c_max
    442442                ENDIF
    443443
  • palm/trunk/SOURCE/init_rankine.f90

    r77 r107  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Initial profiles are reset to constant profiles
    77!
    88! Former revisions:
     
    4848    kc1 = nzb
    4949    kc2 = nzt+1
     50
     51!
     52!-- Reset initial profiles to constant profiles
     53    IF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0 )  THEN
     54       DO  i = nxl-1, nxr+1
     55          DO  j = nys-1, nyn+1
     56             pt(:,j,i) = pt_init
     57             u(:,j,i)  = u_init
     58             v(:,j,i)  = v_init
     59          ENDDO
     60       ENDDO
     61    ENDIF
    5062
    5163!
Note: See TracChangeset for help on using the changeset viewer.