2 variables to ASCII file

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
Received on Wed Nov 28 2007 - 12:38:31 MST