| 1 |
# /bin/csh -f |
| 2 |
# create limb darkening removed Ic record and create image set for Ic_noLD and M for range of time |
| 3 |
# specified by first two arguments. |
| 4 |
|
| 5 |
# set echo |
| 6 |
|
| 7 |
set wantlow = $1 |
| 8 |
set wanthigh = $2 |
| 9 |
|
| 10 |
# from here make script that can make _nrt images in standard place |
| 11 |
|
| 12 |
set CADENCE = 15 |
| 13 |
|
| 14 |
set wantlow_t = `time_convert time=$wantlow` |
| 15 |
set wanthigh_t = `time_convert time=$wanthigh` |
| 16 |
|
| 17 |
@ daylow = $wantlow_t / 86400 |
| 18 |
@ dayhigh = $wanthigh_t / 86400 |
| 19 |
|
| 20 |
set end_movie_t = 0 |
| 21 |
|
| 22 |
set day = $daylow |
| 23 |
while ($day <= $dayhigh) |
| 24 |
@ day_t = 86400 * $day |
| 25 |
@ nextday_t = $day_t + 86400 |
| 26 |
@ last_t = $nextday_t - 45 |
| 27 |
if ($last_t > $wanthigh_t) set last_t = $wanthigh_t |
| 28 |
set first_t = $day_t |
| 29 |
if ($first_t < $wantlow_t) set first_t = $wantlow_t |
| 30 |
@ cadence = $CADENCE * 60 |
| 31 |
@ first_mod_t = $first_t - 45 |
| 32 |
@ imglow_t = $first_mod_t / $cadence |
| 33 |
@ imglow_t = $imglow_t + 1 |
| 34 |
@ first_t = $imglow_t * $cadence |
| 35 |
set first = `time_convert s=$first_t zone=TAI` |
| 36 |
@ last_mod_t = $last_t + $cadence |
| 37 |
@ imghigh_t = $last_mod_t / $cadence |
| 38 |
@ n_images = $imghigh_t - $imglow_t |
| 39 |
if ($n_images < 0) set n_images = 0 |
| 40 |
if ($n_images == 0) then |
| 41 |
@ day = $day + 1 |
| 42 |
continue |
| 43 |
endif |
| 44 |
set end_movie_t = $last_t |
| 45 |
@ n_minutes = $n_images * $CADENCE |
| 46 |
set QRY = '['$first'/'$n_minutes'm@'$CADENCE'm]' |
| 47 |
set yyyymmdd = `echo $first | sed -e 's/_.*//' -e 's/\./Q/' -e 's/\./X/'` |
| 48 |
set YEAR = `echo $yyyymmdd | sed -e 's/Q.*//'` |
| 49 |
set MON = `echo $yyyymmdd | sed -e 's/^.*Q//' -e 's/X.*//'` |
| 50 |
set DAY = `echo $yyyymmdd | sed -e 's/^.*X//'` |
| 51 |
|
| 52 |
set IMGPATH = /home/jsoc/hmi/images/$YEAR |
| 53 |
if (!(-e $IMGPATH)) mkdir $IMGPATH |
| 54 |
set IMGPATH = $IMGPATH/$MON |
| 55 |
if (!(-e $IMGPATH)) mkdir $IMGPATH |
| 56 |
set IMGPATH = $IMGPATH/$DAY |
| 57 |
if (!(-e $IMGPATH)) mkdir $IMGPATH |
| 58 |
|
| 59 |
set n_source = `show_info -cq 'hmi.Ic_45s_nrt'$QRY` |
| 60 |
set n_dest = `show_info -cq 'hmi.Ic_noLimbDark_45s'$QRY` |
| 61 |
if ($n_source > 0 && $n_dest < $n_source) then |
| 62 |
hmi_limbdark in='hmi.Ic_45s_nrt'$QRY out=hmi.Ic_noLimbDark_45s -n |
| 63 |
endif |
| 64 |
|
| 65 |
cd $IMGPATH |
| 66 |
|
| 67 |
setenv RGBDEF /home/phil/workflow/scripts/rgb.txt |
| 68 |
|
| 69 |
render_image in='hmi.Ic_45s_nrt'$QRY \ |
| 70 |
pallette=/home/phil/apps/heat.sao \ |
| 71 |
min=20000 \ |
| 72 |
max=75000 \ |
| 73 |
scale=1,4,8,16 \ |
| 74 |
outid=time \ |
| 75 |
type=jpg \ |
| 76 |
outname=Ic \ |
| 77 |
out='| ppmlabel -color white -size {%0.75:5} -x 15 -y {%98} -text "SDO/HMI Quick-Look Continuum: {ID}" | pnmtojpeg -quality=95' |
| 78 |
|
| 79 |
render_image in='hmi.Ic_noLimbDark_45s'$QRY \ |
| 80 |
pallette=/home/phil/apps/heat.sao \ |
| 81 |
min=0.3 \ |
| 82 |
max=1.2 \ |
| 83 |
scale=1,4,8,16 \ |
| 84 |
outid=time \ |
| 85 |
type=jpg \ |
| 86 |
outname=Ic_flat \ |
| 87 |
out='| ppmlabel -color white -size {%0.75:5} -x 15 -y {%98} -text "SDO/HMI Quick-Look Continuum: {ID}" | pnmtojpeg -quality=95' |
| 88 |
|
| 89 |
render_image in='hmi.M_45s_nrt'$QRY \ |
| 90 |
outname=M \ |
| 91 |
pallette=/home/phil/apps/mag.lut \ |
| 92 |
outid=time \ |
| 93 |
-c \ |
| 94 |
min=-500 \ |
| 95 |
max=500 \ |
| 96 |
type=jpg \ |
| 97 |
scale=1,4,8,16 \ |
| 98 |
out='| ppmlabel -color white -size {%0.75:5} -x 15 -y {%98} -text "SDO/HMI Quick-Look Magnetogram: {ID}" | pnmtojpeg -quality=95' |
| 99 |
|
| 100 |
@ day = $day + 1 |
| 101 |
end |
| 102 |
|
| 103 |
if ($end_movie_t > 0) then |
| 104 |
set end_movie = `time_convert s=$end_movie_t zone=TAI` |
| 105 |
/home/phil/workflow/scripts/makeIcmovie $end_movie |
| 106 |
endif |
