ENEN90031: Quantitative Environmental ModellingMATLAB Programming AssignmentDUE DATE: 25-March-2019, 10:00 AMThis assignment contributes 10% to your final subject mark.OverviewThis assignment requires you to develop a MatLab function to model an infiltration basin from 01-Jan-1990to 31-Dec-2015 (if your catchment don’t have continuous data for this period, you can choose any other 15year time period) and to analyse the behaviour of the system over time. Successful completion of thisassignment will assist you in the completion of later assignments within this subject.You will need to complete a function that models the interaction between rainfall, recharge and runoff.Then you need to analyse the results over time. You will also need to complete a script that calls yourfunction and analyses the outputs. The script needs to produce a series of plots. A brief write-up of yourresults is required.TaskYou are required to finish a MatLab program (this is a function with inputs and outputs, not just a script) tomodel an infiltration basin (look in the appendix for a conceptual diagram and other information about thesystem). The overall purpose of this program is to model the water infiltrated to the underlying unconfinedaquifer and to examine how the recharge and runoff from the basin varies over time. Your program shouldundertake the following tasks.a) Model the overall system and produce appropriate plots. We have provided a partly writtenMatLab function, but you will need to complete it. Before you start adding to the function, pleaseread the programming tips on page 2, the description of the model in Appendix 1 and the functionitself. Your first step should be to add some comments to the existing code as you decipher how itworks. Then make a plan of the key steps involved and add this as comments to the function.b) Your function will produce one output structure. The first line of your function looks like this:function [InfiltrationOutput]=infiltration(Rainfall,PETData,basinCharacter)Where;Rainfall is an nx2 matrix with col 1: date in datenum format and col 2: daily precipitation in mm/dPETData is an nx2 matrix with col 1: date in datenum format and col 2: daily pet in mm/dbasinCharacter is a structure with different basin characteristics such as basinArea: Basin area in m2? maxStorageVolume: Maximum storage capacity in m3 catchmentArea: catchment area contributing runoff from to the basis (assume this isimpervious with 100% runoff) m2 RechargeMax: Maximum recharge rate from the basin in mm/dc) The output structure has already been initialised within the provided function. It includes membersthat are all vectors and are: recharge: the rate at which water is recharged to the aquifer (m3/d) overflowRate: the rate at which water overflows from the basin (m3/d) volume: the volume of water stored in the basin at the end of each time step (m3)d) You should also write a script to load the data, call the modified function and to analyse the modeloutputs by generating the following graphs and values: timeseries graphs of one month of data for the three model outputs (choose a period with somerainfall and overflow events) The first column of the PETData has dates in MatLab datenumformat. Annual timeseries of overflow from the basin Bar graphs showing the total number of days the basin overflows per month (Select any one yearfor this) Your graphs should have appropriate axis labels, axis scales, titles, and legends. Your script should also check that your infiltration.m function calculations are conserving thewater balance over the simulation period (we are assuming zero losses from the system and allthe relationships are linear in nature). This is an example of one way to test whether your modeloutputs make sense. Testing is important part of quality checking. You might think of other waysof checking your model. In your report please include a table with mean annual recharge, mean annual overflow andmean storage volume for the simulation period.Note: Make sure you choose the PET and Precipitation of the same time period for running yourinfiltration basin model and all the inputs and outputs are in correct units.A suggested programming processTo assist you with writing your function (or script), the following tips may be useful. The first line must define the function name, inputs and outputs (not needed for scripts). This isalready done for you. After you’ve written the first line, write comments explaining each sub-section of your program.This will help you plan and structure each aspect of the program. Given some of the code isalready written, start by writing comments explaining what it does. You might also find it usefulto fix up the formatting at this time. The comments should be clear enough so that someone else could understand the overallpurpose of your program. They should also clearly state the programs major assumptions andlimitations; the inputs required and the outputs created; and an explanation of each majorstep. After you’ve written the comments, start writing MatLab commands for the first sub-section of your program. Save the program, and then use the MatLab debugger check that the first subsectionis working as you expect. Repeat this process until the whole program isfinished. The first sub-section of your program should initialise the outputs. This can be done by settingthe output variables, say “output1”, to output1=NaN(rows,1). Again this has been done for youhere, so you can learn from this example. Calling the figure() function before a new plot creates a new figure window and stops an old plotfrom being over-written. This might be useful in your script.You should review the plotting section of help to see what graph types are available and selectthose that suit your purposes.Report/SubmissionYou will need to submit the following1. Your MatLab function that models the infiltration basin and the script that produces the graphs. Youneed to use comments and formatting (indenting if blocks and for loops, etc) to make your functionand script clear to read. You will need to do some formatting of the code supplied! You can indentlines by highlighting them as a block and using the tab key (or shift-tab to reduce indent).2. Your MatLab script that calls the function and organises and plots the data. Also make sure yourscript is clear to read.3. Write a report of up to 3 pages (plus appendices) showing the results requested above and brieflydescribing them. Please note, this does not need to be a formally structured report and only the firstthree pages will be marked (plus your MatLab code in appendices). Your report and MatLab filesshould be submitted online via the subject LMS link for MATLAB file. Follow the instructions there.Your script and function should be included in your report document as Appendix 1 and Appendix 2.Please copy your code from the MatLab editor and paste it by right-clicking and selecting the optionto “Keep Source Formatting”. The MatLab code is not included in your page limit.Your report should include: your plots for your infiltration basin as specifiedabove; tables as specified above; and a brief description of any additional error checking you have added to the function and the variousresults you have produced.All plots should include titles, axis labels and legends (where there are multiple series plotted). Axislimits should be set appropriately. All figures (plots) and tables should have appropriate captions.AssessmentLate Penalty: assignments submitted late will be penalised 10% per day or part thereof; assignmentssubmitted more than 5 days late will not bemarked.Over-length reports: we will stop marking at the end of page 3 of the main body of your report.Your assignment will be assessed as follows:? 20% for the function structure; 10% for clear commenting and indentation; 10% for the numerical outputs from the program;10% for handling of battery full and empty special cases; 30% for the quality of the plots produced, especially the labelling andformatting;10% for function testing; and 10% for the quality of yourscript.Appendix 1Infiltration basins are impounding facilities which temporarily store surface runoff and allow it to infiltrateinto the subsurface.In this problem the infiltration basin is considered as a leaky bucket. The leakage from the bottom of thebucket will be going to the groundwater systems as recharge. In this case, total water (i.e., (Catchment Area+ Basin Area ) x Rainfall - Basin Area x PET) will be available for recharge, overflow or change in volume. Thebasin was assumed to be leaking at a constant rate which is independent of the total water available. Theoverflow will occur only after the recharge is satisfied and the storage is full.We can use a conservation equation as the core of the infiltration model. This conservation equation saysthat the rate of change of volume in the basin is equal to the rate of inflow minus the rate of outflow.�V = Volumet = TimeI = Total inflowRe = RechargeRechargeRunoff InflowMax. Soil moistureholding capacityRu = RunoffFor more details on infiltration basin design follow the link belowhttp://www.water.wa.gov.au/__data/assets/pdf_file/0017/1691/99302.pdfAppendix 2: Downloading Climate DataTo obtain your data for your MatLab Assignment you should:1. Go to http://www.bom.gov.au/climate/data/stations/2. Select the Weather Station Number tab3. Enter you station number in the Number box4. Go to step 25. Select the weather element(Rainfall-total)6. Select the reporting frequency (daily)7. Go to step 38. You can download the data from the links on the right-hand side under Data availability. Yourstation should be the first on the list.The following table shows student ID numbers and station numbers – use this station for step 3.Student ID Station ID Student ID Station ID Student ID Station ID1057975 86071 761462 83083 1007728 84016893750 90015 923060 81049 941470 76047767423 82039 947512 88164 760354 850961006796 86038 754869 76064 1006877 88109971039 87113 905458 88110 1006880 800231039422 85279 918276 84070 1066032 80091963755 87031 1019430 86038 864506 79028813187 85072 807198 88023 882719 88043902586 76031 882020 86282 907975 86104834770 89002 999937 86077 290102 86371696273 85279 733295 82042 884827 86127834204 80015 979055 86383 952258 89085954875 85072 966358 86354 887294 78077934932 82076 1002132 81123 884308 82076925588 86351 1022348 83083 868569 86351974872 89002 990881 84143 953648 90171883380 86361 1033239 83024 982409 901731048165 86375 1050847 85283 948610 852801031359 83084 860336 82138 825653 84016983244 90176 887721 90180 878825 81123880132 83085 1035151 83084 876738 86371899282 88110 1010244 90175 905076 901841021802 87168 886616 86373 854774 84142967920 80015 241334 89085 1029294 82139976341 86361 1006886 88162 724028 84084961280 90184 584954 81124 977846 852961056911 84144 757354 86383 878562 88051852928 88164 928538 84143 1035690 90015876401 81125 864689 86375 901122 84142973193 79105 1006823 77094949946 86351 918040 86371 本团队核心人员组成主要包括硅谷工程师、BAT一线工程师，精通德英语！我们主要业务范围是代做编程大作业、课程设计等等。我们的方向领域：window编程 数值算法 AI人工智能 金融统计 计量分析 大数据 网络编程 WEB编程 通讯编程 游戏编程多媒体linux 外挂编程 程序API图像处理 嵌入式/单片机 数据库编程 控制台 进程与线程 网络安全 汇编语言 硬件编程 软件设计 工程标准规等。其中代写编程、代写程序、代写留学生程序作业语言或工具包括但不限于以下范围:C/C++/C#代写Java代写IT代写Python代写辅导编程作业Matlab代写Haskell代写Processing代写Linux环境搭建Rust代写Data Structure Assginment 数据结构代写MIPS代写Machine Learning 作业 代写Oracle/SQL/PostgreSQL/Pig 数据库代写/代做/辅导Web开发、网站开发、网站作业ASP.NET网站开发Finance Insurace Statistics统计、回归、迭代Prolog代写Computer Computational method代做因为专业，所以值得信赖。如有需要，请加QQ：99515681 或邮箱：99515681@qq.com 微信：codehelp
网友评论