calcGef: Irradiation and irradiance on the generator plane.

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/calcGef.R

Description

This function obtains the global, diffuse and direct irradiation and irradiance on the generator plane from the values of daily or intradaily global irradiation on the horizontal plane. It makes use of the functions calcG0, fTheta, fInclin. Besides, it can calculate the shadows effect with the calcShd function.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
calcGef(lat,
        modeTrk='fixed',
        modeRad='prom',
        dataRad,
        prev, prom, mapa, bd,bdI,
        sample='hour',
        keep.night=TRUE,
        sunGeometry='michalsky',
        corr, f,
        betaLim=90, beta=abs(lat)-10, alfa=0,
        iS=2, alb=0.2, horizBright=TRUE, HCPV=FALSE,
        modeShd='',
        struct=list(),
        distances=data.frame())

Arguments

lat

numeric, latitude (degrees) of the point of the Earth where calculations are needed. It is positive for locations above the Equator.

modeTrk

character, to be chosen from 'fixed', 'two' or 'horiz'. When modeTrk='fixed' the surface is fixed (inclination and azimuth angles are constant). The performance of a two-axis tracker is calculated with modeTrk='two', and modeTrk='horiz' is the option for an horizontal N-S tracker. Its default value is modeTrk='fixed'

modeRad, dataRad, prom, mapa, bd, bdI

Information about the source data of the global irradiation. See calcG0 for details.

prev

A G0 object (or something that can be coerced to a G0 object). It is used when modeRad='prev'. It is deprecated and dataRad should be used instead.

sample, keep.night

See calcSol for details.

sunGeometry

character, method for the sun geometry calculations. See calcSol, fSolD and fSolI.

corr, f

See calcG0 for details.

beta

numeric, inclination angle of the surface (degrees). It is only needed when modeTrk='fixed'.

betaLim

numeric, maximum value of the inclination angle for a tracking surface. Its default value is 90 (no limitation))

alfa

numeric, azimuth angle of the surface (degrees). It is positive to the West. It is only needed when modeTrk='fixed'. Its default value is alfa=0

iS

integer, degree of dirtiness. Its value must be included in the set (1,2,3,4). iS=1 corresponds to a clean surface while iS=4 is the selection for a dirty surface. Its default value is 2.

alb

numeric, albedo reflection coefficient. Its default value is 0.2

modeShd, struct, distances

See calcShd for details.

horizBright

logical, if TRUE, the horizon brightness correction proposed by Reind et al. is used.

HCPV

logical, if TRUE the diffuse and albedo components of the effective irradiance are set to zero. HCPV is the acronym of High Concentration PV system.

Value

A Gef object.

Author(s)

Oscar Perpiñán Lamigueiro.

References

See Also

calcG0, fTheta, fInclin, calcShd.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
###12 Average days.

G0dm=c(2.766,3.491,4.494,5.912,6.989,7.742,7.919,7.027,5.369,3.562,2.814,2.179)*1000;
Ta=c(10, 14.1, 15.6, 17.2, 19.3, 21.2, 28.4, 29.9, 24.3, 18.2, 17.2, 15.2)

##Fixed surface, default values of inclination and azimuth.

gef<-calcGef(lat=37.2, modeRad='prom', dataRad=list(G0dm=G0dm, Ta=Ta))
print(gef)
xyplot(gef)

##Two-axis surface, no limitation angle.

gef2<-calcGef(lat=37.2, modeRad='prom', dataRad=list(G0dm=G0dm, Ta=Ta), modeTrk='two')
print(gef2)
xyplot(gef2)

##Fixed surface
gefAguiar <- calcGef(lat=41, modeRad='aguiar', dataRad=G0dm)

##Two-axis tracker, using the previous result.
##'gefAguiar' is internally coerced to a 'G0' object.

gefAguiar2 <- calcGef(lat=41, modeRad='prev', dataRad=gefAguiar, modeTrk='two')
print(gefAguiar2)
xyplot(gefAguiar2)

###Shadows between two-axis trackers, again using the gefAguiar result.

struct=list(W=23.11, L=9.8, Nrow=2, Ncol=8)
distances=data.frame(Lew=40, Lns=30, H=0)

gefShd<-calcGef(lat=41, modeRad='prev',
                dataRad=gefAguiar, modeTrk='two',
                modeShd=c('area', 'prom'),
                struct=struct, distances=distances)
print(gefShd)
##The Gef0, Bef0 and Def0 values are the same as those contained in the
##                gefAguiar2 object

Example output

Loading required package: zoo

Attaching package: 'zoo'

The following objects are masked from 'package:base':

    as.Date, as.Date.numeric

