Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Plotting Time Series Diagram using RIP4.7

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

midhunmachary

New member
Hello,
I am currently working on fog. I run WRF model for 24 hours and plotted visibility during fog hours using RIP 4.7. Now I want to extract the predicted visibility over a particular location and to draw a time series for 24 hours.

Please help me to draw a time series diagram of visibility over the location (28.56N, 77.11E) using RIP. What modifications should I make in the script given below.

Thanks in advance,
Midhun M
Department of Atmospheric Science
Central University of Rajasthan, India

&userin
idotitle=1,titlecolor='def.foreground',
ptimes=1,-24,1
ptimeunits='h',tacc=120,timezone=+5.5,iusdaylightrule=1,
iinittime=1,ivalidtime=1,inearesth=0,
flmin=.09, frmax=.92, fbmin=.10, ftmax=.85,
ntextq=0,ntextcd=0,fcoffset=0.0,
idescriptive=1,icgmsplit=0,maxfld=10,itrajcalc=0,imakev5d=0
&end
&trajcalc
rtim=15,ctim=6,dtfile=3600.,dttraj=600.,vctraj='s',
xjtraj=95,90,85,80,75,70,65,80.6,80.6,80.6,80.6,80.6,80.6,
yitraj=50,55,60,65,70,75,80,77,77,77,77,77,77,
zktraj=.9,.9,.9,.9,.9,.9,.9,.99,.9,.8,.7,.6,.5,
ihydrometeor=0
&end
===========================================================================
---------------------- Plot Specification Table ---------------------
===========================================================================
feld=visb; ptyp=hc; cmth=fill; nohl=L; cbeg=0; cend=5.0;>
cosq=0,light.blue,0.2,light.violet,0.5,light.purple,1.0,orange,2.0,light.gray,3.0,>
aqua,4.0,white
#feld=visd; ptyp=hc; cmth=cont; nohl=L; cbeg=0; cend=5.0; linw=2
feld=map; ptyp=hb; outy=CO
feld=map; ptyp=hb; outy=topind.ascii
feld=tic; ptyp=hb
feld=bullet; titl=D; crsa=28.58lat,77.20lon; ptyp=hb
feld=bullet; titl=M; crsa=19.00lat,73.00lon; ptyp=hb
feld=bullet; titl=K; crsa=22.34lat,88.24lon; ptyp=hb
feld=bullet; titl=C; crsa=13.00lat,80.25lon; ptyp=hb
===========================================================================
 
Hi,
I first would like to apologize for the delay in response. We've had members of our group on travel lately and have gotten a bit behind with responses.

Unfortunately RIP doesn't have the capability to do time-series plots.You can use it to output numerical values at specific points, but you would have to do some significant modifying of the code to enable time-series plotting. There is a time-series output capability when running WRF, but you would have to know, in advance, the points you're interested in, and the specific fields you'd like to see. If your runs weren't significantly large and time-consuming, you could look into that. Again, unfortunately, this option doesn't have a built-in visibility feature, so you would have to try to make some modifications to do so. If you are interested, you could read about that option here:
http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.1/users_guide_chap5.html#output_diags

Finally, I was given a fortran program that someone wrote to plot time series directly from WRF output (it doesn't use RIP). I will attach it here, but the person hasn't used it in many years and it's not something we support. You're welcome to take a look and see if you are capable of using it for what you need.
 

Attachments

  • time_series.tar
    4.9 MB · Views: 80
Top