Changeset 13537 for NEMO/trunk/src/OCE
- Timestamp:
- 2020-09-29T12:09:41+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/BDY/bdyini.F90
r13286 r13537 1111 1111 CASE( 'N' ) 1112 1112 IF( nbdyind == -1 ) THEN ! Automatic boundary definition: if nbdysegX = -1 1113 nbdyind = jpjglo - 2 ! set boundary to whole side of model domain.1113 nbdyind = Nj0glo - 2 ! set boundary to whole side of model domain. 1114 1114 nbdybeg = 2 1115 nbdyend = jpiglo -11115 nbdyend = Ni0glo -1 1116 1116 ENDIF 1117 1117 nbdysegn = nbdysegn + 1 1118 1118 npckgn(nbdysegn) = kb_bdy ! Save bdy package number 1119 jpjnob(nbdysegn) = nbdyind 1120 jpindt(nbdysegn) = nbdybeg 1121 jpinft(nbdysegn) = nbdyend 1119 jpjnob(nbdysegn) = nbdyind + nn_hls 1120 jpindt(nbdysegn) = nbdybeg + nn_hls 1121 jpinft(nbdysegn) = nbdyend + nn_hls 1122 1122 ! 1123 1123 CASE( 'S' ) … … 1125 1125 nbdyind = 2 ! set boundary to whole side of model domain. 1126 1126 nbdybeg = 2 1127 nbdyend = jpiglo - 11127 nbdyend = Ni0glo - 1 1128 1128 ENDIF 1129 1129 nbdysegs = nbdysegs + 1 1130 1130 npckgs(nbdysegs) = kb_bdy ! Save bdy package number 1131 jpjsob(nbdysegs) = nbdyind 1132 jpisdt(nbdysegs) = nbdybeg 1133 jpisft(nbdysegs) = nbdyend 1131 jpjsob(nbdysegs) = nbdyind + nn_hls 1132 jpisdt(nbdysegs) = nbdybeg + nn_hls 1133 jpisft(nbdysegs) = nbdyend + nn_hls 1134 1134 ! 1135 1135 CASE( 'E' ) 1136 1136 IF( nbdyind == -1 ) THEN ! Automatic boundary definition: if nbdysegX = -1 1137 nbdyind = jpiglo - 2 ! set boundary to whole side of model domain.1137 nbdyind = Ni0glo - 2 ! set boundary to whole side of model domain. 1138 1138 nbdybeg = 2 1139 nbdyend = jpjglo - 11139 nbdyend = Nj0glo - 1 1140 1140 ENDIF 1141 1141 nbdysege = nbdysege + 1 1142 1142 npckge(nbdysege) = kb_bdy ! Save bdy package number 1143 jpieob(nbdysege) = nbdyind 1144 jpjedt(nbdysege) = nbdybeg 1145 jpjeft(nbdysege) = nbdyend 1143 jpieob(nbdysege) = nbdyind + nn_hls 1144 jpjedt(nbdysege) = nbdybeg + nn_hls 1145 jpjeft(nbdysege) = nbdyend + nn_hls 1146 1146 ! 1147 1147 CASE( 'W' ) … … 1149 1149 nbdyind = 2 ! set boundary to whole side of model domain. 1150 1150 nbdybeg = 2 1151 nbdyend = jpjglo - 11151 nbdyend = Nj0glo - 1 1152 1152 ENDIF 1153 1153 nbdysegw = nbdysegw + 1 1154 1154 npckgw(nbdysegw) = kb_bdy ! Save bdy package number 1155 jpiwob(nbdysegw) = nbdyind 1156 jpjwdt(nbdysegw) = nbdybeg 1157 jpjwft(nbdysegw) = nbdyend 1155 jpiwob(nbdysegw) = nbdyind + nn_hls 1156 jpjwdt(nbdysegw) = nbdybeg + nn_hls 1157 jpjwft(nbdysegw) = nbdyend + nn_hls 1158 1158 ! 1159 1159 CASE DEFAULT ; CALL ctl_stop( 'ctypebdy must be N, S, E or W' ) … … 1192 1192 IF(lwp) WRITE(numout,*) 'Number of north segments : ', nbdysegn 1193 1193 IF(lwp) WRITE(numout,*) 'Number of south segments : ', nbdysegs 1194 ! 1194 1195 ! 1. Check bounds 1195 1196 !---------------- … … 1229 1230 IF (jpjwft(ib).gt.jpjglo) CALL ctl_stop( 'End index out of domain' ) 1230 1231 ENDDO 1231 !1232 1232 ! 1233 1233 ! 2. Look for segment crossings
Note: See TracChangeset
for help on using the changeset viewer.