Re: sptial distribution

Hi Adam,

It works. Thanks!

Chao

On Tue, 2006-07-25 at 13:08 -0600, Adam Phillips wrote:
> Hi Chao,
>
> You are setting res1_at_gsnSpreadColors = True , which automatically
> chooses the colors for you. Set that to False and then
> res1_at_cnFillColors = (/ 0,0,5,20,35,50,65,75,85,96/)
> should work...
> Adam
>
> Chao Luo wrote:
> > Hi Adam,
> >
> > I'd like to fill the white color in the min level for each panel, and
> > try something like this:
> >
> > gsn_define_colormap(wks_id,"BlAqGrYeOrRe") ; choose colormap
> > res_at_cnFillColors = (/ 0,0,5,20,35,50,65,75,85,96/)
> > plot(0) = gsn_csm_contour_map_ce(wks_id,so2,res)
> >
> > But it doesn't work. Followings are my script and plot. Would you please
> > take a moment to look at the script and plot?
> >
> > http://dust.ess.uci.edu/chaoluo/temp/tst.ncl
> > http://dust.ess.uci.edu/chaoluo/temp/fig1.ps
> >
> > Thanks,
> > Chao
> >
> >
> > On Tue, 2006-07-25 at 12:17 -0600, Adam Phillips wrote:
> >
> >>Hi Chao,
> >>
> >>When asking for advice it is best if you thoroughly describe what is
> >>going on, what you would like happen, etc. What color is being plotted
> >>as the first two colors in your label bar? An example .gif image of the
> >>resulting plot will also help, along with your updated resource list.
> >>Adam
> >>
> >>Chao Luo wrote:
> >>
> >>>Hi Adam,
> >>>
> >>>Thanks much! It works. One more question is how set white color for the
> >>>min level. I want to fill white color for the minimum level, and set:
> >>>
> >>>gsn_define_colormap(wks_id,"BlAqGrYeOrRe") ; choose colormap
> >>>res_at_cnFillColors = (/ 0,0,5,20,35,50,65,75,85,96/)
> >>>plot(0) = gsn_csm_contour_map_ce(wks_id,so2,res)
> >>>
> >>>But it doesn't work.
> >>>
> >>>Thanks,
> >>>Chao
> >>>
> >>>On Tue, 2006-07-25 at 10:49 -0600, Adam Phillips wrote:
> >>>
> >>>
> >>>>Hi Chao,
> >>>>
> >>>>The fastest varying dimension in fortran is the leftmost dimension of
> >>>>the array. In your case ilon... The fastest varying dimension in ncl is
> >>>>the rightmost dimension... The fastest varying dimension maps to the
> >>>>fastest varying dimension, so upon read-in into NCL, longitude should be
> >>>>the rightmost dimension.
> >>>>
> >>>>data = asciiread("./annl_avg_emi95.dat" \
> >>>> , (/spe_nbr,lat_nbr,lon_nbr/),"float")
> >>>>
> >>>>Adam
> >>>>
> >>>>Chao Luo wrote:
> >>>>
> >>>>
> >>>>>Hi Adam,
> >>>>>
> >>>>>Thanks for your detailed information. I changed the script. The min(so4)
> >>>>>and max(so4) are still not correct. It seems the read data from asc file
> >>>>>is not correct. The array in asc file is (72x46x48) write out by FORTRAN
> >>>>>format: write(99,'(8e12.4)') (((data
> >>>>>(ilon,jlat,kspe),ilon=1,72),jlat=1,46),kspe=1,48)
> >>>>>
> >>>>>The statement of :
> >>>>> data = asciiread("./annl_avg_emi95.dat" \
> >>>>> , (/lon_nbr,lat_nbr,spe_nbr/),"float")
> >>>>>should read right thing?
> >>>>>
> >>>>>Thanks,
> >>>>>
> >>>>>Chao
> >>>>>
> >>>>>
> >>>>>On Mon, 2006-07-24 at 17:33 -0600, Adam Phillips wrote:
> >>>>>
> >>>>>
> >>>>>
> >>>>>>Hi Chao,
> >>>>>>
> >>>>>>Just so you know, you don't need to do 2 do loops here:
> >>>>>>
> >>>>>>do ilon=0,lon_nbr-1
> >>>>>>do ilat=0,lat_nbr-1
> >>>>>>so4(ilon,ilat)=data(ilon,ilat,34-1)*1.e9*0.0409*96. ;;;;ug SO4/m3
> >>>>>>no3(ilon,ilat)=data(ilon,ilat,33-1)*1.e9*0.0409*62. ;;;;ug No3/m3
> >>>>>>nh4(ilon,ilat)=data(ilon,ilat,32-1)*1.e9*0.0409*18. ;;;;ug nh4/m3
> >>>>>>
> >>>>>>so2(ilon,ilat)=data(ilon,ilat,26-1)*1.e9*0.0409*64. ;;;;ug SO2/m3
> >>>>>>hno3(ilon,ilat)=data(ilon,ilat,6-1)*1.e9*0.0409*63. ;;;;ug hno3/m3
> >>>>>>nh3(ilon,ilat)=data(ilon,ilat,30-1)*1.e9*0.0409*17. ;;;;ug nh4/m3
> >>>>>>end do
> >>>>>>end do
> >>>>>>
> >>>>>>should be
> >>>>>>
> >>>>>>so4 = data(:,:,34-1)*1.e9*0.0409*96. ;;;;ug SO4/m3
> >>>>>>
> >>>>>>and so on. You are also creating a new array complete with attributes
> >>>>>>and coordinate variables by doing this:
> >>>>>>so4 =data(:,:,33)
> >>>>>>
> >>>>>>all of data's attributes and coordinate variables are being transferred
> >>>>>>to so4..thus, so4 does not need to be preallocated, so you can remove
> >>>>>>this line:
> >>>>>>so4 =new((/lon_nbr,lat_nbr/),float)
> >>>>>>
> >>>>>>1) if the min/max of so4 is not correct while the min/max of the other
> >>>>>>arrays (no3,so2) are then it is likely that you didn't read the file in
> >>>>>>correctly or that your calculation show in the lines above is not correct.
> >>>>>>
> >>>>>>2) Re: the distribution. NCL expects you to pass arrays dimensioned
> >>>>>>(lat,lon) into the *_map_* plotting functions (ex.
> >>>>>>gsn_csm_contour_map_ce)... You are passing in arrays dimensioned
> >>>>>>(lon,lat), and are likely getting errors related to this. You can fix
> >>>>>>this by dimension reordering:
> >>>>>>plot(0) = gsn_csm_contour_map_ce(wks_id,so2(lat|:,lon:),res1)
> >>>>>>
> >>>>>>3) NCL will always show 1 color > than the max contour level and 1 color
> >>>>>>< the min contour level. In your case set this for plot(0):
> >>>>>> res1_at_cnMinLevelValF = 2.0 ; instead of 0
> >>>>>> res1_at_cnMaxLevelValF = 18.0 ; instead of 20
> >>>>>>
> >>>>>>Finally, latGlobeFo and lonGlobeFo set up your lat and longitude arrays
> >>>>>>for you. You are overwriting them with:
> >>>>>>lon = -177.5 + ispan(0,lon_nbr-1,1) * 5.
> >>>>>>lat(0)=-89.0
> >>>>>>lat(lat_nbr-1)=89.0
> >>>>>>do ilat=1, lat_nbr-2
> >>>>>>lat(ilat)=-86.00 + 4.*(ilat-1)
> >>>>>>end do
> >>>>>>
> >>>>>>Good luck,
> >>>>>>Adam
> >>>>>>
> >>>>>>Chao Luo wrote:
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>>Hi,
> >>>>>>>
> >>>>>>>I am try to make the contour plots by reading asc data of model
> >>>>>>>simulation by previous ncl-talk suggestions. I set 3D variable to couple
> >>>>>>>of 2D variable, and then plot the 2D variable. It seems there are 3
> >>>>>>>types of errors. First, pint out min and max of so4 are not correct.
> >>>>>>>Second, the distribution is not correct, and third one is that there are
> >>>>>>>extra colors less than min and higher than max of levels.
> >>>>>>>I attached script and plot here. Any help and suggestions are very
> >>>>>>>appreciated!
> >>>>>>>
> >>>>>>>Chao
> >>>>>>>http://dust.ess.uci.edu/chaoluo/temp/species_conc_8p_fig1.ncl
> >>>>>>>http://dust.ess.uci.edu/chaoluo/temp/fig1.ps
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>_______________________________________________
> >>>>>>>ncl-talk mailing list
> >>>>>>>ncl-talk_at_ucar.edu
> >>>>>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>>>>>
> >>>>>_______________________________________________
> >>>>>ncl-talk mailing list
> >>>>>ncl-talk_at_ucar.edu
> >>>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>>>
> >>>_______________________________________________
> >>>ncl-talk mailing list
> >>>ncl-talk_at_ucar.edu
> >>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk_at_ucar.edu
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 25 2006 - 15:07:40 MDT