Re: Plotting multiple files...

Mary and Dave,

I appreciate the help on this. It really work the document. I
appreciate the help. Thanks so much...

Cheers,
Erick Rafael Rivera Lemus
Software Technician
Institute of Marine and Coastal Sciences
Coastal Ocean Observation Laboratory

Rutgers, The State University of New Jersey
71 Dudley Road New Brunswick, NJ 08901

Office: 848-932-3261
Fax: 732-932-8578
Cell: 617-447-5106

On Sep 12, 2011, at 3:31 PM, Mary Haley wrote:

> The extra numbers are the old numbering scheme for PNG files that
> was used in NCL versions 5.2.1 and earlier. The first one was xxxx.
> 000001.png, the second one, xxxx.000002.png, etc.
>
> If you upgrade to NCL V6.0.0, then if you just have one frame, NCL
> will call it "xxxx.png".
>
> Dave Allured just posted a link to the FAQ on how to rename these
> internally in your script:
>
> http://www.ncl.ucar.edu/FAQ/#o_formats_005
>
> --Mary
>
> On Sep 12, 2011, at 1:18 PM, Brandon Fisel wrote:
>
>> Erick,
>>
>> You could try saving the .png to another directory called 'mydir',
>> otherwise I'm not sure why the extra numbers are being added on.
>> Try this:
>>
>> pltName = "mydir/" + sfx@fBase
>>
>>
>> -------------------------------------------------------------
>> Brandon Fisel
>>
>>
>> On Mon, Sep 12, 2011 at 13:37, Erick Rivera Lemus <lemus@marine.rutgers.edu
>> > wrote:
>> Hi Brandon,
>>
>> Your suggestion worked like a charm...I appreciate the help on this.
>>
>> The only thing that I notice is whenever I print out the image the
>> name of the image is as follows:
>>
>>
>> OI_MARA_2011_08_23_1000.000001.png
>>
>> Do you have any idea on how to get rid of that 0000001 in between.
>> Not sure why it does that.
>>
>> Thanks for the help
>>
>> Cheers,
>> Erick Rafael Rivera Lemus
>> Software Technician
>> Institute of Marine and Coastal Sciences
>> Coastal Ocean Observation Laboratory
>>
>> Rutgers, The State University of New Jersey
>> 71 Dudley Road New Brunswick, NJ 08901
>>
>> Office: 848-932-3261
>> Fax: 732-932-8578
>> Cell: 617-447-5106
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> On Sep 12, 2011, at 12:08 PM, Brandon Fisel wrote:
>>
>>> Hi Erick,
>>>
>>> gsn_csm_vector_map expects u,v to be 2D but since you are
>>> concatenating u,v they are 3D variables.
>>> Try this:
>>>
>>> do n=0, dimsizes(time) - 1
>>> sfx = get_file_suffix(flista(n),1)
>>> pltName = sfx@fBase
>>> print (pltName)
>>> ;exit
>>> wks = gsn_open_wks(pltType,pltName) ; open a ps file
>>> ; gsn_define_colormap(wks,"BlAqGrYeOrRevi200") ; choose color map
>>> gsn_define_colormap(wks,"temp1") ; choose a color map
>>> i = NhlNewColor(wks,0.7,0.7,0.7) ; add gray to
>>> colormap
>>>
>>> res = True ; plot mods desired
>>> res@gsnSpreadColors = True ; Use full color map
>>> res@gsnSpreadColorStart = 16
>>> res@gsnSpreadColorEnd = -3 ; exlude gray at
>>> color 200
>>> res@mpLandFillColor = "gray" ; set land to be gray
>>> res@mpMinLonF = minLon ; select a
>>> subregion
>>> res@mpMaxLonF = maxLon
>>> res@mpMinLatF = minLat
>>> res@mpMaxLatF = maxLat
>>> res@mpDataBaseVersion = "HighRes"
>>> ; res@mpDataBaseVersion = "MediumRes"
>>> res@pmTickMarkDisplayMode = "Always"
>>> res@pmTickMarkZone = 3
>>>
>>>
>>> res@mpProjection = "mercator" ; projection
>>> res@mpLimitMode = "Corners" ; method to zoom
>>> res@mpLeftCornerLatF = minLat; min(lat)
>>> res@mpLeftCornerLonF = minLon; min(lon)
>>> res@mpRightCornerLatF = maxLat; max(lat)
>>> res@mpRightCornerLonF = maxLon; max(lon)
>>> ; res@mpFillDrawOrder = "PostDraw"
>>>
>>> ;res@tfDoNDCOverlay = True ; do not
>>> transform data
>>>
>>> res@gsnAddCyclic = False
>>> res@gsnMaximize = True ; Maximize plot in frame
>>> res@gsnSpreadColors = True ; span full colormap
>>> res@lbOrientation = "Vertical"
>>> res@lbLabelAutoStride = "True"
>>> res@lbTitleString = "Surface Velocity cm/s"
>>> res@lbTitlePosition = "Right"
>>> res@lbTitleDirection = "Down"
>>> res@lbTitleFontHeightF = "0.015"
>>> res@lbLabelAlignment = "InteriorEdges"
>>> res@vcMonoLineArrowColor = False ; color arrows based on
>>> magnitude
>>> res@vcMonoFillArrowFillColor = False ; Color according to
>>> magnitude
>>> ;res@vcFillArrowEdgeColor = "black"
>>> res@vcLevelSelectionMode = "ManualLevels" ;
>>> res@vcLevelSpacingF = 2.5;
>>> res@vcMinLevelValF = 0;
>>> res@vcMaxLevelValF = 25;
>>>
>>> ; the full domain. If you wish to just consider the domain you are
>>> plotting,
>>> ; you must manually set those levels.
>>>
>>> res@tiMainString = "Standard Radar Network ~C~ Sea Surface
>>> Currents Normal Vectors~C~ " + date(n)
>>>
>>> ;res@gsnCenterString = "Sea Surface Currents"
>>> res@gsnRightString = ""
>>> res@tiXAxisString = " Longitude"
>>> res@tiYAxisString = " Latitude"
>>>
>>> ; res@vcMinFracLengthF = 0.075 ; define
>>> fraction vector
>>> res@vcRefMagnitudeF = 50.0 ; define vector
>>> ref mag
>>> res@vcRefLengthF = 0.08; 0.055 Shorter Vector 0.125
>>> Longer Vector ; define length of vec ref
>>> res@vcMinDistanceF = 0.01 ;0.015 0.0075
>>> 0.000 ; remove out vectors Higher number less are shown
>>> ; 0.015 is the original number. 0.0075 more vectors. 0.000
>>> double the amount of vectors
>>> res@vcRefAnnoOrthogonalPosF = -0.20 ; move ref vector
>>> ;res@vcRefAnnoArrowLineColor = "black" ; change ref
>>> vector color
>>> ;res@vcRefAnnoArrowUseVecColor = True ; don't use vec
>>> color for ref
>>>
>>> res@vcGlyphStyle = "CurlyVector" ; turn on curley
>>> vectors
>>> ;res@vcLineArrowColor = "black" ; change vector
>>> color
>>> res@vcLineArrowThicknessF = 1.0 ; change vector
>>> thickness
>>> ; res@vcVectorDrawOrder = "PostDraw" ; draw vectors
>>> last
>>> res@mpFillDrawOrder = "PostDraw"
>>>
>>> ;res@vcMonoLineArrowColor = True
>>>
>>> plot =gsn_csm_vector_map(wks,u(n,:,:),v(n,:,:),res)
>>>
>>> end do
>>>
>>>
>>> -------------------------------------------------------------
>>> Brandon Fisel
>>>
>>>
>>> On Mon, Sep 12, 2011 at 10:24, Erick Rivera Lemus <lemus@marine.rutgers.edu
>>> > wrote:
>>> Hello,
>>>
>>> I was going thru some of the archiving help list but I couldn't
>>> find something specific to what I want to do.
>>>
>>> I am reading in 24 hrs of hourly data and I will like to plot each
>>> hour as a figure (png works but preferably PS for higher quality
>>> resolution) named with a time stamp. I get this error as well.
>>>
>>> fatal:Number of dimensions in parameter (1) of
>>> (gsn_csm_vector_map) is (3), (2) dimensions were expected
>>> fatal:Execute: Error occurred at or near line 134 in file
>>> mult_plots.ncl
>>>
>>>
>>> I have attach my code.Any input would be greatly appreciate it.
>>>
>>>
>>>
>>> Cheers,
>>> Erick Rafael Rivera Lemus
>>> Software Technician
>>> Institute of Marine and Coastal Sciences
>>> Coastal Ocean Observation Laboratory
>>>
>>> Rutgers, The State University of New Jersey
>>> 71 Dudley Road New Brunswick, NJ 08901
>>>
>>> Office: 848-932-3261
>>> Fax: 732-932-8578
>>> Cell: 617-447-5106
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
Received on Tue Sep 13 19:29:13 2011