def save_Anomaly_as_npz(fy, ly, var, product_n = 3):
a = subroutine.read_meta_data('var')
formal_name = subroutine.celldata(a, 'var', var, 'formal_name')
dim = subroutine.celldata(a, 'var', var, 'dim')
_,title_name,_=subroutine.product_n_to_name(product_n)
tmpyear = 2000
for month in range(1,13):
_,strmonth=subroutine.strym(2000,month)
Clim=subroutine.load_npz('Ave_of_' + formal_name + '_m'+strmonth+'_' + str(fy) + '-' + str(ly) + '_'+title_name)
print 'month = ', month
start = time.time()
for year in range(fy,ly+1):
stryear,strmonth=subroutine.strym(year,month)
j=year-fy
if dim == '3D':
data=subroutine.get_data(year,month,var,0,title_name)
elif dim == '2D':
data=subroutine.get_data(year,month,var,1,title_name)
else:
raise Exception('your dim is not valid!')
Anomaly=data-Clim
DIR = subroutine.get_DIR(year, month, formal_name, product_n)
subroutine.check_and_make_DIR(DIR)
subroutine.save_npz(Anomaly,DIR+'Anomaly_from_Mean_Annual_Cycle_Year'+str(fy)+'-'+str(ly),data_dir_flg='NO')
print 'elapsed time:', time.time() - start
def save_Vertical_section_of_data_as_npz(self, year, month, data, dname, fname, var, product_n): # 色んなデータの鉛直断面を保存することになるだろう。 # だから、dname,fnameは自分で設定できるようにしておく。 # dnameはD2.D2Data[~].dir_name、 # fnameはD2.D2Data[~].fnameが望ましいかと。 import subroutine Data, _, _ = self.Get_VerticalSection_from_data(data, var, product_n) DIR = subroutine.get_DIR(year, month, dname, product_n) subroutine.check_and_make_DIR(DIR) subroutine.save_npz(Data, DIR + fname + '_VerticalSection_of_' + self.AreaName, data_dir_flg = 'NO')
def load_Anomaly_of_npz(year, month, var, fy = 1990, ly = 2011, product_n = 3):
a = subroutine.read_meta_data('var')
formal_name = subroutine.celldata(a, 'var', var, 'formal_name')
DIR = subroutine.get_DIR(year, month, formal_name, product_n)
data = subroutine.load_npz(DIR + 'Anomaly_from_Mean_Annual_Cycle_Year' + str(fy) + '-' + str(ly), data_dir_flg = 'NO')
return data
def load_data_of_npz(self, year, month, product_n): import subroutine DIR = subroutine.get_DIR(year, month, self.dir_name, product_n) data = subroutine.load_npz(DIR + self.fname, data_dir_flg = 'NO') return data
def save_data_as_npz(self, year, month, data, product_n): import subroutine DIR = subroutine.get_DIR(year, month, self.dir_name, product_n) subroutine.check_and_make_DIR(DIR) subroutine.save_npz(data, DIR + self.fname, data_dir_flg = 'NO')
def load_Area_Trimmed_data_of_npz(self, year, month, dname, fname, product_n): import subroutine DIR = subroutine.get_DIR(year, month, dname, product_n) data = subroutine.load_npz(DIR + fname + '_trimmed_at_' + self.AreaName, data_dir_flg = 'NO') return data
def save_Area_Trimmed_data_as_npz(self, data, year, month, dname, fname, product_n): import subroutine DIR = subroutine.get_DIR(year, month, dname, product_n) subroutine.check_and_make_DIR(DIR) subroutine.save_npz(data, DIR + fname + '_trimmed_at_' + self.AreaName, data_dir_flg = 'NO')
def load_UV_of_VerticalSection_of_npz(self, year, month, product_n): import subroutine DIR = subroutine.get_DIR(year, month, 'Current', product_n) Data = subroutine.load_npz(DIR + 'UV_of_VerticalSection_of_' + self.AreaName, data_dir_flg = 'NO') return Data
def save_UV_of_VerticalSection_as_npz(self, year, month, Data, product_n): import subroutine DIR = subroutine.get_DIR(year, month, 'Current', product_n) subroutine.check_and_make_DIR(DIR) subroutine.save_npz(Data, DIR + 'UV_of_VerticalSection_of_' + self.AreaName, data_dir_flg = 'NO')
def load_Vertical_section_of_data_of_npz(self, year, month, dname, fname, product_n): import subroutine DIR = subroutine.get_DIR(year, month, dname, product_n) data = subroutine.load_npz(DIR + fname + '_VerticalSection_of_' + self.AreaName, data_dir_flg = 'NO') return data
