https://teslamotorsclub.com/tmc/threads/model-3-specification-speculation.69796/

Model3.py

# -*- coding: utf-8 -*-
"""
Created on Sat Apr 30 12:02:50 2016

@author: darin
"""

import requests

import os, sys

import urllib

def openfile(file):
    with open(file,'r') as openfile:
        readfile=openfile.read()
    openfile.close()
    return readfile
    
def adrag(speed, cda):
    return 0.5 * speed * speed * 1.225 * cda

txtflag = 0
homedir = os.path.dirname(os.path.abspath(sys.argv[0]))
    
for root, dirnames, filenames in os.walk(homedir):
    for thefile in filenames:
        if '.txt' in thefile:
            txtflag = 1
            
if txtflag == 0:
    linklist = []
    
    r = requests.get('https://www.epa.gov/vehicle-and-fuel-emissions-testing/dynamometer-drive-schedules')
    pref = 'https://www.epa.gov/'
    
    for i in r.text.split('\n'):
        for j in i.split('><'):
            if '.txt' in j:
                linklist.append(j.split('"')[1])
    for link in linklist:
        urllib.urlretrieve(pref + link, link.split('/')[-1])
        

#.23 Cd = .525m^2 CdA
#Model 3 55 weight = 1815, height = 55", width = 71", Cd =.21, CdA=.479
#Model 3 55 weight = 1905, height = 55", width = 71", Cd =.21
#Model 3 55 weight = 2177, height = 56", width = 77", Cd =.24
#Compared to S70, EPA for the 55 will be ~230 miles and ~300 miles for the 75.

vlist = [('Model_3_55',1815,.525), ('Model_3_75',1905,.525), ('Model_S_70D',2180, .604), ('Model_S_90D',2279,.604)]
vreslist = []

for car in vlist:
    lstring = car[0] + '_list = []'
    exec(lstring)


for root, dirnames, filenames in os.walk(homedir):
    for thefile in filenames:
        if '.txt' in thefile:
            tfile = openfile(root + '/' + thefile)
            for car in vlist:
                clistname = car[0] + '_list'
                ptot = 0.0
                vw = car[1]
                cda = car[2]
                for sec in tfile.split('\n'):
                    fair = 0.0
                    froll = 0.0
                    ftot = 0.0
                    try:
                        vel = float(sec.split('\t')[1]) * 0.44704
                        fair = adrag(vel, cda)
                        if vel != 0.0:
                            froll = float(vw) * 0.0089 #from Roadster 3 article
                        ftot = froll + fair
                        ptot = ptot + (ftot * vel) / 3600.0
                        
                    except:
                        pass
                if ptot > 0:
                    acommand = clistname + ".append(('" + thefile + "'," + str(int(ptot)) +'))'
                    exec(acommand)

#S70 unadjusted HWFET = 135, FTP75 = 133.7, GGE =33.7kWh
s70hwycol = 33700/135.0
s70ftpcol = 33700/133.7

zippedresult = eval('[i for i in zip(' + str([eval(i[0]+'_list') for i in vlist])[1:][:-1] + ')]')

for i in zippedresult:
    print i
    if i[0][0] == 'hwycol.txt':
        print '55 Range is about ' + str((float(i[2][1])/float(i[0][1]))*245.0*(55.0/70.0))
        print '75 Range is about ' + str((float(i[2][1])/float(i[1][1]))*245.0*(75.0/70.0))
    if i[0][0] == 'ftpcol.txt':
        print '55 Range is about ' + str((float(i[2][1])/float(i[0][1]))*245.0*(55.0/70.0))
        print '75 Range is about ' + str((float(i[2][1])/float(i[1][1]))*245.0*(75.0/70.0))
    print '\n'

'''
#guesstimate for background power consumption watts, 50W
for i in zip(Model_3_55_list,Model_S_70D_list):
    if i[0][1] != 0:
        if i[0][0] == 'ftpcol.txt':
            whmile = (s70ftpcol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70ftpcol/whmile) * 243.0 * (55.0/70.0)), i[0][1], i[1][1]
        if i[0][0] == 'hwycol.txt':
            whmile = (s70hwycol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70hwycol/whmile) * 246.0 * (55.0/70.0)), i[0][1], i[1][1]

            
for i in zip(Model_3_75_list,Model_S_70D_list):
    if i[0][1] != 0:
        if i[0][0] == 'ftpcol.txt':
            whmile = (s70ftpcol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70ftpcol/whmile) * 243.0 * (75.0/70.0)), i[0][1], i[1][1]
        if i[0][0] == 'hwycol.txt':
            whmile = (s70hwycol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70hwycol/whmile) * 246.0 * (75.0/70.0)), i[0][1], i[1][1]
            
for i in zip(Model_S_90D_list,Model_S_70D_list):
    if i[0][1] != 0:
        if i[0][0] == 'ftpcol.txt':
            whmile = (s70ftpcol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70ftpcol/whmile) * 243.0 * (90.0/70.0)), i[0][1], i[1][1]
        if i[0][0] == 'hwycol.txt':
            whmile = (s70hwycol - 50.0)*(float(i[0][1])/float(i[1][1])) + 50.0
            print i[0][0], 'wh/mile ' + str(whmile), 'Range ' + str((s70hwycol/whmile) * 246.0 * (90.0/70.0)), i[0][1], i[1][1]
'''