Design day weather

Author: chipbarnaby

src/ANCREC.CPP

@@ -93,8 +93,8 @@ record::record( BP _b, TI i, SI noZ/*=0*/)  	// construct record i of basAnc b,
 	   				b->n++ : additional record in static array.
 	   CAUTION can be called from statSetup (til eliminated). */
 
-	_b->validate("record::record");			// abort if base bad -- we depend on it for size
-	if (!noZ)						// can be suppressed for data-init (static) records
+	_b->validate("record::record");		// abort if base bad -- we depend on it for size
+	if (!noZ)							// can be suppressed for data-init (static) records
 		memset( (char *)this+SZVFTP, 0, _b->eSz-SZVFTP);	/* zero all members of DERIVED CLASS record
 						           except virtFcnTblPtr in 1st 2 or 4 bytes (ancrec.h define) */
 	b = _b;
@@ -685,9 +685,9 @@ SI FC basAnc::ancNext( USI &an, BP * _b)		// (registered) record anchor iterator
 //---------------------------------------------------------------------------------------------------------------------------
 RC basAnc::validate(	// validate an anchor: check self-consistency of anchor and its entry 0
 
-	char *fcnName,		// calling fcn name, for err msgs.
-	int erOp /*=ABT*/,	// reporting control. ABT-->PABT, WRN-->PWRN here.
-	SI noStat /*=0*/ )	// 0: caller's action OK on static-storage anchor. nz: not ok.
+	const char* fcnName,	// calling fcn name, for err msgs.
+	int erOp /*=ABT*/,		// reporting control. ABT-->PABT, WRN-->PWRN here.
+	SI noStat /*=0*/ )		// 0: caller's action OK on static-storage anchor. nz: not ok.
 
 // if error, issues message, returns non-RCOK if not ABT
 {

src/ANCREC.H

@@ -108,7 +108,7 @@ class basAnc    	// base class for record anchors: basAnc<recordName>
     static BP FC anc4n( USI an, int erOp=ABT);   				// access anchor by anchor #
     static RC FC findAnchorByNm( char *what, BP *b);
     static SI FC ancNext( USI &an, BP *_b);						// iterate anchors
-    RC validate( char *fcnName, int erOp=ABT, SI noStat=0);				// check for valid anchor
+    RC validate( const char* fcnName, int erOp=ABT, SI noStat=0);		// check for valid anchor
     RC findRecByNm1( const char* _name, TI *_i, record **_r);    		// find record by 1st match on name
     RC findRecByNmU( const char* _name, TI *_i, record **_r);  			// find record by unique name match
     RC findRecByNmO( const char* _name, TI ownTi, TI *_i, record **_r);	// find record by name and owner subscript

src/CGWTHR.CPP

@@ -54,24 +54,33 @@ void FC cgWthrClean( 		// cgwthr overall init/cleanup routine
 		cgWfDone();		// close any open weather file
 
 // init weather file info record
-	Wfile.wf_Init();		// mbr function in wfpak.cpp sets record to say no file open
-	Wthr.wd_Init();			// mbr function in wfpak.cpp clears record
-	WthrNxHr.wd_Init();		// mbr function in wfpak.cpp clears record
+	Wfile.wf_Init();		// sets record to say no file open
+	Wthr.wd_Init();			// clears record
+	WthrNxHr.wd_Init();		// clears record
 
 	// add any other init/cleanup found necessary or desirable
 
 }		// cgWthrClean
 //---------------------------------------------------------------------------
-RC TOPRAT::tp_WfInit()  	// Initialize weather file.  Displays any error messages.
-// also uses: Top.tp_wfName [,.isDT]
+RC TOPRAT::tp_WthrInit()  	// Initialize weather data/ display any error messages.
+// also uses: tp_wfName, tp_TDVfName, [,.isDT]
 // returns non-RCOK if error, message already issued.
 {
 	RC rc = Wfile.wf_Open( tp_wfName, tp_TDVfName);	// open wthr file and option TDV file
 													// init WFILE object, msg if error.
 	if (!rc)
-		rc = Wfile.wf_FillWDYEAR( WRN);
+		rc = Wfile.wf_FillWDYEAR( WRN);		// read/cache entire weather file
+											//   TODO: handle partial years
+
+	// design conditions
+	//  set any derived values
+	//  WHY here: location dependence, call after locinit
+	DESCOND* pDC;
+	RLUP( DcR, pDC)
+		rc |= pDC->dc_RunInit();
+
 	return rc;
-}		// TOPRAT::tp_WfInit
+}		// TOPRAT::tp_WthrInit
 //---------------------------------------------------------------------------
 RC TOPRAT::tp_WthrBegHour()		// start-hour weather stuff: read file, set up public Wthr and Top members.
 // tp_WthrBegSubhr must also be called, after this function.
