| 1 | greg | 1.1 | #!/bin/csh -f | 
| 2 | greg | 2.2 | # RCSid: $Id$ | 
| 3 | greg | 1.1 | # | 
| 4 |  |  | #  Generate source description for window with venetian blinds. | 
| 5 |  |  | # | 
| 6 |  |  | #       1/11/88         Greg Ward | 
| 7 |  |  | # | 
| 8 |  |  | #  Usage: genwindow shellfile | 
| 9 |  |  | # | 
| 10 |  |  | #  Takes as input the following shell variables from file $1: | 
| 11 |  |  | # | 
| 12 |  |  | #       worient         Window orientation (degrees from south) | 
| 13 |  |  | #       wwidth          Window width | 
| 14 |  |  | #       wheight         Window height | 
| 15 |  |  | #       wtrans          Window transmission | 
| 16 |  |  | #       bdepth          Blind depth | 
| 17 |  |  | #       bspac           Blind vertical spacing | 
| 18 |  |  | #       gap             Gap between blinds and window | 
| 19 |  |  | #       brcurv          Blind curvature radius (negative upward, zero none) | 
| 20 |  |  | #       bangle          Blind inclination (degrees from horizontal, up and out) | 
| 21 |  |  | #       material        Blind material type (metal or plastic) | 
| 22 |  |  | #       ucolor          Blind up side color ("red green blue") | 
| 23 |  |  | #       uspec           Blind up side specularity (0 to 1) | 
| 24 |  |  | #       urough          Blind up side roughness (0 to .2) | 
| 25 |  |  | #       dcolor          Blind down side color (opt) | 
| 26 |  |  | #       dspec           Blind down side specularity (opt) | 
| 27 |  |  | #       drough          Blind down side roughness (opt) | 
| 28 |  |  | #       lat             Site latitude (degrees) | 
| 29 |  |  | #       long            Site longitude (degrees) | 
| 30 |  |  | #       mer             Site standard meridian (degrees) | 
| 31 |  |  | #       hour            Hour (standard 24hr time) | 
| 32 |  |  | #       day             Day (of month) | 
| 33 |  |  | #       month           Month (of year) | 
| 34 |  |  | #       grefl           Ground plane reflectance | 
| 35 |  |  | #       sky             Sky conditions (sunny, clear, cloudy) | 
| 36 |  |  | #       nsources        Number of sources to divide window | 
| 37 |  |  | # | 
| 38 |  |  | #  Creates the following output files: | 
| 39 |  |  | # | 
| 40 |  |  | #       stdout          Window description for Radiance | 
| 41 |  |  | #       $1.d$$          Window output distribution | 
| 42 |  |  | # | 
| 43 |  |  |  | 
| 44 |  |  | if ( $#argv != 1 ) then | 
| 45 |  |  | echo Usage: $0 input | 
| 46 |  |  | exit 1 | 
| 47 |  |  | endif | 
| 48 |  |  |  | 
| 49 |  |  | #       Set input and output files | 
| 50 |  |  | set input = $1 | 
| 51 |  |  | set distrib = $1.d$$ | 
| 52 |  |  | set tmpdir = /tmp | 
| 53 |  |  | set descrip = $tmpdir/gw$$.des | 
| 54 |  |  | set distoct = $tmpdir/gw$$.oct | 
| 55 |  |  | set skydesc = $tmpdir/gw$$.sky | 
| 56 |  |  | set remove = ( $distoct $skydesc $descrip ) | 
| 57 |  |  | set removerr = ( $remove $distrib ) | 
| 58 |  |  | onintr error | 
| 59 |  |  |  | 
| 60 |  |  | #       Set default values | 
| 61 |  |  | set worient = 0. | 
| 62 |  |  | set wtrans = ( .96 .96 .96 ) | 
| 63 |  |  | set bdepth = 0.025 | 
| 64 |  |  | set bspac = 0.0155 | 
| 65 |  |  | set gap = 0. | 
| 66 |  |  | set brcurv = 0. | 
| 67 |  |  | set bangle = 0. | 
| 68 |  |  | set material = plastic | 
| 69 |  |  | set ucolor = ( 0.5 0.5 0.5 ) | 
| 70 |  |  | set uspec = 0. | 
| 71 |  |  | set urough = 0.05 | 
| 72 |  |  | set lat = 37.8 | 
| 73 |  |  | set long = 122 | 
| 74 |  |  | set mer = 120 | 
| 75 |  |  | set hour = 12 | 
| 76 |  |  | set day = 21 | 
| 77 |  |  | set month = 10 | 
| 78 |  |  | set grefl = .2 | 
| 79 |  |  | set sky = sunny | 
| 80 |  |  | set nsources = 6 | 
| 81 |  |  |  | 
| 82 |  |  | #       Get input | 
| 83 |  |  | source $input | 
| 84 |  |  | if ( $status ) goto error | 
| 85 |  |  |  | 
| 86 |  |  | #       Create window | 
| 87 |  |  | cat > $descrip <<_EOF_ | 
| 88 |  |  |  | 
| 89 |  |  | void glass clear_glass | 
| 90 |  |  | 0 | 
| 91 |  |  | 0 | 
| 92 |  |  | 3 $wtrans | 
| 93 |  |  |  | 
| 94 |  |  | clear_glass polygon window | 
| 95 |  |  | 0 | 
| 96 |  |  | 0 | 
| 97 |  |  | 12 | 
| 98 |  |  | 0       0       0 | 
| 99 |  |  | 0       0       $wheight | 
| 100 |  |  | $wwidth 0       $wheight | 
| 101 |  |  | $wwidth 0       0 | 
| 102 |  |  |  | 
| 103 |  |  | void $material blind_upmat | 
| 104 |  |  | 0 | 
| 105 |  |  | 0 | 
| 106 |  |  | 5 $ucolor $uspec $urough | 
| 107 |  |  | _EOF_ | 
| 108 |  |  | if ( $status ) goto error | 
| 109 |  |  |  | 
| 110 |  |  | #       Blinds | 
| 111 |  |  | genblinds blind_upmat blind_up $bdepth $wwidth $wheight \ | 
| 112 |  |  | `ev "ceil($wheight/$bspac)"` $bangle -r $brcurv \ | 
| 113 |  |  | | xform -t $gap -$wwidth 1e-4 -rz 90 >> $descrip | 
| 114 |  |  | if ( $status ) goto error | 
| 115 |  |  | if ( $?dcolor ) then | 
| 116 |  |  | cat >> $descrip <<_EOF_ | 
| 117 |  |  |  | 
| 118 |  |  | void $material blind_dnmat | 
| 119 |  |  | 0 | 
| 120 |  |  | 0 | 
| 121 |  |  | 5 $dcolor $dspec $drough | 
| 122 |  |  | _EOF_ | 
| 123 |  |  | genblinds blind_dnmat blind_down $bdepth $wwidth $wheight \ | 
| 124 |  |  | `ev "ceil($wheight/$bspac)"` $bangle -r $brcurv \ | 
| 125 |  |  | | xform -t $gap -$wwidth 0 -rz 90 >> $descrip | 
| 126 |  |  | if ( $status ) goto error | 
| 127 |  |  | endif | 
| 128 |  |  |  | 
| 129 |  |  | #       Make sky | 
| 130 |  |  | switch ($sky) | 
| 131 |  |  | case sun*: | 
| 132 |  |  | set skysw = +s | 
| 133 |  |  | breaksw | 
| 134 |  |  | case clear: | 
| 135 |  |  | set skysw = -s | 
| 136 |  |  | breaksw | 
| 137 |  |  | case cloud*: | 
| 138 |  |  | set skysw = -c | 
| 139 |  |  | breaksw | 
| 140 |  |  | endsw | 
| 141 |  |  | dosky: | 
| 142 |  |  | gensky $month $day $hour $skysw -a $lat -o $long -m $mer -g $grefl \ | 
| 143 |  |  | | xform -rz `ev "-($worient)"` > $skydesc | 
| 144 |  |  | if ( $skysw == +s ) then | 
| 145 |  |  | if ( `sed -n 13p $skydesc | rcalc -e '$1=if($3,-1,1)'` < 0 ) then | 
| 146 |  |  | set skysw = -s | 
| 147 |  |  | goto dosky | 
| 148 |  |  | endif | 
| 149 |  |  | endif | 
| 150 |  |  | cat >> $skydesc <<_EOF_ | 
| 151 |  |  |  | 
| 152 |  |  | skyfunc glow skyglow | 
| 153 |  |  | 0 | 
| 154 |  |  | 0 | 
| 155 |  |  | 4 1 1 1 0 | 
| 156 |  |  |  | 
| 157 |  |  | skyglow source sky | 
| 158 |  |  | 0 | 
| 159 |  |  | 0 | 
| 160 |  |  | 4 0 -1 0 180 | 
| 161 |  |  | _EOF_ | 
| 162 |  |  |  | 
| 163 |  |  | #       Make distribution | 
| 164 |  |  | oconv $skydesc $descrip > $distoct | 
| 165 |  |  | if ( $status ) goto error | 
| 166 |  |  | echo 2 5 85 9 0 340 18 > $distrib | 
| 167 |  |  | makedist -h -d -x1 0 1 0 -x2 0 0 1 -x3 1 0 0 -alpha 5-85:10 -beta 0-340:20 \ | 
| 168 |  |  | -tw $bspac -th $bspac -td `ev "2*($bdepth+$gap)"` \ | 
| 169 |  |  | -tc `ev "$wwidth/2" "$bdepth+$gap" "$wheight/2"` \ | 
| 170 |  |  | -ab 2 -aa .1 -ad 64 -as 64 \ | 
| 171 |  |  | -x 16 -y 16 $distoct >> $distrib | 
| 172 |  |  | if ( $status ) goto error | 
| 173 |  |  |  | 
| 174 |  |  | set wsgrid = (`ev "floor(sqrt($nsources*$wwidth/$wheight)+.5)" "floor(sqrt($nsources*$wheight/$wwidth)+.5)"`) | 
| 175 |  |  | @ nsources = $wsgrid[1] * $wsgrid[2] | 
| 176 |  |  | #       Print header | 
| 177 |  |  | cat <<_EOF_ | 
| 178 |  |  | # | 
| 179 |  |  | #  Window with venetian blinds | 
| 180 |  |  | #  Created from $input `date` | 
| 181 |  |  | # | 
| 182 |  |  | #       Window orientation (degrees from south):        $worient | 
| 183 |  |  | #       Window height:                                  $wheight | 
| 184 |  |  | #       Window width:                                   $wwidth | 
| 185 |  |  | #       Window transmission:                            $wtrans | 
| 186 |  |  | #       Blind depth:                                    $bdepth | 
| 187 |  |  | #       Blind spacing:                                  $bspac | 
| 188 |  |  | #       Gap to window:                                  $gap | 
| 189 |  |  | #       Blind curvature radius:                         $brcurv | 
| 190 |  |  | #       Blind inclination (degrees altitude):           $bangle | 
| 191 |  |  | #       Blind material:                                 $material | 
| 192 |  |  | #       Up side color:                                  $ucolor | 
| 193 |  |  | #       Up side specularity:                            $uspec | 
| 194 |  |  | #       Up side roughness:                              $urough | 
| 195 |  |  | _EOF_ | 
| 196 |  |  | if ( $?dcolor ) then | 
| 197 |  |  | cat <<_EOF_ | 
| 198 |  |  | #       Down side color:                                $dcolor | 
| 199 |  |  | #       Down side specularity:                          $dspec | 
| 200 |  |  | #       Down side roughness:                            $drough | 
| 201 |  |  | _EOF_ | 
| 202 |  |  | else | 
| 203 |  |  | echo \#\        Down side same as up | 
| 204 |  |  | endif | 
| 205 |  |  | cat <<_EOF_ | 
| 206 |  |  | #       Latitude (degrees):                             $lat | 
| 207 |  |  | #       Longitude (degrees):                            $long | 
| 208 |  |  | #       Standard Meridian (degrees):                    $mer | 
| 209 |  |  | #       Month Day Hour:                                 $month $day $hour | 
| 210 |  |  | #       Ground plane reflectance:                       $grefl | 
| 211 |  |  | #       Sky condition:                                  $sky | 
| 212 |  |  | #       Number of window sources:                       $nsources | 
| 213 |  |  | # | 
| 214 |  |  | _EOF_ | 
| 215 |  |  |  | 
| 216 |  |  | #       Send sources | 
| 217 |  |  | xform -e -rz $worient <<_EOF_ | 
| 218 |  |  |  | 
| 219 |  |  | void brightdata wdistrib | 
| 220 |  |  | 10 noop $distrib source.cal src_theta src_phi -rx 90 -ry -90 -mx | 
| 221 |  |  | 0 | 
| 222 |  |  | 0 | 
| 223 |  |  |  | 
| 224 |  |  | wdistrib illum willum | 
| 225 |  |  | 0 | 
| 226 |  |  | 0 | 
| 227 |  |  | 3 1 1 1 | 
| 228 |  |  |  | 
| 229 |  |  | !gensurf willum wsource "$wwidth*t" $bdepth "$wheight*s" $wsgrid[2] $wsgrid[1] | 
| 230 |  |  | _EOF_ | 
| 231 |  |  | if ( $status ) goto error | 
| 232 |  |  |  | 
| 233 |  |  | #       Send window | 
| 234 |  |  | xform -rz $worient $descrip | 
| 235 |  |  | if ( $status ) goto error | 
| 236 |  |  |  | 
| 237 |  |  | #       All done, print and exit | 
| 238 |  |  | rm -f $remove | 
| 239 |  |  | exit 0 | 
| 240 |  |  |  | 
| 241 |  |  | #       Error exit | 
| 242 |  |  | error: | 
| 243 |  |  | rm -f $removerr | 
| 244 |  |  | exit 1 |