<div dir="ltr">Dear Mary,<div><br></div><div>Thanks so much for your email and useful information.</div><div>Is there any example of extracting flight data using NCL codes? I need to do extract the data using post-processing tools. I need to do interpolation over time and layers(height). </div><div><br></div><div>I appreciate your helps.</div><div><br></div><div>Best Regards,</div><div>Negin</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 17, 2015 at 1:36 PM, Mary Barth <span dir="ltr">&lt;<a href="mailto:barthm@ucar.edu" target="_blank">barthm@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Hi Negin,<br>
    <br>
    For doing the same thing as you, I have found the model level by
    using the WRF height array (&quot;z&quot;) that you can get from wrf_getvar
    routine.<br>
    <br>
    It is also possible to create an output file from WRF-chem of
    variables along a flight track. Here is some more info:<br>
    
    <blockquote type="cite">
      
    </blockquote>
    
    “Tracking”  output:  vertical  profiles  of  prescribed 
    meteorological  and  chemical  <br>
    species  at  a  set  of  prescribed  times  and  horizontal 
    coordinates  are  written  to  a  <br>
    special  output  file,  wrfout_track_d&lt;nn&gt;.    The  namelist 
    variable  track_loc_in  is  a  <br>
    count  of  the  track  locations  and  must  be  set  to  a 
    positive  value  otherwise  the  default  <br>
    setting  of  zero  will  result  in  no  track  output  of  any 
    variables.    Times  and  locations  <br>
    must  be  specified  in  the  file  wrfinput_track.txt(  see 
    example  below).  Meteorological  <br>
    variables  z,  p,  t,  u,  v,  w,  alt,  qcloud,  qrain,  qice, 
    qsnow,  qgraup,  and  qvapor  are  output  <br>
    if  track_loc_in  is  non-­zero.    Chemical  species 
    concentrations  may  also  be  output  if  <br>
    both  namelist  variables  track_chem_num  and  track_chem_name 
    are  set.    The  total  <br>
    number  of  chemical  species  to  output  must  be  &lt;=  100.<br>
    <br>
    As  an  example  the  following  namelist  settings  will  output 
    the  default  meteorological  <br>
    variables  and  co  and  o3  species  concentrations  at  the  two 
    times  and  locations  <br>
    specified  in  the  wrfinput_track.txt  file  :<br>
    <br>
    &amp;domains<br>
    track_loc_in  =  2,<br>
    /<br>
    &amp;chem<br>
    track_chem_num  =  2,<br>
    track_chem_name  =  ‘co’,  ‘o3’,<br>
    /<br>
    <br>
    The  following  two  lines  comprise  the  contents  of  the  ascii 
    input  file  <br>
    wrfinput_track.txt:<br>
    <br>
    2010-­08‐10_00:12:00  41.450  -­87.300<br>
    2010‐08‐10_00:36:00  41.510  -­87.390<br>
    <br>
    This  will  result  in  the  indicated  variables  being  output 
    to  wrfout_track_d01  at  the  <br>
    times  00:12:00  UTC  and  00:36:00  UTC  on  August  10,  2010  at 
    the  grid  points  nearest  <br>
    to  the  points  (41.450, -­87.300)  and  (41.510, -­87.390).   
    Note  that  the  tracking  tool  <br>
    does  not  interpolate  in  time.  The  indicated  output  times 
    need  to  be  given  in  <br>
    multiples  of  model  time  steps  else  no  output  is  produced. 
    <br>
    The  exact  Fortran  format  for  the  lines  in  the  file 
    wrfinput_track.txt  is  <br>
    (A19,1X,F7.3,1X,F7.3)  and  the  horizontal  coordinates  are 
    ordered  for  latitude,  <br>
    longitude  with  standard  WRF  conventions  wherein  south 
    latitudes  and  west  <br>
    longitudes  are  negative.<br>
    <br>
    This information is from
    <a href="http://www.acd.ucar.edu/wrf-chem/MOZCART_UsersGuide.pdf" target="_blank">http://www.acd.ucar.edu/wrf-chem/MOZCART_UsersGuide.pdf</a> <br>
    <br>
    -- Mary Barth<br>
    <br>
    <div>On 2/17/15 12:00 PM,
      <a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a> wrote:<br>
    </div>
    <blockquote type="cite">
      <pre>----------------------------------------------------------------------

Message: 1
Date: Mon, 16 Feb 2015 23:45:34 -0600
From: Negin S <a href="mailto:negin513@gmail.com" target="_blank">&lt;negin513@gmail.com&gt;</a>
Subject: [ncl-talk] flight pathway data
To: Ncl Talk <a href="mailto:ncl-talk@ucar.edu" target="_blank">&lt;ncl-talk@ucar.edu&gt;</a>
Message-ID:
        <a href="mailto:CAJpVDhqVPP-5sXWZ2saDg4vpE4OxcWK8LcN7nwgk_gtRN99X7Q@mail.gmail.com" target="_blank">&lt;CAJpVDhqVPP-5sXWZ2saDg4vpE4OxcWK8LcN7nwgk_gtRN99X7Q@mail.gmail.com&gt;</a>
Content-Type: text/plain; charset=&quot;utf-8&quot;

Hi Everyone!

Recently I am working on extracting data along a flight pathway.

I used the following for figuring out grid index of flight data across its
pathway:

nm = getind_latlon2d (lat2d,lon2d, lat, lon)

do k=0,dimsizes(lat)-1
     n = nm(k,0)
     m = nm(k,1)
end

but I also need some help calculating the layer number based on flight
altitude and time step based on flight time.
I could not find any example for extracting data along flight pathway.

Basically what I wanted to do is to extract data along flight pathway from
model output(i.e. wrf-chem output) and write all the variables from model
in a txt file.

I highly appreciate all your helps.

Best Regards,
Negin
</pre>
    </blockquote>
    <br>
  </div>

<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>