예제 #1
0
def make_nprof_timeseries(fy,ly,Area_id):
	import AQC
	import numpy as np
	import subroutine
	N=np.ones([ly-fy+1,12])
	for year in range(fy,ly+1):
		for month in range(1,13):
			n=AQC.specified_NPROF(year,month,Area_id)
			N[year-fy,month-1]=len(n)

	ym=subroutine.ym_timeseries(fy,ly)
	N=subroutine.to_oneline(N)
	return ym,N
예제 #2
0
def specified_LD(year,month,dt,Area_id): # 領域で絞り込んだプロファイルのみ層厚計算
	import subroutine
	import AQC
	import numpy as np
	import LD
	temp,salt,pres,lon,lat=AQC.get_data(year,month)
	N_PROF=lon.size
	n=AQC.specified_NPROF(year,month,Area_id)
	MLD=np.zeros(n.size)
	ILD=np.zeros(n.size)
	BLT=np.zeros(n.size)
	for i in range(0,n.size):
		MLD[i],ILD[i],BLT[i]=LD.easy_cal(pres[n[i],:],temp[n[i],:],salt[n[i],:],dt)
		MLD[i]=LD.to_5m(MLD[i])
		ILD[i]=LD.to_5m(ILD[i])
	return MLD,ILD,BLT