Steps for running the fortran model (https://opensource.ncsa.illinois.edu/bitbucket/projects/VBD/repos/mir_model_f/browse)
Original Code and setup on idea.
Steps in run process (cron)
- The 3 routines that are run on the original server
- Get forecast
/home/nan/WNV_WX/Not_smooth/run.get_forecast
gets csv of forcast
copies (renames) to ILCD_<year><month><day> in same directory
- then copies all IL* files to ../MIR_Model, ../MIR_cmb, and MIR_indiv
Get prism data
/home/nan/WNV_WX/MIR_model/run.cp_prism
get prism data
/home/mtimlin/estSTCD/estCDwestnile/estCD.csh
- set d1 to 19810101 and d2 to `date +"%Y%m%d"` the get 4 files:
curl -o estCD.dir/estCDp --data 'params={"state":"IL","sdate":"'$d1'","edate":"'$d2'","grid":"21","elems":[{"name":"pcpn","area_reduce":"climdiv_mean"}]}' http://data.rcc-acis.org/GridData
curl -o estCD.dir/estCDn --data 'params={"state":"IL","sdate":"'$d1'","edate":"'$d2'","grid":"21","elems":[{"name":"mint","area_reduce":"climdiv_mean"}]}' http://data.rcc-acis.org/GridData
curl -o estCD.dir/estCDx --data 'params={"state":"IL","sdate":"'$d1'","edate":"'$d2'","grid":"21","elems":[{"name":"maxt","area_reduce":"climdiv_mean"}]}' http://data.rcc-acis.org/GridData
curl -o estCD.dir/estCDt --data 'params={"state":"IL","sdate":"'$d1'","edate":"'$d2'","grid":"21","elems":[{"name":"avgt","area_reduce":"climdiv_mean"}]}' http://data.rcc-acis.org/GridData
- set d1 to 19810101 and d2 to `date +"%Y%m%d"` the get 4 files:
copy prism data
from timlin directory cd /home/mtimlin/estSTCD/estCDwestnile/estCD.dir/110x where x=[1,9] where x is the climate division
- to WNV_WX/Not_smooth/110x
Run Model
/home/nan/WNV_WX/MIR_model/run.cmp_n_nweek
- run the model
- ./cmp_n_nweek 11 0x 2017 where x=[1,9]
- rename plots
- copy plots (rename) (110x_2017_plot → 110x_YRcur_plot) to webapp directory
- run the model
Step through fortran code
Lines Function Important variables Description Important variables names 38-185 initialize variables 188-191 characters to int ayear→iyear CD->icd 194 ? c=',' 196-202 create 3 20x360 arrays initialized as: t,p = 999.9 d =0.0 209-212 read week_days_2016.csv, input days of week for year wmonth(i),wday(i),wweek(i),wyear,wdate,wdow 216-231 create output files normals MIR plot, (move to ./output) 241-423 compute normals