Open
Conversation
Author: Jonathan R. Commit of an updated Beam Deposition module. The following changes were added: * Included ability to run beamdep.py as a stand alone script. The input is the data from shot 118890.1560. The values are hard coded in a debug version of the Core class and Inp class. These debug versions are classes defined within beamdep.py. The profiles are found in a debug file in ./inputs/ generated using an np.savetxt command of the ion and electron densities and temperatures. * Deposition profiles have always been significantly off. When integrating over dPdr for all 3 beam types, we consistently underestimated the total power. This was addressed by multiplying by the integral of H(r) until we can figure out why this is necessary. * Debug dictionary added to be able to plot various quantities of interest to see how changes to the code affect the results.
Date: 09/04/2018
# Conflicts: # .idea/gt3-dev.iml # .idea/misc.xml
Nbeams dev
Date: 2/26/2019 Time: 20:32 More work
Date: 2/27/2019 Time: 12:09pm Including debug functions in radial transport to verify source terms, fluxes, etc. are reasonable. Can be triggered by passing debugFlag=True to radialtransport. The neutrals calculation and related calculations (e.g., Sion) can also be deactivated by passing neutFlag=False to radialtransport.
Date: 2/27/2019 Time: 12:09pm Including debug functions in radial transport to verify source terms, fluxes, etc. are reasonable. Can be triggered by passing debugFlag=True to radialtransport. The neutrals calculation and related calculations (e.g., Sion) can also be deactivated by passing neutFlag=False to radialtransport. gt3.py now includes modules that let you interactively create input files from commandline. The input files are saved to the input folder as togt3_shotnum_timeid. See the code comments. The main entrance to the interactive mode is gt3Prep. You should provide the shot id, timeid, and run id, as well as the filename. The "quiet" parameter disables interactive file generation and tells the code to look for an input file. If no input file is found, program exits. This will be helpful when/if we want to run gt3 on many shots at once and we want the code to not bug us.
Radial transport dev. Merging in Nick's changes to radial transport.
Date: 8/27/2019 Time: 11:06pm
GT3 Radial Transport Branch 8/27 commit
Date: 05/18/2020 Time: 11:28AM Minor error handling modifications.
Consolidated a few RadialTransport calculations Added new GT3TestBase class to TestBase module
Added plotting routines to RadialTransport Added mode flag to getGT3Test method
…than the main GT3 code if needed. Updated and added some Core plotting routines. Removed some unnecessary input variables Deleted Neutrals.Functions that were not necessary Attached core/inp objects to neutrals object for referencing
… Beams.py Added additional plotting routines to Core Added some explicit float to values to ensure they are calculated as floats if import division is not set. Remvoed NeutPatcy in favor or implementing method in Core Removed Test python files Removed Chi class Added IOL Comparison utility
Added an interpolator for d3d data that are a bit wonky.
…t object. Signed-off-by: Jonathan R <veto1024@gmail.com>
…ethods for calculating gamma/Q as instance variables. Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
…updated using update functions. Particle/Heat flux differential methods are updated to take in a splineVerify keyword to verify the splines are fitting correctly. Signed-off-by: Jonathan R <veto1024@gmail.com>
…lues of the asymmetry coefficients. Signed-off-by: Jonathan R <veto1024@gmail.com>
…ry coefficients. Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
…ectorialBase was not computing correct total velocities. Signed-off-by: Jonathan R <veto1024@gmail.com>
…rom pertrubation theory Signed-off-by: Jonathan R <veto1024@gmail.com>
…rom pertrubation theory Signed-off-by: Jonathan R <veto1024@gmail.com>
Improved ability to handle multiple x-points in core Added ability to take D rotation date as frequencies instead of velocities Added GT3LineString class, which supers(LineString) and adds functionality to quickly plot linestrings against our wall Signed-off-by: Jonathan R <veto1024@gmail.com>
…etermination of separatrix, strike points, etc. This will be helpful when wanting NeutPy to be able to get these flux surface data without running its own routine. This allows just GT3's algorithms to be updated and simply called from NeutPy whenever we need to improve the algorithm for various shot geometries. This commit is still broken! Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
Updates to ProfileClasses.py: Added set_to_zeros method to set profiles to zeroes. Added private spline parameters and ability to update on the fly. Added flag to indicate a profile has been overwritten manually. Added a multiplier argument to set_raw_data when setting data that needs to be converted between, e.g., kev/ev/J. Added method to attenuate a profile if necessary. Added method to overwrite a profile with a splien interpolation if the data are messy but otherwise fine. Necessary when derivatives of data are taken. Added TwoDProfileWithHM to programmatically add heat mapping capabilities to 2D data. Psi updates: Added kwargs to allow kwargs to be passed to control how the x-point/psi normalization procedure occurs. Misc. Added kwargs to be able to pass through to Psi when it becomes necessary to force an x-point to normalize to. When neutral data are updated, any NaNs will raise an exception. Cooling rate in imp_radiation was still being calculated even if neutrals were turned off, confounding any neutral sensitivity study. Moved chi calculation private and allow the ability to programmatically change the viscous heating asymmetry coefficients. Added methods to catch flags in rtrans that will overwrite T/n profiles with splined profiles. Fixed input reading process so that velocities coming in as frequencies are properly handled. Added neutpy configuration to testing folder Added new unit tests Signed-off-by: Jonathan R <veto1024@gmail.com>
…n radial transport. Impurity radiation would incorrectly turn off if neutrals calculation was turned off. Signed-off-by: Jonathan R <veto1024@gmail.com>
Added a couple of collision frequencies to rtrans. Signed-off-by: Jonathan R <veto1024@gmail.com>
Added a couple of collision frequencies to rtrans. Signed-off-by: Jonathan R <veto1024@gmail.com>
Changed docstring in OneDProfiles,TwoDProfiles to be properties. Better error handling in OneDProfiles,TwoDProfiles. Added ability to easily clear Profile L, derivative, FSA when modifying data. Added additional testing of ProfileClasses.py Signed-off-by: Jonathan R <veto1024@gmail.com>
Signed-off-by: Jonathan R <veto1024@gmail.com>
…ive. Unsure this actually matters, though. Exposed epsilon from the IOL calculation as a OneDProfile. Added ability to directly scale input profiles/floats to run sensitivity tests. PlotBase can now take in kwargs for the x/y-labels. Signed-off-by: Jonathan R <veto1024@gmail.com>
…ive. Unsure this actually matters, though. Exposed epsilon from the IOL calculation as a OneDProfile. Added ability to directly scale input profiles/floats to run sensitivity tests. PlotBase can now take in kwargs for the x/y-labels. Signed-off-by: Jonathan R <veto1024@gmail.com>
# Conflicts: # GT3/utilities/PlotBase.py
…ow IOL/non-IOL comparisons. Added module-level constants for plot markers and sizing. Will want to propogate these changes further through the code. q95 was being calculated incorrectly (used rho when should have used r in the spline) Neutrals module can more easily be rerun if iterating on neutpy outputs. Added ability in RadialTransport to iterate on neutpy output. Initial SOL n,T are generated by constant D/chi. Once calculated in RadialTransport, these fluxes are fed back as sources to NeutPy. Added noIOL values to various calculations in RadialTransport for comparisons. Added method to change viscous heating asymmetries. Added collisionality plotter in RadialTranposrt, although probably should be a OneDProfile. Added a few kwargs for graphing chi to make it more presentable. Added ability to scale and translate values upon input in ReadInFile for sensitivity studies, etc. Signed-off-by: Jonathan R <veto1024@gmail.com>
wdeshazer
approved these changes
Oct 16, 2024
Member
There was a problem hiding this comment.
@veto1024 Would you mind putting nbeams /home/jonathan/Dropbox/GTEDGE/MyPrograms/GTEDGE/lib/beams/NBeamsMDS/NBeams/bin/Release/nbeams' in a shared location or if you are hosting, resubmit this change with the appropriate URL.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
4 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Bringing all the changes to GT3 to the master branch since it now works pretty fluidly with NeutPy, and it looks like most the important bugs are worked out.