Ignore:
Timestamp:
Mar 20, 2007 8:33:14 AM (17 years ago)
Author:
raasch
Message:

preliminary changes for radiation conditions

File:
1 edited

Legend:

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

    r72 r73  
    88! and pt_reference are checked,
    99! output of precipitation amount/rate and roughnes length + check
    10 ! possible negative humidities are avoided in initial profile
     10! possible negative humidities are avoided in initial profile,
     11! dirichlet/neumann changed to dirichlet/radiation, etc.
    1112!
    1213! Former revisions:
     
    705706!-- Check boundary conditions and set internal variables:
    706707!-- Lateral boundary conditions
    707     IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/neumann'  .AND. &
    708          bc_lr /= 'neumann/dirichlet' )  THEN
     708    IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
     709         bc_lr /= 'radiation/dirichlet' )  THEN
    709710       IF ( myid == 0 )  THEN
    710711          PRINT*, '+++ check_parameters:'
     
    713714       CALL local_stop
    714715    ENDIF
    715     IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/neumann'  .AND. &
    716          bc_ns /= 'neumann/dirichlet' )  THEN
     716    IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
     717         bc_ns /= 'radiation/dirichlet' )  THEN
    717718       IF ( myid == 0 )  THEN
    718719          PRINT*, '+++ check_parameters:'
     
    22842285    ENDIF
    22852286
    2286     IF ( outflow_l )  THEN
     2287    IF ( bc_lr == 'radiation/dirichlet' )  THEN
    22872288       dist_nxr    = nx - inflow_disturbance_begin
    22882289       dist_nxl(1) = nx - inflow_disturbance_end
    2289     ELSEIF ( outflow_r )  THEN
     2290    ELSEIF ( bc_lr == 'dirichlet/radiation' )  THEN
    22902291       dist_nxl    = inflow_disturbance_begin
    22912292       dist_nxr(1) = inflow_disturbance_end
    2292     ELSEIF ( outflow_s )  THEN
     2293    ENDIF
     2294    IF ( bc_ns == 'dirichlet/radiation' )  THEN
    22932295       dist_nyn    = ny - inflow_disturbance_begin
    22942296       dist_nys(1) = ny - inflow_disturbance_end
    2295     ELSEIF ( outflow_n )  THEN
     2297    ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
    22962298       dist_nys    = inflow_disturbance_begin
    22972299       dist_nyn(1) = inflow_disturbance_end
Note: See TracChangeset for help on using the changeset viewer.