/home/eoseret/qaas_runs_CPU_9468/171-152-3172/intel/CloverLeafFC/build/CloverLeafFC/CloverLeaf_ref/kernels/generate_chunk_kernel.f90: 85 - 161
--------------------------------------------------------------------------------

85:     !$OMP PARALLEL SHARED(x_cent,y_cent)
86:     !$OMP DO
87:     DO k=y_min-2,y_max+2
88:     !$OMP SIMD
89:       DO j=x_min-2,x_max+2
90:         energy0(j,k)=state_energy(1)
91:       ENDDO
92:     ENDDO
93:     !$OMP END DO
94:     !$OMP DO
95:     DO k=y_min-2,y_max+2
96:     !$OMP SIMD
97:       DO j=x_min-2,x_max+2
98:         density0(j,k)=state_density(1)
99:       ENDDO
100:     ENDDO
101:     !$OMP END DO
102:     !$OMP DO
103:     DO k=y_min-2,y_max+2
104:     !$OMP SIMD
105:       DO j=x_min-2,x_max+2
106:         xvel0(j,k)=state_xvel(1)
107:       ENDDO
108:     ENDDO
109:     !$OMP END DO
110:     !$OMP DO
111:     DO k=y_min-2,y_max+2
112:     !$OMP SIMD
113:       DO j=x_min-2,x_max+2
114:         yvel0(j,k)=state_yvel(1)
115:       ENDDO
116:     ENDDO
117:     !$OMP END DO
118: 
119:     DO state=2,number_of_states
120: 
121:       ! Could the velocity setting be thread unsafe?
122:       x_cent=state_xmin(state)
123:       y_cent=state_ymin(state)
124: 
125:       !$OMP DO PRIVATE(radius,jt,kt)
126:       DO k=y_min-2,y_max+2
127:       !$OMP SIMD
128:         DO j=x_min-2,x_max+2
129:           IF(state_geometry(state).EQ.g_rect ) THEN
130:             IF(vertexx(j+1).GE.state_xmin(state).AND.vertexx(j).LT.state_xmax(state)) THEN
131:               IF(vertexy(k+1).GE.state_ymin(state).AND.vertexy(k).LT.state_ymax(state)) THEN
132:                 energy0(j,k)=state_energy(state)
133:                 density0(j,k)=state_density(state)
134:                 DO kt=k,k+1
135:                   DO jt=j,j+1
136:                     xvel0(jt,kt)=state_xvel(state)
137:                     yvel0(jt,kt)=state_yvel(state)
138:                   ENDDO
139:                 ENDDO
140:               ENDIF
141:             ENDIF
142:           ELSEIF(state_geometry(state).EQ.g_circ ) THEN
143:             radius=SQRT((cellx(j)-x_cent)*(cellx(j)-x_cent)+(celly(k)-y_cent)*(celly(k)-y_cent))
144:             IF(radius.LE.state_radius(state))THEN
145:               energy0(j,k)=state_energy(state)
146:               density0(j,k)=state_density(state)
147:               DO kt=k,k+1
148:                 DO jt=j,j+1
149:                   xvel0(jt,kt)=state_xvel(state)
150:                   yvel0(jt,kt)=state_yvel(state)
151:                 ENDDO
152:               ENDDO
153:             ENDIF
154:           ELSEIF(state_geometry(state).EQ.g_point) THEN
155:             IF(vertexx(j).EQ.x_cent .AND. vertexy(k).EQ.y_cent) THEN
156:               energy0(j,k)=state_energy(state)
157:               density0(j,k)=state_density(state)
158:               DO kt=k,k+1
159:                 DO jt=j,j+1
160:                   xvel0(jt,kt)=state_xvel(state)
161:                   yvel0(jt,kt)=state_yvel(state)
