strange values ​​extracted with grads

All other topics about postprocessing model data (GrADS and other software), about other numerical weather prediction software (including WRF-NMM and WRF-ARW discussion unrelated to UEMS/WRF EMS), and general meteorology talk go in this forum.
meteoadriatic
Posts: 1584
Joined: Wed Aug 19, 2009 10:05 am

Re: strange values ​​extracted with grads

Post by meteoadriatic » Thu Jan 02, 2014 11:01 am

It has nothing to do with wrf version but with your scripts.

Antonix
Posts: 260
Joined: Fri Oct 16, 2009 8:53 am

Re: strange values ​​extracted with grads

Post by Antonix » Thu Jan 02, 2014 11:21 am

Also I think it's a problem in your script.
them have modified from the working version???

just one question:
when you change version of wrf
you create a new domain or copy it from the old version???

create always a new domain
It does not happen "always" but often the copy of the domain due to some problem in different version.

dominic
Posts: 336
Joined: Thu Jun 14, 2012 7:19 am

Re: strange values ​​extracted with grads

Post by dominic » Thu Jan 02, 2014 11:29 am

Guys I'm doing a run on Italy with both core nmm that arw at 9km, then I will also plot with grads and I pull the data.

I post everything here around 14.00 so you see what can be the problem

dominic
Posts: 336
Joined: Thu Jun 14, 2012 7:19 am

Re: strange values ​​extracted with grads

Post by dominic » Thu Jan 02, 2014 11:30 am

Antonix wrote:Also I think it's a problem in your script.
them have modified from the working version???

just one question:
when you change version of wrf
you create a new domain or copy it from the old version???

create always a new domain
It does not happen "always" but often the copy of the domain due to some problem in different version.

hello Antonix

I've used both scripts and both give the same data extracted. with regard to the domains I have tried to recreate them, using and not using Modis, but nothing.

I will post something clearer to study with you a possible solution

Antonix
Posts: 260
Joined: Fri Oct 16, 2009 8:53 am

Re: strange values ​​extracted with grads

Post by Antonix » Thu Jan 02, 2014 11:46 am

if you can plot the data as a map, and you can not extract them with txt, then the problem is in your script.
I had a similar problem and it was the fault of the localization domain.
Make a simple domain (no nesting, no modis). not change physical settings.
test to see if it works.
(ems_domain --newdom aaa --core arw
go in to domain folder and type
ems_domain --localize).
if you have other problems is the fault of your script.
ps. what is the size of your grib file?
use the right precautions for large files?

dominic
Posts: 336
Joined: Thu Jun 14, 2012 7:19 am

Re: strange values ​​extracted with grads

Post by dominic » Thu Jan 02, 2014 12:00 pm

Antonix, Ivan this is my Output on Run 06z core NMM for city of Napoli

Code: Select all

