The URL can be used to link to this page

Home My WebLink About03_EFDC_Model_OverviewEFDC Model TrainingDrew Ackerman, Cardno ENTRIXSeptember X, 2013 EFDC Model Background•Solves the three-dimensional, vertically hydrostatic, free surface, turbulent averaged equations of motions for a variable density fluid•Sigma coordinates and Cartesian or curvilinear, orthogonal grid•Assumes incompressible flow and hydrostatic pressure distribution with dynamically coupled salinity and temperature transport EFDC Model Background, Continued•Simulates:•Wetting and drying•Simulation of controlled flow structures•Vegetation resistance•Wave-current boundary layers and wave-induced currents•Embedded single port buoyant jet module EFDC Model Linkage•Provides dynamics for sediment, water quality and toxics fate and transport•Sediment simulations–Cohesive/non-cohesive, deposition/resuspension, sediment load and bed loads•Water quality–Eutrophication model, sediment diagenesis•Toxics–Trace metals/organic hydrocarbons and sediment interactions EFDC Versions•Tetra Tech version•EPA version•Original EFDC model by Hamrick•Last updated in 2002•Sandia Labs•Research version•Dynamic Solutions•Commercial version with GUI•Additional capabilities and support•Last updated in 2013 EFDC Application•Original model•EFDC (EPA) version•Expanded model•Dynamic Solutions (without GUI)•DOS window•Text Editor Model Configuration•Model developed to support hydrodynamic and water quality simulations•Stage, temperature, salinity•Nutrient transformations, dissolved oxygen kinetics, and eutrophication•Focus of this training is temperature and salinity Model Files File Types•Model control file•Describes model and sets up its output•Model grid files•Defines model grid orientation, elevations•Model boundary/input files•Characterizes what is entering/leaving model domain Model FilesEFDC.inp Control primary simulation control fileSHOW.inpFile controlling screen printWQ3DWC.inpWater quality simulation control file EFDC Model Model Grid FilesCELL.inp Grid cell typeCELLLT.inp Grid cell typeCORNERS.inpCorner coordinates for each cellDXDY.inp Grid cell dimensionsLXLY.inp Grid cell orientation Initial Condition FilesBEDBDN.inp Sediment bed bulk densityBEDDDN.inp Sediment bed porosityBEDLAY.inp Sediment bed thicknessSALT.inp Salinity concentrationSEDB.inp Cohesive sediment in sediment bedSEDW.inp Cohesive sediment in water columnTEMP.inpTemperature Boundary/Input FilesCWQSERxx.inpWater quality input filesGWMAP.inpDefines groundwater seepage zonesGWSEEP.inpDefines groundwater flow ratePSER.inp Time series tide boundary conditionQSER.inp Time series flowSDSER.inp Time series cohesive sedimentSNSER.inp Time series non-cohesive sedimentSSER.inp Time series salinityTSER.inpTime series temperature Boundary/Input Files, ContinuedWQALGG.inpAlgal dynamicsWQBNENFLX.inpTime series benthic fluxWQBENMAP.inpMapping info for spatial benthic fluxesWQSETL.inpWater quality settling ratesWQWCMAP.inpWater quality kinetic zonesWSER.inp Time series wind General File Comments•Discuss relevant files to hydrodynamic simulation and model grid•Free format•Sensitive to integer/real numbers•View in Notepad++ Model Grid Model Grid Overview•Need detailed information describing model cells•Cell orientation and where it is in the world•Current grid has one open boundary to east defining oceanic conditions Grid Development•Three options•GRIDGEN•Dynamic Solutions •Delft3d GEFDC•Original model grid generation program•Define boundary points•Define model grid type•Specify grid relaxation parameters•Cartesian, curvilinear •DOS and text based•Outputs image of grid in DXF format Dynamic Solutions•More flexible grid generation•Cartesian, curvilinear •Use geographic files to help define grid boundary•Can import other grids to EFDC grids•Deflt RGFGrid, Grid95, SEAGRID Delft3d RGFGrid•Intuitive GUI •Good online support•Active user community•Online training videos LTPR EFDC Model Grid•Lower Tar Pamlico River •6 sigma coordinate layer models•593 horizontal cells LTPREFDC Model Grid Grid Files •Model consists of a series of quadrilateral cells•Orthogonal and curvilinear•Files define cells•Location•Orientation•ShapeEFDC Grid FilesCELL.inp Grid cell typeCELLLT.inp Grid cell typeCORNERS.inpCoordinates for each cellDXDY.inp Grid cell dimensionsLXLY.inp Grid cell orientation CELL.inp i1234567891011121234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000911000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000099999999999999999999100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000959959955555555555559099999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999995555555555555555555809555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555570999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999555556000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000009555555555500000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000999999999940000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000093000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001314151617181234567890123456789012345678901234567890123456789012345678903500000000000000000000000000000000000000000000000000000000000034000000000000000000000000000000000000000000000000000000000000330000000000000000000000000000000000000009999900000000000000003200000000000000000000000000000000000000095559000000000000000031000000000000000000000000000000999999999959990000000000000000300000000000000000000000000000999555555555590000000000000000002900000000000000000000000000009555555555555900000000000000000028000000000000000000000000000095555555555599000000000000000000270000000000000000000000000000995559999999900000000000000000002600000000000000000000000000000995555900000000000000000000000025000000000000000000000000000000955559000000000000000000000000240000000000000000000000000000009555590000000000000000000000002300000000000000000000000000000095555900000000000000000000000022000000000000000000000000000000955559000000000000000000000000210000000000000000000000000000009555599000000000000000000000002000000000000000000000000000000095555590000000000000000000000019000000000000000000000000000000955555900000000000000000000000180000000000000000000000000000009555559900000000000000000000001700000000000000000000000000000095555559900000000000000000000016000000000000000000000000000000955555559900000000000000000000150000000000009990000099999999999555555559990000000000000000001400000009999995999999955555555555555555555900000000000000000013009999995555555555555555555555555555555559000000000000000000129995555555555555555555555555555555555555590000000000000000001155555555555555555555555555555555555555555900000000000000000010555555555555555555555555555555555555555559000000000000000000955555555555555555555555555555555555555555900000000000000000085555555555555555555555555595555555555555590000000000000000007555555555555555555555555555555555555555559000000000000000000655555995555995555555559995595999999999999900000000000000000055555999999999999559999909999990000000000000000000000000000004955590000000000999900000000000000000000000000000000000000000399599000000000000000000000000000000000000000000000000000000020999000000000000000000000000000000000000000000000000000000001000000000000000000000000000000000000000000000000000000000000 Grid Numbering CORNERS.inp DXDY.inp LXLY.inp Initial Conditions •Specify model conditions at start of the simulation•Focus on salinity and temperatureInitial Condition FilesBEDBDN.inp Sediment bed bulk densityBEDDDN.inp Sediment bed porosityBEDLAY.inp Sediment bed thicknessSALT.inp Salinity concentrationSEDB.inp Cohesive sediment in sediment bedSEDW.inp Cohesive sediment in water columnTEMP.inpTemperature Model Layers•Sigma (σ-stretched) layers•Layers expand and contract•Like an accordion•6 layers defined for the model•Lowest number on bottom123456 SALT.inp TEMP.inp Boundary Files Boundary/Input FilesASER.inpAtmospheric condition fileCWQSERxx.inpWater quality input filesGWMAP.inpDefines groundwater seepage zonesGWSEEP.inpDefines groundwater flow ratePSER.inp Time series tide boundary conditionQSER.inp Time series flowSDSER.inp Time series cohesive sedimentSNSER.inp Time series non-coheisve sedimentSSER.inp Time series salinityTSER.inpTime series temperatureWQALGG.inpAlgal dynamicsWQBNENFLX.inpTime series benthic fluxWQBENMAP.inpMapping info for spatial benthic fluxesWQSETL.inpWater quality settling ratesWQWCMAP.inpWater quality kinetic zonesWSER.inp Time series wind Atmospheric Conditions Atmospheric Conditions•Describe atmospheric conditions•Processes affected•Heating•Dissolved oxygen•Rainfall•Wet deposition of nutrients•Photosynthesis Atmospheric condition input file: ASER.inp•Define environmental conditions•Atmospheric pressure•Dry air temperature•Relative humidity•Rainfall•Evaporation•Solar radiation•Fraction cloud cover ASER.inp Wind input file: WSER.inp•Wind conditions•Speed•Direction WSER.inp Tidal Boundary Pressure time series: PSER.inp•Define seaward tidal condition at the “open boundary”•Three options•Define elevations (current model configuration for calibration and validation years)•Define model tidal harmonics•Use regression models based on stage at Washington (used to set up 2007 and 2008) Define Elevations•Current model configuration•Can capture storm surges•Need good local monitoring information for simulation period•Need time series for entire simulation period•Utilize monitoring data to define absolute or relative water surface elevation PSER.inpMPSERTCPSERTAPSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changingthe input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion Define Harmonics•Modeled data from other sources•JTides, ACRIC model•Estimate tidal harmonics for different constituents•M2, S2, N2, K2, O1, K1, Q1, M4, M6•http://www.unc.edu/ims/ccats/tides/tides.htm•Extrapolate observed stage at Washington to open boundary condition at the estuary mouth •Based on lag time and distance (Xu et al, 2008)•Can model any period once defined•More computationally efficient Model Card 14/15 Model Card 17/18 Surface Flow Inputs QSER.inp•Flow time series•Define flows into/out of model•River/stream flows•Point sources•Return flows•19 identified sources in the model QSER FlowsNS = 1, Greenville + GUC WTP IntakeNS = 2, Chicod Creek flow (cfs)NS = 3, Grindle Creek flow (cfs)NS = 4, Tranters Creek flow (cfs)NS = 5, Greenville WWTP Flow (cfs)NS = 6, Washington WWTP Flow (cfs)NS = 7, Belhaven WWTP Flow (cfs)NS = 8, PCPS WWTP Flow (cfs)NS = 9, ChocotowingNS = 10, BlountsNS = 11, DurhamNS = 12, SouthNS = 13, GooseNS = 14, BathNS = 15, Pungo CreekNS = 16, PantegoNS = 17, Pungo CanalNS = 18, GUC WTP Raw Water Withdrawal Rate (cms)NS = 19, GUC WTP Return Rate (cms) Sources QSER.inpMQSERTCQSERTAQSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion Groundwater Groundwater Flows•Different flows into model grid•Define zones and constant inflow rates•4 zones with 2 inflow rates (m/d) GWMAP.inp GWSEEP.inp Temperature and Salinity SSER.inpMCSERTCCSERTACSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion TSER.inpMCSERTCCSERTACSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion EFDC Input Files EFDC Control File•Main control file: EFDC.inp•Controls•Grid definitions•Inputs•Time steps•Output•Model calibration parameters•Arranged in “cards” Model Execution Simulations on Citrix Server•Run EFDC model on appserver.ncwater.org server•Same log in information as with OASIS modeling Virtual DesktopLogging onto the Citrix server allows for virtual desktop access to the model Virtual Desktop Model SimulationRun model from icon on desktop•Runs base case model•D:\Tar_users\dackerman\EFDCBaseCase Scenario Development Example Scenarios and Required Input Changes•Alteration of conditions for a pre-developed year•Development of a new simulation year•Simulation of sea level rise•Simulation of potential impacts of climate change Alteration of Conditions for a Pre-developed Year •Potential model years include 2001, 2003, 2007, 2008•Key parameters/factors could be altered•Flow•Temperature•Salinity•Files to change•QSER.INP •TSER.INP •SSER.INP Development of a New Simulation Year•Key parameters/factors•Flow, temperature, salinity, boundary conditions, atmospheric conditions, initial conditions•Files to change•QSER.INP, TSER.INP, SSER.INP, ASER.INP, WSER.INP, PSER.INP, TEMP.INP, SALT.INP, DXDY.INP•It will be important to process all data to develop accurate representation of new conditions•This will be a significant effort to undertake Simulation of Sea Level Rise•Key parameters/factors•Boundary conditions, initial conditions•Files to change•TSER.INP, SSER.INP, TEMP.INP, SALT.INP, DXDY.INP, PSER.INP, EFDC.INP•Requires changing seawater elevation but also salinity and temperature Simulation of Potential Impacts of Climate Change•Potentially hotter/cooler or dryer/wetter conditions •Key parameters/factors•Flow, temperature, salinity, boundary conditions, atmospheric conditions, initial conditions•List files to change•QSER.INP, TSER.INP, SSER.INP, ASER.INP, WSER.INP, PSER.INP, TEMP.INP, SALT.INP, DXDY.INP, EFDC.INP•Climate change will impact terrestrial and oceanic conditions•There will be a large degree of uncertainty and care should be taken with these simulations Model Scenario Runs Scenario Simulations on ServerCopy EFDCBaseCase directory into a new directoryModify input files in new directoryRun simulation in new directory Model Output EFDC Output•EFDC example output file•Copy post-processing files (PostProc.exe) into directory where simulation was run•Copy post-processing Excel file (PostProcessing.xlsx) onto your computer Your local computerCitrix Server Data from EFDC Output•Post-process results on server•Will speed up download times (output files are large)•Copy “PostProc.exe” into scenario folder on the server•From \EFDC_PostProcessing Data from EFDC Output•Extract data from salinity and temperature output files•16 stations specified in EFDC.INP file–Card 87•SALTSxx.out and TEMTSxx.out•Run the executable file to write a file that averages the top, middle, and bottom two layers (temperature and salinity)•Double click on PostProc.exe•To generate TempSalt.csv Example model outputadfdfafafaf>werwreqwrewreafdafasfdasfdas>adffasfdasfafdsa Calculate Percentiles•Steps to calculate percentiles1.Double click “PostProc.exe”2.Open PostProcessing.xlsx on your computer3.Copy TempSalt.csv from the scenario folder to your computer4.Open TempSalt.csv in EXCEL5.Select columns B through J and copy6.Paste into Cell B1 in the “Data Import” tab7.Results are shown in “Temp Salt Percentiles” PostProcessing.xlsxOpen TempSalt.csv in EXCEL and paste those columns into Columns “B” through “J”Percentiles for the top two layers, middle two layers, and bottom two layers are calculated in the “Temp Salt Percentiles” tab Temp Salt Percentiles Additional Model Post Processing Advanced Model Post Processing•Requires model user to write post processing code•FORTRAN, C++, BASIC, Python, Matlab•Will require more knowledge of model output structure than has been presented here Issues with direct linkage to the OASIS model Constraints to Linking the LTPR Model to the OASIS Model•LTPR model has been set up for four specific years •Headwater flows, open boundary condition, atmospheric conditions, withdrawals, and discharges•OASIS model is a time-series model developed for a longer period•Changes to the LTPR model must be hard wired into input files •Cannot read output files from OASIS directly•Running additional years with the LTPR model will require reconfiguring the model boundary and forcing files Questions?