ggg_pp = gpmdpr_pp_b * res_bin
ggg_z = gpmdpr_z_b * res_bin
ggg_top = dpr_top_b * res_bin
#name = '/automount/ags/velibor/gpmdata/dumpdata/gpm_dpr_rado_all/' \
# 'dprrado_'+sc+'/dprrado_'+ sc + str(gpm_zeit)
name = '/automount/ags/velibor/gpmdata/dumpdata/gpm_dpr_rado_all_int25/' \
'dprrado_'+sc+'/dprrado_'+ sc + str(gpm_zeit)
w2h(name, gpm_x, gpm_y, ggg_pp, ggg_bbh, ggg_bbw, ggg_typ,
ggg_phase, ggg_pp, ggg_z, ggg_top, rrr, rrr2)
print '/automount/ags/velibor/gpmdata/dumpdataV7/SKILL_' + sc + str(
gpm_zeit)
del (gprof_lat, gprof_lon, gpmdpr_pp, res_bin, rrr, ggg_pp,
ggg_bbh, ggg_bbw, ggg_phase, ggg_typ, rwdata, x, y, gpm_x,
gpm_y, gpm_time, xy, grid_gpm_xy, grid_xy, rrr2, rn,
rwattrs, result, pfad, pfad_radolan, ht, m, d, ye, mt,
year, dpr_pp, dpr_bbh, dpr_bbw, dpr_phase, dpr_typ,
dpr_pp_b, dpr_bbh_b, dpr_bbw_b, dpr_phase_b, dpr_typ_b,
ggg_z, gpmdpr_pp_b, gpmdpr_z_b, r_pro2, pfad2,
pfad_radolan2, rwattrs2, rwdata2, rw_filename,
rw_filename2, result2, name, dpr_top, dpr_top_b, ggg_top)
except:
print 'Datei ' + pfad_radolan + ' nicht gefunden!'
from pcc import melde_dich
melde_dich('Das Program ist fertig!')
SS['M'],SS['F'],SS['C'],SS['HR'],SS['POD'],
SS['FAR'],SS['BID'],SS['HSS'],SS['bias'],
SS['RMSE'],np.nanmean(ggg),np.nanmean(rrr),
np.nanmedian(ggg),np.nanmedian(rrr)])
del(text, slope, intercept, r_value, p_value, std_err, line,
width, height, ell,maske,SS)
except:
pass
#plt.show()
del(gprof_lat, gprof_lon, gprof_pp, res_bin, rrr, ggg, rwdata, x, y,
gpm_x, gpm_y, gpm_time, xy,grid_gpm_xy, grid_xy, mask, rn,
rwattrs, result, pfad, pfad_radolan, ht, m, d, ye ,mt, year )
from pcc import melde_dich
melde_dich('Das Program pcc_alldproverpass_stats_rr.py ist fertig!')
f.close()
#meinftp.retrlines('LIST')
daten=[] #Initialisierung einer Liste (leere Liste)
meinftp.dir(daten.append)
print
print 'Ort und Name der lokalen Datei: ' + directory_local + dpr_overpass
print
file = open(directory_local+dpr_overpass, 'wb')
print 'Download: ftp-Server: ' + directory +dpr_overpass
meinftp.retrbinary('RETR '+dpr_overpass, file.write)
print
print 'Die lokale Datei ' + directory_local+dpr_overpass +' wird geschlossen.'
file.close()
print meinftp.quit()
print
print 'Die FTP-Verbindung wurde von mir getrennt.'
melde_dich('Das Program ist fertig!')
import glob
import wradlib as wrl
from osgeo import osr
from pcc import get_time_of_gpm
from pcc import cut_the_swath
from satlib import write2hdf as w2h
import time
sekunden = 60*60*5
time.sleep(sekunden)
print 'start nach ',str(sekunden), ' Sekunden.'
from pcc import melde_dich
melde_dich('Das Program REF STARTET!')
# Ref.Threshold nach RADOLAN_Goudenhoofdt_2016
TH_ref = 0.1
scc = ['NS', 'HS', 'MS']
pfad = ('/automount/ags/velibor/gpmdata/dpr/*.HDF5')
pfad_gpm = sorted(glob.glob(pfad))
print 'Es sind ', len(pfad_gpm), ' vorhanden!'
for iii in scc:
sc = iii
#len(pfad_gpm)