#!/bin/bash # ----------------------------INPUT--------------------------------------- # Seting year and month (initial and final) for running year='2018' monthI='01' monthF='12' # Folder of global meteorological model data WPS_DATA='/home/nobre/Build_WRF/DATA' # Folder of WPS files WPS_FOLDER='/home/nobre/Build_WRF/WPS' #WPS_FOLDER='/home/nobre/Build_WRF/AERMOD_robao' # WRF folder WRF_FOLDER='/home/nobre/Build_WRF/WRF/test/em_real' #WRF_FOLDER='/home/nobre/Build_WRF/AERMOD_robao' # Number of processors NOP='1' # Seting runing days jan=(1 10 20 28) #FALTA AUTOMATIZAR O DOMÍNIO EM TODOS OS NAMELISTS # ==================================Namelist.wps============================================================== # &geogrid parent_id1='1' parent_id2='1' parent_id3='1' parent_grid_ratio1='1' parent_grid_ratio2='3' parent_grid_ratio3='3' i_parent_start1='1' i_parent_start2='20' i_parent_start3='6' j_parent_start1='1' j_parent_start2='7' j_parent_start3='6' e_we1='29' e_we2='10' e_we3='19' e_sn1='18' e_sn2='10' e_sn3='13' geog_data_res1="'30s'" geog_data_res2="'30s'" geog_data_res3="'30s'" dx1='30000' dy1='30000' map_proj1="'mercator'" ref_lat1='-27.388' ref_lon1='-50.524' truelat11='-27.388' truelat21='0' stand_lon1='-50.524' geog_data_path='/home/nobre/Build_WRF/WPS_GEOG', ref_x1='14.5' ref_y1='9.0' # ===================================Namelist.imput========================================================== # &domains time_step1='180' time_step_fract_num1='0' time_step_fract_den1='1' max_dom1='3' e_vert1='33' e_vert2='33' e_vert3='33' p_top_requested1='5000' num_metgrid_levels1='32' num_metgrid_soil_levels1='4' dx1='30000' dx2='10000' dx3='10000' dy1='30000' dy2='10000' dy3='10000' grid_id1='1' grid_id2='2' grid_id3='3' parent_time_step_ratio1='1' parent_time_step_ratio2='3' parent_time_step_ratio3='3' feedback1='1' smooth_option1='0' # &physics mp_physics1='10' mp_physics2='10' mp_physics3='10' ra_lw_physics1='1' ra_lw_physics2='1' ra_lw_physics3='1' ra_sw_physics1='1' ra_sw_physics2='1' ra_sw_physics3='1' radt1='30' radt2='30' radt3='30' sf_sfclay_physics1='1' sf_sfclay_physics2='1' sf_sfclay_physics3='1' sf_surface_physics1='2' sf_surface_physics2='2' sf_surface_physics3='2' bl_pbl_physics1='1' bl_pbl_physics2='1' bl_pbl_physics3='1' bldt1='0' bldt2='0' bldt3='0' cu_physics1='1' cu_physics2='1' cu_physics3='1' cudt1='5' cudt2='5' cudt3='5' isfflx1='1' ifsnow1='0' icloud1='1' surface_input_source1='1' num_soil_layers1='4' sf_urban_physics1='0' sf_urban_physics2='0' sf_urban_physics3='0' maxiens1='1' maxens1='3' maxens21='3' maxens31='16' ensdim1='144' NUM_LAND_CAT1='24' # &fdda # &dynamics hybrid_opt1='1' w_damping1='0' diff_opt1='1' diff_opt2='1' diff_opt3='1' km_opt1='4' km_opt2='4' km_opt3='4' diff_6th_opt1='0' diff_6th_opt2='0' diff_6th_opt3='0' diff_6th_factor1='0.12' diff_6th_factor2='0.12' diff_6th_factor3='0.12' base_temp1='290.' damp_opt1='0' zdamp1='5000' zdamp2='5000' zdamp3='5000' dampcoef1='0.2' dampcoef2='0.2' dampcoef3='0.2' khdif1='0' khdif2='0' khdif3='0' kvdif1='0' kvdif2='0' kvdif3='0' non_hydrostatic1='.true.' non_hydrostatic2='.true.' non_hydrostatic3='.true.' moist_adv_opt1='1' moist_adv_opt2='1' moist_adv_opt3='1' scalar_adv_opt1='1' scalar_adv_opt2='1' scalar_adv_opt3='1' # &bdy_control spec_bdy_width1='5' spec_zone1='1' relax_zone1='4' specified1='.true.' specified2='.false.' specified3='.false.' nested1='.false.' nested2='.true.' nested3='.true.' # &grib2 # &namelist_quilt nio_tasks_per_group1='0' nio_groups1='1' # =========================================================================================================== echo "" echo "" echo "######################## Starting program runningWRFv2 ################################" echo " ___ ______ ______ _ | | | ___| | ___ | / \ _ | | | | | | | | / | \ | \_/| | |___ | |___ | |\|\| / /_\ \ | __/ |_____| |_____| |___\_\ /_/ \_\ |_| " echo " Laboratório de Controle da Qualidade do Ar " echo "" echo " Universidade Federal de Santa Catarina - Brasil" echo "" echo " email: leonardo.hoinaski@ufsc.br" echo "" echo " last update: 07-10-2019 " echo echo "#######################################################################################" echo "" echo "Initializing in $year - $monthI - $monthF" echo "" #==============================LOOP OVER EACH MONTH======================================= echo "=======================================================================================" k=0 for i in $(seq $monthI $monthF); do echo "" echo "############################ Month number: $i ############################" echo h=$((i + 1)) rm -rf $WPS_DATA/MetInUse mkdir $WPS_DATA/MetInUse # -------------------------LINKING MET FILES------------------------ if [ $i -lt 10 ]; then ln -sf $WPS_DATA/gfsanl_4_${year}0${i}* $WPS_DATA/MetInUse/. fi if [ $i -lt 9 ]; then ln -sf $WPS_DATA/gfsanl_4_${year}0${h}* $WPS_DATA/MetInUse/. fi if [ $i -gt 9 ]; then ln -sf $WPS_DATA/gfsanl_4_${year}${i}* $WPS_DATA/MetInUse/. fi if [ $i -gt 8 ]; then ln -sf $WPS_DATA/gfsanl_4_${year}${h}* $WPS_DATA/MetInUse/. fi if [ $i -eq 12 ]; then rm -rf $WPS_DATA/MetInUse/gfsanl_4_${year}${h}* fi echo "Link of met file is done" echo "" # -------------------------LINKING MET IS DONE------------------------ for j in $(seq 0 2); do k=$((j + 1)) echo " ======== > Day ${jan[$j]} to ${jan[$k]} ============================" cd $WPS_FOLDER if [ $i -lt 10 ]; then echo "month smaller than 10" if [ $j -lt 1 ]; then echo "day smaller than 10" echo "start_date = '${year}-0${i}-0${jan[$j]}_00:00:00', '${year}-0${i}-0${jan[$j]}_00:00:00', '${year}-0${i}-0${jan[$j]}_00:00:00'," echo "end_date = '${year}-0${i}-${jan[$k]}_00:00:00', '${year}-0${i}-${jan[$k]}_00:00:00', '${year}-0${i}-${jan[$k]}_00:00:00'," cat > namelist.wps < namelist.wps < namelist.wps < namelist.wps <& log.geogrid echo "Linking Vtable" ln -sf $WPS_FOLDER/ungrib/Variable_Tables/Vtable.GFS Vtable echo "Runnig ungrib.exe" $WPS_FOLDER/ungrib.exe >& log.ungrib echo "Runnig metgrid.exe" $WPS_FOLDER/metgrid.exe >& log.metgrid echo "" cd $WRF_FOLDER echo " ** Writing WRF namelist.input file ** " echo "" rm -rf namelist.input cat > namelist.input < Day 28 to 01 ============================' cd $WPS_FOLDER if [ $i -lt 9 ]; then echo "Month smaller than 10" echo "start_date = '${year}-0${i}-28_00:00:00', '${year}-0${i}-28_00:00:00', '${year}-0${i}-28_00:00:00', " echo "end_date = '${year}-0${h}-01_00:00:00', '${year}-0${h}-01_00:00:00', '${year}-0${h}-01_00:00:00', " # Wrinting WPS file from day 31 to 01 cat > namelist.wps < namelist.wps < namelist.wps < namelist.wps <& log.geogrid echo "Linking Vtable" ln -sf $WPS_FOLDER/ungrib/Variable_Tables/Vtable.GFS Vtable echo "Runnig ungrib.exe" $WPS_FOLDER/ungrib.exe >& log.ungrib echo "Runnig metgrid.exe" $WPS_FOLDER/metgrid.exe >& log.metgrid cd $WRF_FOLDER rm -rf namelist.input if [ $i -lt 12 ]; then echo " ** Writing WRF namelist.input files ** " echo "" cat > namelist.input < namelist.input <