请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 235|回复: 1

[其他] 求助:将1维txt整合为2维txt的问题

[复制链接]

新浪微博达人勋

发表于 2019-7-8 16:49:41 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
先生成了一个月一个txt文件,总共3年,每年10个月的,把它们整合到一个2维txt里面。
树立了很多遍,总是提示循环有错误,不知道到底哪里错了,求教大神。
错误提示如图:
fatal:syntax error:line 78 in file 8.ncl before or near do end do
fatal:error in statement
fatal:syntax error:line 79 in file 8.ncl before or near do end do
fatal:error in statement

fatal:syntax error:line 132 in file 8.ncl before or near \n end
fatal:error in statement

语句如下:

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
lonL = 80
lonR = 120
latN = 50
latS = 10
nn=lonR-lonL+1
mm=latN-latS+1
len_grid=mm*nn
;print(mm)
;print(nn)
;print(len_grid)
;len_t=1
yearnum=3
monthnum=8
sta_year=2016
sta_month=3
anualDaynum=245
Daynum=anualDaynum*yearnum
rawdata=new((/len_grid,Daynum/),"float")
outdir="/mnt/d/fx/"
indir="/mnt/d/fx/"
;datapath="/mnt/d/fx/"
startj=0
do iyear=0,yearnum-1
do imonth=0,monthnum-1
  yeartmp=sta_year+iyear
  monthtmp=sta_month+imonth
  years=sprinti("%0.2i",yeartmp)
  months=sprinti("%0.2i",monthtmp)
  yearstr=years+months;actually it is yearmonthstr
  print(yearstr)
;end do
;end do
;yearstr="201603"
datapath=indir+"hgt500mb"+yearstr+"mn1D.txt"
print(datapath)
if (monthtmp.eq.3) then
endj=startj+30
len_t=31
end if
if (monthtmp.eq.4) then
endj=startj+29
len_t=30
end if
if (monthtmp.eq.5) then
endj=startj+30
len_t=31
end if
if (monthtmp.eq.6) then
endj=startj+29
len_t=30
end if
if (monthtmp.eq.7) then
endj=startj+30
len_t=31
end if
if (monthtmp.eq.8) then
endj=startj+30
len_t=31
end if
if (monthtmp.eq.9) then
endj=startj+29
len_t=30
end if
if (monthtmp.eq.10) then
endj=startj+30
len_t=31
  
rawdata(:,startj:endj)=asciiread(datapath,(/len_grid,len_t/),"float")
startj=endj+1
delete(datapath)
delete(len_t)
end do
end do
;printVarSummary(rawdata)
;printMinMax(rawdata,True)  
rawdata!0="grid"
rawdata!1="day";time in columns
printVarSummary(rawdata)
outdata=rawdata(day|:,grid|:);time in lows
printVarSummary(outdata)
s=dimsizes(outdata)
print(s(0))
print(s(1))
;print(outdata(:,1))
ncol=s(1)
  ;fmtf="501f9.4"
  fmtf   = ncol + "f12.4";s(1) number of columns  
  opt        =True
  opt@fout=outdir+"hgt500mbnm2D.txt"
  write_matrix(outdata,fmtf,opt)
  
  ;ncol=s(0)
  ;fmtf="501f9.4"
  ;fmtf   = ncol + "f12.4";s(1) number of columns  
  ;opt        =True
  ;opt@fout=outdir+"hgt500mb_201404Dailyavg_nm2D.txt"
  ;write_matrix(rawdata,fmtf,opt)
;**************************************************************************;
if (1.eq.0) then
wks=gsn_open_wks("eps",picdir+"hgt925")
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
gsn_define_colormap(wks,"GMT_polar")               
res                 = True                    ; plot mods desired
;res@gsnDraw         = False                   ; don't draw yet
;res@gsnFrame        = False                   ; don't advance frame yet
res@gsnAddCyclic    =False
;res@cnLevelSelectionMode = "ExplicitLevels"   ; set explicit contour levels
;res@cnLevels        = (/232.7,234.2,238,240,244,248.4,252,258.3,276,286.5,292.1,300,306/)
res@tiMainOn        = False
res@cnFillOn        = True                    ; turn on color
; res@cnFillPalette   = "gui_default"           ; set color map
res@lbLabelBarOn    = False                   ; will draw a panel label bar instead
res@cnLinesOn       = False                   ; turn off contour lines
res@cnLineLabelsOn  = False                   ; tuen off line labels
res@lbLabelBarOn=True
res@mpMaxLatF       =latN
res@mpMinLatF       =latS
res@mpMaxLonF       =lonR
res@mpMinLonF       =lonL
;plot= gsn_csm_contour_map_ce(wks,hgt0({lv_ISBL0|ilev},{lat_0|-90:90},{lon_0|0:359}), res)
;plot= gsn_csm_contour_map_ce(wks,hgt, res)
end if
end



error.png
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2019-7-8 17:30:13 | 显示全部楼层
不好意思了大家,问题已解决,少了个end if
密码修改失败请联系qq:937062711
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表