การใช้งาน Weather Research and Forecasting (WRF) model เบื้องต้นมีขั้นตอนดังนี้
ต้องลง WRF เรียบร้อยแล้ว
ต้องลง WPS เรียบร้อยแล้ว
ต้องดาวน์โหลด Geography data for WPS geogrid.exe เรียบร้อยแล้ว
ต้องลง NCL เรียบร้อยแล้ว (สำหรับหัวข้อ 2.2)
1. Prepare global model output
1.1 Download global model output
สามารถดาวน์โหลด Global model output ฟรีได้จากเว็บไซต์ เช่น https://www2.mmm.ucar.edu/wrf/users/download/free_data.html
และสามารถเลือกดาวน์โหลดข้อมูลจาก Global Forecasting System (GFS) โดยตรงได้จาก ที่นี่
ข้อมูล GFS ถูกสำรองไว้บนเว็บไซต์นี้ เพียง 10 วันเท่านั้น
Example data
สามารถคัดลอกข้อมูล Global model output สำหรับตัวอย่างนี้ ได้โดยใช้คำสั่ง
cp -r XXX . # Copy the prepared example data to the current directory
1.2 Link ./WPS/Vtable
ในการนำเข้าข้อมูล จำเป็นต้องสร้าง Variable table เพื่อใช้ในการจับคู่ชื่อ ประเภทและหน่วยของตัวแปรในชุดข้อมูลนำเข้ากับตัวแปรที่ WRF model ใช้
WPS มี Variable table สำหรับนำเข้าข้อมูลจาก Global model บางตัวอยู่ในโฟลเดอร์ ./WPS/ungrib/Variable_Tables/
ก่อนการแปลงข้อมูลด้วย ungrib.exe จำเป็นต้องลิ้งค์หรือคัดลอกตารางที่เหมาะสมออกมายังตำแหน่งเดียวกันกับ ungrib.exe โดยใช้ชื่อ Vtable
cd WPS ln -s ./ungrib/Variable_Tables/Vtable.GFS ./Vtable # For GFS data, link Vtable.GFS to the same location as ungrib.exe, naming it Vtable
1.3 Run ./WPS/link_grib.csh
ungrib.exe จะอ่านเฉพาะชื่อไฟล์ GRIBFILE.AAA, GRIBFILE.AAB, …
สคริปต์ link_grib.csh ทำหน้าที่ลิ้งค์และแปลงชื่อไฟล์ในโฟลเดอร์ที่มึชุดข้อมูลมายังโฟลเดอร์ WPS ให้โดยอัตโนมัติ
./link_grib.csh Path_to_your_data # Path to the folder containing the data *** {USER INPUT} ***
Common mistake
Setting path to the wrong folder
The folder contains some unrelated files or files from other models
2. Run WPS
2.1 Edit ./WPS/namelist.wps
แก้ไขไฟล์ namelist.wps ในโฟลเดอร์ WPS ดังนี้
# Edit namelist.wps # as follows while leave others unchanged start_date = '2022-05-09_00:00:00', '2022-05-09_00:00:00', end_date = '2022-05-09_12:00:00', '2022-05-09_12:00:00', i_parent_start = 1, 43, j_parent_start = 1, 53, e_we = 120, 121, e_sn = 150, 151, dx = 36000, dy = 36000, map_proj = 'mercator', ref_lat = 13.208, ref_lon = 101, truelat1 = 13.208, truelat2 = 0, stand_lon = 100.861, geog_data_path = '../WPS_GEOG' # *** {USER INPUT} ***
สำหรับการตั้งค่า geog_data_path ในบรรทัดที่ 17 ดู Geography data --- for WPS geogrid.exe
2.2 Check domains (Optional, require NCL)
สามารถตรวจสอบ Domain ที่ถูกตั้งค่าไว้ใน namelist.wps ได้โดยใช้ plotgrids_new.ncl ซึ่งอยู่ในโฟลเดอร์ ./WPS/util/ โดยมีขั้นตอนดังนี้
คัดลอก plotgrids_new.ncl มาอยู่ในโฟลเดอร์เดียวกันกับ namelist.wps
เลือกรูปแบบการแสดงผลโดยแก้ไขตัวแปร type ใน plotgrids_new.ncl
สำหรับ X11: เลือก type = “x11“ (ค่าตั้งต้น)
สำหรับ PDF: เลือก type = ”pdf”
ใช้คำสั่ง ncl plotgrids_new.ncl
สำหรับ X11 ต้องเปิดใช้งาน X11 Forwarding หน้าต่างแสดงผลถึงขึ้นมาได้
สัญลักษณ์ ; คือ การ comment ใน NCL
cd WPS cp ./util/plotgrids_new.ncl . # Copy the script to the same directory as namelist.wps nano plotgrids_new.ncl # Comment/uncomment the variable "type" # For X11 (default), X11 Forwarding must be enabled. ncl plotgrids_new.ncl # Run the script
2.3 WPS job submission script
#!/bin/bash #SBATCH -p devel #SBATCH -N 1 --ntasks-per-node=5 #SBATCH -t 02:00:00 #SBATCH -J WPS #SBATCH -A thaisc module purge module load HDF5/1.10.5-iimpi-2019b module load JasPer/1.900.1-intel-2019b module load libpng/1.6.37-GCCcore-10.3.0 export JASPERINC=$EBROOTJASPER/include export JASPERLIB=$EBROOTJASPER/lib export NETCDF=/tarafs/data/project/thaisc/skanoksi/Libraries/NET$ export LD_LIBRARY_PATH=${NETCDF}/lib:${LD_LIBRARY_PATH} export PATH=${NETCDF}/bin:${PATH} ulimit -s unlimited srun -n4 --exclusive ./geogrid.exe & srun -n1 --exclusive ./ungrib.exe & wait srun -n4 --exclusive ./metgrid.exe
สามารถรวม Job submission script ได้
2.4 Submit job
3. Run WRF
3.1 Edit ./em_real/namelist.input
3.2 WRF job submission script
#!/bin/bash #SBATCH -p devel #SBATCH -N 1 --ntasks-per-node=40 #SBATCH -t 02:00:00 #SBATCH -J WRF #SBATCH -A thaisc module purge module load HDF5/1.10.5-iimpi-2019b module load JasPer/1.900.1-intel-2019b module load libpng/1.6.37-GCCcore-10.3.0 export JASPERINC=$EBROOTJASPER/include export JASPERLIB=$EBROOTJASPER/lib export NETCDF=/tarafs/data/project/thaisc/skanoksi/Libraries/NET$export LD_LIBRARY_PATH=${NETCDF}/lib:${LD_LIBRARY_PATH} export PATH=${NETCDF}/bin:${PATH} ulimit -s unlimited srun ./real.exe & wait srun ./wrf.exe
3.3 Submit job
Contact Us
ThaiSC support service : thaisc-support@nstda.or.th