Description: "Weather data reader"

- Home > Modelica

## Buildings.BoundaryConditions.WeatherData

Weather data reader

### Information

This package contains models to read weather data. It also contains an expandable connector called `Bus` that is used in the library to provide weather data to the different models.

Extends from Modelica.Icons.VariantsPackage (Icon for package containing variants).

### Package Content

Name | Description
--- | ---
Bus | Data bus that stores weather data
ReaderTMY3 | Reader for TMY3 weather data
Examples | Collection of models that illustrate model use and test models
BaseClasses | Package with base classes for Buildings.BoundaryConditions.WeatherData

## Buildings.BoundaryConditions.WeatherData.Bus

Data bus that stores weather data

### Information

This component is an expandable connector that is used to implement a bus that contains the weather data.

Extends from Modelica.Icons.SignalBus (Icon for signal bus).

### Modelica definition

expandable connector Bus "Data bus that stores weather data" extends Modelica.Icons.SignalBus; end Bus;

## Buildings.BoundaryConditions.WeatherData.ReaderTMY3

Reader for TMY3 weather data

### Information

This component reads TMY3 weather data (Wilcox and Marion, 2008) or user specified weather data. The weather data format is the Typical Meteorological Year (TMY3) as obtained from the EnergyPlus web site at http://apps1.eere.energy.gov/buildings/energyplus/cfm/weather_data.cfm. These data, which are in the EnergyPlus format, need to be converted as described in the next paragraph.

#### Adding new weather data

To add new weather data, proceed as follows:

1. Download the weather data file with the `epw` extension from http://apps1.eere.energy.gov/buildings/energyplus/cfm/weather_data.cfm.
2. Add the file to `Buildings/Resources/weatherdata` (or to any directory for which you have write permission).
3. On a console window, type

         cd Buildings/Resources/weatherdata
         java -jar ../bin/ConvertWeatherData.jar inputFile.epw

   This will generate the weather data file `inputFile.mos`, which can be read by the model Buildings.BoundaryConditions.WeatherData.ReaderTMY3.

#### Location data that are read automatically from the weather data file

The following location data are automatically read from the weather file:

- The latitude of the weather station, `lat`,
- the longitude of the weather station, `lon`, and
- the time zone relative to Greenwich Mean Time, `timZone`.

#### Wet bulb temperature

By default, the data bus contains the wet bulb temperature. This introduces a nonlinear equation. However, we have not observed an increase in computing time because of this equation. To disable the computation of the wet bulb temperature, set `computeWetBulbTemperature=false`.

#### Using constant or user-defined input signals for weather data

This model has the option of using a constant value, using the data from the weather file, or using data from an input connector for the following variables:

- The atmospheric pressure,
- the ceiling height,
- the total sky cover pressure,
- the opaque sky cover pressure,
- the dry bulb temperature,
- the dew point temperature,
- the sky black body temperature,
- the relative humidity,
- the wind direction,
- the wind speed,
- the global horizontal radiation, direct normal and diffuse horizontal radiation, and
- the infrared horizontal radiation.

By default, all data are obtained from the weather data file, except for the atmospheric pressure, which is set to the parameter `pAtm=101325` Pascals.

The parameter `*Sou` configures the source of the data. For the atmospheric pressure, temperatures, relative humidity, wind speed and wind direction, the enumeration Buildings.BoundaryConditions.Types.DataSource is used as follows:

Parameter `*Sou` | Data used to compute weather data.
--- | ---
File | Use data from file.
Parameter | Use value specified by the parameter.
Input | Use value from the input connector.

Because global, diffuse and direct radiation are related to each other, the parameter `HSou` is treated differently. It is set to a value of the enumeration Buildings.BoundaryConditions.Types.RadiationDataSource, and allows the following configurations:

Parameter `HSou` | Data used to compute weather data.
--- | ---
File | Use data from file.
Input_HGloHor_HDifHor | Use global horizontal and diffuse horizontal radiation from input connector.
Input_HDirNor_HDifHor | Use direct normal and diffuse horizontal radiation from input connector.
Input_HDirNor_HGloHor | Use direct normal and global horizontal radiation from input connector.

Notes

