dim_avg

;************************************************
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
   fn = "/storage/sdc1/nlcheng/api/"+"20101230.nc";
   a = addfile(fn,"r")

m=86
n=125
do time = 1,23
 
NO2=a->NO2(time,0,m,n)
NO=a->NO(time,0,m,n)
SO2=a->SO2(time,0,m,n)

ASO4J=a->ASO4J(time,0,m,n)
ASO4I=a->ASO4I(time,0,m,n)
ANH4J=a->ANH4J(time,0,m,n)
ANH4I=a->ANH4I(time,0,m,n)
ANO3J=a->ANO3J(time,0,m,n)
ANO3I=a->ANO3I(time,0,m,n)
AALKJ=a->AALKJ(time,0,m,n)
AXYL1J=a->AXYL1J(time,0,m,n)
AXYL2J=a->AXYL2J(time,0,m,n)
AXYL3J=a->AXYL3J(time,0,m,n)
ATOL1J=a->ATOL1J(time,0,m,n)
ATOL2J=a->ATOL2J(time,0,m,n)
ATOL3J=a->ATOL3J(time,0,m,n)
ABNZ1J=a->ABNZ1J(time,0,m,n)
ABNZ2J=a->ABNZ2J(time,0,m,n)
ABNZ3J=a->ABNZ3J(time,0,m,n)
ATRP1J=a->ATRP1J(time,0,m,n)
ATRP2J=a->ATRP2J(time,0,m,n)
AISO1J=a->AISO1J(time,0,m,n)
AISO2J=a->AISO2J(time,0,m,n)
ASQTJ=a->ASQTJ(time,0,m,n)
AORGCJ=a->AORGCJ(time,0,m,n)
AORGPAJ=a->AORGPAJ(time,0,m,n)
AORGPAI=a->AORGPAI(time,0,m,n)
AECJ=a->AECJ(time,0,m,n)
AECI=a->AECI(time,0,m,n)
A25J=a->A25J(time,0,m,n)
A25I=a->A25I(time,0,m,n)
ACORS=a->ACORS(time,0,m,n)
ASOIL=a->ASOIL(time,0,m,n)
ANAJ=a->ANAJ(time,0,m,n)
ANAI=a->ANAI(time,0,m,n)
ACLJ=a->ACLJ(time,0,m,n)
ACLI=a->ACLI(time,0,m,n)
ANAK=a->ANAK(time,0,m,n)
ACLK=a->ACLK(time,0,m,n)
ASO4K=a->ASO4K(time,0,m,n)
ANH4K=a->ANH4K(time,0,m,n)
ANO3K=a->ANO3K(time,0,m,n)
AISO3J=a->AISO3J(time,0,m,n)
AOLGAJ=a->AOLGAJ(time,0,m,n)
AOLGBJ=a->AOLGBJ(time,0,m,n)

PM10=(ASO4J+ASO4I+ANH4J +ANH4I+ANO3J+ANO3I+AALKJ+AXYL1J +AXYL2J+AXYL3J+ ATOL1J+ATOL2J+ATOL3J+ABNZ1J+ABNZ2J+ABNZ3J+ATRP1J+ATRP2J+AISO1J+AISO2J+ ASQTJ+AORGCJ+AORGPAJ+AORGPAI+AECJ+AECI+A25J+A25I+ ACORS+ASOIL+ ANAJ+ANAI+ACLJ+ACLI+ANAK+ACLK+ASO4K+ANH4K+ANO3K)/1000
 
SO2 = (64/22.4)*SO2
NO2 = (46/22.4)*NO2
NO = (30/22.4)* NO
NOx =NO2+NO
print("SO2: "+SO2+" NOx:"+NOx+" PM10:"+PM10)
end do
print("-------------------------------------------")
 
;SO2 = (64/22.4)*(dim_avg( SO2(TSTEP|1:23, LAY|0 ROW|m, COL|n)))
;NO2 = (46/22.4)*(dim_avg( NO2(TSTEP|1:23, LAY|0 ROW|m, COL|n)))
;NO = (30/22.4)*(dim_avg( NO(TSTEP|1:23, LAY|0 ROW|m, COL|n)))
;NOx =NO2+NO
;PM10 =dim_avg( PM10(TSTEP|1:23, LAY|0 ROW|m, COL|n)))
 
do time =1,23
SO2 = (64/22.4)*(dim_avg( SO2(time|1:23)))
NO2 = (46/22.4)*(dim_avg( NO2(time|1:23)))
NO = (30/22.4)*(dim_avg( NO(time|1:23)))
NOx =NO2+NO
PM10 =dim_avg( PM10(time|1:23))
end do

print("SO2: "+SO2+" NOx:"+NOx+" PM10:"+PM10)

end

I want to calculate the average value SO2,NOx,and PM10 from 1:00 to 23:00 on 30th December. The 2010123.nc file contains 25 hours from 0:00 to 24:00. And i want to use
SO2 = (64/22.4)*(dim_avg( SO2(TSTEP|1:23, LAY|0 ROW|m, COL|n))) and
SO2 = (64/22.4)*(dim_avg( SO2(time|1:23))) but it turns errors. and the line 97 is SO2 = (64/22.4)*(dim_avg( SO2(time|1:23)))
i do not know how to change it .hope for your help! thank you very much

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk