Changeset 2177 in Sophya for trunk/SophyaPI/ProgPI/skymapmodule.pic
- Timestamp:
- Aug 11, 2002, 4:53:54 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaPI/ProgPI/skymapmodule.pic
r2175 r2177 1 #---------------------- 1 #---------------------- Create Maps and Save 2 2 shell rm -f h4*.fits t4*.ppf h8*.fits t8*.ppf clh*.ppf clt*.ppf almh*.ppf almt*.ppf 3 3 delobjs * … … 39 39 disp h4 zoom/4 40 40 disp t4 zoom/4 41 zone 42 43 #---------------------- 44 delobjs * 41 42 #---------------------- test Map2Double 43 delobjs * 44 zone 45 45 openfits h4.fits 46 46 … … 50 50 51 51 delobjs * 52 zone 52 53 openppf t4.ppf 53 54 map2double t4 … … 55 56 saveppf t8 t8.ppf 56 57 57 #---------------------- 58 delobjs * 58 #---------------------- test Map2Float 59 delobjs * 60 zone 61 openfits h8.fits 62 openppf t8.ppf 63 64 map2float h8 65 map2float t8 66 typemap h8 67 typemap t8 68 69 #---------------------- test TypeMap 70 delobjs * 71 zone 59 72 openfits h4.fits 60 73 openfits h8.fits … … 66 79 typemap t8 67 80 68 #---------------------- 69 delobjs * 81 #---------------------- test Map2Map 82 delobjs * 83 zone 70 84 openfits h8.fits 71 85 openppf t8.ppf … … 73 87 map2map t8 h 74 88 75 #---------------------- 76 delobjs * 89 #---------------------- test MapMult 90 delobjs * 91 zone 77 92 openfits h8.fits 78 93 openppf t8.ppf … … 80 95 mapmult t8 1000 81 96 82 #---------------------- 83 delobjs * 97 #---------------------- test MapProj 98 delobjs * 99 zone 84 100 openfits h8.fits 85 101 maproj h8 h8p … … 92 108 93 109 delobjs * 110 zone 94 111 openppf t8.ppf 95 112 maproj t8 t8p … … 101 118 saveppf t8p100 t8p100.ppf 102 119 103 #---------------------- 104 delobjs * 120 #---------------------- test Map2Cl 121 delobjs * 122 zone 105 123 openfits h8p64.fits 106 124 map2cl h8p64 clh64 … … 112 130 113 131 delobjs * 132 zone 114 133 openppf t8p100.ppf 115 134 map2cl t8p100 clt100 191 … … 117 136 n/plot clt100.n*(n+1)*val%n ! ! crossmarker5 118 137 119 #---------------------- 120 delobjs * 138 #---------------------- test Map2Alm 139 delobjs * 140 zone 121 141 openfits h8p64.fits 122 142 map2alm h8p64 almh64 … … 125 145 126 146 delobjs * 147 zone 127 148 openppf t8p100.ppf 128 149 map2alm t8p100 almt100 191 … … 130 151 disp almt100 131 152 132 #---------------------- 133 delobjs * 153 #---------------------- test Alm2Cl 154 delobjs * 155 zone 156 openppf clh64.ppf 157 openppf almh64.ppf 158 alm2cl almh64 clh64fralm 159 n/plot clh64.n*(n+1)*val%n ! ! crossmarker5 160 n/plot clh64fralm.n*(n+1)*val%n ! ! "red same circlemarker5" 161 c++exec clh64fralm -= clh64; 162 n/plot clh64fralm.n*(n+1)*val%n ! ! crossmarker5 163 164 delobjs * 165 zone 166 openppf clt100.ppf 167 openppf almt100.ppf 168 alm2cl almt100 clt100fralm 169 n/plot clt100.n*(n+1)*val%n ! ! crossmarker5 170 n/plot clt100fralm.n*(n+1)*val%n ! ! "red same circlemarker5" 171 c++exec clt100fralm -= clt100; 172 n/plot clt100fralm.n*(n+1)*val%n ! ! crossmarker5 173 174 #---------------------- test Cl2Map et SetTypeMap 175 delobjs * 176 zone 134 177 openppf clh64.ppf 135 178 settypemap h … … 139 182 140 183 delobjs * 184 zone 141 185 openppf clt100.ppf 142 186 settypemap t … … 145 189 disp t8p100frcl 146 190 147 #---------------------- 148 delobjs * 191 #---------------------- test Alm2Map et SetTypeMap 192 delobjs * 193 zone 149 194 openppf almh64.ppf 150 195 settypemap h … … 154 199 155 200 delobjs * 201 zone 156 202 openppf almt100.ppf 157 203 settypemap t … … 160 206 disp t8p100fralm 161 207 162 #---------------------- 163 delobjs * 208 #---------------------- test Cl2llCl 209 delobjs * 210 zone 211 openppf clh64.ppf 212 cl2llcl clh64 llclh64 213 n/plot clh64.n*(n+1)*val%n ! ! crossmarker5 214 n/plot llclh64.val%n ! ! "red same circlemarker5" 215 216 delobjs * 217 zone 218 openppf clt100.ppf 219 cl2llcl clt100 llclt100 220 n/plot clt100.n*(n+1)*val%n ! ! crossmarker5 221 n/plot llclt100.val%n ! ! "red same circlemarker5" 222 223 #---------------------- test ClMean 224 delobjs * 225 zone 226 openppf clh64.ppf 227 cl2llcl clh64 llclh64 228 n/plot llclh64.val%n ! ! circlemarker5 229 clmean llclh64 230 clmean llclh64 ! 231 clmean llclh64 100,9999 mean 232 echo $mean 233 234 delobjs * 235 zone 236 openppf clt100.ppf 237 cl2llcl clt100 llclt100 238 n/plot llclt100.val%n ! ! circlemarker5 239 clmean llclt100 240 clmean llclt100 ! 241 clmean llclt100 100,9999 mean 242 echo $mean 243 244 #---------------------- test ClMult 245 delobjs * 246 openppf clh64.ppf 247 zone 1 2 248 n/plot clh64.val%n ! ! circlemarker5 249 clmult clh64 1000. 250 n/plot clh64.val%n ! ! circlemarker5 251 252 delobjs * 253 openppf clt100.ppf 254 zone 1 2 255 n/plot clt100.val%n ! ! circlemarker5 256 clmult clt100 1000. 257 n/plot clt100.val%n ! ! circlemarker5 258 259 #---------------------- test ClOper 260 delobjs * 261 zone 262 openppf clh64.ppf 263 cp clh64 dum 264 clop dum + clh64 - clh64 * clh64 / clh64 265 n/plot clh64.val%n ! ! crossmarker5 266 n/plot dum.val%n ! ! "same red circlemarker5" 267 268 delobjs * 269 zone 270 openppf clt100.ppf 271 cp clt100 dum 272 clop dum + clt100 - clt100 * clt100 / clt100 273 n/plot clt100.val%n ! ! crossmarker5 274 n/plot dum.val%n ! ! "same red circlemarker5" 275 276 delobjs * 277 zone 278 openppf clh64.ppf 279 openppf clt100.ppf 280 clop clh64 - clt100 281 n/plot clh64.val%n ! ! crossmarker5 282 283 #---------------------- test ClRebin 284 delobjs * 285 zone 286 openppf clh64.ppf 287 cl2llcl clh64 llclh64 288 clrebin llclh64 clntu 10,0 289 290 n/plot clntu.n%l ! ! crossmarker5 291 292 n/plot llclh64.val%n ! ! "crossmarker3" 293 n/plot clntu.clmean%l ! ! "same circlemarker5 red" 294 n/plot clntu.cllin%l ! ! "same boxmarker5 blue" 295 n/plot clntu.clpar%l ! ! "same trianglemarker5 violet" 296 297 n/plot llclh64.val%n ! ! "crossmarker3" 298 nt2d clntu l clmean 0 sclmean 1 " " "same circlemarker5 red thinline" 299 300 n/plot llclh64.val%n ! ! "crossmarker3" 301 nt2d clntu l cllin 0 scllin 1 " " "same boxmarker5 blue thinline" 302 303 n/plot llclh64.val%n ! ! "crossmarker3" 304 nt2d clntu l clpar 0 sclpar 1 " " "same trianglemarker5 violet thinline" 305 306 #---------------------- test CrMaskMap et SetTypeMap 307 delobjs * 308 zone 164 309 settypemap h 165 310 crmapmask h8m 256 -20,20 100,130 0,1 … … 168 313 169 314 delobjs * 315 zone 170 316 settypemap t 171 317 crmapmask t8m 400 -20,20 100,130 0,1 … … 173 319 disp t8m zoom/4 174 320 175 #---------------------- 176 delobjs * 321 #---------------------- test CrMaskFrMap 322 delobjs * 323 zone 177 324 openfits h8.fits 178 325 crmaskfrmap h8fm 256 h8 -1.e-30,1.e-30 0,1 … … 181 328 182 329 delobjs * 330 zone 183 331 openppf t8.ppf 184 332 crmaskfrmap t8fm 400 t8 -1.e-30,1.e-30 0,1 … … 186 334 disp t8fm zoom/4 187 335 188 #---------------------- 189 delobjs * 336 #---------------------- test MaskMap 337 delobjs * 338 zone 190 339 openfits h8.fits 191 340 openfits h8m.fits … … 194 343 195 344 delobjs * 345 zone 196 346 openppf t8.ppf 197 347 openppf t8m.ppf … … 200 350 201 351 delobjs * 352 zone 202 353 openfits h8.fits 203 354 openppf t8m.ppf … … 206 357 207 358 delobjs * 359 zone 208 360 openppf t8.ppf 209 361 openfits h8m.fits … … 211 363 disp t8 zoom/4 212 364 213 #---------------------- 214 delobjs * 365 #---------------------- test MapCover 366 delobjs * 367 zone 215 368 openfits h8m.fits 216 369 mapcover h8m 0.9,1 couvh8m … … 218 371 219 372 delobjs * 373 zone 220 374 openppf t8m.ppf 221 375 mapcover t8m 0.9,1 couvt8m 222 376 echo $couvt8m 223 377 224 #---------------------- 225 delobjs * 378 #---------------------- test Map2Local 379 delobjs * 380 zone 226 381 openfits h8.fits 227 382 map2local h8 h8loc 200,300 20,30 0,90 ! … … 231 386 232 387 delobjs * 388 zone 233 389 openppf t8.ppf 234 390 map2local t8 t8loc 200,300 20,30 0,90 ! … … 237 393 disp t8loc 238 394 239 #---------------------- 240 delobjs * 241 openfits h8.fits 242 rename h8 h8save 243 openfits h8.fits 244 mapop h8 + h8save - h8save * h8save / h8save 395 #---------------------- test MapOper 396 delobjs * 397 zone 398 openfits h8.fits 399 cp h8 h8save 400 mapop h8save + h8 - h8 * h8 / h8 401 disp h8save zoom/4 402 mapop h8save - h8 403 disp h8save zoom/4 404 405 delobjs * 406 zone 407 openppf t8.ppf 408 cp t8 t8save 409 mapop t8save + t8 - t8 * t8 / t8 410 disp t8save zoom/4 411 mapop t8save - t8 412 disp t8save zoom/4 413 414 delobjs * 415 zone 416 openfits h8.fits 417 openppf t8.ppf 418 mapop h8 - t8 245 419 disp h8 zoom/4 246 420 247 openppf t8.ppf 248 rename t8 t8save 249 openppf t8.ppf 250 mapop t8 + t8save - t8save * t8save / t8save 251 disp t8 zoom/4 252 253 mapop h8 - t8save 254 disp h8 zoom/4 255 mapop t8 - h8save 256 disp t8 zoom/4 257 258 #---------------------- 259 delobjs * 421 #---------------------- test MapStat 422 delobjs * 423 zone 260 424 openfits h8.fits 261 425 openfits h8fm.fits … … 271 435 mapstat t8 t8fm mean sig 272 436 echo mean=$mean sig=$sig 437 438 #---------------------- test VarOper 439 varop 1000 + 100 result 440 echo result is $result 441 varop 1000 - 100 result 442 echo result is $result 443 varop 1000 * 100 result 444 echo result is $result 445 varop 1000 / 100 result 446 echo result is $result 447 varop 1000 / 0 result 448 echo result is $result 449 450 varop 1000 sqrt ! result 451 echo result is $result 452 varop -1000 sqrt ! result 453 echo result is $result 454 varop 10 exp ! result 455 echo result is $result 456 varop 1000 log ! result 457 echo result is $result 458 varop -1000 log ! result 459 echo result is $result 460 varop 1000 log10 ! result 461 echo result is $result 462 varop -1000 log10 ! result 463 echo result is $result 464 465 varop 1000 pow 2 result 466 echo result is $result 467 varop 1000 pow 1000 result 468 echo result is $result
Note:
See TracChangeset
for help on using the changeset viewer.