การทำค่าเฉลี่ยรายเดือนสำหรับผลคำนวณการไหลจาก MIKE 21 FM ที่ได้ผลเป็นไฟล์นามสกุล dfsu

ขั้นตอนการทำข้อมูลเฉลี่ยรายเดือน สำหรับผลคำนวณของแบบ MIKE 21 Flexible Mesh ซึ่งจะได้ผลคำนวณเป็นไฟล์นามสกุล DFSU มีประโยชน์เช่น หากต้องการหาทิศทางการไหลเฉลี่ยรายเดือนหรือรายปีในพื้นที่ผลคำนวณ 2 มิติ

ก่อนจะทำตามขั้นตอนด้านล่าง จะต้องมีไฟล์ผลคำนวณจากแบบจำลอง MIKE 21 FM เรียบร้อยแล้ว โดยเป็นผลคำนวณที่ครอบคลุมช่วงเวลาที่ต้องการหาค่าสถิติ

มีขั้นตอนดังนี้

1 สร้าง Data Extraction FM จากการเปิดโปรแกรม MIKE Zero แล้ว New File แล้วเลือกดังรูป

2 เพิ่ม Outputs ตัวแรก เลือกเป็น Area แล้วเลือกที่เก็บไฟล์ โดยตัวอย่างต้องการเฉพาะเดือน มค. หากผลคำนวณเดิมเป็นรายชั่วโมง เริ่มจาก 1 มค. ทั้งเดือน 31 วันจะมี 744 timestep

3 เมื่อตั้งค่าเสร็จแล้วให้ Save แล้วรัน หลังรันจะได้ไฟล์ที่ Extract ออกมาตามที่ตั้งค่าไว้

4 ให้เปิดไฟล์ ที่ตัดมาสำหรับทั้งเดือน มค. (ที่ต้องการหาค่าเฉลี่ย) โดยคลิกขวาแล้วเลือก Open with MIKE Zero Data Manager

5 เมื่อเปิดไฟล์ได้แล้วให้คำนวณค่าสถิติทีละ parameter ตัวอย่างรูปด้านล่างเริ่มทำสำหรับ Surface elevation

6 เสร็จแล้วให้ save เก็บเป็นชื่อตาม parameter แต่ละตัวแล้วปิด stat ไป

7 ทำการคำนวณ Stat สำหรับ parameter อื่นที่เหลือคือ U และ V เช่นเดิม แล้ว Save เพื่อให้รู้ว่าเป็น U และ V ตามลำดับ

8 เมื่อได้ไฟล์ Stat สำหรับเดือนที่ต้องการ ของ parameter ทั้งหมดคือ Surface, U และ V แล้วให้เปิดไฟล์ Stat Surface ขึ้นมาด้วย data manager แล้ว Save as เป็น Ave_HUV ซึ่งเป็นชื่อไฟล์ผลสุดท้ายที่ต้องการ

9 เบื้องต้น ไฟล์นี้จะมี 3 item อยู่แล้ว โดยเป็นค่า minimum, maximum และ average ของ Surface elevation ดังนั้น เราจะแก้ไขให้ minimum และ maximum เป็น average U และ average V ตามลำดับ โดยให้ใช้เมนู View/Items…

10 เลือกแก้ไข Item แรก ดังรูป

11 แก้ไข Item ที่สอง ดังรูป

12 เมื่อแก้ไข Item ทั้ง 2 แล้วให้คลิก OK จากนั้นเลือก item Average : U แล้วใช้เมนู Edit/Calculator 

13 จากนั้นคลิกที่ Add variable แล้วเลือกไฟล์ Stat U โดยให้เลือก Item ให้ถูกต้อง ดังรูป

14 เชคซ้ำอีกครั้งตรง Item ของ Variables จะต้องตรงกันคือ Average : U กับ Statistical mean : U velocity จากนั้นให้พิมพ์ Var3 ตรง Expression แล้วคลิก calculate

15 ทำเช่นเดียวกันกับ Item ที่เป็น Average : V

16 เมื่อทำเสร็จแล้วให้ Save แล้วปิด แล้วเปิดด้วย data viwer (ดับเบิลคลิกไฟล์ dfsu) แล้วตรวจสอบโดยการให้แสดง Vector โดยคลิกขวาที่รูปแล้วเลือก Vector ปรับแต่งความยาว vector โดยคลิกขวาเลือก data option 

การเตรียม 2D distributed rainfall ด้วยข้อมูลฝนรายสถานี

ขั้นตอนด้านล่าง แนะนำวิธีการเตรียมข้อมูลฝน 2 มิติที่เปลี่ยนแปลงตามเวลา