1. In HVAC systems, when the fan is off, changes in atmospheric pressure can cause small air flow rates in the duct system due to change in pressure and hence in the mass of air that is stored in air volumes (such as in fluid junctions or in the room model). This may increase computing time. Therefore, the default value for the atmospheric pressure is set to a constant. Furthermore, if the initial pressure of air volumes are different from the atmospheric pressure, then fast pressure transients can happen in the first few seconds of the simulation. This can cause numerical problems for the solver. To avoid this problem, set the atmospheric pressure to the same value as the medium default pressure, which is typically set to the parameter `Medium.p_default`. For medium models for moist air and dry air, the default is `Medium.p_default=101325` Pascals.

2. Different units apply depending on whether data are obtained from a file, or from a parameter or an input connector:

   - When using TMY3 data from a file (e.g. `USA_IL_Chicago-OHare.Intl.AP.725300_TMY3.mos`), the units must be the same as the original TMY3 file used by EnergyPlus (e.g. `USA_IL_Chicago-OHare.Intl.AP.725300_TMY3.epw`). The TMY3 data used by EnergyPlus are in both SI units and non-SI units. If `Resources/bin/ConvertWeatherData.jar` is used to convert the `.epw` file to an `.mos` file, the units of the TMY3 data are preserved and the file can be directly used by this data reader. The data reader will automatically convert units to the SI units used by Modelica. For example, the dry bulb temperature `TDryBul` in TMY3 is in degree Celsius. The data reader will automatically convert the data to Kelvin. The wind direction `winDir` in TMY3 is degrees and will be automatically converted to radians.
   - When using data from a parameter or from an input connector, the data must be in the SI units used by Modelica. For instance, the unit must be `Pa` for pressure, `K` for temperature, `W/m2` for solar radiations and `rad` for wind direction.
3. The ReaderTMY3 should only be used with TMY3 data. It contains a time shift for solar radiation data that is explained below. This time shift needs to be removed if the user may want to use the ReaderTMY3 for other weather data types.

#### Implementation

##### Start and end data for annual weather data files

The TMY3 weather data, as well as the EnergyPlus weather data, start at 1:00 AM on January 1, and provide hourly data until midnight on December 31. Thus, the first entry for temperatures, humidity, wind speed etc. are values at 1:00 AM and not at midnight. Furthermore, the TMY3 weather data files can have values at midnight of December 31 that may be significantly different from the values at 1:00 AM on January 1. Since annual simulations require weather data that start at 0:00 on January 1, data need to be provided for this hour. Due to the possibly large change in weatherdata between 1:00 AM on January 1 and midnight at December 31, the weather data files in the Buildings library do not use the data entry from midnight at December 31 as the value for t=0. Rather, the value from 1:00 AM on January 1 is duplicated and used for 0:00 on January 1. To maintain a data record with 8760 hours, the weather data record from midnight at December 31 is deleted. These changes in the weather data file are done in the Java program that converts EnergyPlus weather data file to Modelica weather data files, and which is described below.

##### Time shift for solar radiation data

To read weather data from the TMY3 weather data file, there are two data readers in this model. One data reader obtains all data except solar radiation, and the other data reader reads only the solar radiation data, shifted by 30 minutes. The reason for this time shift is as follows: The TMY3 weather data file contains for solar radiation the "...radiation received on a horizontal surface during the 60-minute period ending at the timestamp." Thus, as the figure below shows, a more accurate interpolation is obtained if time is shifted by 30 minutes prior to reading the weather data.

#### References

- Wilcox S. and W. Marion. Users Manual for TMY3 Data Sets. Technical Report, NREL/TP-581-43156, revised May 2008.

### Parameters

