Changeset 236 for trunk/SRC/Grid/checkperio.pro
- Timestamp:
- 03/21/07 18:18:44 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Grid/checkperio.pro
r232 r236 47 47 48 48 IF file_test(file) EQ 0 THEN BEGIN 49 print, 'file '+file+' not found'49 ras = report( 'file '+file+' not found') 50 50 RETURN 51 51 ENDIF … … 78 78 part2 = reverse(shift(part2, -1)) 79 79 IF array_equal(part1, part2) NE 1 THEN BEGIN 80 print, vinq.name+' bad y periodicity (1)'80 ras = report( vinq.name+' bad y periodicity (1)') 81 81 IF keyword_set(write) THEN $ 82 82 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1], offset = [1, jpj-1] … … 88 88 part2 = reverse(part2) 89 89 IF array_equal(part1, part2) NE 1 THEN BEGIN 90 print, vinq.name+' bad y periodicity (1)'90 ras = report( vinq.name+' bad y periodicity (1)') 91 91 IF keyword_set(write) THEN $ 92 92 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1], offset = [1, jpj-1] 93 93 ENDIF 94 94 END 95 ELSE: print, vinq.name+' nothing to check'95 ELSE:ras = report( vinq.name+' nothing to check') 96 96 ENDCASE 97 97 ; y periodicity (2) … … 102 102 part2 = reverse(part2) 103 103 IF array_equal(part1, part2) NE 1 THEN BEGIN 104 print, vinq.name+' bad y periodicity (2)'104 ras = report( vinq.name+' bad y periodicity (2)') 105 105 IF keyword_set(write) THEN $ 106 106 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1], offset = [(jpi-1)/2+2, jpj-2] … … 109 109 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 110 110 END 111 ELSE: print, vinq.name+' nothing to check'111 ELSE:ras = report( vinq.name+' nothing to check') 112 112 ENDCASE 113 113 ; x periodicity (1) … … 115 115 ncdf_varget, cdfid, vid, part2, count = [1, jpj], offset = [jpi-2, 0] 116 116 IF array_equal(part1, part2) NE 1 THEN BEGIN 117 print, vinq.name+' bad x periodicity (1)'117 ras = report( vinq.name+' bad x periodicity (1)') 118 118 IF keyword_set(write) THEN $ 119 119 ncdf_varput, cdfid, vid, part2, count = [1, jpj], offset = [0, 0] … … 123 123 ncdf_varget, cdfid, vid, part2, count = [1, jpj], offset = [1, 0] 124 124 IF array_equal(part1, part2) NE 1 THEN BEGIN 125 print, vinq.name+' bad x periodicity (2)'125 ras = report( vinq.name+' bad x periodicity (2)') 126 126 IF keyword_set(write) THEN $ 127 127 ncdf_varput, cdfid, vid, part2, count = [1, jpj], offset = [jpi-1, 0] … … 141 141 part2 = reverse(shift(part2, -1, 0, 0), 1) 142 142 IF array_equal(part1, part2) NE 1 THEN BEGIN 143 print, vinq.name+' bad y periodicity (1)'143 ras = report( vinq.name+' bad y periodicity (1)') 144 144 IF keyword_set(write) THEN $ 145 145 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk], offset = [1, jpj-1, 0] … … 152 152 part2 = reverse(part2, 1) 153 153 IF array_equal(part1, part2) NE 1 THEN BEGIN 154 print, vinq.name+' bad y periodicity (1)'154 ras = report( vinq.name+' bad y periodicity (1)') 155 155 IF keyword_set(write) THEN $ 156 156 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk], offset = [1, jpj-1, 0] 157 157 ENDIF 158 158 END 159 ELSE: print, vinq.name+' nothing to check'159 ELSE:ras = report( vinq.name+' nothing to check') 160 160 ENDCASE 161 161 ; y periodicity (2) … … 166 166 part2 = reverse(part2, 1) 167 167 IF array_equal(part1, part2) NE 1 THEN BEGIN 168 print, vinq.name+' bad y periodicity (2)'168 ras = report( vinq.name+' bad y periodicity (2)') 169 169 IF keyword_set(write) THEN $ 170 170 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1, jpk], offset = [(jpi-1)/2+2, jpj-2, 0] … … 173 173 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 174 174 END 175 ELSE: print, vinq.name+' nothing to check'175 ELSE:ras = report( vinq.name+' nothing to check') 176 176 ENDCASE 177 177 ; x periodicity (1) … … 179 179 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk], offset = [jpi-2, 0, 0] 180 180 IF array_equal(part1, part2) NE 1 THEN BEGIN 181 print, vinq.name+' bad x periodicity (1)'181 ras = report( vinq.name+' bad x periodicity (1)') 182 182 IF keyword_set(write) THEN $ 183 183 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk], offset = [0, 0, 0] … … 187 187 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk], offset = [1, 0, 0] 188 188 IF array_equal(part1, part2) NE 1 THEN BEGIN 189 print, vinq.name+' bad x periodicity (2)'189 ras = report( vinq.name+' bad x periodicity (2)') 190 190 IF keyword_set(write) THEN $ 191 191 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk], offset = [jpi-1, 0, 0] … … 206 206 part2 = reverse(shift(part2, -1, 0, 0, 0), 1) 207 207 IF array_equal(part1, part2) NE 1 THEN BEGIN 208 print, vinq.name+' bad y periodicity (1)'208 ras = report( vinq.name+' bad y periodicity (1)') 209 209 IF keyword_set(write) THEN $ 210 210 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk, jpt], offset = [1, jpj-1, 0, 0] … … 217 217 part2 = reverse(part2, 1) 218 218 IF array_equal(part1, part2) NE 1 THEN BEGIN 219 print, vinq.name+' bad y periodicity (1)'219 ras = report( vinq.name+' bad y periodicity (1)') 220 220 IF keyword_set(write) THEN $ 221 221 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk, jpt], offset = [1, jpj-1, 0, 0] 222 222 ENDIF 223 223 END 224 ELSE: print, vinq.name+' nothing to check'224 ELSE:ras = report( vinq.name+' nothing to check') 225 225 ENDCASE 226 226 ; y periodicity (2) … … 231 231 part2 = reverse(part2, 1) 232 232 IF array_equal(part1, part2) NE 1 THEN BEGIN 233 print, vinq.name+' bad y periodicity (2)'233 ras = report( vinq.name+' bad y periodicity (2)') 234 234 IF keyword_set(write) THEN $ 235 235 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1, jpk, jpt], offset = [(jpi-1)/2+2, jpj-2, 0, 0] … … 238 238 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 239 239 END 240 ELSE: print, vinq.name+' nothing to check'240 ELSE:ras = report( vinq.name+' nothing to check') 241 241 ENDCASE 242 242 ; x periodicity (1) … … 244 244 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [jpi-2, 0, 0, 0] 245 245 IF array_equal(part1, part2) NE 1 THEN BEGIN 246 print, vinq.name+' bad x periodicity (1)'246 ras = report( vinq.name+' bad x periodicity (1)') 247 247 IF keyword_set(write) THEN $ 248 248 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [0, 0, 0, 0] … … 252 252 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [1, 0, 0, 0] 253 253 IF array_equal(part1, part2) NE 1 THEN BEGIN 254 print, vinq.name+' bad x periodicity (2)'254 ras = report( vinq.name+' bad x periodicity (2)') 255 255 IF keyword_set(write) THEN $ 256 256 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [jpi-1, 0, 0, 0] 257 257 ENDIF 258 258 END 259 ELSE: print, vinq.name+' nothing to check'259 ELSE:ras = report( vinq.name+' nothing to check') 260 260 ENDCASE 261 261 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.