Skip to content

push fixes while debugging OF PR #3025 #3

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
ptrbortolotti wants to merge 1 commit into
base: main
Choose a base branch
from
Open

push fixes while debugging OF PR #3025 #3

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
61 changes: 31 additions & 30 deletions openfast_toolbox/fastfarm/FASTFarmCaseCreation.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -507,7 +507,7 @@ def _checkInputs(self):
mod_wake_str = ['','polar', 'curled', 'cartesian']
print(f'WARNING: One or more temporal or spatial resolution for low- and high-res domains were not given.')
print(f' Estimated values for {mod_wake_str[self.mod_wake]} wake model shown below.')
self._determine_resolutions_from_dummy_amrwind_grid()
# self._determine_resolutions_from_dummy_amrwind_grid()

# Check the temporal and spatial resolutions if provided
if self.dt_low != None and self.dt_high!= None:
Expand Down Expand Up @@ -1289,7 +1289,7 @@ def checkIfExists(f):
self.turbfilename = value

elif key == 'libdisconfilepath':
if not value.endswith('.so'):
if not value.endswith('.so') and not value.endswith('.dylib'):
raise ValueError(f'The libdiscon file should end in "*.so"')
if os.path.isabs(value):
self.libdisconfilepath = value
Expand Down Expand Up @@ -2148,41 +2148,42 @@ def FF_setup(self, outlistFF=None, **kwargs):


# Planes to save in FAST.Farm. We want the planes through the original farm, so let's get the position of the turbines at wdir=0
alignedTurbs = self.allCases.where(self.allCases['inflow_deg']==0, drop=True).isel(case=0)
if self.inflowStr == 'TurbSim':
# Turbine location in TurbSim reference frame
xWT = alignedTurbs['Tx'].values + self.xoffset_turbsOrigin2TSOrigin
yWT = alignedTurbs['Ty'].values + self.yoffset_turbsOrigin2TSOrigin
elif self.inflowStr == 'LES':
# Turbine location in LES reference frame
xWT = alignedTurbs['Tx'].values
yWT = alignedTurbs['Ty'].values

offset=10
planes_xy = [self.zhub+self.zbot]
planes_yz = np.unique(np.round(xWT+offset, 2))
planes_xz = np.unique(np.round(yWT, 2))

# Number of planes must be at most 9
self.planes_xy = planes_xy[0:9]
self.planes_yz = planes_yz[0:9]
self.planes_xz = planes_xz[0:9]
if self.allCases['inflow_deg']==0:
alignedTurbs = self.allCases.where(self.allCases['inflow_deg']==0, drop=True).isel(case=0)
if self.inflowStr == 'TurbSim':
# Turbine location in TurbSim reference frame
xWT = alignedTurbs['Tx'].values + self.xoffset_turbsOrigin2TSOrigin
yWT = alignedTurbs['Ty'].values + self.yoffset_turbsOrigin2TSOrigin
elif self.inflowStr == 'LES':
# Turbine location in LES reference frame
xWT = alignedTurbs['Tx'].values
yWT = alignedTurbs['Ty'].values

offset=10
planes_xy = [self.zhub+self.zbot]
planes_yz = np.unique(np.round(xWT+offset, 2))
planes_xz = np.unique(np.round(yWT, 2))

# Number of planes must be at most 9
self.planes_xy = planes_xy[0:9]
self.planes_yz = planes_yz[0:9]
self.planes_xz = planes_xz[0:9]

if self.inflowStr == 'LES':
self._FF_setup_LES(**kwargs)

elif self.inflowStr == 'TurbSim':
# We need to make sure the TurbSim boxes have been executed. Let's check the last line of the logfile
highboxlog_path = os.path.join(self.path, self.condDirList[0], self.caseDirList[0], 'Seed_0', 'TurbSim', 'log.hight1.seed0.txt')
if not os.path.isfile(highboxlog_path):
raise ValueError(f'All TurbSim boxes need to be completed before this step can be done.')
# # We need to make sure the TurbSim boxes have been executed. Let's check the last line of the logfile
# highboxlog_path = os.path.join(self.path, self.condDirList[0], self.caseDirList[0], 'Seed_0', 'TurbSim', 'log.hight1.seed0.txt')
# if not os.path.isfile(highboxlog_path):
# raise ValueError(f'All TurbSim boxes need to be completed before this step can be done.')

with open(highboxlog_path) as f:
last = None
for last in (line for line in f if line.rstrip('\n')): pass
# with open(highboxlog_path) as f:
# last = None
# for last in (line for line in f if line.rstrip('\n')): pass

if last is None or 'TurbSim terminated normally' not in last:
raise ValueError(f'All TurbSim boxes need to be completed before this step can be done.')
# if last is None or 'TurbSim terminated normally' not in last:
# raise ValueError(f'All TurbSim boxes need to be completed before this step can be done.')

self._FF_setup_TS(**kwargs)

Expand Down
4 changes: 4 additions & 0 deletions openfast_toolbox/fastfarm/TurbSimCaseCreation.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -124,6 +124,10 @@ def _discretization_lowres(self):
if self.ds_low is None:
self.ds_low = self.Cmeander*self.D*self.URef/150

# Spatial resolution for high-res
if self.ds_high is None:
self.ds_high = self.cmax

self.dy = np.floor(self.ds_low/self.ds_high)*self.ds_high
self.dz = np.floor(self.ds_low/self.ds_high)*self.ds_high
self.dt = self.dt_low
Expand Down