Changeset 73 for palm/trunk/SOURCE/init_pegrid.f90
 Timestamp:
 Mar 20, 2007 8:33:14 AM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/init_pegrid.f90
r4 r73 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! uxrp, vynp eliminated, 7 ! dirichlet/neumann changed to dirichlet/radiation, etc. 7 8 ! 8 9 ! Former revisions: … … 495 496 ! 496 497 ! For noncyclic boundaries extend array u (v) by one gridpoint 497 IF ( bc_lr /= 'cyclic' ) uxrp = 1498 IF ( bc_ns /= 'cyclic' ) vynp = 1498 ! IF ( bc_lr /= 'cyclic' ) uxrp = 1 499 ! IF ( bc_ns /= 'cyclic' ) vynp = 1 499 500 500 501 ! … … 743 744 ! by one gridpoint on the left/rightmost (northest/southest) processor 744 745 IF ( pleft == MPI_PROC_NULL ) THEN 745 IF ( bc_lr == 'dirichlet/ neumann' ) THEN746 IF ( bc_lr == 'dirichlet/radiation' ) THEN 746 747 inflow_l = .TRUE. 747 ELSEIF ( bc_lr == ' neumann/dirichlet' ) THEN748 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 748 749 outflow_l = .TRUE. 749 750 ENDIF … … 751 752 752 753 IF ( pright == MPI_PROC_NULL ) THEN 753 IF ( bc_lr == 'dirichlet/ neumann' ) THEN754 IF ( bc_lr == 'dirichlet/radiation' ) THEN 754 755 outflow_r = .TRUE. 755 ELSEIF ( bc_lr == ' neumann/dirichlet' ) THEN756 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 756 757 inflow_r = .TRUE. 757 758 ENDIF 758 uxrp = 1759 ! uxrp = 1 759 760 ENDIF 760 761 761 762 IF ( psouth == MPI_PROC_NULL ) THEN 762 IF ( bc_ns == 'dirichlet/ neumann' ) THEN763 IF ( bc_ns == 'dirichlet/radiation' ) THEN 763 764 outflow_s = .TRUE. 764 ELSEIF ( bc_ns == ' neumann/dirichlet' ) THEN765 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 765 766 inflow_s = .TRUE. 766 767 ENDIF … … 768 769 769 770 IF ( pnorth == MPI_PROC_NULL ) THEN 770 IF ( bc_ns == 'dirichlet/ neumann' ) THEN771 IF ( bc_ns == 'dirichlet/radiation' ) THEN 771 772 inflow_n = .TRUE. 772 ELSEIF ( bc_ns == ' neumann/dirichlet' ) THEN773 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 773 774 outflow_n = .TRUE. 774 775 ENDIF 775 vynp = 1776 ! vynp = 1 776 777 ENDIF 777 778 … … 787 788 ENDIF 788 789 #else 789 IF ( bc_lr == 'dirichlet/ neumann' ) THEN790 IF ( bc_lr == 'dirichlet/radiation' ) THEN 790 791 inflow_l = .TRUE. 791 792 outflow_r = .TRUE. 792 uxrp = 1793 ELSEIF ( bc_lr == ' neumann/dirichlet' ) THEN793 ! uxrp = 1 794 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 794 795 outflow_l = .TRUE. 795 796 inflow_r = .TRUE. 796 797 ENDIF 797 798 798 IF ( bc_ns == 'dirichlet/ neumann' ) THEN799 IF ( bc_ns == 'dirichlet/radiation' ) THEN 799 800 inflow_n = .TRUE. 800 801 outflow_s = .TRUE. 801 ELSEIF ( bc_ns == ' neumann/dirichlet' ) THEN802 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 802 803 outflow_n = .TRUE. 803 804 inflow_s = .TRUE. 804 vynp = 1805 ! vynp = 1 805 806 ENDIF 806 807 #endif
Note: See TracChangeset
for help on using the changeset viewer.