@@ -432,7 +441,7 @@ RC WDHR::wd_WfReader( 	// read an hour's weather data and make adjustments
 			if (Top.tp_AuszWthrSource() == TOPRAT_COOLDSCOND)
 			{	// DSCOND design day: overwrite/adjust weather file values with generated
 				int iDC = Top.tp_coolDsCond[ Top.tp_dsDayI-1];
-				const DESCOND& DC = DCiB[ iDC];
+				const DESCOND& DC = DcR[ iDC];
 				wd_FillFromDESCOND( DC, iHrST);
 			}
 		}
@@ -517,7 +526,7 @@ RC WDHR::wd_WfReader( 	// read an hour's weather data and make adjustments
 			rc = pWF->wf_Read( jDayST, iHrST, this, WRN|wfOp);	// Read hour's data from weather file
 			if (Top.tp_AuszWthrSource() == TOPRAT_COOLDSCOND)
 			{	int iDC = Top.tp_coolDsCond[ Top.tp_dsDayI-1];
-				const DESCOND& DC = DCiB[ iDC];
+				const DESCOND& DC = DcR[ iDC];
 				DC.dc_GenerateTemps( iHrST, this);
 
 			}
@@ -611,11 +620,12 @@ void WDHR::wd_FillFromDESCOND(		// overwrite/adjust hourly data for design condi
 	wd_taDbAvg14 += dbAvgDiff / 14.f;
 	wd_taDbAvg31 += dbAvgDiff / 31.f;
 
-
 	wd_wndSpd = dc.dc_wndSpd;
 	wd_wndDir = 0.f;
 
-	if (wd_sunup)
+	// derive irradiance
+	//   note: tauB/tauD both 0 -> no solar
+	if (wd_sunup && dc.dc_tauB > 0.f && dc.dc_tauD > 0.f)
 		wd_glrad = slASHRAETauModel(
 			kPiOver2 - wd_slAlt,
 			dc.dc_tauB, dc.dc_tauD,
@@ -636,17 +646,107 @@ void WDHR::wd_FillFromDESCOND(		// overwrite/adjust hourly data for design condi
 /////////////////////////////////////////////////////////////////////////////
 RC DESCOND::dc_CkF()		// check after input
 {	RC rc = RCOK;
-	if (dc_MCWB > dc_DB)
-		rc |= oer( "dcMCWB (%g) must be <= dcDB (%g)", dc_MCWB, dc_DB);
 	return rc;
 }		// DESCOND::dc_CkF
 //-----------------------------------------------------------------------------
-DOY DESCOND::dc_GetDOY() const
+#if defined( _DEBUG)
+static void DCTauTest( int doy)
 {
-	return dc_day;
+float ebnlist[] = { 150.f, 200.f, 250.f, 300.f, -1.f };
+float edflist[] = { 20.f, 30.f, 40.f };
+
+	DESCOND dc( &DcR, 1);
+	dc.dc_doy = doy;
+	strcpy( dc.name, "Test");
+
+	for (int ib=0; ebnlist[ ib]>=0.f; ib++)
+	{	dc.dc_ebnSlrNoon = ebnlist[ ib];
+		for (int id=0; edflist[ id]>=0.f; id++)
+		{	dc.dc_edhSlrNoon = edflist[ id];
+			dc.dc_CheckFixSolar( 0);
+		}
+	}
+}		// DCTauTest
+#endif
+//------------------------------------------------------------------------------
+RC DESCOND::dc_RunInit()		// init for run
+// call after all input-time expressions have been evaluated
+// call after locInit (re latitude dependency)
+{
+#if defined( _DEBUG)
+static int tested = 0;
+	if (!tested)
+	{	DCTauTest( dc_doy);
+		tested++;
+	}
+#endif
 
-}	// DESCOND::dc_GetDOY
+	RC rc = RCOK;
+
+	if (dc_MCWB > dc_DB)
+		rc |= oer( "dcMCWB (%g) must be <= dcDB (%g)", dc_MCWB, dc_DB);
 
+	// if any solar values given, check and cross-derive
+	//   else leave all 0 = no solar
+	if (dc_tauB + dc_tauD + dc_ebnSlrNoon + dc_edhSlrNoon > 0.f)
+	{	int nTau = IsSetCount( DESCOND_TAUB, DESCOND_TAUD, 0);
+		int nE   = IsSetCount( DESCOND_EBNSLRNOON, DESCOND_EDHSLRNOON, 0);
+		if (nTau == 0)
+			rc |= dc_CheckFixSolar( 0);		// derive tauB / tauD
+		else if (nE == 0)
+			rc |= dc_CheckFixSolar( 1);		// derive Ebn / Edh
+		else
+			rc |= oer( "must give either dcTauB/dcTauD and"
+				 " not dcEbnSlrNoon/dcEdhSlrNoon or vice-versa");
+	}
+	else
+		oInfo( "no irradiance due to dcTauB/dcTauD/dcEbnSlrNoon/dcEdhSlrNoon all 0");
+
+	return rc;
+}	// DESCOND::dc_RunInit
+//-----------------------------------------------------------------------------
+RC DESCOND::dc_CheckFixSolar(
+	int options)	// option bits
+					//  0: derive tau from solar noon
+					//  1: derive solar noon from tau
+{
+	RC rc = RCOK;
+
+	// solar noon angles
+	slday( dc_doy, SLTMSOLAR, 1);
+	float sunAlt, sunAzm;
+	slaltazm( 12.f, &sunAlt, &sunAzm);
+	float sunZen = kPiOver2 - sunAlt;
+
+	if (options & 1)
+	{	rc |=   limitCheck( DESCOND_TAUB, .08, 1.)
+		      | limitCheck( DESCOND_TAUD, 1.2, 3.2);
+		if (!rc)
+			slASHRAETauModel( sunZen, dc_tauB, dc_tauD,
+						dc_ebnSlrNoon, dc_edhSlrNoon);
+	}
+	else	
+	{	rc |=   limitCheck( DESCOND_EBNSLRNOON, 0., 370.)
+		      | limitCheck( DESCOND_EDHSLRNOON, 0., 110.);
+		if (!rc)
+		{	bool ret = slASHRAETauModelInv( sunZen, dc_ebnSlrNoon, dc_edhSlrNoon,
+					dc_tauB, dc_tauD);
+			if (!ret)
+				rc = oer( "failed to derive dcTauB and dcTauD."
+					"  Check dcEbnSlrNoon and dcEdhSlrNoon input values.");
+#if defined( _DEBUG)
+			else
+			{	float ebn, edh;
+				float egl = slASHRAETauModel( sunZen, dc_tauB, dc_tauD, ebn, edh);
+				if (frDiff( ebn, dc_ebnSlrNoon) > .01f
+				  || frDiff( edh, dc_edhSlrNoon) > .01f)
+					oWarn( "ebn/edh mismatch");
+			}
+#endif
+		}
+	}
+	return rc;
+}	// DESCOND::dc_CheckFixSolar
 //-----------------------------------------------------------------------------
 void DESCOND::dc_GenerateTemps(
 	int iHr,				// hour of day (0 - 23)

src/CNAUSZ.CPP

@@ -105,7 +105,7 @@ const char* TOPRAT::tp_AuszDoing() const	// current autosizing activity
 		else
 		{	int auszDoy;
 			if (iAWS == TOPRAT_COOLDSCOND)
-				auszDoy = DCiB[ tp_coolDsCond[ iDs]].dc_GetDOY();
+				auszDoy = DcR[ tp_coolDsCond[ iDs]].dc_GetDOY();
 			else if (iAWS == TOPRAT_COOLDSDAY)
 				auszDoy = tp_coolDsDay[ iDs];
 			if (auszDoy > 0)
@@ -232,7 +232,7 @@ x	tp_dsDay = 1 + (tp_dsDayI != 0);			// 1st desDay is heating (1), others coolin
 		auszMon = tp_date.month;
 	}
 	else if (tp_AuszWthrSource() == TOPRAT_COOLDSCOND)
-	{	const DESCOND* pDC = DCiB.GetAtSafe( tp_coolDsCond[ tp_dsDayI-1]);
+	{	const DESCOND* pDC = DcR.GetAtSafe( tp_coolDsCond[ tp_dsDayI-1]);
 		jDay = pDC->dc_GetDOY();
 		tddyi( tp_date, jDay);
 		auszMon = tp_date.month;

src/CNCULT.CPP

@@ -117,7 +117,7 @@ makAncRI2(XiB,"export");		// EXPORT input info
 makAncCOL(RcoliB);				// REPORTCOL input info for UDT reports
 makAncCOL2(XcoliB,"exportCol"); // EXPORTCOL input info for UDT exports
 makAncIMPF(ImpfiB);				// IMPORTFILE input info, 2-94
-// holidays
+
 makAncHDAY(HdayiB);				// HDAY input info
 makAncDESCOND(DCiB);			// DESCOND input info
 
@@ -2433,15 +2433,19 @@ static CULT dcT[] =
 //----------        ----- ---------------  -------  -- ------  ------  ----  ----------  ---- ----
 CULT( "*",          STAR,  0,               0,       0, 0,      0,     0,      0.f,    N,   dcStarCkf),
 
-CULT( "dcDay",	    DAT,  DESCOND_DAY,	    RQD,      0, VEOI,   TYDOY,  0,    200,    N,   N),
+CULT( "dcDay",	    DAT,  DESCOND_DOY,	    RQD,      0, VEOI,   TYDOY,  0,    200,    N,   N),
 
 CULT( "dcDB",       DAT,  DESCOND_DB,       RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
 CULT( "dcMCDBR",    DAT,  DESCOND_MCDBR,    RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
 CULT( "dcMCWB",     DAT,  DESCOND_MCWB,     RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
 CULT( "dcMCWBR",	DAT,  DESCOND_MCWBR,    RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
 CULT( "dcWindSpeed",DAT,  DESCOND_WNDSPD,   RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
-CULT( "dcTauB",     DAT,  DESCOND_TAUB,     RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
-CULT( "dcTauD",     DAT,  DESCOND_TAUD,     RQD,      0, VEOI,   TYFL,   0,    0.f,       N,   N),
+// Taub/Taud and Ebn/Ebh interact, see dc_TopDC
+CULT( "dcTauB",     DAT,  DESCOND_TAUB,      0,       0, VEOI,   TYFL,   0,    0.f,       N,   N),
+CULT( "dcTauD",     DAT,  DESCOND_TAUD,      0,       0, VEOI,   TYFL,   0,    0.f,       N,   N),
+CULT( "dcEbnSlrNoon",DAT, DESCOND_EBNSLRNOON,0,       0, VEOI,   TYFL,   0,    0.f,       N,   N),
+CULT( "dcEdhSlrNoon",DAT, DESCOND_EDHSLRNOON,0,       0, VEOI,   TYFL,   0,    0.f,       N,   N),
+
 
 CULT( "endDesCond", ENDER, 0,           0,        0, 0,      0,      0,    0.f,       N,   N),
 CULT()

src/CNCULT2.CPP

@@ -815,9 +815,8 @@ const char* TOPRAT::When( IVLCH ivl) const		// date / time doc for error message
 	return s;
 }			// TOPRAT::When
 //-----------------------------------------------------------------------------
-LOCAL RC topDC()		// finalize DESCONDs
-{	RC rc = RCOK;
-	// input records used at runtime, nothing to do
+LOCAL RC topDC()		// copy input DESCONDs to run records
+{	RC rc = DcR.RunDup( DCiB);
 	return rc;
 }	// topDC
 //-----------------------------------------------------------------------------

src/CNGUTS.CPP

@@ -208,7 +208,7 @@ void FC runFazDataFree()
 
 	// Note: entries are put into most of these in input code in cncult2.cpp:topCkf and its subfcns.
 
-	ZrB.free(); 	// Zones runtime info.  basAnc::free is in lib\ancrec.cpp.
+	ZrB.free(); 	// Zones runtime info
 	ZnresB.free();	// zone simulation results
 	RSysResR.free();	// RSYS results
 	IzxR.free();	// interzone transfers
@@ -925,9 +925,10 @@ RC FC cgRddInit(	// Perform initialization common to main simulation run and eac
 						// inits Locsolar (comments below)
 						// EF: return RCBAD on any non-RCOK return.
 
-// Set up weather file
-// do each phase, allowing name to change. Messages errors.
-	EF( Top.tp_WfInit() )	// also uses Top .tp_wfName, .skyModel [,.isDT].  EF: return RCBAD if error.
+// Set up weather file, TDV file, and DESCONDs
+//   Do each phase, allowing file or location changes.
+//   Messages errors.
+	EF( Top.tp_WthrInit() )	// also uses Top .tp_wfName, .skyModel [,.isDT]
 
 // Open import files used in this run, if any
 	EF( impfStart() )					// impf.cpp. related code: cncult4.cpp, cuparse.cpp.
@@ -2503,7 +2504,7 @@ LOCAL void FC NEAR binResInit( int isAusz)	// initialize & open binary results (
 					/* ***** TEMPORARY weather file info for initial DL test graphs ***** 9-4-94
 					   Make names match values passed to setHourInfo() in doIvlAccum()!
 					   Expect will depend on wf type.
-					   Coordinate via tp_WfInit (xxx called in this fcn about 100 lines up) ***** */
+					   Coordinate via tp_WthrInit(xxx called in this fcn about 100 lines up) ***** */
 #ifdef TWOWV // bfr.h 9-94
 					2, "Est Direct", "Est Diffuse" ) )	// # dl weather vbls and their names. Estimated from CEC wthr file for now.
 #else
@@ -2555,15 +2556,15 @@ static SLLOCDAT* NEAR Locsolar = NULL;	/* Ptr to solar info structure (slpak.h)
 				           Daily info set by slpak:slday from cgwthr.cpp and cgsolar.cpp at least, 1-95. */
 //------------------------------------------------------------------------
 LOCAL RC FC NEAR locInit()		// location-initializer for cne
-// inits lib\slpak.cpp using location info in Top members (gotten weather file and/or user input in cncultx.cpp)
+// inits slpak.cpp using location info in Top members (gotten weather file and/or user input in cncultx.cpp)
 // returns RCOK if ok, else message already issued.
 {
-/* init lib\slpak.cpp: alloc SLLOCDAT (slpak.h) and init for location.
-   More of Locsolar is set each day by slday(). It is used by most other slpak calls.
-   slpak remembers its location internally. */
+// init slpak.cpp: alloc SLLOCDAT (slpak.h) and init for location.
+//   More of Locsolar is set each day by slday(). It is used by most other slpak calls.
+//   slpak remembers its location internally.
 
 	Locsolar = slinit( RAD( Top.latitude),	// latitude (to radians)
-					   RAD( Top.longitude), 	// longitude
+					   RAD( Top.longitude), // longitude
 					   Top.timeZone, 		// time zone
 					   Top.elevation );		// site altitude (ft)
 	return RCOK;

src/CNGUTS.H

@@ -113,8 +113,7 @@ extern anc<WFDATA> WthrR;		// basAnc for hour's weather data record. Runtime onl
 extern WFDATA Wthr;				// the one static WFDATA record containing unpacked & adjusted data for hour 1-94
 extern anc<WFDATA> WthrNxHrR;	// basAnc for next hour's weather data record ("weatherNextHour")
 extern WFDATA WthrNxHr;			// static record for next hour's unpacked & adjusted data for cgwthr.cpp read-ahead
-
-extern anc<DESCOND> DCiB;		// design conditions (input records only)
+extern anc<DESCOND> DcR;		// design conditions
 
 extern anc<ZNR> ZrB;			// zones runtime info
 extern anc<ZNRES> ZnresB;      	// Zones simulation results info