Re: 2 variables to ASCII file

Hi. Thank you for this tip.
unfortunately I am still getting an error that says NCL does not
recognize "table_attach_row" command.
Attached are both my code and my errors (log file) where the problem
is on line 52. any suggestions

On Nov 29, 2007 8:49 AM, Dennis Shea <shea_at_ucar.edu> wrote:
> Try
>
> https://www.ncl.ucar.edu/Document/Functions/Contributed/table_attach_rows.shtml
>
> uv = table_attach_row (u, v, 0)
> printVarSummary(uv)
>
> write_matrix (uv, fmtx, opt)
> delete(uv)
>
> ---
> Good luck
>
>
> Erik Noble wrote:
> > Dear NCL,
> > My code below currently takes both zonal and meridional wind from a model
> > and intends to write the data out to a file using the matrix ASCII command.
> >
> > I have been sucessful writing one variable out as a matrix ascii file. How
> > can I write both "sets of variables out to one file so that when I use a
> > another program to plot the ascii data, the user "sees" both U and V in the
> > file?
> >
> > Thank you,
> > Erik
> >
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> > ;load "WRFUserARW.ncl"
> >
> > begin
> > ;
> > ; The WRF ARW input file.
> > ; This needs to have a ".nc" appended, so just do it.
> > diri1 ="/Volumes/Data_and_Models/Model-Output/Athena/"
> > ifile1 ="WRF-SOP3_Athena_3_1_2.nc"
> > title = "U_V_700_r2_"
> > a=addfile(diri1+ifile1,"r")
> > ; What times and how many time steps are in the data set?
> > times = wrf_user_list_times(a) ; get times in the file
> > ntimes = dimsizes(times) ; number of times in the file
> >
> > ; The specific pressure levels that we want the data interpolated to.
> > pressure_levels = (/ 850., 700., 500., 300./) ; pressure levels to plot
> > nlevels = dimsizes(pressure_levels) ; number of pressure
> > levels
> >
> > ;do it = 0,ntimes-1 ; TIME LOOP
> > do it =248, ntimes-1, 8 ; TIME LOOP
> >
> > ; First get the variables we will need
> >
> > p = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical
> > coordinate
> > v = wrf_user_getvar(a, "va",it) ; grid point variable
> > u = wrf_user_getvar(a, "ua",it) ; grid point variable
> > printVarSummary(v)
> > do level = 0,nlevels-1 ; LOOP OVER LEVELS
> >
> > pressure = pressure_levels(level)
> > v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
> > u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
> >
> > if ( pressure .eq. 700 ) then
> > N = 109 ; number of rows
> > M = 139 ; number of columns
> > fWidth = 10 ; specify the format width
> > fDec = 3 ; specify the number to the right of decimal point
> > fmtx = M + "f" + fWidth + "." + fDec ; fmtx="50f10.3"
> >
> > opt = True
> > opt_at_fout = title+times(it)
> >
> > write_matrix (v_plane, fmtx, opt)
> > end if
> > end do
> > ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> >
> > end do ; END OF TIME LOOP
> >
> > end
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk_at_ucar.edu
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk