Steps for running the fortran model (https://opensource.ncsa.illinois.edu/bitbucket/projects/VBD/repos/mir_model_f/browse)
/home/nan/WNV_WX/Not_smooth/run.get_forecast
gets csv of forcast
copies (renames) to ILCD_<year><month><day> in same directory
Get prism data
/home/nan/WNV_WX/MIR_model/run.cp_prism
get prism data
/home/mtimlin/estSTCD/estCDwestnile/estCD.csh
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
copy prism data
from timlin directory cd /home/mtimlin/estSTCD/estCDwestnile/estCD.dir/110x where x=[1,9] where x is the climate division
Run Model
/home/nan/WNV_WX/MIR_model/run.cmp_n_nweek
Lines | Function | Important variables Description | Important variables names | Issues |
---|---|---|---|---|
38-185 | initialize variables | |||
188-191 | characters to int | ayear→iyear CD->icd | ||
194 | ? | c=',' | ||
196-202 | create 3 20x360 arrays | initialized as: t=temperature p=precipitation d=? | 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 | |||
243-247 | set some time period length values | bmon = 12 | ||
250-293 | calculate normal for temp and precip | running to the end of the file without indexes every updating |