Type | Name | Default | Description
--- | --- | --- | ---
Boolean | computeWetBulbTemperature | true | If true, then this model computes the wet bulb temperature
String | filNam | "" | Name of weather data file
Data source
DataSource | pAtmSou | Buildings.BoundaryConditions... | Atmospheric pressure
Pressure | pAtm | 101325 | Atmospheric pressure (used if pAtmSou=Parameter) [Pa]
DataSource | ceiHeiSou | Buildings.BoundaryConditions... | Ceiling height
Real | ceiHei | 20000 | Ceiling height (used if ceiHei=Parameter) [m]
DataSource | totSkyCovSou | Buildings.BoundaryConditions... | Total sky cover
Real | totSkyCov | 0.5 | Total sky cover (used if totSkyCov=Parameter). Use 0 <= totSkyCov <= 1 [1]
DataSource | opaSkyCovSou | Buildings.BoundaryConditions... | Opaque sky cover
Real | opaSkyCov | 0.5 | Opaque sky cover (used if opaSkyCov=Parameter). Use 0 <= opaSkyCov <= 1 [1]
DataSource | TDryBulSou | Buildings.BoundaryConditions... | Dry bulb temperature
Temperature | TDryBul | 293.15 | Dry bulb temperature (used if TDryBul=Parameter) [K]
DataSource | TDewPoiSou | Buildings.BoundaryConditions... | Dew point temperature
Temperature | TDewPoi | 283.15 | Dew point temperature (used if TDewPoi=Parameter) [K]
DataSource | TBlaSkySou | Buildings.BoundaryConditions... | Black-body sky temperature
Temperature | TBlaSky | 273.15 | Black-body sky temperature (used if TBlaSkySou=Parameter) [K]
DataSource | relHumSou | Buildings.BoundaryConditions... | Relative humidity
Real | relHum | 0.5 | Relative humidity (used if relHum=Parameter) [1]
DataSource | winSpeSou | Buildings.BoundaryConditions... | Wind speed
Velocity | winSpe | 1 | Wind speed (used if winSpe=Parameter) [m/s]
DataSource | winDirSou | Buildings.BoundaryConditions... | Wind direction
Angle | winDir | 1.0 | Wind direction (used if winDir=Parameter) [rad]
DataSource | HInfHorSou | Buildings.BoundaryConditions... | Infrared horizontal radiation
HeatFlux | HInfHor | 0.0 | Infrared horizontal radiation (used if HInfHorSou=Parameter) [W/m2]
RadiationDataSource | HSou | Buildings.BoundaryConditions... | Global, diffuse, and direct normal radiation
Sky temperature
SkyTemperatureCalculation | calTSky | Buildings.BoundaryConditions... | Computation of black-body sky temperature

### Connectors

Type | Name | Description
--- | --- | ---
input RealInput | pAtm_in | Input pressure [Pa]
input RealInput | ceiHei_in | Input ceiling height [m]
input RealInput | totSkyCov_in | Input total sky cover [1]
input RealInput | opaSkyCov_in | Input opaque sky cover [1]
input RealInput | TDryBul_in | Input dry bulb temperature [K]
input RealInput | TDewPoi_in | Input dew point temperature [K]
input RealInput | TBlaSky_in | Black-body sky temperature [K]
input RealInput | relHum_in | Input relative humidity [1]
input RealInput | winSpe_in | Input wind speed [m/s]
input RealInput | winDir_in | Input wind direction [rad]
input RealInput | HInfHor_in | Input infrared horizontal radiation [W/m2]
input RealInput | HGloHor_in | Input global horizontal radiation [W/m2]
input RealInput | HDifHor_in | Input diffuse horizontal radiation [W/m2]
input RealInput | HDirNor_in | Input direct normal radiation [W/m2]
Bus | weaBus | Weather data bus

### Modelica definition

