Re: sptial distribution

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
Received on Tue Jul 25 2006 - 12:42:22 MDT