gendaymtx returns -1.#J -1.#J -1.#J

Here are the solid angles (in steradians) for the full Klems rows:

0.023909417
0.0237731163
0.0234119493
0.0273807841
0.0293333813
0.0349581625
0.0478687897
0.051801754
0.135517335

Note that we start at the zenith in the Klems subdivision, not at the horizon the way Tregenza's goes.

-Greg

P.S. I used the command:

  cnt 145 | rcalc -f klems_ang.cal -e 'kbin=recno;$1=omega' | uniq

···

From: Mostapha Sadeghipour <[email protected]>
Date: October 3, 2013 10:21:55 AM PDT

Hi Greg,

Thank you for checking. My bad! I use my own written epw2wea convertor so the user don't need to download epw2wea file separately. That was a bug indeed. Sorry about that.

One more question about conversion of the results of gendaymtx. Where can I find solid angle subtended for each patch of Reinhart sky? I already have the values for Tregenza sky as (0.0435, 0.0416, 0.0474, 0.0407, 0.0429, 0.0445, 0.0455, 0.0344) for the 8 rows. Is it fine to use the same numbers for every two rows?

Thanks again,
Mostapha

On Thu, Oct 3, 2013 at 12:21 PM, Gregory J. Ward <[email protected]> wrote:
Hi Mostapha,

Thanks for sending your files. For some reason, your WEA file has the month and day columns swapped, which is the source of your error. What did you use to translate the EPW file? The latest HEAD includes Christoph Reinharts epw2wea program, which does the correct translation. With the correct input, gendaymtx seems to produce appropriate output.

Best,
-Greg

From: Greg Ward <[email protected]>
Date: October 3, 2013 8:37:15 AM PDT

Hi Mostapha,

This sounds like a bug in gendaymtx. Some versions of the standard library print out "-1.#J" for NaN (Not a Number), which can result from a divide-by-zero or similar illegal math operation. If you send me your weather file in a separate message (not to the mailing list), I will see if I can reproduce the problem. There may be a corner case or two that gendaymtx still doesn't deal with gracefully.

Cheers,
-Greg

From: Mostapha Sadeghipour <[email protected]>
Date: October 2, 2013 1:37:49 PM PDT

Hi all,

I tried to read the results of the gendaymtx with a custom script and I noticed that there are multiple values of -1.#J -1.#J -1.#J for RGB. Does anybody know why this is happening?

I'm trying to get the total solar radiation for sky patches (-O1). I'm running RADIANCE on Windows and here is the line:

gendaymtx -m 1 -O1 .\USA_AK_Anchorage.Intl.AP.702730_TMY3.wea> .\USA_AK_Anchorage.Intl.AP.702730_TMY3.mtx

Thanks in advance,
Mostapha