Loading required package: lattice
Loading required package: latticeExtra
Loading required package: RColorBrewer
Time Zone set to UTC.

Object of class  Gef 

Source of meteorological information: prom- 

Latitude of source:  37.2 degrees
Latitude for calculations:  37.2 degrees

Monthly averages:
               Bod       Bnd       Gd       Dd       Bd     Gefd     Defd
Jan 2018  8.734722  4.889585 4.475137 1.205450 3.239099 4.207434 1.123637
Feb 2018  9.605390  4.987463 4.900181 1.459503 3.402073 4.610792 1.359468
Mar 2018 10.275636  5.208261 5.431596 1.771373 3.610527 5.118680 1.649543
Apr 2018 10.528948  6.442369 6.296541 1.929253 4.301912 5.932536 1.796901
May 2018 10.437208  7.636610 6.793157 1.933014 4.782856 6.380243 1.798179
Jun 2018 10.294239  9.107754 7.176597 1.761494 5.329488 6.728800 1.638456
Jul 2018 10.302746 10.015944 7.503074 1.539231 5.876272 7.049591 1.435190
Aug 2018 10.392720  8.601781 7.280359 1.555253 5.647399 6.864813 1.453270
Sep 2018 10.236603  6.648405 6.316107 1.659148 4.597587 5.964057 1.549835
Oct 2018  9.668861  4.500495 4.727579 1.545645 3.142544 4.452040 1.438902
Nov 2018  8.807771  4.596611 4.375751 1.249438 3.095195 4.117991 1.164572
Dec 2018  8.177756  3.430839 3.473967 1.129231 2.320640 3.270820 1.051404
             Befd
Jan 2018 3.062799
Feb 2018 3.224822
Mar 2018 3.435019
Apr 2018 4.090753
May 2018 4.529007
Jun 2018 5.031569
Jul 2018 5.554282
Aug 2018 5.358196
Sep 2018 4.373462
Oct 2018 2.986096
Nov 2018 2.932056
Dec 2018 2.202874

Yearly values:
          Bod      Bnd       Gd       Dd       Bd     Gefd     Defd     Befd
2018 3572.657 2316.292 2092.386 569.9012 1502.183 1969.056 531.0119 1424.107
-----------------
Mode of tracking:  fixed 
    Inclination:  27.2 
    Orientation:  0 
Object of class  Gef 

Source of meteorological information: prom- 

Latitude of source:  37.2 degrees
Latitude for calculations:  37.2 degrees

Monthly averages:
              Bod       Bnd        Gd       Dd       Bd      Gefd     Defd
Jan 2018 14.13536  4.889585  6.488819 1.441834 4.889585  6.315669 1.385165
Feb 2018 15.42754  4.987463  6.831049 1.675463 4.987463  6.636480 1.602692
Mar 2018 16.59044  5.208261  7.369418 1.992999 5.208261  7.144765 1.898037
Apr 2018 17.65060  6.442369  8.871764 2.259771 6.442369  8.606524 2.152656
May 2018 18.71581  7.636610 10.193612 2.390555 7.636610  9.900628 2.280694
Jun 2018 19.87688  9.107754 11.611074 2.328790 9.107754 11.297920 2.229562
Jul 2018 18.51450 10.015944 11.786489 2.034929 9.568613 11.480652 1.953898
Aug 2018 17.33412  8.601781 10.747814 1.957717 8.601781 10.463678 1.878950
Sep 2018 16.24384  6.648405  8.788987 1.957080 6.648405  8.542775 1.873084
Oct 2018 15.16121  4.500495  6.375531 1.714157 4.500495  6.183453 1.634005
Nov 2018 13.99889  4.596611  6.206646 1.454352 4.596611  6.035891 1.394325
Dec 2018 12.70444  3.430839  4.816029 1.254509 3.430839  4.676497 1.198636
             Befd
Jan 2018 4.791793
Feb 2018 4.887714
Mar 2018 5.104096
Apr 2018 6.313522
May 2018 7.483878
Jun 2018 8.925599
Jul 2018 9.377241
Aug 2018 8.429746
Sep 2018 6.515437
Oct 2018 4.410485
Nov 2018 4.504678
Dec 2018 3.362223

Yearly values:
          Bod      Bnd       Gd       Dd       Bd    Gefd     Defd     Befd
2018 5972.909 2316.292 3046.733 683.3005 2302.425 2961.44 653.4751 2256.376
-----------------
Mode of tracking:  two 
    Inclination limit: 90 
Warning message:
In checkG0Ta(file, maxmin = TRUE) :
  Ambient temperature data is not available. A new column with a constant value has been added.
