pro iug_load_mf_rish_pam_nc, downloadonly=downloadonly, $
trange=trange, $
verbose=verbose
if ~keyword_set(verbose) then verbose=2
if ~size(fns,/type) then begin
file_names = file_dailynames( $
file_format='YYYY/YYYYMMDD',trange=trange,times=times,/unique)+'_pam.nc'
source = file_retrieve(/struct)
source.verbose=verbose
source.local_data_dir = root_data_dir() + 'iugonet/rish/misc/pam/mf/nc/'
source.remote_data_dir = 'http://database.rish.kyoto-u.ac.jp/arch/iugonet/data/mf/pameungpeuk/nc/ver1_0_1/'
local_paths=file_retrieve(file_names,_extra=source)
local_paths_all = ~(~size(local_paths_all,/type)) ? $
[local_paths_all, local_paths] : local_paths
if ~(~size(local_paths_all,/type)) then local_paths=local_paths_all
endif else file_names=fns
if (not keyword_set(downloadonly)) then downloadonly=0
if (downloadonly eq 0) then begin
pam_time=0
zon_wind=0
mer_wind=0
ver_wind=0
height = fltarr(36)
for j=0L,n_elements(local_paths)-1 do begin
file= local_paths[j]
if file_test(/regular,file) then dprint,'Loading pameungpeuk file: ',file $
else begin
dprint,'pameungpeuk file ',file,' not found. Skipping'
continue
endelse
cdfid = ncdf_open(file,/NOWRITE)
glob = ncdf_inquire( cdfid )
print,'Dimensions', glob.ndims
for i=0L,glob.ndims-1 do begin
ncdf_diminq, cdfid, i, name,size
if i EQ glob.recdim then $
print,' ', name, size, '(Unlimited dim)' $
else $
print,' ', name, size
endfor
print
print, 'Variables'
for m=0L,glob.nvars-1 do begin
info = ncdf_varinq(cdfid, m)
FmtStr = '(A," (",A," ) Dimension Ids = [ ", 10(I0," "),$)'
print, FORMAT=FmtStr, info.name,info.datatype, info.dim[*]
print, ']'
for l=0L,info.natts-1 do begin
attname = ncdf_attname(cdfid,m,l)
ncdf_attget,cdfid,m,attname,attvalue
print,' Attribute ', attname, '=', string(attvalue)
endfor
endfor
ncdf_varget, cdfid, 'time', time
ncdf_varget, cdfid, 'range', range
ncdf_varget, cdfid, 'uwind', uwind
ncdf_varget, cdfid, 'vwind', vwind
ncdf_varget, cdfid, 'wwind', wwind
uwind_pam=fltarr(n_elements(time),n_elements(range))
vwind_pam=fltarr(n_elements(time),n_elements(range))
wwind_pam=fltarr(n_elements(time),n_elements(range))
height = range/1000
for i=0L, n_elements(time)-1 do begin
uwind_pam[i,*]=uwind[0,*,i]
vwind_pam[i,*]=vwind[0,*,i]
wwind_pam[i,*]=wwind[0,*,i]
a = uwind_pam[i,*]
wbad = where(a eq -9999,nbad)
if nbad gt 0 then a[wbad] = !values.f_nan
uwind_pam[i,*] =a
b = vwind_pam[i,*]
wbad = where(b eq -9999,nbad)
if nbad gt 0 then b[wbad] = !values.f_nan
vwind_pam[i,*] =b
c = wwind_pam[i,*]
wbad = where(c eq -9999,nbad)
if nbad gt 0 then c[wbad] = !values.f_nan
wwind_pam[i,*] =c
endfor
append_array, pam_time, time
append_array, zon_wind, uwind_pam
append_array, mer_wind, vwind_pam
append_array, ver_wind, wwind_pam
ncdf_close,cdfid
endfor
acknowledgstring = 'Note: If you would like to use following data for scientific purpose, please read and follow the DATA USE POLICY '$
+'(http://database.rish.kyoto-u.ac.jp/arch/iugonet/data_policy/Data_Use_Policy_e.html '$
+'The distribution of MF radar data has been partly supported by the IUGONET (Inter-university Upper '$
+ 'atmosphere Global Observation NETwork) project (http://www.iugonet.org/) funded '$
+ 'by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan.'
if size(zon_wind,/type) eq 4 then begin
dlimit=create_struct('data_att',create_struct('acknowledgment',acknowledgstring,'PI_NAME', 'T. Tsuda'))
store_data,'iug_mf_pam_uwnd',data={x:pam_time, y:zon_wind, v:height},dlimit=dlimit
new_vars=tnames('iug_mf_pam_uwnd')
if new_vars[0] ne '' then begin
options,'iug_mf_pam_uwnd',ytitle='MF-pam!Cheight!C[km]',ztitle='uwnd!C[m/s]'
endif
store_data,'iug_mf_pam_vwnd',data={x:pam_time, y:mer_wind, v:height},dlimit=dlimit
new_vars=tnames('iug_mf_pam_vwnd')
if new_vars[0] ne '' then begin
options,'iug_mf_pam_vwnd',ytitle='MF-pam!Cheight!C[km]',ztitle='vwnd!C[m/s]'
endif
store_data,'iug_mf_pam_wwnd',data={x:pam_time, y:ver_wind, v:height},dlimit=dlimit
new_vars=tnames('iug_mf_pam_wwnd')
if new_vars[0] ne '' then begin
options,'iug_mf_pam_wwnd',ytitle='MF-pam!Cheight!C[km]',ztitle='wwnd!C[m/s]'
endif
new_vars=tnames('iug_mf_pam_*')
if new_vars[0] ne '' then begin
options, ['iug_mf_pam_uwnd','iug_mf_pam_vwnd','iug_mf_pam_wwnd'], 'spec', 1
low_en=-100
high_en=100
low_v=-20
high_v=20
tclip, 'iug_mf_pam_uwnd',low_en,high_en,/overwrite
tclip, 'iug_mf_pam_vwnd',low_en,high_en,/overwrite
tclip, 'iug_mf_pam_wwnd',low_v,high_v,/overwrite
DT=1800
tdegap, 'iug_mf_pam_uwnd',dt=DT,/overwrite
tdegap, 'iug_mf_pam_vwnd',dt=DT,/overwrite
tdegap, 'iug_mf_pam_wwnd',dt=DT,/overwrite
endif
endif
pam_time=0
zon_wind=0
mer_wind=0
ver_wind=0
endif
new_vars=tnames('iug_mf_pam_*')
if new_vars[0] ne '' then begin
print,'******************************
print, 'Data loading is successful!!'
print,'******************************
endif
print, '****************************************************************
print, 'Acknowledgement'
print, '****************************************************************
print, 'Note: If you would like to use following data for scientific purpose,
print, 'please read and follow the DATA USE POLICY'
print, '(http://database.rish.kyoto-u.ac.jp/arch/iugonet/data_policy/Data_Use_Policy_e.html'
print, 'The distribution of MF radar data has been partly supported by the IUGONET'
print, '(Inter-university Upper atmosphere Global Observation NETwork) project'
print, '(http://www.iugonet.org/) funded by the Ministry of Education, Culture, Sports, Science'
print, 'and Technology (MEXT), Japan.'
end