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
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