[ncl-talk] vertically distributing WRF-Chem emissions

[Tabish Ansari]

Hi

I had surface emissions files for WRF-Chem which I'm now trying to
redistribute across 8 vertical levels by creating a new file. I tried
working through this script but got an error:

begin

a = addfile("wrfchemi_00z_d01_pow_Oct.nc","r")
b = addfile("dummy.nc","c")
SO2 =  a->E_SO2(:,0,:,:) ;;;;;EXTRACTING THE SURFACE EMISSIONS
E_SO2 = new((/12,8,120,147/),float)
;;NOW PARTITIONING INTO 8 LEVELS
b->E_SO2(:,0,:,:) = 0.0 * SO2
b->E_SO2(:,1,:,:) = 0.3 * SO2
b->E_SO2(:,2,:,:) = 0.3 * SO2
b->E_SO2(:,3,:,:) = 0.1 * SO2
b->E_SO2(:,4,:,:) = 0.1 * SO2
b->E_SO2(:,5,:,:) = 0.1 * SO2
b->E_SO2(:,6,:,:) = 0.05 * SO2
b->E_SO2(:,7,:,:) = 0.05 * SO2

end


please help.

Thanks

Tabish

Tabish U Ansari
PhD student, Lancaster Environment Center
Lancaster Univeristy
Bailrigg, Lancaster,
LA1 4YW, United Kingdom
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161021/1f09118a/attachment.html