Changeset 722 for palm/trunk/SOURCE/init_pegrid.f90
- Timestamp:
- Apr 11, 2011 6:21:09 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_pegrid.f90
r710 r722 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! 6 ! Bugfix: bc_lr/ns_cyc/dirrad/raddir replaced by bc_lr/ns, because variables 7 ! are not yet set here; grid_level set to 0 7 8 ! 8 9 ! ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!! … … 193 194 ! 194 195 !-- If necessary, set horizontal boundary conditions to non-cyclic 195 IF ( .NOT. bc_lr_cyc) cyclic(1) = .FALSE.196 IF ( .NOT. bc_ns_cyc) cyclic(2) = .FALSE.196 IF ( bc_lr /= 'cyclic' ) cyclic(1) = .FALSE. 197 IF ( bc_ns /= 'cyclic' ) cyclic(2) = .FALSE. 197 198 198 199 ! … … 977 978 ENDIF 978 979 979 grid_level = maximum_grid_level 980 ! 981 !-- Default level 0 tells exchange_horiz that all ghost planes have to be 982 !-- exchanged. grid_level is adjusted in poismg, where only one ghost plane 983 !-- is required. 984 grid_level = 0 980 985 981 986 #if defined( __parallel ) … … 1067 1072 !-- horizontal boundary conditions. 1068 1073 IF ( pleft == MPI_PROC_NULL ) THEN 1069 IF ( bc_lr _dirrad) THEN1074 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1070 1075 inflow_l = .TRUE. 1071 ELSEIF ( bc_lr _raddir) THEN1076 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1072 1077 outflow_l = .TRUE. 1073 1078 ENDIF … … 1075 1080 1076 1081 IF ( pright == MPI_PROC_NULL ) THEN 1077 IF ( bc_lr _dirrad) THEN1082 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1078 1083 outflow_r = .TRUE. 1079 ELSEIF ( bc_lr _raddir) THEN1084 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1080 1085 inflow_r = .TRUE. 1081 1086 ENDIF … … 1083 1088 1084 1089 IF ( psouth == MPI_PROC_NULL ) THEN 1085 IF ( bc_ns _dirrad) THEN1090 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1086 1091 outflow_s = .TRUE. 1087 ELSEIF ( bc_ns _raddir) THEN1092 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1088 1093 inflow_s = .TRUE. 1089 1094 ENDIF … … 1091 1096 1092 1097 IF ( pnorth == MPI_PROC_NULL ) THEN 1093 IF ( bc_ns _dirrad) THEN1098 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1094 1099 inflow_n = .TRUE. 1095 ELSEIF ( bc_ns _raddir) THEN1100 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1096 1101 outflow_n = .TRUE. 1097 1102 ENDIF … … 1123 1128 1124 1129 #else 1125 IF ( bc_lr _dirrad) THEN1130 IF ( bc_lr == 'dirichlet/radiation' ) THEN 1126 1131 inflow_l = .TRUE. 1127 1132 outflow_r = .TRUE. 1128 ELSEIF ( bc_lr _raddir) THEN1133 ELSEIF ( bc_lr == 'radiation/dirichlet' ) THEN 1129 1134 outflow_l = .TRUE. 1130 1135 inflow_r = .TRUE. 1131 1136 ENDIF 1132 1137 1133 IF ( bc_ns _dirrad) THEN1138 IF ( bc_ns == 'dirichlet/radiation' ) THEN 1134 1139 inflow_n = .TRUE. 1135 1140 outflow_s = .TRUE. 1136 ELSEIF ( bc_ns _raddir) THEN1141 ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN 1137 1142 outflow_n = .TRUE. 1138 1143 inflow_s = .TRUE.
Note: See TracChangeset
for help on using the changeset viewer.