Transcript 강의4_ELCOM-CAEDYM 사용법

Prepared by S.W. Chung
Dept. of Environmental. Eng., Chungbuk National University
E-mail: [email protected]

MODEL SET UP

PRE-PROCESS

OUTPUT CONTROL

MODEL RUN

POST-PROCESS

DOWNLOAD

APPLICATION EXAMPLES
2

모델 다운로드 (http://www.cwr.uwa.edu.au)

실행파일 폴더 만들기 c:\cwrsoft\bin_2013

실행파일 폴더 Path 설정
①
②
③
④
⑤
⑥
⑦
Window 시작에서 “제어판” 열기
“시스템” 클릭
“고급” 탭 클릭
“환경변수” 클릭
사용자 변수에서 “Path” 더블클릭
변수값 텍스트 박스에 c:\cwrsoft\bin_2013; 입력
“확인” 클릭
3

실행파일 폴더 Path 설정 확인
Window 시작에서 “실행” 클릭
실행 윈도창에 “cmd” 입력
4

실행파일 폴더 Path 설정 확인
> echo %path%
5

예제파일 복사 및 모델링 작업 폴더 구성
6

모델링 작업 폴더 구성


Pre는 Pre-process에 필요한 파일과 bathymetry 폴더 포함해야 함
Run은 모델 실행에 필요한 파일과 “dbconv”, “infiles”, “ncfiles”,
“txtfiles”, ”unffiles” 를 포함해야 함
7
 입력 파일 준비
bathymetry.dat : 저수지 지형과 수치해석을 위한 Grid 정보 파일
2. bc.dat: 경계조건(유량, 수질, 기상)을 지정하는 셀 지정 파일
3. levee.dat: levee 경계조건을 지정하는 파일
4. update.dat: 경계조건을 업데이트 하는 내용 지정하는 파일
5. run_pre.dat: pre_elcom.exe 파일 실행 위한 입력파일 설정 파일
6. sparsedata.unf and usedata.unf : pre_elcom.exe 실행 후 생성되는 파일
7. datablock file (*.db): 출력결과 지정 파일
8. temporal boundary condition files (*.dat): 경계조건 입력 파일
9. initial conditions file(*.ini): 초기조건 입력파일
10. Chemical parameters file (*.chm): 이화학적 수질매개변수 파일
11. Biological parameters file (*.bio): 생물학적 수질매개변수 파일
12. Sediment parameters file(*.sed): 유사와 관련된 매개변수 파일
1.
8

Test Run : RoundLake

cmd 창 열기: 시작-프로그램 및 파일 검색에 cmd 입력

폴더 이동 : >cd c:\cwrsoft\bin_2013\RoundLake

폴더 내용 확인: > dir
9

Pre-Process (pre_elcom.exe) 실행화면
> pre_elcom run_pre.dat
두 개 파일이 /run/infiles 에 생성됨
sparsedata.unf
usedata.unf
10
 ELCOM-CAEDYM 모의 흐름도
11
bathymetry.dat
bc.dat
run_pre.dat
Preprocessing
pre_elcom
sparsedata.unf
usedata.unf
12
bathymetry.dat
bc.dat
Temporal boundary conditions files
- meteorological data
- inflow / outflow data
run_pre.dat
Initial condition files
pre_elcom
datablock.db
sparsedata.unf
run_elcom.dat
Preprocessing
usedata.unf
13
bathymetry.dat
bc.dat
Temporal boundary conditions files
- meteorological data
- inflow / outflow data
run_pre.dat
Initial condition files
pre_elcom
datablock.db
sparsedata.unf
run_elcom.dat
Preprocessing
usedata.unf
ELCOM
Model
unffiles
14
bathymetry.dat
bc.dat
Temporal boundary conditions files
- meteorological data
- inflow / outflow data
run_pre.dat
Initial condition files
pre_elcom
datablock.db
sparsedata.unf
run_elcom.dat
Preprocessing
usedata.unf
ELCOM
Model
unffiles
dbconv
run_dbconv.dat
NetCDF (*.nc) files
Postprocessing,
visualization
15
 지형자료 작성 절차
16
 Creation of Bathymetry File (convertXYZtoELCOMbathy.m)
function [] = convertXYZtoELCOMbathy(bathFile,xyzFile,outlineFile,dx,dy,dz,land_value,open_value)
% function [] = convertXYZtoELCOMbathy(bathFile,,outlineFile,dx,dy,dz,land_value,open_value)
%
% Generate an ELCOM bathymetry file form xyz data
% Two input files are required. One containg xyz data of bathymetry and one
% with xyz data of the lake outline.
%
% Inputs:
%
bathFile : filename of output bathymetry file
%
xyzFile : file containing three colums of xyz bathymetry data
%
outlineFile : file containing three colums of xyz data for perimiter of lake
% dx
: dx value to use
% dy
: dy value to use
% dz
: dz value to use
% land_value : (optional) land_value for bath file (default is 999.9)
% open_value : (optional) land_value for bath file (default is 888.8)
17
 Creation of Bathymetry File
>> convertXYZtoELCOMbathy(‘dc_test.dat’, ‘20all-1.xyz’, ‘dc_outline.prn’, 150, 150, 1, 999.9, 888.8)
>>modifyELCOMbathy(‘dc_test.dat’)
18
 Creation of Bathymetry File
Contains
• Location information (latitude,
longitude, altitude)
• Grid spacing (horizontal and
vertical)
• Bathymetry matrix (x – rows,
positive down, y – columns
positive right)
19
 Defining Flow Direction
y
Inflow
x
흐름방향은 셀의 중앙을 기준으로 Positive (xp, yp)와 Negative (xm, ym) 판단
• flow_xp : Flow coming from the x-positive direction when viewed from the cell center
• flow_xm : Flow coming from the x-minus direction
• flow_yp : Flow coming from the y-positive direction
• flow_ym : Flow coming from the y-minus direction
20
 Defining Flow Direction
y
Inflow
point-of-view
flow_xm
x
흐름방향은 셀의 중앙을 기준으로 Positive (xp, yp)와 Negative (xm, ym) 판단
• flow_xp : Flow coming from the x-positive direction when viewed from the cell center
• flow_xm : Flow coming from the x-minus direction
• flow_yp : Flow coming from the y-positive direction
• flow_ym : Flow coming from the y-minus direction
21
 Creation of Boundary Condition File
DC_BC_150.dat
BC identifier
(i, j, k) location
Keyword
• open – open boundary
• flow – inflow/outflow
• land – land cell with no-slip
• slip – land cell with slip condition
• section – variable surface forcing
• bubbler – for destratification systems
22
 Run Pre-Processor

> pre_elcom.exe run_pre.dat
23
 datablock.xml
24
 datablock.xml <= ELCOM User Manual 참조










GlobalTime
globalVariable
PROFILE
CURTAIN
SHEET
AggregateSheet
OUTFLOW
GENERAL3D
ALL3D
DRIFTER
25
 Input Files
26
 Run ELCOM

> elcd.exe
run_elcom.dat
27
 Run ELCOM-CAEDYM


> createCDsim caedym.con CAEDYM.nc
> elcd run_elcom.dat
run_elcom.dat
28
 Creation of ncfiles

> dbconv.exe run_dbconv_profile1.dat
Convert unformatted files (*.unf) to NETCDF format files (*.nc)
run_dbconv_profile1.dat
29
 Creation of ncfiles using DOS batch command

> run_dbconv.bat
Convert unformatted files (*.unf) to NETCDF format files (*.nc)
run_dbconv.bat
30
 Using Matlab for Profiles Comparison

>> profile_NTU_1
31
 Using Matlab for Curtain Plot


>> data = nldnc (‘curtain2.nc’)
>>[H]= plotELCOMcurtain(data, ‘TEMPERATURE’, 900)
32
 Make a Movie

>> callMakeMovieDaecheong
33
 Using ELCOM-CAEDYM GUI : Modeller, ARMS-Lite
34
http://www.cwr.uwa.edu.au
35
35
DOWNLOAD
CWR Software Policy
36
DOWNLOAD
Commercial Membership Packages
37
DOWNLOAD
Download
1. ELCOM User’s Manual
2. ELCOM Science Manual
3. CAEDYM User’s Guide
4. CAEDYM Science Manual
5. Other Publications – http://www.cwr.uwa.edu.au/research/publications.php
38
 Example Applications
•
•
•
•
•
•
•
•
•
Actively used in more than 60 countries
 Turbidity
 Eutrophication
 Nutrient and carbon budgets
 Fisheries, aquaculture and food web studies
 Public health risk assessments (metals, pathogens)
Upper Yarra Reservoir, Melbourne, AU
Lake Burrangorang, Sydney, AU
North Pine Dam, Queensland, AU
Rivers (eg Murray River, Tiete River)
Lake Como and Venice Lagoon, Italy
Lake Kinneret, Israel
Lake Constance, Germany
Brownlee Reservoir, USA
39
 Temperature Structure – Mundaring, WA
40
 Inflow Dynamics – North Pine Dam
41
 Tracer Study – Lake Como
42
 State-Of-The-Art: Case Studies
•
Pelagic Food Web: Lake Kinneret, Israel
–
Nutrient cycling
–
Species succession
–
Phytoplankton-Zooplankton interaction
–
Bacterial biomass
•
Carbon Budgeting: Lake Burragorang, AU
•
Microbial pollution: Sugarloaf Reservoir, AU
•
Heavy Metal Cycling & Bioaccumulation: Lake Coeur d’Alene, USA
•
Sediment Diagenesis: Cockburn Sound, AU
•
Benthic Macroalgae: Lake Erie, USA
43
 3D Movie – Lake Kinneret
44
 Lake Kinneret: Carbon Fluxes
45
 Burragorang Reservoir: 1997 Flood Event
46
 North Pine Reservoir: Pathogen
47
 Daecheong Reservoir: Turbidity Flow
A2
R4
R3
R2
A1
R1
 curtain1: R1-R2-R3-R4, curtain2: A1-R4-A2
48
 Daecheong Reservoir: Algal Bloom
2001 July ~ August Microcystis bloom event
Temperature
Chl-a (CYANO)
Tracer from Tributary
49
 소옥천 영향분석: PO4-P
L:\CWR\Daecheong\DC_20
01\Caedym3.1\run8_1\run\
ncfiles\DC_2001_PO4_Sheet.
m
50
 용담댐 저수지 수치추적자
[TRACER]
51
51
 용담댐 저수지 탁수모의
52
52
 용담댐 저수지 조류모의
53
53