Object of class  Gef 

Source of meteorological information: bd-aguiar 

Latitude of source:  41 degrees
Latitude for calculations:  41 degrees

Monthly averages:
              Bod      Bnd        Gd       Dd       Bd      Gefd     Defd
Jan 2018 13.03684 5.728808  7.296037 1.430154 5.696204  7.113004 1.380318
Feb 2018 14.66766 5.207308  7.126519 1.742846 5.202936  6.925320 1.668397
Mar 2018 16.17329 5.150032  7.412249 2.111674 5.120133  7.181894 2.010001
Apr 2018 17.64733 6.316460  8.958839 2.457335 6.316460  8.683089 2.338579
May 2018 19.35631 7.249819 10.208367 2.800317 7.227588  9.894338 2.663035
Jun 2018 19.86122 8.074723 11.202331 2.940279 8.074723 10.865678 2.798886
Jul 2018 19.71148 9.113853 11.976680 2.709651 9.069337 11.643334 2.593204
Aug 2018 17.71146 7.817770 10.633483 2.624763 7.806599 10.326138 2.508322
Sep 2018 16.80571 6.682766  9.003549 2.175523 6.631300  8.744328 2.079072
Oct 2018 14.77668 4.330195  6.377514 1.888522 4.318021  6.176764 1.796580
Nov 2018 13.46521 5.346747  6.927312 1.450989 5.310629  6.748435 1.397603
Dec 2018 12.70720 4.172781  5.666588 1.352989 4.172781  5.514191 1.299567
             Befd
Jan 2018 5.582280
Feb 2018 5.098877
Mar 2018 5.017730
Apr 2018 6.190130
May 2018 7.083036
Jun 2018 7.913229
Jul 2018 8.887951
Aug 2018 7.650467
Sep 2018 6.498674
Oct 2018 4.231661
Nov 2018 5.204416
Dec 2018 4.089325

Yearly values:
         Bod      Bnd       Gd       Dd       Bd     Gefd     Defd     Befd
2018 5961.75 2288.886 3129.002 781.9837 2281.406 3038.494 746.9211 2235.778
-----------------
Mode of tracking:  two 
    Inclination limit: 90 
Object of class  Gef 

Source of meteorological information: bd-aguiar 

Latitude of source:  41 degrees
Latitude for calculations:  41 degrees

Monthly averages:
             Gef0d    Def0d    Bef0d        Gd       Dd       Bd      Gefd
Jan 2018  7.113004 1.380318 5.582280  7.296037 1.430154 5.696204  6.859068
Feb 2018  6.925320 1.668397 5.098877  7.126519 1.742846 5.202936  6.681436
Mar 2018  7.181894 2.010001 5.017730  7.412249 2.111674 5.120133  6.908321
Apr 2018  8.683089 2.338579 6.190130  8.958839 2.457335 6.316460  8.401855
May 2018  9.894338 2.663035 7.083036 10.208367 2.800317 7.227588  9.707779
Jun 2018 10.865678 2.798886 7.913229 11.202331 2.940279 8.074723 10.747494
Jul 2018 11.643334 2.593204 8.887951 11.976680 2.709651 9.069337 11.441005
Aug 2018 10.326138 2.508322 7.650467 10.633483 2.624763 7.806599 10.065549
Sep 2018  8.744328 2.079072 6.498674  9.003549 2.175523 6.631300  8.330185
Oct 2018  6.176764 1.796580 4.231661  6.377514 1.888522 4.318021  6.000919
Nov 2018  6.748435 1.397603 5.204416  6.927312 1.450989 5.310629  6.515058
Dec 2018  5.514191 1.299567 4.089325  5.666588 1.352989 4.172781  5.292413
             Defd     Befd
Jan 2018 1.340229 5.368433
Feb 2018 1.627854 4.895537
Mar 2018 1.962436 4.791722
Apr 2018 2.290722 5.956752
May 2018 2.631752 6.927759
Jun 2018 2.779525 7.814405
Jul 2018 2.561465 8.717361
Aug 2018 2.465291 7.432908
Sep 2018 2.013053 6.150551
Oct 2018 1.765422 4.086974
Nov 2018 1.359627 5.009015
Dec 2018 1.260091 3.907024

Yearly values:
        Gef0d    Def0d    Bef0d       Gd       Dd       Bd     Gefd    Defd
2018 3038.494 746.9211 2235.778 3129.002 781.9837 2281.406 2951.445 732.455
         Befd
2018 2163.194
-----------------
Mode of tracking:  two 
    Inclination limit: 90 

solaR documentation built on May 2, 2019, 6:07 p.m.

Related to calcGef in solaR...