Dear All,
I would like to plot cyclone tracks for the period JJAS data of each
years (for ex: June 2001,july 2001,August 2001,September 2001 like as well
as for the years 2004,2005).Would you please give me any suggestions and
code to do this.
Please find the below summary of time variable.
===========
[nagarjuna@localhost waste]$ ncl test_time_ibtracs.ncl
Copyright (C) 1995-2011 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.0.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
Variable: time
Type: double
Total Size: 7478008 bytes
934751 values
Number of Dimensions: 2
Dimensions and sizes: [storm | 6823] x [time | 137]
Coordinates:
Number Of Attributes: 3
long_name : Modified Julian Day
units : days since 1858-11-17 00:00:00
_FillValue : 9.969209999999999e+36
===========
============================================
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;; Carl Schreck (carl@cicsnc.org) ;; February 2012
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;; Description: Draw a blank map and put genesis locations on it
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
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
;;***************************************************************
;; These are some parameters that could be useful to have up top
;;***************************************************************
plotType = "ps"
plotName = "ibtracs"
plotDpi = 150 ; only applicable if plotType = "png"
;;***************************************************************
;; data path
;;***************************************************************
ibtDir = "./"
ibtPath = "Allstorms.ibtracs_wmo.v03r04.nc
<http://allstorms.ibtracs_wmo.v03r04.nc/>"
;;***************************************************************
;; region
;; -1=All, 0=North Atlantic, 1=South Atlantic, 2=West Pacific,
;; 3=East Pacific, 4=South Pacific, 5=North Indian, 6=South Indian
;;***************************************************************
basin = -1
;; region
minLon = 50 ; 30
maxLon = 120
minLat = 00
maxLat = 30
centerLon = 180 ; graphics
;;***************************************************************
;; time period
;;***************************************************************
year = 2007
timeUnits = "days since 1800-01-01 00:00:00"
minTime = ut_inv_calendar( year, 01, 01, 00, 0, 0, timeUnits, 0 )
maxTime = ut_inv_calendar( year, 12, 31, 18, 0, 0, timeUnits, 0 )
;;***************************************************************
;; Read the data
;;***************************************************************
inFile = addfile( ibtDir+ibtPath, "r" )
stormYear = inFile->season
stormBasin = inFile->genesis_basin
if( basin.ge.0 ) then
stormInd = ind( ( year.eq.stormYear ).and.( basin.eq.stormBasin ) )
else
stormInd = ind( year.eq.stormYear )
end if
stormLat = short2flt( inFile->lat_wmo(stormInd,:) )
stormLon = short2flt( inFile->lon_wmo(stormInd,:) )
stormLon = where( stormLon.lt.0, stormLon+360, stormLon )
printMinMax(stormLat, True)
printMinMax(stormLon, True)
;;***************************************************************
;; Customize base plot
;;***************************************************************
res = True
res@gsnDraw = False
res@gsnFrame = False
res@tiMainString = year
;; ...set the spacing of the lat/lon labels
res@gsnMajorLonSpacing = 20
res@gsnMinorLonSpacing = 10
res@gsnMajorLatSpacing = 20
res@gsnMinorLatSpacing = 10
;; ...map gridlines
res@mpGridAndLimbOn = True
res@mpGridLatSpacingF = 10
res@mpGridLonSpacingF = 10
res@mpGridMaskMode = "MaskLand"
res@mpGridLineColor = "blue"
;; ...set the bounds of a map plot
res@mpMinLatF = minLat
res@mpMaxLatF = maxLat
res@mpMinLonF = minLon
res@mpMaxLonF = maxLon
res@mpCenterLonF = centerLon
;; ...set map resources
res@mpFillOn = False
res@mpGeophysicalLineThicknessF = 2
res@mpGeophysicalLineColor = "darkgreen"
res@mpNationalLineColor = "darkgreen"
res@mpUSStateLineColor = "darkgreen"
res@mpOutlineBoundarySets = "AllBoundaries"
; res@mpDataBaseVersion = "MediumRes"
; res@mpDataSetName = "Earth..1"
;; Set up an x-y line graph
xyRes = True
xyRes@gsnDraw = False
xyRes@gsnFrame = False
xyRes@xyMarkLineMode = "MarkLines"
xyRes@xyMonoDashPattern = True
xyRes@xyLineColor = "red"
xyRes@xyLineThicknessF = 2
xyRes@xyMarker = 1
xyRes@xyMarkerColor = "black"
;; ...allows png or gif to work
if( ( plotType.eq."png" ).or.( plotType.eq."gif" ) ) then
;plotTypeLocal = "eps"
plotTypeLocal = "ps"
else
plotTypeLocal = plotType
end if
;; ...open the workstation
wks = gsn_open_wks( plotTypeLocal, plotName )
gsn_define_colormap( wks, "default" )
map = gsn_csm_map( wks, res )
plot = gsn_csm_xy( wks, stormLon, stormLat, xyRes )
overlay( map, plot )
draw( map )
frame( wks )
;; convert the image if necessary
if( ( plotType.eq."png" ).or.( plotType.eq."gif" ) ) then
system( "convert -trim +repage -density " + plotDpi + " " \\
+ plotName + "."+plotTypeLocal+" " + plotName + "." + plotType )
system( "rm -f " + plotName + "."+plotTypeLocal )
end if
end
============================================================
But above script is used for Annual data(Jan to Dec).I need to
separate a season
wise .Does anybody have any idea how to do this type of work.
Thanks in advance for any help,
Regards,
D.Nagarjuna Rao.
Cc
Bcc
Recipients
Replies will be sent from dnrao.au@gmail.com
Loading...
Attach a file
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 3 05:50:38 2013
This archive was generated by hypermail 2.1.8 : Mon Jun 24 2013 - 11:46:47 MDT