ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/dayfact.csh
(Generate patch)

Comparing ray/src/util/dayfact.csh (file contents):
Revision 2.1 by greg, Tue Nov 12 17:19:31 1991 UTC vs.
Revision 2.7 by greg, Thu Oct 6 05:49:44 2005 UTC

# Line 1 | Line 1
1   #!/bin/csh -f
2 < # SCCSid "$SunId$ LBL"
2 > # RCSid: $Id$
3   #
4   # Interactive script to calculate daylight factors
5   #
6   set nofile="none"
7 < set genskyf=$nofile
7 > set illumpic=$nofile
8   set octree=$nofile
9   set dfpict=$nofile
10   set ilpict=$nofile
11 + set dspict=$nofile
12 + set nodaysav=1
13   set fcopts=($*)
14 + set designlvl=500
15   set wporig=(0 0 0)
16   set wpsize=(1 1)
17 < set rtargs=(-ab 1 -ad 256 -as 128 -aa .15 -av .1 .1 .1)
17 > set rtargs=(-ab 1 -ad 256 -as 128 -aa .15 -av .3 .3 .3)
18   set maxres=128
19  
20   alias readvar 'echo -n Enter \!:1 "[$\!:1]: ";set ans="$<";if("$ans" != "")set \!:1="$ans"'
# Line 19 | Line 22 | alias readvar 'echo -n Enter \!:1 "[$\!:1]: ";set ans=
22   cat <<_EOF_
23                          DAYLIGHT FACTOR CALCULATION
24  
25 < This script calculates daylight factors and/or illuminance levels on a
26 < rectangular workplane and produces a contour plot from the result.  The
27 < input is a Radiance scene description (and octree) and the output is one
28 < or more color Radiance pictures.
25 > This script calculates daylight factors, illuminance levels, and/or
26 > energy savings due to daylight on a rectangular workplane and produces
27 > a contour plot from the result.  The input is a Radiance scene description
28 > (and octree) and the output is one or more color Radiance pictures.
29  
30   _EOF_
31 + echo "Have you already calculated an illuminance picture using dayfact?"
32 + readvar illumpic
33 + if ( $illumpic != $nofile ) then
34 +        if ( ! -r $illumpic ) then
35 +                echo "Cannot read $illumpic"
36 +                exit 1
37 +        endif
38 +        set title=$illumpic:r
39 +        set gotillumpic
40 +        goto getgenskyf
41 + endif
42 +
43   readvar octree
44 < if ( "$octree" == "$nofile" || ! -f "$octree" ) then
44 > if ( $octree == $nofile || ! -f $octree ) then
45          echo "You must first create an octree with"
46          echo "oconv before running this script."
47          exit 1
48   endif
49   set title="$octree:r"
35 echo "In what scene file is the gensky command located?"
36 readvar genskyf
37 if ( "$genskyf" == "$nofile" || ! -r "$genskyf" ) then
38        echo "You will not be able to compute daylight"
39        echo "factors since there is no gensky file."
40 else
41        set title=$title\ `sed -n 's/^.*gensky  *\([0-9][0-9]*  *[0-9][0-9]*  *[0-9][0-9.]*\).*$/\1/p' $genskyf`
42        set extamb=`xform -e $genskyf|sed -n 's/^# Ground ambient level: //p'`
43 endif
50   echo -n "Is the z-axis your zenith direction? "
51   if ( "$<" !~ [yY]* ) then
52          echo "I'm sorry, you cannot use this script"
# Line 57 | Line 63 | set rtargs=($rtargs -ar `getinfo -d<$octree|rcalc -e '
63   echo "What calculation options do you want to give to rtrace?"
64   echo "(It is very important to set the -a* options correctly.)"
65   readvar rtargs
66 + echo "Do you want to save the illuminance picture for later runs?"
67 + readvar illumpic
68 + ############
69 + getgenskyf:
70 + set genskyf=$nofile
71 + echo "In what scene file is the gensky command located?"
72 + readvar genskyf
73 + if ( $genskyf == $nofile ) then
74 +        echo "You will not be able to compute daylight factors"
75 +        echo "or energy savings since there is no gensky file."
76 + else
77 +        xform -e $genskyf > /tmp/gsf$$
78 +        grep '^# gensky ' /tmp/gsf$$
79 +        if ( $status ) then
80 +                echo "The file $genskyf does not contain a gensky command\!"
81 +                rm -f /tmp/gsf$$
82 +                goto getgenskyf
83 +        endif
84 +        set title=$title\ `sed -n 's/^# gensky  *\([0-9][0-9]*  *[0-9][0-9]*  *[0-9][0-9.]*\).*$/\1/p' /tmp/gsf$$`
85 +        set extamb=`sed -n 's/^# Ground ambient level: //p' /tmp/gsf$$`
86 +        grep -s '^# gensky .* -c' /tmp/gsf$$
87 +        set nodaysav=$status
88 +        rm -f /tmp/gsf$$
89 +        if ( $nodaysav ) then
90 +                echo "The gensky command was not done for an overcast sky"
91 +                echo "(-c option), so energy savings cannot be calculated."
92 +                echo -n "Continue anyway? "
93 +                if ( "$<" =~ [nN]* ) then
94 +                        exit 0
95 +                endif
96 +        endif
97 + endif
98   echo "Illuminance contour picture if you want one"
99   readvar ilpict
100   if ( $?extamb ) then
101          echo "Daylight factor contour picture if you want one"
102          readvar dfpict
103   endif
104 < if ( "$ilpict" == "$nofile" && "$dfpict" == "$nofile" ) then
104 > if ( ! $nodaysav ) then
105 >        echo "Energy savings contour picture if you want one"
106 >        readvar dspict
107 >        if ( $dspict != $nofile ) then
108 >                echo "Workplane design level (lux)"
109 >                readvar designlvl
110 >        endif
111 > endif
112 > if ( $ilpict == $nofile && $dfpict == $nofile && $dspict == $nofile ) then
113          echo "Since you don't want any output, I guess we're done."
114          exit 0
115   endif
116   echo "Title for output picture"
117   readvar title
118 < set sctemp=/usr/tmp/sc$$.csh
118 > set sctemp=/tmp/sc$$.csh
119   cat <<'_EOF_' > $sctemp
120 < set iltemp=/usr/tmp/il$$.pic
121 < set tltemp=/usr/tmp/tl$$.pic
122 < set tempfiles=($iltemp $sctemp $tltemp)
120 > if ( $illumpic != $nofile ) then
121 >        set iltemp=""
122 > else
123 >        set iltemp=/tmp/il$$.pic
124 >        set illumpic=$iltemp
125 > endif
126 > set tltemp=/tmp/tl$$.pic
127 > set dstemp=/tmp/ds$$.pic
128 > set temp1=/tmp/tfa$$
129 > set tempfiles=($iltemp $sctemp $tltemp $dstemp $temp1)
130   echo "Your dayfact job is finished."
131   echo "Please check for error messages below."
132   echo ""
133   set echo
134 < cnt $wpres[2] $wpres[1] \
134 > if ( ! $?gotillumpic ) then
135 >        cnt $wpres[2] $wpres[1] \
136          | rcalc -e '$1=($2+.5)/'"$wpres[1]*$wpsize[1]+$wporig[1]" \
137                  -e '$2=(1-($1+.5)/'"$wpres[2])*$wpsize[2]+$wporig[2]" \
138                  -e '$3='"$wporig[3]" -e '$4=0;$5=0;$6=1' \
139 <        | rtrace $rtargs -h+ -I+ -ov -faf $octree \
140 <        | pvalue -r -x $wpres[1] -y $wpres[2] -df \
141 <        | pfilt -h 20 -n 0 -x 350 -y 350 -p 1 -r 1 > $iltemp
142 < set maxval=`getinfo < $iltemp | rcalc -i 'EXPOSURE=${e}' -e '$1=3/e'`
143 < psign -h 50 " $title " | pfilt -1 -x /2 -y /2 > $tltemp
139 >        | rtrace $rtargs -h+ -I+ -ov -fac -x $wpres[1] -y $wpres[2] $octree \
140 >        > $temp1
141 >        pfilt -h 20 -n 0 -x 300 -y 300 -p 1 -r 1 $temp1 > $illumpic
142 > endif
143 > set maxval=`getinfo < $illumpic | rcalc -i 'EXPOSURE=${e}' -e '$1=3/e'`
144 > if ("$maxval" == "") set maxval=3
145 > psign -h 42 " $title " | pfilt -1 -x /2 -y /2 > $tltemp
146   '_EOF_'
147 < if ( "$ilpict" != "$nofile" ) then
148 <        echo 'falsecolor -cb -l Lux -s "$maxval*470" \\
149 <                $fcopts -m 470 -ip $iltemp \\
147 > if ( $ilpict != $nofile ) then
148 >        echo 'falsecolor -cb -l Lux -s "$maxval*179" \\
149 >                $fcopts -m 179 -ip $illumpic \\
150                  | pcompos -a 1 - $tltemp > $ilpict' >> $sctemp
151   endif
152 < if ( "$dfpict" != "$nofile" ) then
153 <        echo 'falsecolor -cb -l DF -s 100 \\
154 <                $fcopts -m "100/PI/$extamb" -ip $iltemp \\
152 > if ( $dfpict != $nofile ) then
153 >        echo 'falsecolor -cb -l DF -s 16 \\
154 >                $fcopts -m "100/PI/$extamb" -ip $illumpic \\
155                  | pcompos -a 1 - $tltemp > $dfpict' >> $sctemp
156 + endif
157 + if ( $dspict != $nofile ) then
158 +        echo 'pcomb -e "lo=1-$designlvl/20000*3.1416*$extamb/li(1)" \\
159 +                -o $illumpic | falsecolor -cb -l "%Save" -s 100 \\
160 +                $fcopts -m 100 -p $illumpic \\
161 +                | pcompos -a 1 - $tltemp > $dspict' >> $sctemp
162   endif
163   echo 'rm -f $tempfiles' >> $sctemp
164   (source $sctemp) |& mail `whoami` &

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines