*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:05:04 MDT
This archive was generated by hypermail 2.2.0 : Fri Oct 09 2009 - 08:29:22 MDT