<div dir="ltr"><div><div><div>Deal NCL Talk team<br><br></div>I got a little strange about my solar irradiation output as shown in attachment file using the code below.<br><br> rad = 4.0*atan(1)/180.<br>   cza = acos(cz)/rad<br>   cza@long_name = &quot;zenith angle in degrees&quot;<br><br>     xlat = in-&gt;XLAT<br>  <br> nt = 5     ; example specific time<br> nl = 5<br> ml = 4<br><br> ;  print(&quot; XLAT            cz             cza        &quot;)<br>    print(xlat(:,nl,ml)+&quot;     &quot;+cz(:,nl,ml)+&quot;      &quot;+cza(:,nl,ml))<br><br></div><br><br>Appreciate for your help<br><br></div>Soares<br><div><div><div><br></div></div></div></div>