Re: mass change when regridding rectilinear grid to a curvilinear grid

Dear Bin Zhang,

I had to reject your post because it had a 25 megabyte attachment. This is too large for the ncl-talk mail list.

I have your files, and will take a look soon.

-Mary

> .
>
> From: "zhangbinrc" <zhangbinrc@163.com>
> Subject: mass change when regridding rectilinear grid to a curvilinear grid
> Date: October 22, 2013 1:53:29 AM MDT
> To: <ncl-talk@ucar.edu>
>
>
> Hi all,
>
> I want to regrid global ammonia emissions (0.5 x 0.667 degree) to a lambert confromal domain (27km x 27km) covering China.
>
> The total ammonia mass is about ten times of that before regridding. Does anyone knows why ?
>
> Here are the script:
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> begin
>
> fil = addfile("./geo_em.d01.nc","r")
> lat2d = fil->XLAT_M(0,:,:)
> lon2d = fil->XLONG_M(0,:,:)
> nlat2d = dimsizes(lat2d(:,0))
> nlon2d = dimsizes(lon2d(0,:))
>
>
> nh3 = new((/nlat2d,nlon2d/),"float")
> fil_nh3 = addfile("./Streets_NH3_FF_2000.geos5.05x0666.nc","r")
> s_nh3 = fil_nh3->NH3
>
> printVarSummary(s_nh3)
> nh3_d = rgrid2rcm(s_nh3&lat,s_nh3&lon,s_nh3,lat2d,lon2d,0)
> nh3 = doubletofloat(nh3_d)
>
> ; printVarSummary(nh3)
> print(dim_sum(dim_sum(nh3))/10^6)
> print(dim_sum(dim_sum(s_nh3))/10^6)
>
> end
>
>
>
>
> 张斌 博士生/Bin Zhang Ph.D. Candidate
> 清华大学地球系统科学研究中心/ Center for Earth System Science, Tsinghua
> Email:zhangbin06@mails.tsinghua.edu.cn
> <geo_em.d01.nc><Streets_NH3_FF_2000.geos5.05x0666.nc><regrid_nh3.ncl>
>
>
>
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 22 08:19:10 2013