Calculation of Perez model

I was looking into Perez model as implemented in SSC in cmod_pvwattsv0.cpp and comparing it to the technical handbook.

According to the module and other resources, the airmass should be



where b is the incidence angle.

Somehow in  this line  the equation for airmass is:
AIRMASS = 1.0 / (CZ + 0.15 * pow(93.9 - ZENITH, -1.253) );

and CZ in  this line  is:
CZ = cos(zen);

where zen is the sun's zenith in radiations and ZENITH the sun's zenith in degrees.

Shouldn't the equation for airmass replace CZ with cos(inc)? or am I missing something?


Thanks in advance.

Best,
Mustafa Alotbah

Hi Paul,

Thank you, I figured that out after I had published the post, but even in lib_irradproc.cpp,  this line  defines the airmass, which I think has the same issue, it uses CZ (Cos Zenith) instead of COSINC which is only defined afterwards.


Best regards,
Mustafa.

Hi Paul,

another thing,  in this line  B is calculated as 
B = D * F1 * ZC / ZH;

ZC being "a" in the handbook and ZH assumed to be "b", however, this does not seem to be correct, according to the handbook:



which is B in the code, and b is defined as 



while, ZH is the maxiumum of 85 degrees and CZ again, instead of ZC/COSINC, check these lines out: 

1967  CZ = cos(zen);
1968  ZH = (CZ > 0.0871557) ? CZ : 0.0871557;  /* Maximum of 85 deg */



Best regards,
Mustafa