the line 68 is end do
the script is below
I am not sure whether the red part of the script is true! hope for your help
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/ncl/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
load "/usr/local/ncl/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
z= asciiread("./site-i-j.txt",(/6336,3/),"float")
I = z(:,0) ;
i = z(:,1);
j = z(:,2);
do I =1,6336
i(I)=z(I,1);
j(I)=z(I,1);
m=j(I)
n=i(I)
do ntime = 1,29
ifile = ntime
mtime =ntime+10
fn = "/storage/sdc1/nlcheng/cmaq/site/METCRO2D/04/"+ifile+".nc";
a = addfile(fn,"r")
fon = "201104-"+ ntime+"-"+I+"-"+n+"-"+m
PBL = a->PBL(:,0,m,n)
TEMP2= a->TEMP2(:,0,m,n)
WSPD10 = a->WSPD10(:,0,m,n)
WDIR10=a->WDIR10(:,0,m,n)
print(PBL)
I = 25 ; number of elements per variable
K = 4 ; number of variables (1D arrays)
p = new ( (/I,K/), typeof(PBL) )
p(:,0) = (/ PBL /)
p(:,1) = (/ TEMP2 /)
p(:,2) = (/ WSPD10 /)
p(:,3) = (/ WDIR10 /)
fmtx = "4f11.5"
opt = True
opt@row = True
opt@title = "T(ÊÀ½çʱ), PBL(m) ,TEMP2(K), WSPD10(m/s) ,WDIR10(DEGREES)"
opt@fout = fon
write_matrix (p,fmtx, opt)
end do
end do
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Tue Oct 18 2011 - 09:42:15 MDT