Changeset 1964 for palm/trunk/SOURCE/init_pegrid.f90
 Timestamp:
 Jul 14, 2016 3:35:18 PM (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/init_pegrid.f90
r1933 r1964 19 19 ! Current revisions: 20 20 !  21 ! 21 ! Bugfix: erroneous setting of nest_bound_l/r/s/n = .TRUE. for vertical nesting mode removed. 22 22 ! 23 23 ! Former revisions: … … 1099 1099 ! 1100 1100 ! Setting of flags for inflow/outflow/nesting conditions. 1101 IF ( nesting_mode == 'vertical' .AND. nest_domain ) THEN 1102 IF ( nxl == 0 ) THEN 1101 IF ( pleft == MPI_PROC_NULL ) THEN 1102 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1103 inflow_l = .TRUE. 1104 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1105 outflow_l = .TRUE. 1106 ELSEIF ( bc_lr == 'nested' ) THEN 1103 1107 nest_bound_l = .TRUE. 1104 1108 ENDIF 1105 ELSE1106 IF ( pleft == MPI_PROC_NULL ) THEN1107 IF ( bc_lr == 'dirichlet/radiation' ) THEN1108 inflow_l = .TRUE.1109 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN1110 outflow_l = .TRUE.1111 ELSEIF ( bc_lr == 'nested' ) THEN1112 nest_bound_l = .TRUE.1113 ENDIF1114 ENDIF1115 1109 ENDIF 1116 1110 1117 IF ( nesting_mode == 'vertical' .AND. nest_domain ) THEN 1118 IF ( nxr == nx ) THEN 1111 IF ( pright == MPI_PROC_NULL ) THEN 1112 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1113 outflow_r = .TRUE. 1114 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1115 inflow_r = .TRUE. 1116 ELSEIF ( bc_lr == 'nested' ) THEN 1119 1117 nest_bound_r = .TRUE. 1120 1118 ENDIF 1121 ELSE 1122 IF ( pright == MPI_PROC_NULL ) THEN 1123 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1124 outflow_r = .TRUE. 1125 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1126 inflow_r = .TRUE. 1127 ELSEIF ( bc_lr == 'nested' ) THEN 1128 nest_bound_r = .TRUE. 1129 ENDIF 1130 ENDIF 1131 ENDIF 1132 1133 IF ( nesting_mode == 'vertical' .AND. nest_domain ) THEN 1134 IF ( nys == 0 ) THEN 1119 ENDIF 1120 1121 IF ( psouth == MPI_PROC_NULL ) THEN 1122 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1123 outflow_s = .TRUE. 1124 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1125 inflow_s = .TRUE. 1126 ELSEIF ( bc_ns == 'nested' ) THEN 1135 1127 nest_bound_s = .TRUE. 1136 1128 ENDIF 1137 ELSE 1138 IF ( psouth == MPI_PROC_NULL ) THEN 1139 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1140 outflow_s = .TRUE. 1141 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1142 inflow_s = .TRUE. 1143 ELSEIF ( bc_ns == 'nested' ) THEN 1144 nest_bound_s = .TRUE. 1145 ENDIF 1146 ENDIF 1147 ENDIF 1148 1149 IF ( nesting_mode == 'vertical' .AND. nest_domain ) THEN 1150 IF ( nyn == ny ) THEN 1129 ENDIF 1130 1131 IF ( pnorth == MPI_PROC_NULL ) THEN 1132 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1133 inflow_n = .TRUE. 1134 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1135 outflow_n = .TRUE. 1136 ELSEIF ( bc_ns == 'nested' ) THEN 1151 1137 nest_bound_n = .TRUE. 1152 ENDIF1153 ELSE1154 IF ( pnorth == MPI_PROC_NULL ) THEN1155 IF ( bc_ns == 'dirichlet/radiation' ) THEN1156 inflow_n = .TRUE.1157 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN1158 outflow_n = .TRUE.1159 ELSEIF ( bc_ns == 'nested' ) THEN1160 nest_bound_n = .TRUE.1161 ENDIF1162 1138 ENDIF 1163 1139 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.