02JAN201406|napoli|2|11.9464|10.8526|11.7457|8.00147|79.4635|0.00803626|0|0|0|4.78907|7.58418|263.524|5.2297|1019.76|82.229|7229.21|1916.31|0.00531171|
02JAN201406|napoli|3|12.4012|11.3872|11.9235|8.16202|78.0645|1.91109|0|0|0|5.08395|7.68722|228.623|2.02987|1019.96|82.229|7229.21|1990.51|0.635216|
02JAN201406|napoli|4|12.9651|11.8204|12.4758|7.80751|73.4175|1.91109|0|0|0|8.66672|12.8768|216.615|3.83249|1019.63|82.229|7497.19|2022.24|0.793386|
02JAN201406|napoli|5|13.3847|12.419|12.805|7.81817|71.6375|3.57712|0|0|0|13.3169|15.7322|214.66|9.24088|1019.48|82.229|7585.71|1965.01|2.06884|
02JAN201406|napoli|6|13.6249|12.8151|13.0591|8.25545|72.6055|3.99371|0|0|0|15.4418|17.3192|211.037|21.6854|1019.21|82.229|7473.83|1905.01|1.92928|
02JAN201406|napoli|7|13.7398|13.0654|13.1765|8.84796|75.0327|6.17222|0|0|0|17.1885|19.2635|205.581|31.4639|1018.99|82.229|7185.71|1910.86|2.33227|
02JAN201406|napoli|8|13.6813|13.0973|13.1996|9.58174|79.1444|27.3614|0|0|0|17.3167|19.9686|198.598|46.7376|1018.92|82.229|6874.24|1942.64|4.0593|
02JAN201406|napoli|9|13.6075|12.9981|13.1977|10.3909|83.8509|99.9759|0|0|0|16.8519|20.4618|189.077|68.0098|1018.7|82.229|6557.54|2005.77|3.07461|
02JAN201406|napoli|10|13.5081|12.8959|13.1321|10.924|87.3219|100|0.0140669|0|0|16.7367|20.3969|179.53|88.2777|1018.62|82.229|6282.9|2049.58|5.31472|
02JAN201406|napoli|11|13.3329|12.7037|13.1112|11.4164|91.0366|100|0.0665943|0|0|16.6286|21.3694|171.419|108.662|1018.66|82.229|6087.55|2075.86|8.19654|
02JAN201406|napoli|12|13.4297|12.8153|13.0137|11.7453|92.181|100|0.209077|0|0|16.9928|23.6507|170.584|123.825|1018.58|82.229|5976.07|2139.11|11.6249|
02JAN201406|napoli|13|13.5683|12.9738|13.1515|11.7904|91.449|100|0.889452|0|0|17.2883|24.1582|180.491|124.929|1018.58|82.229|6003.84|2197.81|2.15812|
02JAN201406|napoli|14|13.6877|13.1579|13.2972|11.8696|91.0871|100|1.07614|0|0|16.8526|23.1828|192.87|136.924|1018.64|82.229|5063.57|2212.29|0.142955|
02JAN201406|napoli|15|13.753|13.2996|13.3725|11.9049|90.8649|67.3963|0.782905|0|0|16.6148|22.5228|202.98|136.3|1018.77|82.229|5079.29|2210.81|4.5465|
02JAN201406|napoli|16|13.7832|13.3674|13.4183|11.9198|90.732|70.3944|0.670379|0|0|16.7358|22.6482|208.943|139.141|1018.95|82.229|5088.93|2204.29|4.97492|
02JAN201406|napoli|17|13.8705|13.41|13.5192|11.9224|90.1099|100|0.864568|0|0|17.3197|24.2975|214.144|155.209|1019.27|82.229|6206.25|2189.39|3.91536|
02JAN201406|napoli|18|13.9894|13.5245|13.6687|11.8357|88.7133|100|1.15443|0.0273707|0|17.5608|25.2809|216.068|163.431|1019.38|82.229|6317.73|2176.49|14.9672|
02JAN201406|napoli|19|14.0415|13.6719|13.7365|11.7287|87.7089|100|1.01876|0.0216544|0|17.1869|23.1668|212.938|154.173|1019.38|82.229|6406.25|2179.41|12.8979|
02JAN201406|napoli|20|14.0444|13.7297|13.7514|11.7693|87.8934|100|0.642746|0.000511485|0|16.3955|22.4283|208.839|161.557|1019.42|82.229|6406.25|2179.7|12.4189|
02JAN201406|napoli|21|14.043|13.7413|13.7668|11.739|87.6916|100|0.525636|0.000476412|0|16.102|22.2375|202.706|157.41|1019.4|82.229|6406.25|2180.93|11.8336|
02JAN201406|napoli|22|14.0376|13.7307|13.7693|11.8504|88.3494|100|0.436236|0|0|15.9281|22.0384|193.581|159.621|1019.34|82.229|6322.55|2182.88|11.5152|
02JAN201406|napoli|23|14.0145|13.7064|13.7567|11.8753|88.6937|100|0.36939|9.48417e-05|0|16.376|22.1665|183.882|139.214|1019.32|82.229|6310.5|2196.04|12.5448|
02JAN201406|napoli|24|14.0174|13.6928|13.7485|12.0499|89.7298|96.4968|0.347983|0.000117028|0|16.9007|23.2412|175.862|157.35|1019.47|82.229|6199.02|2215.99|13.5117|
02JAN201406|napoli|25|14.0292|13.6907|13.7407|12.2021|90.6046|100|0.474118|0.0002|0|17.836|24.5988|171.543|181.724|1019.64|82.229|6170.44|2250.43|10.799|

the column of precipitation is the tenth column, the eleventh column is that of the snow

'set lat 40.86'
'set lon 14.28'

Antonix
Posts: 260
Joined: Fri Oct 16, 2009 8:53 am

Re: strange values ​​extracted with grads

Post by Antonix » Thu Jan 02, 2014 12:08 pm

ok
the error is in your script.
I repeat: you can plot a map of precipitation??

dominic
Posts: 336
Joined: Thu Jun 14, 2012 7:19 am

Re: strange values ​​extracted with grads

Post by dominic » Thu Jan 02, 2014 12:14 pm

Antonix wrote:ok
the error is in your script.
I repeat: you can plot a map of precipitation??

this is the plot

http://postimg.org/image/bkqidws4j/

Antonix
Posts: 260
Joined: Fri Oct 16, 2009 8:53 am

Re: strange values ​​extracted with grads

Post by Antonix » Thu Jan 02, 2014 12:26 pm

okay...summarizing
1) your model works
2) if you are using screen viewing grads work well.
3) the map makes
4) in the txt extracts all the values ​​except rain.
I think the problem is in your script.
or in some mystery...for example one mystery is:

why in your script uses

'd apcpsfc'
regel = sublin(result,9)
res_apcp=subwrd(regel,4)

and not

'd ACPCPaccsfc'
regel = sublin(result,9)
res_apcp=subwrd(regel,4)

???

dominic
Posts: 336
Joined: Thu Jun 14, 2012 7:19 am

Re: strange values ​​extracted with grads

Post by dominic » Thu Jan 02, 2014 12:47 pm

Antonix wrote:okay...summarizing
1) your model works
2) if you are using screen viewing grads work well.
3) the map makes
4) in the txt extracts all the values ​​except rain.
I think the problem is in your script.
or in some mystery...for example one mystery is:

why in your script uses

'd apcpsfc'
regel = sublin(result,9)
res_apcp=subwrd(regel,4)

and not

'd ACPCPaccsfc'
regel = sublin(result,9)
res_apcp=subwrd(regel,4)

???

Antonix In fact, it is quite a mystery :/

Yes I put ACPCPaccsfc as a variable but nothing!

Now I try to use the script by Cory and I let you know the result

Post Reply