Re: NCL Help

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