*Dear Ncl users,**
*
*The error i found is not coming now, by some changes i made. Mainly i used
latGlobeF& lonGlobeF function only.*
*
*
*The new script is like this:*
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
;
; Read a netCDF file using addfile.
;
a = addfile("TROUT.01-02-00.nc","r")
wks=gsn_open_wks("ps","BCconcentration")
;
;Add Coordinate variables
;
nlat = 221
lat = latGlobeF(nlat, "lat", "latitude", "degrees_north")
mlon = 331
lon = lonGlobeF(mlon, "lon", "longitude", "degrees_east")
bcc1 = a->BC1(0,0,:,:)
bcc2 = a->BC2(0,0,:,:)
TotBC = bcc1+bcc2
TotBC!0 = "lat"
TotBC!1 = "lon"
TotBC&lon = lon
TotBC&lat = lat
printVarSummary(TotBC)
print(TotBC&lat)
print(TotBC&lon)
gsn_define_colormap(wks,"gui_default") ; select color map
res = True ; plot mods desired
res_at_gsnSpreadColors = True ; make cn span entire color
map
res_at_mpFillOn = False ; turn off gray continents
res_at_cnFillOn = True ; turn on color fill
res_at_mpMinLonF = 48 ; set min lon
res_at_mpMaxLonF = 100 ; set max lon
res_at_mpMinLatF = 0 ; set min lat
res_at_mpMaxLatF = 30 ; set max lat
plot = gsn_csm_contour_map_ce(wks,TotBC,res) ; create plot
end
*But the problem is still there with fspan.....!!!!!*
*
*
*Cant we use fspan for setting lat and lon?*
*
*
*
*
*The regional setting made was:(which still shows error)*
;nlat = 221
;mlon = 331
;lat = fspan(0,30,nlat)
;lon = fspan(48,100,mlon)
;lat!0 = "lat"
;lat&lat = lat
;lat_at_units = "degress_north"
;lat_at_long_name = "latitude"
;lon!0 = "lon"
;lon&lon = lon
;lon_at_units = "degress_east"
;lon_at_long_name = "longitude"
*Thanks.....!!!!*
-- A.R.Ragi or simply A.R MTech Atmospheric Science CUSAT/CDAC +919960462099 Sent from Maharashtra, India On Thu, Oct 8, 2009 at 10:35, A.R Ragi <ar.ragi_at_gmail.com> wrote: > *Dear Saji , David and all,** > * > *I printed TotBC and TotBC&lat and TotBC&lon.* > *The TotBC output is like:* > > Variable: TotBC > Type: float > Total Size: 292604 bytes > 73151 values > Number of Dimensions: 2 > Dimensions and sizes: [lat | 221] x [lon | 331] > Coordinates: > lat: [-59.875..20] > lon: [60..120] > (0) check_for_y_lat_coord: Warning: Data either does not contain a > valid latitude coordinate array or doesn't contain one at all. > (0) A valid latitude coordinate array should have a 'units' attribute > equal to one of the following values: > (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' > 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north' > (0) check_for_lon_coord: Warning: Data either does not contain a valid > longitude coordinate array or doesn't contain one at all. > (0) A valid longitude coordinate array should have a 'units' attribute > equal to one of the following values: > (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' > 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east' > > > > The TotBC&lat and TotBC&lon is also printing. > * > * > *But the error is still remaining.* > * > * > *Thanks for help.* > -- > A.R.Ragi or simply A.R > MTech Atmospheric Science > CUSAT/CDAC > +919960462099 > > > Sent from Maharashtra, India > > ---------- Forwarded message ---------- > From: A.R Ragi <ar.ragi_at_gmail.com> > Date: Thu, Oct 8, 2009 at 10:22 > Subject: NCL Help > To: dbrown_at_ucar.edu > > > > -- > A.R.Ragi or simply A.R > MTech Atmospheric Science > CUSAT/CDAC > +919960462099 > > > Sent from Maharashtra, India > > ---------- Forwarded message ---------- > From: A.R Ragi <ar.ragi_at_gmail.com> > Date: Thu, Oct 8, 2009 at 10:19 > Subject: NCL Help > To: ncl-talk_at_ucar.edu > > > > Sent from Maharashtra, India > > ---------- Forwarded message ---------- > From: A.R Ragi <ar.ragi_at_gmail.com> > Date: Thu, Oct 8, 2009 at 10:17 > Subject: Re: NCL Help > To: Saji <saji_at_apcc21.net> > > > *Dear Ncl Users,* > * > * > *I tried by adding the steps you recommended . The errors showing now are: > * > > > > Variable: bcc1 > Type: float > Total Size: 292604 bytes > 73151 values > Number of Dimensions: 2 > Dimensions and sizes: [ROW | 221] x [COL | 331] > Coordinates: > Number Of Attributes: 3 > long_name : BC1 > units : microgram/m3 > var_desc : Species Concentration > fatal:No coordinate variable exists for dimension (ROW) in variable (bcc1) > fatal:Execute: Error occurred at or near line 53 in file bcoc.ncl > > > > > *Again I changed the following steps and tried:* > > ;print(bcc1&ROW) > ;print(bcc1&COL) > > > *So i got the same errors that lat and lon is not there.* > > > Variable: bcc1 > Type: float > Total Size: 292604 bytes > 73151 values > Number of Dimensions: 2 > Dimensions and sizes: [ROW | 221] x [COL | 331] > Coordinates: > Number Of Attributes: 3 > long_name : BC1 > units : microgram/m3 > var_desc : Species Concentration > (0) check_for_y_lat_coord: Warning: Data either does not contain a > valid latitude coordinate array or doesn't contain one at all. > (0) A valid latitude coordinate array should have a 'units' attribute > equal to one of the following values: > (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' > 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north' > (0) check_for_lon_coord: Warning: Data either does not contain a valid > longitude coordinate array or doesn't contain one at all. > (0) A valid longitude coordinate array should have a 'units' attribute > equal to one of the following values: > (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' > 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east' > > > > *Is it needed to add lat and lon values in both bcc1 and bcc2 then can > find TotBC?* > * > * > *Or anything else?* > * > * > * > * > *Thanks in advance....!!!!!* > > -- > A.R.Ragi or simply A.R > MTech Atmospheric Science > CUSAT/CDAC > +919960462099 > > > Sent from Maharashtra, India > > On Wed, Oct 7, 2009 at 17:57, Saji <saji_at_apcc21.net> wrote: > >> Ragi, >> >> > mlon = 150 >> > lon = lonGlobeF(mlon, "lon", "longitude", "degrees_east") >> > nlat = 100 >> > lat = latGlobeF(nlat, "lat", "latitude", "degrees_north") >> > TotBC!0 = "lat" >> > TotBC!1 = "lon" >> >> First of all, you need to also do the following >> >> TotBC&lon = lon >> TotBC&lat = lat >> >> Secondly, I see a mismatch between your description of the size of >> lat and lon co-ordinates from the output from ncl_filedump. According to >> ncl_filedump: >> > dimensions: >> > TSTEP = 4 // unlimited >> > DATE-TIME = 2 >> > LAY = 41 >> > VAR = 31 >> > ROW = 221 >> > COL = 331 >> >> lat should have a size of 221 and lon of 331. In contrast you use >> 100 and 150 resp. Maybe you should print out the following and check >> if it looks like a latitude/longitude description. >> >> printVarSummary(bcc1) >> print(bcc1&ROW) >> print(bcc1&COL) >> >> You may get a clue on how to define your latitude/longitude co-ordinate >> values based on this information... >> >> >> saji >> >> >> >> >> * A.R Ragi <ar.ragi_at_gmail.com> [2009-10-07 17:00:20 +0530]: >> >> > Dear NCL users, >> > >> > >> > I was trying to write the black carbon and organic carbon concentration >> > script from the STEM model output. >> > >> > The output contains BC1, BC2,OC1, OC2. >> > >> > But it doesn't contains the Latitude and longitude. Last time asked >> about >> > the doubt in ncl talk and so from the reply I tried to add lat and lon >> > using both fspan and latGlobeF and lonGlobeF. >> > >> > But still saying error remains. >> > >> > I'm sure that the error remains with my code itself. >> > >> > Can you please check the code and help me. >> > >> > The file description and my code is following. >> > >> > >> > >> > The description of the file is: >> > >> > >> > Variable: f (file variable) >> > >> > filename: TROUT.01-02-00 >> > path: TROUT.01-02-00.nc >> > file global attributes: >> > IOAPI_VERSION : $Id: @(#) ioapi library version 3.0 $ >> > EXEC_ID : 3.0 >> > FTYPE : 1 >> > CDATE : 2009252 >> > CTIME : 61455 >> > WDATE : 2009252 >> > WTIME : 61455 >> > SDATE : 2006002 >> > STIME : 0 >> > TSTEP : 60000 >> > NTHIK : 1 >> > NCOLS : 331 >> > NROWS : 221 >> > NLAYS : 41 >> > NVARS : 31 >> > GDTYP : 3 >> > P_ALP : 17.5 >> > P_BET : 17.5 >> > P_GAM : 0 >> > XCENT : 75 >> > YCENT : 17.49999618530273 >> > XORIG : 28.14615058898926 >> > YORIG : -20.71478080749512 >> > XCELL : 27000 >> > YCELL : 27000 >> > VGTYP : 2 >> > VGTOP : 1 >> > VGLVLS : <ARRAY of 42 elements> >> > GDNAM : WRF >> > UPNAM : >> > VAR-LIST : SO2 VSO2 SO2LPS NOX >> > NOY HC1 HC2 HC3 >> BIOMCO >> > Hg CO CO_TRACER1 CO_TRACER2 >> > CO_TRACER3 CO_TRACER4 CO_TRACER5 CO_TRACER6 >> CO_TRACER7 >> > CO_TRACER8 CO_TRACER9 CO_TRACER10 CO_TRACER11 >> SULF >> > BC1 BC2 OC1 OC2 >> SSF >> > SSC DUST1 DUST2 >> > FILEDESC : Generated by program AQMS >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > >> > HISTORY : >> > dimensions: >> > TSTEP = 4 // unlimited >> > DATE-TIME = 2 >> > LAY = 41 >> > VAR = 31 >> > ROW = 221 >> > COL = 331 >> > variables: >> > integer TFLAG ( TSTEP, VAR, DATE-TIME ) >> > units : <YYYYDDD,HHMMSS> >> > long_name : TFLAG >> > var_desc : Timestep-valid flags: (1) YYYYDDD or (2) HHMMSS >> > >> > float SO2 ( TSTEP, LAY, ROW, COL ) >> > long_name : SO2 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float VSO2 ( TSTEP, LAY, ROW, COL ) >> > long_name : VSO2 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float SO2LPS ( TSTEP, LAY, ROW, COL ) >> > long_name : SO2LPS >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float NOX ( TSTEP, LAY, ROW, COL ) >> > long_name : NOX >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float NOY ( TSTEP, LAY, ROW, COL ) >> > long_name : NOY >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float HC1 ( TSTEP, LAY, ROW, COL ) >> > long_name : HC1 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float HC2 ( TSTEP, LAY, ROW, COL ) >> > long_name : HC2 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float HC3 ( TSTEP, LAY, ROW, COL ) >> > long_name : HC3 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float BIOMCO ( TSTEP, LAY, ROW, COL ) >> > long_name : BIOMCO >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float Hg ( TSTEP, LAY, ROW, COL ) >> > long_name : Hg >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO ( TSTEP, LAY, ROW, COL ) >> > long_name : CO >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER1 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER1 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER2 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER2 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER3 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER3 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER4 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER4 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER5 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER5 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER6 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER6 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER7 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER7 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER8 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER8 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER9 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER9 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER10 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER10 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float CO_TRACER11 ( TSTEP, LAY, ROW, COL ) >> > long_name : CO_TRACER11 >> > units : ppbv >> > var_desc : Species Concentration >> > >> > float SULF ( TSTEP, LAY, ROW, COL ) >> > long_name : SULF >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float BC1 ( TSTEP, LAY, ROW, COL ) >> > long_name : BC1 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float BC2 ( TSTEP, LAY, ROW, COL ) >> > long_name : BC2 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float OC1 ( TSTEP, LAY, ROW, COL ) >> > long_name : OC1 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float OC2 ( TSTEP, LAY, ROW, COL ) >> > long_name : OC2 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float SSF ( TSTEP, LAY, ROW, COL ) >> > long_name : SSF >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float SSC ( TSTEP, LAY, ROW, COL ) >> > long_name : SSC >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float DUST1 ( TSTEP, LAY, ROW, COL ) >> > long_name : DUST1 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > float DUST2 ( TSTEP, LAY, ROW, COL ) >> > long_name : DUST2 >> > units : microgram/m3 >> > var_desc : Species Concentration >> > >> > >> > >> > > >> > > The code i attempted is: >> > > >> > > >> > > >> > > >> > > ;Add Coordinate variables >> > > >> > > >> > > >> > > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" >> > > >> > > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" >> > > >> > > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" >> > > >> > > >> > > >> > > begin >> > > >> > > ; >> > > >> > > ; Read a netCDF file using addfile. >> > > >> > > ; >> > > >> > > >> > > >> > > a = addfile("TROUT.01-02-00.nc","r") >> > > >> > > >> > > >> > > wks=gsn_open_wks("ps","BCconcentration") >> > > >> > > >> > > >> > > ; >> > > >> > > ;Add Coordinate variables >> > > >> > > ; >> > > >> > > >> > > >> > > ;nlat = 100 >> > > >> > > ;mlon = 150 >> > > >> > > >> > > >> > > ;lat = fspan(-59.875,20,nlat) >> > > >> > > ;lon = fspan(60,120,mlon) >> > > >> > > >> > > >> > > ;lat!0 = "lat" >> > > >> > > ;lat&lat = lat >> > > >> > > ;lat_at_units = "degress_north" >> > > >> > > ;lat_at_long_name = "latitude" >> > > >> > > >> > > >> > > ;lon!0 = "lon" >> > > >> > > ;lon&lon = lon >> > > >> > > ;lon_at_units = "degress_east" >> > > >> > > ;lon_at_long_name = "longitude" >> > > >> > > >> > > >> > > mlon = 150 >> > > >> > > lon = lonGlobeF(mlon, "lon", "longitude", "degrees_east") >> > > >> > > >> > > >> > > nlat = 100 >> > > >> > > lat = latGlobeF(nlat, "lat", "latitude", "degrees_north") >> > > >> > > >> > > >> > > >> > > >> > > bcc1 = a->BC1(0,0,:,:) >> > > >> > > bcc2 = a->BC2(0,0,:,:) >> > > >> > > >> > > >> > > TotBC = bcc1+bcc2 >> > > >> > > >> > > >> > > TotBC!0 = "lat" >> > > >> > > TotBC!1 = "lon" >> > > >> > > >> > > >> > > gsn_define_colormap(wks,"gui_default") >> > > >> > > >> > > >> > > res = False >> > > >> > > res_at_gsnSpreadColors = True >> > > >> > > >> > > >> > > res_at_mpFillOn = False >> > > >> > > res_at_cnFillOn = True >> > > >> > > res_at_mpMinLonF = -140 >> > > >> > > res_at_mpMaxLonF = -65 >> > > >> > > res_at_mpMinLatF = 10. >> > > >> > > >> > > >> > > plot = gsn_csm_contour_map_ce(wks,TotBC,res) ; create plot >> > > >> > > >> > > >> > > >> > > >> > > end >> > > >> > > --- >> > > I'm waiting for your valuable reply. Please reply when you get time. >> > > Thanks in advance. >> > > >> > > >> > > >> > > -- >> > > A.R.Ragi or simply A.R >> > > MTech Atmospheric Science >> > > CUSAT/CDAC >> > > +919960462099 >> > > >> > > >> > > Sent from Maharashtra, India >> > > >> >> > _______________________________________________ >> > ncl-talk mailing list >> > List instructions, subscriber options, unsubscribe: >> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk >> >> >> -- >> Saji N. Hameed >> >> APEC Climate Center >> 1463 U-dong, Haeundae-gu, +82 51 745 3951 >> BUSAN 612-020, KOREA saji_at_apcc21.net >> Fax: +82-51-745-3999 >> >> >> >> > > > >
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 07 2009 - 23:51:07 MDT
This archive was generated by hypermail 2.2.0 : Fri Oct 09 2009 - 08:29:22 MDT