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

preliminary changes for radiation conditions

File:
1 edited

Legend:

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

    r4 r73  
    44! Actual revisions:
    55! -----------------
    6 !
     6! uxrp, vynp eliminated,
     7! dirichlet/neumann changed to dirichlet/radiation, etc.
    78!
    89! Former revisions:
     
    495496!
    496497!-- For non-cyclic boundaries extend array u (v) by one gridpoint
    497     IF ( bc_lr /= 'cyclic' )  uxrp = 1
    498     IF ( bc_ns /= 'cyclic' )  vynp = 1
     498!    IF ( bc_lr /= 'cyclic' )  uxrp = 1
     499!    IF ( bc_ns /= 'cyclic' )  vynp = 1
    499500
    500501!
     
    743744!-- by one gridpoint on the left/rightmost (northest/southest) processor
    744745    IF ( pleft == MPI_PROC_NULL )  THEN
    745        IF ( bc_lr == 'dirichlet/neumann' )  THEN
     746       IF ( bc_lr == 'dirichlet/radiation' )  THEN
    746747          inflow_l  = .TRUE.
    747        ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
     748       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
    748749          outflow_l = .TRUE.
    749750       ENDIF
     
    751752
    752753    IF ( pright == MPI_PROC_NULL )  THEN
    753        IF ( bc_lr == 'dirichlet/neumann' )  THEN
     754       IF ( bc_lr == 'dirichlet/radiation' )  THEN
    754755          outflow_r = .TRUE.
    755        ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
     756       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
    756757          inflow_r  = .TRUE.
    757758       ENDIF
    758        uxrp      = 1
     759!       uxrp      = 1
    759760    ENDIF
    760761
    761762    IF ( psouth == MPI_PROC_NULL )  THEN
    762        IF ( bc_ns == 'dirichlet/neumann' )  THEN
     763       IF ( bc_ns == 'dirichlet/radiation' )  THEN
    763764          outflow_s = .TRUE.
    764        ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
     765       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
    765766          inflow_s  = .TRUE.
    766767       ENDIF
     
    768769
    769770    IF ( pnorth == MPI_PROC_NULL )  THEN
    770        IF ( bc_ns == 'dirichlet/neumann' )  THEN
     771       IF ( bc_ns == 'dirichlet/radiation' )  THEN
    771772          inflow_n  = .TRUE.
    772        ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
     773       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
    773774          outflow_n = .TRUE.
    774775       ENDIF
    775        vynp      = 1
     776!       vynp      = 1
    776777    ENDIF
    777778
     
    787788    ENDIF
    788789#else
    789     IF ( bc_lr == 'dirichlet/neumann' )  THEN
     790    IF ( bc_lr == 'dirichlet/radiation' )  THEN
    790791       inflow_l  = .TRUE.
    791792       outflow_r = .TRUE.
    792        uxrp      = 1
    793     ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
     793!       uxrp      = 1
     794    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
    794795       outflow_l = .TRUE.
    795796       inflow_r  = .TRUE.
    796797    ENDIF
    797798
    798     IF ( bc_ns == 'dirichlet/neumann' )  THEN
     799    IF ( bc_ns == 'dirichlet/radiation' )  THEN
    799800       inflow_n  = .TRUE.
    800801       outflow_s = .TRUE.
    801     ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
     802    ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
    802803       outflow_n = .TRUE.
    803804       inflow_s  = .TRUE.
    804        vynp      = 1
     805!       vynp      = 1
    805806    ENDIF
    806807#endif
Note: See TracChangeset for help on using the changeset viewer.