Re: Plotting multiple files...

From: Brandon Fisel <bjfisel_at_nyahnyahspammersnyahnyah>
Date: Mon Sep 12 2011 - 13:18:49 MDT

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 Mon Sep 12 13:19:17 2011

This archive was generated by hypermail 2.1.8 : Mon Sep 12 2011 - 13:37:01 MDT