IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jan 31, 2007, 1:09:27 PM (19 years ago)
Author:
Paul Price
Message:

Fixing conflicts.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/psdemo/simdata/sim.pro

    r11489 r11496  
    55
    66$flatlevel = 10000
    7 $biaslevel = 200
    8 $skylevel  = 20000
     7$biaslevel = 500
     8$skyrate  = 200
     9$darkrate = 100
    910
    1011macro cleanup
    … …  
    2122  exec mkdir -p proc
    2223  mkbiases
     24  mkdarks
    2325  mkflats
    2426  mkobjects
    … …  
    2628end
    2729
    28 # create a bias detrun:
    29 # exec dettool -pretend -definebyquery -det_type bias -select_exp_type bias -workdir path://SIMTEST/mkbias.0 -dbname simtest
    30 # exec dettool -updatedetrun -det_id 2 -state stop -dbname isp
    31 
    32 # create a flat detrun:
    33 # exec dettool -pretend -definebyquery -inst ISP-Apogee -det_type dark -select_exp_type dark -workdir path://ISP/mkdark.1 -dbname isp -mode verify  -simple
     30### create a bias detrun:
     31# exec dettool -pretend -definebyquery -inst SIMPLE -filter NONE -det_type bias -select_exp_type bias -workdir path://SIMTEST/mkbias.0 -dbname simtest
     32
     33### If something goes wrong:
     34# exec dettool -updatedetrun -det_id 1 -state stop -dbname simtest
     35
     36### create a dark detrun:
     37# exec dettool -pretend -definebyquery -inst SIMPLE -filter NONE -det_type dark -select_exp_type dark -workdir path://SIMTEST/mkdark.0 -dbname simtest
     38
     39### create a flat detrun:
    3440# exec dettool -pretend -definebyquery -inst SIMPLE -filter B -det_type flat -select_exp_type flat -select_filter B -workdir path://SIMTEST/mkflat.0 -dbname simtest
    3541# exec dettool -pretend -definebyquery -inst SIMPLE -filter V -det_type flat -select_exp_type flat -select_filter V -workdir path://SIMTEST/mkflat.1 -dbname simtest
    … …  
    6167end
    6268
     69macro mkdarks
     70  echo "generating 20 dark images"
     71  for i 0 20
     72    $exptime = 3 * ($i + 1)
     73    mkimages.init 1024 1024 5 1000 0 2.0
     74    set dark = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime)*mki_noisei + $darkrate*$exptime + $biaslevel
     75    keyword dark AIRMASS  -wf {1 + 0.1*$i}
     76    keyword dark FILTER   -w  NONE
     77    keyword dark POSANGLE -wf 0.0
     78    keyword dark RA       -wf 0.0
     79    keyword dark DEC      -wf 0.0
     80    keyword dark EXPTIME  -wf $exptime
     81    keyword dark DARKTIME -wf $exptime
     82    keyword dark DATE-OBS -w "2006/10/01"
     83    sprintf time "00:01:%02d" $i
     84    keyword dark UTC-OBS  -w $time
     85    keyword dark CCDSUM   -w "1 1"
     86    keyword dark TELESCOP -w simscope
     87    keyword dark INSTRUME -w SIMPLE
     88    keyword dark OBSTYPE  -w DARK
     89    keyword dark OBJECT   -w DARK
     90    sprintf name "raw/dark.%02d.fits" $i
     91    wd dark $name
     92  end 
     93end
     94
    6395macro mkflats
    6496  echo "generating 20 flat images in 2 filters"
     97  $exptime = 10.0
    6598  for i 0 10
    6699    mkimages.init 1024 1024 5 1000 0 2.0
    67     set flat = blank + $mkiRD*mki_noisei + sqrt($flatlevel)*mki_noisei + $flatlevel + $biaslevel
     100    set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel)*mki_noisei + $darkrate*$exptime + $flatlevel + $biaslevel
    68101    keyword flat AIRMASS  -wf {1 + 0.1*$i}
    69102    keyword flat FILTER   -w  V
    … …  
    71104    keyword flat RA       -wf 0.0
    72105    keyword flat DEC      -wf 0.0
    73     keyword flat EXPTIME  -wf 5.0
    74     keyword flat DARKTIME -wf 0.0
     106    keyword flat EXPTIME  -wf $exptime
     107    keyword flat DARKTIME -wf $exptime
    75108    keyword flat DATE-OBS -w "2006/10/01"
    76109    sprintf time "00:01:%02d" $i
    … …  
    87120  for i 10 20
    88121    mkimages.init 1024 1024 5 1000 0 2.0
    89     set flat = blank + $mkiRD*mki_noisei + sqrt($flatlevel)*mki_noisei + $flatlevel + $biaslevel
     122    set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel)*mki_noisei + $darkrate*$exptime + $flatlevel + $biaslevel
    90123    keyword flat AIRMASS  -wf {1 + 0.1*$i}
    91124    keyword flat FILTER   -w  B
    … …  
    93126    keyword flat RA       -wf 0.0
    94127    keyword flat DEC      -wf 0.0
    95     keyword flat EXPTIME  -wf 5.0
    96     keyword flat DARKTIME -wf 0.0
     128    keyword flat EXPTIME  -wf $exptime
     129    keyword flat DARKTIME -wf $exptime
    97130    keyword flat DATE-OBS -w "2006/10/01"
    98131    sprintf time "00:02:%02d" $i
    … …  
    111144  echo "generating 10 object images in 2 filters"
    112145  for i 0 5
    113     mkimages.init 1024 1024 5 1000 0 2.0
    114     set flat = blank + $mkiRD*mki_noisei + sqrt($skylevel)*mki_noisei + $skylevel + $biaslevel
    115     keyword flat AIRMASS  -wf {1 + 0.1*$i}
    116     keyword flat FILTER   -w  V
    117     keyword flat POSANGLE -wf 0.0
    118     keyword flat RA       -wf 0.0
    119     keyword flat DEC      -wf 0.0
    120     keyword flat EXPTIME  -wf 5.0
    121     keyword flat DARKTIME -wf 0.0
    122     keyword flat DATE-OBS -w "2006/10/01"
     146    $exptime = 3 * ($i + 1)
     147    mkimages.init 1024 1024 5 1000 0 2.0
     148    set obj = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$exptime)*mki_noisei + $darkrate*$exptime + $skyrate*$exptime + $biaslevel
     149    keyword obj AIRMASS  -wf {1 + 0.1*$i}
     150    keyword obj FILTER   -w  V
     151    keyword obj POSANGLE -wf 0.0
     152    keyword obj RA       -wf 0.0
     153    keyword obj DEC      -wf 0.0
     154    keyword obj EXPTIME  -wf $exptime
     155    keyword obj DARKTIME -wf $exptime
     156    keyword obj DATE-OBS -w "2006/10/01"
    123157    sprintf time "00:01:%02d" $i
    124     keyword flat UTC-OBS  -w $time
    125     keyword flat CCDSUM   -w "1 1"
    126     keyword flat TELESCOP -w simscope
    127     keyword flat INSTRUME -w SIMPLE
    128     keyword flat OBSTYPE  -w OBJECT
    129     keyword flat OBJECT   -w OBJECT
     158    keyword obj UTC-OBS  -w $time
     159    keyword obj CCDSUM   -w "1 1"
     160    keyword obj TELESCOP -w simscope
     161    keyword obj INSTRUME -w SIMPLE
     162    keyword obj OBSTYPE  -w OBJECT
     163    keyword obj OBJECT   -w OBJECT
    130164    sprintf name "raw/obj.%02d.fits" $i
    131     wd flat $name
     165    wd obj $name
    132166  end 
    133167
    134168  for i 5 10
    135     mkimages.init 1024 1024 5 1000 0 2.0
    136     set flat = blank + $mkiRD*mki_noisei + sqrt($skylevel)*mki_noisei + $skylevel + $biaslevel
     169    $exptime = 3 * ($i - 4)
     170    mkimages.init 1024 1024 5 1000 0 2.0
     171    set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$exptime)*mki_noisei + $darkrate*$exptime + $skyrate*$exptime + $biaslevel
    137172    keyword flat AIRMASS  -wf {1 + 0.1*$i}
    138173    keyword flat FILTER   -w  B
    … …  
    140175    keyword flat RA       -wf 0.0
    141176    keyword flat DEC      -wf 0.0
    142     keyword flat EXPTIME  -wf 5.0
    143     keyword flat DARKTIME -wf 0.0
     177    keyword flat EXPTIME  -wf $exptime
     178    keyword flat DARKTIME -wf $exptime
    144179    keyword flat DATE-OBS -w "2006/10/01"
    145180    sprintf time "00:02:%02d" $i
Note: See TracChangeset for help on using the changeset viewer.