block ReaderTMY3 "Reader for TMY3 weather data" parameter Boolean computeWetBulbTemperature = true "If true, then this model computes the wet bulb temperature"; //-------------------------------------------------------------- // Atmospheric pressure parameter Buildings.BoundaryConditions.Types.DataSource pAtmSou=Buildings.BoundaryConditions.Types.DataSource.Parameter "Atmospheric pressure"; parameter Modelica.SIunits.Pressure pAtm=101325 "Atmospheric pressure (used if pAtmSou=Parameter)"; Modelica.Blocks.Interfaces.RealInput pAtm_in( final quantity="Pressure", final unit="Pa", displayUnit="Pa") if (pAtmSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input pressure"; //-------------------------------------------------------------- // Ceiling height parameter Buildings.BoundaryConditions.Types.DataSource ceiHeiSou=Buildings.BoundaryConditions.Types.DataSource.File "Ceiling height"; parameter Real ceiHei( final quantity="Height", final unit="m", displayUnit="m") = 20000 "Ceiling height (used if ceiHei=Parameter)"; Modelica.Blocks.Interfaces.RealInput ceiHei_in( final quantity="Height", final unit="m", displayUnit="m") if (ceiHeiSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input ceiling height"; //-------------------------------------------------------------- // Total sky cover parameter Buildings.BoundaryConditions.Types.DataSource totSkyCovSou=Buildings.BoundaryConditions.Types.DataSource.File "Total sky cover"; parameter Real totSkyCov( min=0, max=1, unit="1") = 0.5 "Total sky cover (used if totSkyCov=Parameter). Use 0 <= totSkyCov <= 1"; Modelica.Blocks.Interfaces.RealInput totSkyCov_in( min=0, max=1, unit="1") if (totSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input total sky cover"; // Opaque sky cover parameter Buildings.BoundaryConditions.Types.DataSource opaSkyCovSou=Buildings.BoundaryConditions.Types.DataSource.File "Opaque sky cover"; parameter Real opaSkyCov( min=0, max=1, unit="1") = 0.5 "Opaque sky cover (used if opaSkyCov=Parameter). Use 0 <= opaSkyCov <= 1"; Modelica.Blocks.Interfaces.RealInput opaSkyCov_in( min=0, max=1, unit="1") if (opaSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input opaque sky cover"; //-------------------------------------------------------------- // Dry bulb temperature parameter Buildings.BoundaryConditions.Types.DataSource TDryBulSou=Buildings.BoundaryConditions.Types.DataSource.File "Dry bulb temperature"; parameter Modelica.SIunits.Temperature TDryBul(displayUnit="degC") = 293.15 "Dry bulb temperature (used if TDryBul=Parameter)"; Modelica.Blocks.Interfaces.RealInput TDryBul_in( final quantity="ThermodynamicTemperature", final unit="K", displayUnit="degC") if (TDryBulSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input dry bulb temperature"; //-------------------------------------------------------------- // Dew point temperature parameter Buildings.BoundaryConditions.Types.DataSource TDewPoiSou=Buildings.BoundaryConditions.Types.DataSource.File "Dew point temperature"; parameter Modelica.SIunits.Temperature TDewPoi(displayUnit="degC") = 283.15 "Dew point temperature (used if TDewPoi=Parameter)"; Modelica.Blocks.Interfaces.RealInput TDewPoi_in( final quantity="ThermodynamicTemperature", final unit="K", displayUnit="degC") if (TDewPoiSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input dew point temperature"; //-------------------------------------------------------------- // Black body sky temperature parameter Buildings.BoundaryConditions.Types.DataSource TBlaSkySou=Buildings.BoundaryConditions.Types.DataSource.File "Black-body sky temperature"; parameter Modelica.SIunits.Temperature TBlaSky=273.15 "Black-body sky temperature (used if TBlaSkySou=Parameter)"; Modelica.Blocks.Interfaces.RealInput TBlaSky_in( final quantity="ThermodynamicTemperature", displayUnit="degC", final unit="K") if (TBlaSkySou == Buildings.BoundaryConditions.Types.DataSource.Input) "Black-body sky temperature"; //-------------------------------------------------------------- // Relative humidity parameter Buildings.BoundaryConditions.Types.DataSource relHumSou=Buildings.BoundaryConditions.Types.DataSource.File "Relative humidity"; parameter Real relHum( min=0, max=1, unit="1") = 0.5 "Relative humidity (used if relHum=Parameter)"; Modelica.Blocks.Interfaces.RealInput relHum_in( min=0, max=1, unit="1") if (relHumSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input relative humidity"; //-------------------------------------------------------------- // Wind speed parameter Buildings.BoundaryConditions.Types.DataSource winSpeSou=Buildings.BoundaryConditions.Types.DataSource.File "Wind speed"; parameter Modelica.SIunits.Velocity winSpe(min=0) = 1 "Wind speed (used if winSpe=Parameter)"; Modelica.Blocks.Interfaces.RealInput winSpe_in( final quantity="Velocity", final unit="m/s", min=0) if (winSpeSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input wind speed"; //-------------------------------------------------------------- // Wind direction parameter Buildings.BoundaryConditions.Types.DataSource winDirSou=Buildings.BoundaryConditions.Types.DataSource.File "Wind direction"; parameter Modelica.SIunits.Angle winDir=1.0 "Wind direction (used if winDir=Parameter)"; Modelica.Blocks.Interfaces.RealInput winDir_in( final quantity="Angle", final unit="rad", displayUnit="deg") if (winDirSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input wind direction"; //-------------------------------------------------------------- // Infrared horizontal radiation parameter Buildings.BoundaryConditions.Types.DataSource HInfHorSou=Buildings.BoundaryConditions.Types.DataSource.File "Infrared horizontal radiation"; parameter Modelica.SIunits.HeatFlux HInfHor=0.0 "Infrared horizontal radiation (used if HInfHorSou=Parameter)"; Modelica.Blocks.Interfaces.RealInput HInfHor_in( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") if (HInfHorSou == Buildings.BoundaryConditions.Types.DataSource.Input) "Input infrared horizontal radiation"; parameter Buildings.BoundaryConditions.Types.RadiationDataSource HSou=Buildings.BoundaryConditions.Types.RadiationDataSource.File "Global, diffuse, and direct normal radiation"; //-------------------------------------------------------------- // Global horizontal radiation Modelica.Blocks.Interfaces.RealInput HGloHor_in( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") if (HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HGloHor_HDifHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HGloHor) "Input global horizontal radiation"; //-------------------------------------------------------------- // Diffuse horizontal radiation Modelica.Blocks.Interfaces.RealInput HDifHor_in( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") if (HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HGloHor_HDifHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HDifHor) "Input diffuse horizontal radiation"; //-------------------------------------------------------------- // Direct normal radiation Modelica.Blocks.Interfaces.RealInput HDirNor_in(final quantity="RadiantEnergyFluenceRate", final unit="W/m2") if (HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HDifHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HGloHor) "Input direct normal radiation"; //-------------------------------------------------------------- parameter String filNam="" "Name of weather data file"; final parameter Modelica.SIunits.Angle lon(displayUnit="deg")= Buildings.BoundaryConditions.WeatherData.BaseClasses.getLongitudeTMY3( filNam) "Longitude"; final parameter Modelica.SIunits.Angle lat(displayUnit="deg")= Buildings.BoundaryConditions.WeatherData.BaseClasses.getLatitudeTMY3( filNam) "Latitude"; final parameter Modelica.SIunits.Time timZon(displayUnit="h")= Buildings.BoundaryConditions.WeatherData.BaseClasses.getTimeZoneTMY3(filNam) "Time zone"; Bus weaBus "Weather data bus"; parameter Buildings.BoundaryConditions.Types.SkyTemperatureCalculation calTSky=Buildings.BoundaryConditions.Types.SkyTemperatureCalculation.TemperaturesAndSkyCover "Computation of black-body sky temperature"; constant Real epsCos = 1e-6 "Small value to avoid division by 0"; protected Modelica.Blocks.Tables.CombiTable1Ds datRea( final tableOnFile=true, final tableName="tab1", final fileName=Buildings.BoundaryConditions.WeatherData.BaseClasses.getAbsolutePath(filNam), final smoothness=Modelica.Blocks.Types.Smoothness.ContinuousDerivative, final columns={2,3,4,5,6,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27, 28,29,30}) "Data reader"; Buildings.BoundaryConditions.WeatherData.BaseClasses.CheckTemperature cheTemDryBul "Check dry bulb temperature "; Buildings.BoundaryConditions.WeatherData.BaseClasses.CheckTemperature cheTemDewPoi "Check dew point temperature"; Buildings.BoundaryConditions.WeatherData.BaseClasses.ConvertRelativeHumidity conRelHum "Convert the relative humidity from percentage to [0, 1] "; BaseClasses.CheckPressure chePre "Check the air pressure"; BaseClasses.CheckSkyCover cheTotSkyCov "Check the total sky cover"; BaseClasses.CheckSkyCover cheOpaSkyCov "Check the opaque sky cover"; BaseClasses.CheckRadiation cheGloHorRad "Check the global horizontal radiation"; BaseClasses.CheckRadiation cheDifHorRad "Check the diffuse horizontal radiation"; BaseClasses.CheckRadiation cheDirNorRad "Check the direct normal radiation"; BaseClasses.CheckCeilingHeight cheCeiHei "Check the ceiling height"; BaseClasses.CheckWindSpeed cheWinSpe "Check the wind speed"; BaseClasses.CheckRadiation cheHorRad "Check the horizontal radiation"; BaseClasses.CheckWindDirection cheWinDir "Check the wind direction"; SkyTemperature.BlackBody TBlaSkyCom(final calTSky=calTSky) if not (TBlaSkySou == Buildings.BoundaryConditions.Types.DataSource.Parameter or TBlaSkySou == Buildings.BoundaryConditions.Types.DataSource.Input) "Computation of the sky black-body temperature"; Utilities.Time.ModelTime modTim "Model time"; Modelica.Blocks.Math.Add add "Add 30 minutes to time to shift weather data reader"; Modelica.Blocks.Sources.Constant con30mins(final k=1800) "Constant used to shift weather data reader"; Buildings.BoundaryConditions.WeatherData.BaseClasses.LocalCivilTime locTim( final lon=lon, final timZon=timZon) "Local civil time"; Modelica.Blocks.Tables.CombiTable1Ds datRea1( final tableOnFile=true, final tableName="tab1", final fileName=Buildings.BoundaryConditions.WeatherData.BaseClasses.getAbsolutePath(filNam), final smoothness=Modelica.Blocks.Types.Smoothness.ContinuousDerivative, final columns=8:11) "Data reader"; Buildings.BoundaryConditions.WeatherData.BaseClasses.ConvertTime conTim1 "Convert simulation time to calendar time"; BaseClasses.ConvertTime conTim "Convert simulation time to calendar time"; BaseClasses.EquationOfTime eqnTim "Equation of time"; BaseClasses.SolarTime solTim "Solar time"; // Conditional connectors Modelica.Blocks.Interfaces.RealInput pAtm_in_internal( final quantity="Pressure", final unit="Pa", displayUnit="bar") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput ceiHei_in_internal( final quantity="Height", final unit="m", displayUnit="m") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput totSkyCov_in_internal( final quantity="1", min=0, max=1) "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput opaSkyCov_in_internal( final quantity="1", min=0, max=1) "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput TDryBul_in_internal( final quantity="ThermodynamicTemperature", final unit="K", displayUnit="degC") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput TDewPoi_in_internal( final quantity="ThermodynamicTemperature", final unit="K", displayUnit="degC") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput TBlaSky_in_internal( final quantity="ThermodynamicTemperature", final unit="K", displayUnit="degC") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput relHum_in_internal( final quantity="1", min=0, max=1) "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput winSpe_in_internal( final quantity="Velocity", final unit="m/s") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput winDir_in_internal( final quantity="Angle", final unit="rad", displayUnit="deg") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput HGloHor_in_internal( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput HDifHor_in_internal( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput HDirNor_in_internal( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") "Needed to connect to conditional connector"; Modelica.Blocks.Interfaces.RealInput HInfHor_in_internal( final quantity="RadiantEnergyFluenceRate", final unit="W/m2") "Needed to connect to conditional connector"; Modelica.Blocks.Math.UnitConversions.From_deg conWinDir "Convert the wind direction unit from [deg] to [rad]"; Modelica.Blocks.Math.UnitConversions.From_degC conTDryBul; BaseClasses.ConvertRadiation conHorRad; Modelica.Blocks.Math.UnitConversions.From_degC conTDewPoi "Convert the dew point temperature form [degC] to [K]"; BaseClasses.ConvertRadiation conDirNorRad; BaseClasses.ConvertRadiation conGloHorRad; BaseClasses.ConvertRadiation conDifHorRad; BaseClasses.CheckRelativeHumidity cheRelHum; SolarGeometry.BaseClasses.AltitudeAngle altAng "Solar altitude angle"; SolarGeometry.BaseClasses.ZenithAngle zenAng( final lat = lat) "Zenith angle"; SolarGeometry.BaseClasses.Declination decAng "Declination angle"; SolarGeometry.BaseClasses.SolarHourAngle solHouAng; Modelica.Blocks.Sources.Constant latitude(final k=lat) "Latitude"; Modelica.Blocks.Sources.Constant longitude(final k=lon) "Longitude"; //--------------------------------------------------------------------------- // Optional instanciation of a block that computes the wet bulb temperature. // This block may be needed for evaporative cooling towers. // By default, it is enabled. This introduces a nonlinear equation, but // we have not observed an increase in computing time because of this equation. Buildings.Utilities.Psychrometrics.TWetBul_TDryBulPhi tWetBul_TDryBulXi( redeclare package Medium = Buildings.Media.Air, TDryBul(displayUnit="degC")) if computeWetBulbTemperature; //--------------------------------------------------------------------------- // Conversion blocks for sky cover Modelica.Blocks.Math.Gain conTotSkyCov(final k=0.1) if totSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.File "Convert sky cover from [0...10] to [0...1]"; Modelica.Blocks.Math.Gain conOpaSkyCov(final k=0.1) if opaSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.File "Convert sky cover from [0...10] to [0...1]"; equation //--------------------------------------------------------------------------- // Select atmospheric pressure connector if pAtmSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then pAtm_in_internal = pAtm; elseif pAtmSou == Buildings.BoundaryConditions.Types.DataSource.File then connect(datRea.y[4], pAtm_in_internal); else connect(pAtm_in, pAtm_in_internal); end if; connect(pAtm_in_internal, chePre.PIn); //--------------------------------------------------------------------------- // Select ceiling height connector if ceiHeiSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then ceiHei_in_internal = ceiHei; elseif ceiHeiSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(ceiHei_in, ceiHei_in_internal); else connect(datRea.y[16], ceiHei_in_internal); end if; connect(ceiHei_in_internal, cheCeiHei.ceiHeiIn); //--------------------------------------------------------------------------- // Select total sky cover connector if totSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then totSkyCov_in_internal = totSkyCov; elseif totSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(totSkyCov_in, totSkyCov_in_internal); else connect(conTotSkyCov.u, datRea.y[13]); connect(conTotSkyCov.y, totSkyCov_in_internal); end if; connect(totSkyCov_in_internal, cheTotSkyCov.nIn); //--------------------------------------------------------------------------- // Select opaque sky cover connector if opaSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then opaSkyCov_in_internal = opaSkyCov; elseif opaSkyCovSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(opaSkyCov_in, opaSkyCov_in_internal); else connect(conOpaSkyCov.u, datRea.y[14]); connect(conOpaSkyCov.y, opaSkyCov_in_internal); end if; connect(opaSkyCov_in_internal, cheOpaSkyCov.nIn); //--------------------------------------------------------------------------- // Select dew point temperature connector if TDewPoiSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then TDewPoi_in_internal = TDewPoi; elseif TDewPoiSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(TDewPoi_in, TDewPoi_in_internal); else connect(conTDewPoi.y, TDewPoi_in_internal); end if; connect(TDewPoi_in_internal, cheTemDewPoi.TIn); //--------------------------------------------------------------------------- // Select dry bulb temperature connector if TDryBulSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then TDryBul_in_internal = TDryBul; elseif TDryBulSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(TDryBul_in, TDryBul_in_internal); else connect(conTDryBul.y, TDryBul_in_internal); end if; connect(TDryBul_in_internal, cheTemDryBul.TIn); //--------------------------------------------------------------------------- // Select sky black body temperature connector if TBlaSkySou == Buildings.BoundaryConditions.Types.DataSource.Parameter then TBlaSky_in_internal = TBlaSky; elseif TBlaSkySou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(TBlaSky_in, TBlaSky_in_internal); else connect(TBlaSkyCom.TBlaSky, TBlaSky_in_internal); end if; connect(TBlaSky_in_internal, weaBus.TBlaSky); //--------------------------------------------------------------------------- // Select relative humidity connector if relHumSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then relHum_in_internal = relHum; elseif relHumSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(relHum_in, relHum_in_internal); else connect(conRelHum.relHumOut, relHum_in_internal); end if; connect(relHum_in_internal, cheRelHum.relHumIn); //--------------------------------------------------------------------------- // Select wind speed connector if winSpeSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then winSpe_in_internal = winSpe; elseif winSpeSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(winSpe_in, winSpe_in_internal); else connect(datRea.y[12], winSpe_in_internal); end if; connect(winSpe_in_internal, cheWinSpe.winSpeIn); //--------------------------------------------------------------------------- // Select wind direction connector if winDirSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then winDir_in_internal = winDir; elseif winDirSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(winDir_in, winDir_in_internal); else connect(conWinDir.y, winDir_in_internal); end if; connect(winDir_in_internal, cheWinDir.nIn); //--------------------------------------------------------------------------- // Select global horizontal radiation connector if HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HGloHor_HDifHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HGloHor then connect(HGloHor_in, HGloHor_in_internal) "Get HGloHor using user input file"; elseif HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HDifHor then HDirNor_in_internal*cos(zenAng.zen)+HDifHor_in_internal = HGloHor_in_internal "Calculate the HGloHor using HDirNor and HDifHor according to (A.4.14) and (A.4.15)"; else connect(conGloHorRad.HOut, HGloHor_in_internal) "Get HGloHor using weather data file"; end if; connect(HGloHor_in_internal, cheGloHorRad.HIn); //--------------------------------------------------------------------------- // Select diffuse horizontal radiation connector if HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HGloHor_HDifHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HDifHor then connect(HDifHor_in, HDifHor_in_internal) "Get HDifHor using user input file"; elseif HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HGloHor then HGloHor_in_internal - HDirNor_in_internal*cos(zenAng.zen) = HDifHor_in_internal "Calculate the HGloHor using HDirNor and HDifHor according to (A.4.14) and (A.4.15)"; else connect(conDifHorRad.HOut, HDifHor_in_internal) "Get HDifHor using weather data file"; end if; connect(HDifHor_in_internal, cheDifHorRad.HIn); //--------------------------------------------------------------------------- // Select direct normal radiation connector if HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HGloHor or HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HDirNor_HDifHor then connect(HDirNor_in, HDirNor_in_internal) "Get HDirNor using user input file"; elseif HSou == Buildings.BoundaryConditions.Types.RadiationDataSource.Input_HGloHor_HDifHor then (HGloHor_in_internal -HDifHor_in_internal)/Buildings.Utilities.Math.Functions.smoothMax(x1=cos(zenAng.zen), x2=epsCos, deltaX=0.1*epsCos) = HDirNor_in_internal "Calculate the HDirNor using HGloHor and HDifHor according to (A.4.14) and (A.4.15)"; else connect(conDirNorRad.HOut, HDirNor_in_internal) "Get HDirNor using weather data file"; end if; connect(HDirNor_in_internal, cheDirNorRad.HIn); //--------------------------------------------------------------------------- // Select infrared radiation connector if HInfHorSou == Buildings.BoundaryConditions.Types.DataSource.Parameter then HInfHor_in_internal = HInfHor; elseif HInfHorSou == Buildings.BoundaryConditions.Types.DataSource.Input then connect(HInfHor_in, HInfHor_in_internal); else connect(conHorRad.HOut, HInfHor_in_internal); end if; connect(HInfHor_in_internal, cheHorRad.HIn); connect(chePre.POut, weaBus.pAtm); connect(cheTotSkyCov.nOut, weaBus.nTot); connect(cheOpaSkyCov.nOut, weaBus.nOpa); connect(cheGloHorRad.HOut, weaBus.HGloHor); connect(cheDifHorRad.HOut, weaBus.HDifHor); connect(cheDirNorRad.HOut, weaBus.HDirNor); connect(cheCeiHei.ceiHeiOut, weaBus.celHei); connect(cheWinSpe.winSpeOut, weaBus.winSpe); connect(cheHorRad.HOut, weaBus.radHorIR); connect(cheWinDir.nOut, weaBus.winDir); connect(cheOpaSkyCov.nOut,TBlaSkyCom. nOpa); connect(cheHorRad.HOut,TBlaSkyCom. radHorIR); connect(TBlaSkyCom.TBlaSky, weaBus.TBlaSky); connect(modTim.y, weaBus.cloTim); connect(modTim.y, add.u2); connect(con30mins.y, add.u1); connect(add.y, conTim1.modTim); connect(conTim1.calTim, datRea1.u); connect(modTim.y, locTim.cloTim); connect(modTim.y, conTim.modTim); connect(conTim.calTim, datRea.u); connect(modTim.y, eqnTim.nDay); connect(eqnTim.eqnTim, solTim.equTim); connect(locTim.locTim, solTim.locTim); connect(solTim.solTim, weaBus.solTim); connect(datRea.y[11], conWinDir.u); connect(datRea1.y[1], conHorRad.HIn); connect(cheTemDryBul.TOut,TBlaSkyCom. TDryBul); connect(datRea.y[1], conTDryBul.u); connect(datRea.y[2], conTDewPoi.u); connect(cheTemDewPoi.TOut, weaBus.TDewPoi); connect(TBlaSkyCom.TDewPoi, cheTemDewPoi.TOut); connect(datRea1.y[3], conDirNorRad.HIn); connect(datRea1.y[2], conGloHorRad.HIn); connect(datRea1.y[4], conDifHorRad.HIn); connect(conRelHum.relHumIn, datRea.y[3]); connect(cheRelHum.relHumOut, weaBus.relHum); connect(cheTemDryBul.TOut, weaBus.TDryBul); connect(decAng.decAng, zenAng.decAng); connect(solHouAng.solHouAng, zenAng.solHouAng); connect(solHouAng.solTim, solTim.solTim); connect(decAng.nDay, modTim.y); connect(zenAng.zen, altAng.zen); // Connectors for wet bulb temperature. // These are removed if computeWetBulbTemperature = false connect(chePre.POut, tWetBul_TDryBulXi.p); connect(tWetBul_TDryBulXi.TWetBul, weaBus.TWetBul); connect(cheTemDryBul.TOut, tWetBul_TDryBulXi.TDryBul); connect(cheRelHum.relHumOut, tWetBul_TDryBulXi.phi); connect(altAng.alt, weaBus.solAlt); connect(zenAng.zen, weaBus.solZen); connect(decAng.decAng, weaBus.solDec); connect(solHouAng.solHouAng, weaBus.solHouAng); connect(longitude.y, weaBus.lon); connect(latitude.y, weaBus.lat); end ReaderTMY3;