2D Distributed rainfall (ไฟล์ dfs2)

ด้วยการใช้ข้อมูลฝนจากสถานีตรวจวัดต่างๆในบริเวณนั้น (ไฟล์ dfs0)

โดยใช้เครื่องมือ MIKE SHE Tool box ใน MIKE Zero

โหลดไฟล์ตัวอย่างที่ได้ทดสอบทำด้วยวิธีการด้านล่างนี้จากที่นี่

มีขั้นตอนดังด้านล่าง

1) เตรียมไฟล์ Grid.dfs2 ซึ่งต้องมี EUM Type เป็น Grid Codes และมี Unit เป็น Integer โดยตัวเลขในกริดจะเป็นตัวแทนรหัสสถานีฝนที่ตรงกับ Item name ในไฟล์ข้อมูลฝน ดังนั้น การสร้าง Grid.dfs2 จึงต้องรู้ตำแหน่งสถานีแต่ละสถานีว่าอยู่ตำแหน่งใดในกริด เพื่อนำเลขรหัสสถานีฝนไปกรอกไว้ที่กริดนั้น

และให้เว้นว่างกริดที่เหลือที่ไม่ตรงกับตำแหน่งของสถานีไว้

2) เตรียมไฟล์ Rain.dfs0 ซึ่งเป็นไฟล์ข้อมูลฝนที่สถานีต่าง โดยมี ItemName เป็นรหัสสถานี ตรงกับที่เตรียมไว้ในไฟล์ Grid.dfs2 และให้มี Type เป็น Precipitation Rate และ TS Type เป็น Step Accumulated

3) เปิดโปรแกรม MIKE Zero แล้วสร้างไฟล์ใหม่ เลือกสร้างไฟล์ MIKE SHE / MIKE SHE Toolbox จากนั้นใช้เครื่องมือ File Converter โดยขั้นตอนภายในเป็นการเลือกไฟล์ Rain.dfs0, Grid.dfs2 และ ตั้งชื่อ output file โดยตัวอย่างนี้ให้ชื่อ Rain.dfs2

 4) ผลที่ได้จากเครื่องมือในข้อ 3 เมื่อเปิดตรวจสอบ จะได้ข้อมูลฝนรายสถานีที่กรอกลงในไฟล์ Rain.dfs2 ตามตำแหน่งของสถานีฝน (ตำแหน่งตามรหัสใน Grid.dfs2) โดยมีข้อมูลที่ timestep ต่างๆตรงกับข้อมูลฝนใน Rain.dfs0 (ตัวอย่างนี้ทำไว้เพียง 3 timestep)

5) จากนั้นให้ใช้เครื่องมือช่วยในการ interpolate โดยเลือกเมนู Tools/Interpolation...

6) ตรวจสอบใน Sub-Set ว่าได้เลือกทุก time steps และเลือก Entire set of grid points

7) ในส่วนของ Interpolation แนะนำให้ใช้แบบ Area integrated โดย Number of points ควร >2 (แทนจำนวนสถานีที่จะนำมาใช้ในการ Interpolate และกรอก Area side length ที่เหมาะสม (จำนวนกริด) ที่มั่นใจว่าในทุกตำแหน่งกริดว่างๆที่จะ interpolate ค่า จะสามารถหาค่าจากกริดข้างเคียงที่มีข้อมูลได้ในระยะห่างไม่เกินจำนวนกริดที่กำหนด โดยตัวอย่างนี้คือ 10 กริด

8) ตรวจสอบผลที่ได้จากการ Interpolate ถ้าได้ผลที่เหมาะสมให้ Save แต่ถ้าผลไม่ดี และต้องการ interpolate ใหม่ด้วยค่า number/length อื่นๆ ให้ปิดโดยไม่ save ไฟล์ แล้วเปิดไฟล์ขึ้นมาใหม่แทน เพราะการทำ interpolation จะไม่สามารถ Undo ได้

เท่านี้ก็จะทำให้ได้ไฟล์ข้อมูลฝนที่กระจายในพื้นที่แบบ 2 มิติ สำหรับนำไปใช้ในแบบจำลองที่ต้องการเช่น แบบจำลองน้ำท่วมในพื้นที่ 2 มิติ MIKE 21 และ MIKE 21 FM หรือแบบจำลองอุทกวิทยาน้ำ MIKE SHE ได้แล้ว

เทคนิคนี้ยังสามารถนำไปใช้ได้กับไฟล์ข้อมูลอื่นๆเช่น ข้อมูลการระเหยของน้ำเป็นต้น