from numpy import float32 import os site='AB01' # get data [lati,loni,on,bd]=getemolt_latlon(site)# extracts lat/lon based on site code [lati,loni]=dm2dd(lati,loni)#converts decimal-minutes to decimal degrees ''' dept=[0,5] (obs_dt,obs_temps,obs_salt)=getobs_tempsalt(site, input_time=[dt.datetime(2004,6,7),dt.datetime(2004,6,24)], dep=dept) obs_tempslist=[] for i in range(len(obs_temps)): obs_tempslist.append(f2c(obs_temps[i])) dfsur=pd.DataFrame(obs_tempslist,index=obs_dt) ''' dept=[bd[0]-0.4*bd[0],bd[0]+0.35*bd[0]] (obs_dt,obs_tempb,obs_salt)=getobs_tempsalt(site, input_time=[dt.datetime(2004,6,7),dt.datetime(2004,6,24)], dep=dept) obs_tempblist=[] for i in range(len(obs_tempb)): obs_tempblist.append(f2c(obs_tempb[i])) dfbot=pd.DataFrame(obs_tempblist,index=obs_dt) #dfsur=dfsur.ix[dfsur.index.hour==0] #dfbot=dfbot.ix[dfbot.index.hour==0] # get model modtso=getFVCOM_bottom_tempsalt_netcdf(lati,loni,dt.datetime(2004,6,7),dt.datetime(2004,6,24),layer=range(45),vname='temp') #modtso=modtso.ix[modtso.index.hour==0] for i in range(len(dfbot)): fig=plt.figure() plt.plot(modtso.values[i],range(0,-45,-1),'g*') # plt.plot([float32(dfsur[0][i]),float32(dfbot[0][i])],[0,-44],'r*') plt.plot([float32(dfbot[0][i])],[-44],'r*')
# -*- coding: utf-8 -*- """ Created on Mon Mar 25 13:20:47 2013 Routine to look at model temperature profiles ar a particular eMOLT location where we had both surf & bot sensors @author: jmanning """ import datetime as dt from matplotlib import pyplot as plt from getdata_yw import getobs_tempsalt from getdata import getemolt_latlon from models_yw import getFVCOM_bottom_tempsalt_netcdf from conversions import dm2dd site='AG01' # get data [lati,loni,on,bd]=getemolt_latlon(site)# extracts lat/lon based on site code [lati,loni]=dm2dd(lati,loni)#converts decimal-minutes to decimal degrees dept=[0,5] (obs_dt,obs_temps,obs_salt)=getobs_tempsalt(site, input_time=[dt.datetime(2006,9,10),dt.datetime(2006,9,11)], dep=dept) dept=[bd[0]-0.25*bd[0],bd[0]+0.25*bd[0]] (obs_dt,obs_tempb,obs_salt)=getobs_tempsalt(site, input_time=[dt.datetime(2006,9,10),dt.datetime(2006,9,11)], dep=dept) # get model for k in range(44): modtso=getFVCOM_bottom_tempsalt_netcdf(lati,loni,dt.datetime(2006,9,10),dt.datetime(2006,9,11),layer=k,vname='temp') plt.plot(modtso[0],-k,'g*') plt.plt([obs_temps[0],obs_tempb[0]],[0,-bd[0]],'r*')
import datetime as dt from matplotlib import pyplot as plt from getdata_yw import getobs_tempsalt from getdata import getemolt_latlon from models_yw import getFVCOM_bottom_tempsalt_netcdf from conversions import dm2dd site = 'AG01' # get data [lati, loni, on, bd] = getemolt_latlon(site) # extracts lat/lon based on site code [lati, loni] = dm2dd(lati, loni) #converts decimal-minutes to decimal degrees dept = [0, 5] (obs_dt, obs_temps, obs_salt) = getobs_tempsalt( site, input_time=[dt.datetime(2006, 9, 10), dt.datetime(2006, 9, 11)], dep=dept) dept = [bd[0] - 0.25 * bd[0], bd[0] + 0.25 * bd[0]] (obs_dt, obs_tempb, obs_salt) = getobs_tempsalt( site, input_time=[dt.datetime(2006, 9, 10), dt.datetime(2006, 9, 11)], dep=dept) # get model for k in range(44): modtso = getFVCOM_bottom_tempsalt_netcdf(lati, loni, dt.datetime(2006, 9, 10), dt.datetime(2006, 9, 11),