Re: NCL script for meteograms

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri, 14 Mar 2008 07:23:48 -0600

Merhaba Abdullah

I speculate there are no WRF meteogram scripts. The NCL example
demonstrates how to create a simple meteogram.
               http://www.ncl.ucar.edu/Applications/meteo.shtml

I am not sure what you mean by an "advanced" meteogram.

---
By definition, a meteogram consists of a time series of observations
for different variables at a specific location and, generally, at
a single level.   Most commonly the surface level is visualized.
There are no "meteogram standards", hence, meteograms
take on a variety of forms.
For a WRF type model, a simple approach might be as follows.
[a]
Specify the subscripts corresponding to a grid point closest
to some station (eg: Istanbul, Turkey).
[b]
Read the various variables for that grid point.
[c]
The WRF "Times" variable is awkward to work with for plotting.
NCL offers several functions that will transform this character
variable to numeric forms.
            http://www.ncl.ucar.edu/Document/Functions/wrf.shtml
Specifically:
   
http://www.ncl.ucar.edu/Document/Functions/WRF_contributed/wrf_times_c.shtml
So a thumb nail approach might be as follows.
procedure metogram_WRF (time[*]:numeric,t[*],q[*],u[*],v[*],t[*],s[*],...)
begin
    ; create various plots as per simple example.
end
begin ; MAIN program
    ii =          ; index of desired longitude grid point
    jj =         ;                           latitude
    kk=        ;                           level
    diri = "./"
    fili  = "wrf....nc"
    f = addfile (diri+fili, "r")
    Times = f->Times
    time    =  wrf_times_c( Times, 0 )
     t = f->T(:,kk,jj,ii)   read all the times
    qflx = f->QFX(:,kk,jj,ii)
        :
     meteogram_WRF(time, t,qflx,.....)
    
end
If meteogram_WRF was created as a function
function meteogram_WRF (...)
begin
   ; create various plots as per simple example.
   return(meteogram)
end
 
later
    plot =    meteogram_WRF (...)
    You could overlay more information after returning
     draw(plot)
      frame(wks)
Taking on a task like this can be tedious BUT it is often
the best way to learn a computer language. My 1st
NCL project was to develeop the skewT program.
   http://www.ncl.ucar.edu/Applications/skewt.shtml
Good Luck
Abdullah Kahraman wrote:
> Hello all,
> I am a beginner in NCL. Is there anyone who uses NCL for creating 
> meteograms? There is an example at NCL's website, but actually I am 
> looking for an advanced one to be used with WRF-ARW output, to process 
> my forecasts.
> I will be pleased if one sends such a script.
> Regards,
> Abdullah Kahraman
> Istanbul.
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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
Received on Fri Mar 14 2008 - 07:23:48 MDT

This archive was generated by hypermail 2.2.0 : Mon Mar 17 2008 - 13:59:23 MDT