DATASET FILE TRANSFER

FORM. file_to_transfer.

  PERFORM. make_file.

  PERFORM. dih_send_itts.

ENDFORM.

 

FORM. make_file.

  LOOP AT xtab.
    CLEAR wts_tab.
    wts_tab-line+0(5)    = '1'.
    wts_tab-line+5(10)   = xtab-bukrs.
    wts_tab-line+15(10)  = xtab-byr.
    wts_tab-line+25(20)  = xtab-vbeln_d.
    wts_tab-line+45(20)  = xtab-belnr.
    wts_tab-line+65(10)  = xtab-qty.
    wts_tab-line+75(20)  = xtab-invno.
    wts_tab-line+95(8)   = xtab-inv_dt.
    wts_tab-line+103(20) = xtab-bl_no.
    wts_tab-line+123(20) = xtab-pono.
    wts_tab-line+143(14) = xtab-gi_dt.
    wts_tab-line+157(14) = xtab-cl_dt.
    wts_tab-line+171(8)  = xtab-acc_dt.
    wts_tab-line+189(8)  = xtab-sh_dt.
    wts_tab-line+187(14) = xtab-etd.
    wts_tab-line+201(14) = xtab-eta.
    wts_tab-line+215(3)  = xtab-inco1.
    wts_tab-line+218(5)  = xtab-trtyp.
    wts_tab-line+223(20) = xtab-vessel_name2.
    wts_tab-line+243(3)  = xtab-cnsz.
    wts_tab-line+246(3)  = xtab-unct.
    wts_tab-line+249(10) = xtab-ct_qty.
    wts_tab-line+259(30) = xtab-shln.
    wts_tab-line+289(30) = xtab-vessel_name2.
    wts_tab-line+319(30) = xtab-dtpt.
    wts_tab-line+349(30) = xtab-arpt.
    wts_tab-line+379(30) = xtab-sltp.
    wts_tab-line+409(40) = xtab-slnm.
    wts_tab-line+449(30) = xtab-shtp.
    wts_tab-line+479(40) = xtab-shnm.
    wts_tab-line+519(5)  = xtab-doit+1(5).
    wts_tab-line+524(5)  = xtab-boit+1(5).
    wts_tab-line+529(5)  = xtab-poit+1(5).
    wts_tab-line+534(20) = xtab-btno.
    wts_tab-line+554(10) = xtab-qty.
    wts_tab-line+564(20) = xtab-spart.
    wts_tab-line+584(20) = xtab-bmmd.
    wts_tab-line+604(40) = xtab-mldsc.
    wts_tab-line+644(20) = xtab-bymd.
    wts_tab-line+664(5)  = xtab-prdcd.
    wts_tab-line+669(1)  = xtab-rcdflg.
    wts_tab-line+670(8)  = xtab-tf_dt.
    wts_tab-line+678(1)  = xtab-dropflg.
    wts_tab-line+679(14) = xtab-tsetd.
    wts_tab-line+693(14) = xtab-tseta.
    wts_tab-line+707(30) = xtab-tsvslnam.
    wts_tab-line+737(4)  = xtab-plant.
    wts_tab-line+741(1)  = xtab-bl_flg.
    APPEND wts_tab.
  ENDLOOP.

  IF p_day = 'X'.
    CONCATENATE 'C6K0.Q.ITTS.' sy-datum sy-uzeit 'DSA'
                INTO filename.
  ELSE.
    CONCATENATE 'C6K0.Q.ITTS.' sy-datum sy-uzeit 'HSA'
                INTO filename.
  ENDIF.

  CONCATENATE dir_n filename INTO datafile.

  OPEN DATASET datafile FOR OUTPUT IN TEXT MODE MESSAGE msg_text.
  IF sy-subrc <> 0.
    WRITE : 'File cannot be opened. Reason:', msg_text.
    EXIT.
  ENDIF.

  LOOP AT wts_tab.
    TRANSFER wts_tab TO datafile.
  ENDLOOP.
  CLOSE  DATASET datafile.

ENDFORM.                    " MAKE_FILE

来自 “ ITPUB博客 ” ,链接:http://blog.itpub.net/80479/viewspace-265137/,如需转载,请注明出处,否则将追究法律责任。

转载于:http://blog.itpub.net/80479/viewspace-265137/

import netCDF4 as nc import numpy as np from netCDF4 import Dataset import matplotlib.pyplot as plt from matplotlib.cm import get_cmap from matplotlib.colors import from_levels_and_colors import cartopy.crs as crs import cartopy.feature as cfeature from cartopy.feature import NaturalEarthFeature from wrf import to_np, getvar, interplevel, smooth2d, get_cartopy, cartopy_xlim, cartopy_ylim, latlon_coords, vertcross, smooth2d, CoordPair, GeoBounds,interpline import warnings warnings.filterwarnings('ignore') file = 'D:/transfer/wrfout_d01_2016-03-01_00_00_00' dataset = nc.Dataset(file) latitude = dataset.variables['XLAT'][0][:] longitude = dataset.variables['XLONG'][0][:] tp1 = dataset.variables['RAINC'][1][:][:] co = dataset.variables['co'][1][1][:][:] time = dataset.variables['Times'][:] co2 = dataset.variables['co2'][:] #var = ds.variables['co2'] #print(co2[:]) plt.imshow(co2[ :, :, 98, 78], cmap='hot_r', vmax=400, vmin=350, alpha=0.5) plt.colorbar() #plt.scatter(latitude,longitude, c=co, s=3, cmap='Reds', vmax=1, vmin=0) proj = crs.PlateCarree(central_longitude=180) proj_data = crs.PlateCarree()#LambertCylindrical() #plt.contourf(co[:, :, 98, 78], cmap='hot') fig , ax = plt.subplots(1,1,figsize=(8,8),subplot_kw={'projection':proj}) #plt.imshow(longitude, latitude, co) ax.set_title('CO2 concentration') #ax.set_xlabel('Longitude') #ax.set_ylabel('Latitude') ax.add_feature(cfeature.COASTLINE.with_scale('50m'),lw=0.5) ax.add_feature(cfeature.BORDERS) leftlon, rightlon, lowerlat, upperlat = (90, 110, 4, 31) ######## 调节绘图经纬度范围 Region = [leftlon, rightlon, lowerlat, upperlat] ax.set_extent(Region, crs=proj_data) #经纬度范围,坐标参考系转换 plt.show()
05-28
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值