@@ -131,58 +131,58 @@ program miranda_io
131131 close (readunit)
132132
133133 errors = 0
134- ierr = count ( r1 - w1 /= ishift )
134+ ierr = count ( R1 - W1 /= ishift )
135135 if ( ierr > 0 ) then
136136 print * , ' Error1 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
137137 errors = errors + ierr
138- mloc(1 :3 ) = maxloc (abs (r1 - w1 ))
139- wvalue = w1 (mloc(1 ),mloc(2 ),mloc(3 ))+ ishift
140- rvalue = r1 (mloc(1 ),mloc(2 ),mloc(3 ))
138+ mloc(1 :3 ) = maxloc (abs (R1 - W1 ))
139+ wvalue = W1 (mloc(1 ),mloc(2 ),mloc(3 ))+ ishift
140+ rvalue = R1 (mloc(1 ),mloc(2 ),mloc(3 ))
141141 print * , ' Maximum error at ' , mloc(1 :3 ),' write ' ,wvalue,' ; read ' ,rvalue
142142 endif
143- ierr = count ( r2 - w2 /= ishift )
143+ ierr = count ( R2 - W2 /= ishift )
144144 if ( ierr > 0 ) then
145145 print * , ' Error2 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
146146 errors = errors + ierr
147- mloc(1 :3 ) = maxloc (abs (r2 - w2 ))
148- wvalue = w2 (mloc(1 ),mloc(2 ),mloc(3 ))+ ishift
149- rvalue = r2 (mloc(1 ),mloc(2 ),mloc(3 ))
147+ mloc(1 :3 ) = maxloc (abs (R2 - W2 ))
148+ wvalue = W2 (mloc(1 ),mloc(2 ),mloc(3 ))+ ishift
149+ rvalue = R2 (mloc(1 ),mloc(2 ),mloc(3 ))
150150 print * , ' Maximum error at (' , mloc(1 :3 ),' ); write ' ,wvalue,' ; read ' ,rvalue
151151 endif
152- ierr = count ( r3 (:,:,:,1 )- w3 (:,:,:,1 ) /= ishift )
152+ ierr = count ( R3 (:,:,:,1 )- W3 (:,:,:,1 ) /= ishift )
153153 if ( ierr > 0 ) then
154154 print * , ' Error3 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
155155 errors = errors + ierr
156- mloc(1 :3 ) = maxloc (abs (r3 (:,:,:,1 )- w3 (:,:,:,1 )))
157- wvalue = w3 (mloc(1 ),mloc(2 ),mloc(3 ),1 )+ ishift
158- rvalue = r3 (mloc(1 ),mloc(2 ),mloc(3 ),1 )
156+ mloc(1 :3 ) = maxloc (abs (R3 (:,:,:,1 )- W3 (:,:,:,1 )))
157+ wvalue = W3 (mloc(1 ),mloc(2 ),mloc(3 ),1 )+ ishift
158+ rvalue = R3 (mloc(1 ),mloc(2 ),mloc(3 ),1 )
159159 print * , ' Maximum error at (' , mloc(1 :3 ),' ); write ' ,wvalue,' ; read ' ,rvalue
160160 endif
161- ierr = count ( r4 (:,:,:,1 )- w4 (:,:,:,1 ) /= ishift )
161+ ierr = count ( R4 (:,:,:,1 )- W4 (:,:,:,1 ) /= ishift )
162162 if ( ierr > 0 ) then
163163 print * , ' Error4 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
164164 errors = errors + ierr
165- mloc(1 :3 ) = maxloc (abs (r4 (:,:,:,1 )- w4 (:,:,:,1 )))
166- wvalue = w4 (mloc(1 ),mloc(2 ),mloc(3 ),1 )+ ishift
167- rvalue = r4 (mloc(1 ),mloc(2 ),mloc(3 ),1 )
165+ mloc(1 :3 ) = maxloc (abs (R4 (:,:,:,1 )- W4 (:,:,:,1 )))
166+ wvalue = W4 (mloc(1 ),mloc(2 ),mloc(3 ),1 )+ ishift
167+ rvalue = R4 (mloc(1 ),mloc(2 ),mloc(3 ),1 )
168168 print * , ' Maximum error at (' , mloc(1 :3 ),' ); write ' ,wvalue,' ; read ' ,rvalue
169169 endif
170- ierr = count ( r5 - w5 /= ishift )
170+ ierr = count ( R5 - W5 /= ishift )
171171 if ( ierr > 0 ) then
172172 print * , ' Error5 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
173173 errors = errors + ierr
174- mloc = maxloc (abs (r5 - w5 ))
175- wvalue = w5 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))+ ishift
176- rvalue = r5 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))
174+ mloc = maxloc (abs (R5 - W5 ))
175+ wvalue = W5 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))+ ishift
176+ rvalue = R5 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))
177177 print * , ' Maximum error at (' , mloc,' ); write ' ,wvalue,' ; read ' ,rvalue
178178 endif
179- ierr = count ( r6 - w6 /= ishift )
179+ ierr = count ( R6 - W6 /= ishift )
180180 if ( ierr > 0 ) then
181181 print * , ' Error6 on task ' ,mynod,' at ' ,ierr,' points reading ' ,readunit - 10000
182182 errors = errors + ierr
183- mloc = maxloc (abs (r6 - w6 ))
184- wvalue = w6 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))+ ishift
185- rvalue = r6 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))
183+ mloc = maxloc (abs (R6 - W6 ))
184+ wvalue = W6 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))+ ishift
185+ rvalue = R6 (mloc(1 ),mloc(2 ),mloc(3 ),mloc(4 ))
186186 print * , ' Maximum error at (' , mloc,' ); write ' ,wvalue,' ; read ' ,rvalue
187187 endif
188188
0 commit comments