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