Re: dim_avg

   m=86
   n=125
   tStrt = 1
   tLast = 23
          ; USE ARRAY SYNTAX
   NO2 = a->NO2(tStrt:tLast,0,m,n)
       :
   AOLGBJ = a->AOLGBJ(tStrt:tLast,0,m,n)

   printVarSummary(NO2) ; NO2(*)
   printVarSummary(AOLGBJ) ; AOLGBJ(*)
          :
   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
   printVarSummary(PM10)

The above are one dimensional. You can use 'avg' or 'dim_avg'
PLEASE read the documentation.

   NO2_avg = avg(NO2)
           :
   AOLGBJ_avg = avg(AOLGBJ)

On 11/1/11 9:27 AM, cnl88 wrote:
> ;************************************************
> 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
>
>
> ------------------------------------------------------------------------
> 网易公开课推出哈佛最受欢迎《幸福》课程，帮您修满幸福学分!
> <http://v.163.com/special/positivepsychology/?from=mailfooter>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 2 08:20:36 2011