Re: ِAscii data coordinate adjusting

From: Dave Allured <dave.allured_at_nyahnyahspammersnyahnyah>
Date: Mon Feb 28 2011 - 12:54:36 MST

Ahmed,

The 8 statements from example 2 are added to the program that you
already have. Insert the 8 new statements after "vwnd = new (...)"
in your original program. Then change the name "z" to "vwnd", in 6
lines. For example:

    vwnd!0 = "lat"

Etc.

--Dave

On 2/26/2011 5:06 PM, ahmed lasheen wrote:
> Hello,
> thanks for your reply
> if i use the fspan so i have to get ride of the do loop , but how
> can i fill the
> vwnd array with the data.
>
> On Sun, Feb 27, 2011 at 1:44 AM, Dave Allured <dave.allured@noaa.gov
> <mailto:dave.allured@noaa.gov>> wrote:
>
> You need to generate proper coordinate variables for latitude
> and longitude, and attach them to your data variable "vwnd".
> This method is shown in example 2 on this page:
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/generate_2d_array.shtml
>
> Use *only* the 8 program statements following this comment:
> "Generate 1D lat/lon coord arrays and attach to data." lat1,
> lat2, lon1, lon2 are placed within the two fspan function calls.
>
> Also, it is good practice to set the following variables at the
> top of your program, and never use explicit constants after
> these statements:
>
> nx = 144
> ny = 73
> lat1 = ...
> lat2 = ...
> lon1 = ...
> lon2 = ...
>
> --Dave
>
>
> On 2/26/2011 4:06 PM, ahmed lasheen wrote:
>
> Hello
> I have NC file which contain data about wind , and i wish to
> make
> some calculation on it so i have print
> to ASCII file and then plot it to ensure that i am reading
> it correctly , but when i plot it i found that there is a
> shifting
> in the
> plot(the north is in the south , and the east is in the
> west) .i
> have attached both the original plot from plotting the nc file
> itself and the image of plotting the ASCII file .
> i don not know how to tell ncl that the data is from lat1 to
> lat2
> and lon2 to lon2 specifically.
> Here is the ncl code
> *load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"*
> *load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"*
> *;*************************************************
> *begin*
> *;*************************************************
> *; read in netCDF file*
> *;*************************************************
> * a = asciiread("vdata_ncl",(/10512,1/),"float")*
> *;*************************************************
> *; read in zonal winds*
> *;*************************************************
> *vwnd = new ( (/73,144/), "float") ;; this is fantastic*
> *i=0*
> * do lat=0,72*
> * do lon=0,143*
> * vwnd(lat,lon)=a(i,0)*
> * i=i+1*
> * end do*
> * end do*
> *
> *
> * res = True*
> * res@cnFillOn = True*
> * res@gsnSpreadColors = True*
> * res@tfDoNDCOverlay = True*
> *;*************************************************
> *; create default plot*
> *;*************************************************
> *wks = gsn_open_wks("png","ascii") ; open a ps
> file*
> *gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap*
> *plot2 = gsn_csm_contour_map_ce(wks,vwnd(:,:),res) ; create a
> default plot*
> *end*
>
> I have also attached the ASCII file data.
> thanks in advance
>
> --
> ===============
> Ahmed Lasheen
> Junior researcher at Cairo Numerical Weather Prediction
> Center (CNWPC)
> Egyptian Meteorological Authority(EMA)
> Cairo,Egypt
> ===============
>
>
> nc.000001.png
>
>
>
> ascii.000001.png
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> --
> ===============
> Ahmed Lasheen
> Junior researcher at Cairo Numerical Weather Prediction Center (CNWPC)
> Egyptian Meteorological Authority(EMA)
> Cairo,Egypt
> ===============
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 28 12:54:50 2011

This archive was generated by hypermail 2.1.8 : Wed Mar 02 2011 - 09:18:13 MST