################################################# Katie and I have gone through the sample data you provided us a week ago from your automated collection system. The format of the data file is good - we just need adjustments in the columns included in the file. Katie put together a spreadsheet (attached) tabulating the data we use and the data your system is collecting so far. Column H has notes for you. I had told you that we would not need GSOD data, but Katie has indicated that it should be used instead of CIMIS for the Madera meteorology station. I’m following up with her to find out why. For meteorology input files we actually use C for temperatures even though for measured data and output we use F, so please collect MET data in C. We don’t need vapor pressure data (just dewpoint temperature), and the model isn’t currently equipped to use solar radiation data so you can skip those. At the bottom of the attached list are four stations for which we don’t enter downloaded data directly but rather calculate their values from other downloaded files. Please let me know if you have the capability to do the calculations listed or if we should develop a simple processor to do so. ############################################## http://waterservices.usgs.gov/nwis/iv/?format=json,1.1&sites=11289000,11289500,11274630,11255575&period=P1D¶meterCd=00060 "http://et.water.ca.gov/api/data?appKey=67847b45-a48e-461f-afb6-2a0b72f7635e&targets=56,92,188,148,70,124,7,80,71&startDate="+start_date+"&endDate="+end_date+"&dataItems=day-precip,day-air-tmp-min,day-sol-rad-avg,day-dew-pnt,day-vap-pres-avg,day-wind-spd-avg&unitOfMeasure=E" opener = urllib2.build_opener() req = urllib2.Request('http://localhost:8000', data=json.dumps({'a': 'b'}), headers={'Content-Type': 'application/json'}) response = opener.open(req) TRP,MST,RIP,OCL,ORE,SMN,NEW,SCL,VNS,FFB,SSH,SJS,MOD DMC,SJS,MST,MOD,RIP,OCL,ORE,SMN,SCL,VER,MRB,HSP,ING,WSW,MON,MSR,FFB,SSH,MSG cimisData = cimis_json['Data']['Providers'][0]['Records'] cemisArrays = {} for dayValues in cimisData: print dayValues['Date'] + ',' + dayValues['DayAirTmpMin'] { u'DayPrecip': { u'Qc': u'', u'Unit': u'(in)', u'Value': u'0.01' }, u'DayWindSpdAvg': { u'Qc': u'', u'Unit': u'(MPH)', u'Value': u'2.5' }, u'ZipCodes': u'93622, 93665, 93620', u'Julian': u'303', u'DaySolRadAvg': { u'Qc': u'', u'Unit': u'(Ly/day)', u'Value': u'232' }, u'Standard': u'english', u'DayDewPnt': { u'Qc': u'', u'Unit': u'(F)', u'Value': u'48.9' }, u'Station': u'7', u'DayVapPresAvg': { u'Qc': u'', u'Unit': u'(mBars)', u'Value': u'11.8' }, u'Date': u'2014-10-30', u'Scope': u'daily', u'DayAirTmpMin': { u'Qc': u'', u'Unit': u'(F)', u'Value': u'49.5' } } 'DayPrecip','DayWindSpdAvg','DaySolRadAvg','DayDewPnt','DayVapPresAvg','DayAirTmpMin' cemisStations = [56:'Los Banos',92:'Kesterson',188:'Madera II',148:'Merced',70:'Manteca',124:'Panoche',7:'Firebaugh/Telles',80:'Fresno State',71:'Modesto']