Events

import sys
import os
import matplotlib
matplotlib.use('Agg') # Must be before importing matplotlib.pyplot or pylab!
import matplotlib.pyplot as plt
import matplotlib.pyplot as P
from mpl_toolkits.basemap import Basemap,shiftgrid,maskoceans
from matplotlib.path import Path
from matplotlib.patches import PathPatch
import shapefile
#from netCDF4 import Dataset
import numpy as N
from netCDF4 import Dataset, num2date,date2num
import datetime
from scipy import signal,stats
import matplotlib.colors as colors
#from datetime import datetime
from matplotlib.offsetbox import AnchoredText
#from scipy.interpolate import griddata
#import scipy
import matplotlib as mpl


plt.style.use('seaborn-talk')

params = {      'legend.fontsize': 8,\
                'font.family': 'serif',\
                }
plt.rcParams.update(params)

def ma(a,n=3):
    ret=N.cumsum(a,dtype=float)
    ret[n:]=ret[n:]-ret[:-n]
    return ret[n-1:]/n

file = '../../data/gwl/gwlneudatum.dat'
ta = N.genfromtxt(file,usecols=(0),skip_header=1,dtype='i')
mo = N.genfromtxt(file,usecols=(1),skip_header=1,dtype='i')
ja = N.genfromtxt(file,usecols=(2),skip_header=1,dtype='i')
gw = N.genfromtxt(file,usecols=(3),skip_header=1,dtype='str')

nd = len(gw)

dates = []

for d in range(nd):

    tmp = '%4i-%02i-%02i'%(ja[d],mo[d],ta[d])

    dates.append(tmp)

dates = N.array(dates)

city = {

   'Amman':[35.930359,31.963158],
   'Aqaba':[35.00778,29.52667],
   'WadiMusa':[35.480125,30.321635],

}

datas = ['UERRA','ERA5','W5E5','CHIRPS']

dummy = {}

for c in city:
    for data in datas:

        if(data=='UERRA'): file = '../../data/ClimateExplorer/tp_daily_uerra_19810101-20181231_rm.cdf'
        if(data=='CHIRPS'): file = '../../data/ClimateExplorer/chirps-v2.0.1981-2022.days_p25.cdf'
        if(data=='ERA5'): file = '../../data/ClimateExplorer/era5_tp_daily_af_34-40E_29-34N_su.cdf'
        if(data=='W5E5'): file = '../../data/ClimateExplorer/pr_W5E5v2.0_19810101-20191231.cdf'

        nc = Dataset(file,'r')
        lon = N.array(nc.variables['lon'][:])
        lat = N.array(nc.variables['lat'][:])
        dat = N.array(nc.variables['tp'][:])
        tim = nc.variables['time']

        #print (data,c,dat.shape,N.max(dat),N.min(dat))

        #dat[dat>1000] = N.nan

        tim = num2date(tim[:],units=tim.units,calendar=tim.calendar) 

        jj = []
        mm = []
        dd = []

        for it in tim:

            jj.append(it.year)
            mm.append(it.month)
            dd.append(it.day)

        jj = N.array(jj)
        mm = N.array(mm)
        dd = N.array(dd)

        nd = len(dd)

        nc.close()

        lo,la = N.meshgrid(lon,lat)

        distance = (lo-city[c][0])**2 + (la-city[c][1])**2
        iy,ix = N.where(distance==distance.min())

        name = '%s_%s'%(data,c)

        dummy[name] = []

        datum = 'DATE_%s'%(data)

        dummy[datum] = []

        for d in range(nd):

            dummy[name].append(dat[d,iy,ix])

            date = '%04i-%02i-%02i'%(jj[d],mm[d],dd[d])

            dummy[datum].append(date)

print (dummy.keys())

os.system('rm ./includes/ncep/*.png')
os.system('rm ./includes/uerra/*.png')
os.system('rm ./includes/chirps/*.png')
os.system('rm ./includes/era5/*.png')
os.system('rm ./includes/w5e5/*.png')

for c in city:

    f = open('./includes/%s.md'%c,'w')

    f.write('DATE|UERRA|ERA5|W5E5|GWL|Z500|UERRA|ERA5|W5E5\n')
    f.write('---|---|---|---|---|---|---|---|---\n')

    name1 = 'UERRA_%s'%(c)
    name2 = 'ERA5_%s'%(c) 
    name3 = 'CHIRPS_%s'%(c)
    name4 = 'W5E5_%s'%(c)

    for d in dummy['DATE_ERA5']:

        if(int(d[0:4])>=2000):

             id0 = N.where(dates==d)[0]
             id1 = N.where(N.array(dummy['DATE_UERRA'])==d)[0]
             id2 = N.where(N.array(dummy['DATE_ERA5'])==d)[0]
             id3 = N.where(N.array(dummy['DATE_CHIRPS'])==d)[0]
             id4 = N.where(N.array(dummy['DATE_W5E5'])==d)[0] 

             if(len(id1)==1): value1 = dummy[name1][id1[0]]
             else: value1 = -999.9

             if(value1<0.): value1 = -999.9 

             if(len(id2)==1): value2 = dummy[name2][id2[0]]
             else: value2 = -999.9

             if(len(id3)==1): value3 = dummy[name3][id3[0]]
             else: value3 = -999.9

             if(len(id4)==1): value4 = dummy[name4][id4[0]]
             else: value4 = -999.9

             if((value1>10)|(value2>10)|(value4>10)):

                print (d)

                if(len(id0)==1): os.system('cp ../../data/png/ncep/%s.png ./includes/ncep/'%d)
                os.system('cp ../../data/png/uerra/%s.png ./includes/uerra/'%d) 
                os.system('cp ../../data/png/era5/%s.png ./includes/era5/'%d)
                os.system('cp ../../data/png/chirps/%s.png ./includes/chirps/'%d)
                os.system('cp ../../data/png/w5e5/%s.png ./includes/w5e5/'%d)

                image0 = '![](./ncep/%s.png){: style="height:50px;width:120px"}'%d
                image1 = '![](./uerra/%s.png){: style="height:50px;width:100px"}'%d           
                image2 = '![](./era5/%s.png){: style="height:50px;width:100px"}'%d
                image3 = '![](./chirps/%s.png){: style="height:50px;width:100px"}'%d
                image4 = '![](./w5e5/%s.png){: style="height:50px;width:100px"}'%d

                if(len(id0)==1): gwl = gw[id0[0]] 
                else: gwl = 'N.N.'

                f.write('%s|%05.1f|%05.1f|%05.1f|%s|%s|%s|%s|%s\n'%(d,value1,value2,value4,gwl,image0,image1,image2,image4))

    f.close()