[ncl-talk] OPeNDAP
Ehsan Taghizadeh
ehsantaghizadeh at yahoo.com
Sat Oct 17 10:05:21 MDT 2020
Dear NCL experts,I hope you are doing great.A few months ago I'd asked about reading gfs from OpenDAP via NCL. Actually, I'm still struggling with that and I'll be thankful if I could have your help, please.The problem is that I get zero for the minimum and maximum of some variables, like geopotential height (hgtprs).I've attached a simple script and I'll be thankful if somebody could run it and inform me the result, please. I know one of the experts, Jorge Humberto, could take result via python, and I appreciate him.
Sincerely,Ehsan
On Monday, October 12, 2020, 08:16:05 PM GMT+3:30, Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com> wrote:
Dear Jorge,May I have some lines in python to read GFS from OpenDAP and get the desired values of it and also its values, please?
Sincerely,Ehsan
On Monday, October 12, 2020, 07:34:07 PM GMT+3:30, Jorge Humberto Bravo Méndez <jh.herk at gmail.com> wrote:
I have never used GeoCAT so I can't compare
I like to make my own products, that's why I like python (or other languages I also use a bit of NCL)
I attach an image making a single point graph (located in Mexico) for the variable you want, so I think there is data. In the graph the blue line is GFS, the yellow line is the NAM model
El lun., 12 oct. 2020 a las 10:30, Ehsan Taghizadeh (<ehsantaghizadeh at yahoo.com>) escribió:
Dear Jorge,Thank you for your reply.Actually, that irimo webpage you sent me is from very old. My works doesn't sit in the irimo webpage yet. I've attached one sample of my output.Could you test this page with python, and check output of for example "tmax2m".
https://nomads.ncep.noaa.gov/dods/gfs_0p25_1hr/
Also, do you think using GeoCAT is better or python?
Sincerely,Ehsan
On Monday, October 12, 2020, 6:29:31 PM GMT+3:30, Jorge Humberto Bravo Méndez <jh.herk at gmail.com> wrote:
hello and an apology for the late reply
Thank you for writing to me, due to this pandemic situation the last few months I work at home, and I do not have all the software installed as in my workplace.
I have not installed NCL here, sorry. My main work is in python
I was checking your webpage: https://www.irimo.ir/eng/wd/1070-GFS-Model-Maps.html
Are this your images?
I work in a small (regional) center for weather forecast in Mexico, in one of it's 32 states called Veracruz, I did this website for my wrf local ejecution:
https://cepmpc.github.io/wrfmodel/index_WRF.htmlthe images are created with python
I highly recommend that you use python, I think it is a bit easier to use than NCL, when there is a Python error message showing you where
If there is anything else I can help you with, please write to me
El lun., 5 oct. 2020 a las 2:54, Ehsan Taghizadeh (<ehsantaghizadeh at yahoo.com>) escribió:
Dear Jorge,I hope you are doing great.If you remember, I've asked before about reading GFS from GrADS server (OPeNDAP). You helped me a lot, however, I have a problem with reading variables yet. I'll be thankful if I could have your help again.May I ask you to run the attached script, please? I still get zero for min and max of "hgtprs", "tmax2m", "tmin2m" and "tmp2m".I'll be thankful if I hear from you.
Sincerely,Ehsan TaghizadehForecast Center, IRIMO
On Friday, July 17, 2020, 08:13:49 PM GMT+4:30, Jorge Humberto Bravo Méndez <jh.herk at gmail.com> wrote:
Hello Ehsan!
Oh yes I use python, and ofcourse it can read Opendap servers!
you need install netcdf4 package and use some like
from netcdf4 import Dataset
nc = Dataset('https://server.opendap.com/path/to/you/data')
var = nc['variable_name']
I like python, and I use it a lot. I use software like IDV for a quick view
El vie., 17 jul. 2020 a las 9:52, Ehsan Taghizadeh (<ehsantaghizadeh at yahoo.com>) escribió:
Dear Jorge,Thank you so much for your nice reply. I'll consider UNIDATA-IDV and check its output. However, do you think python could read from opendap servers, too?
Sincerely,Ehsan
On Friday, July 17, 2020, 02:33:42 AM GMT+4:30, Jorge Humberto Bravo Méndez <jh.herk at gmail.com> wrote:
Hello Ehsan!Sorry for the late reply!
I send you a response only to your own email (without ncl-talk) due I not write about ncl
Really I never before try to watch the speed for reading the data for both opendap servers.
Just now I try to see both variables for two servers, using UNIDATA-IDV and in all four variables I do not see zero values, really I have no idea why it is shown in your output
I use the region which is shown in teh images because I'm from Mexico. I recommend you use UNIDATA-IDV for a quick view for variables
I hope this can help you. Regards
El mar., 14 jul. 2020 a las 1:43, Ehsan Taghizadeh via ncl-talk (<ncl-talk at mailman.ucar.edu>) escribió:
Dear JorgeThank you so much for your helpful reply.Actually, I need more help and I'll be thankful if I could have any guidance, again.
It seems the speed of reading from (https://nomads.ncep.noaa.gov:9090/dods/) and (https://thredds.ucar.edu/thredds/idd/forecastModels.html) are different and I got higher speed (but not much) from the later. Is it true?
Also, I tried to read different variables from them. Surprisingly, I took 0 (zero) for min and max for some of the variables, which it couldn't be true.
I've attched the scripts for reading from both. Following are summary of the output for printMinMax.
Reading from nomdas:$ ncl gfsopendapplotncltalk1.nclVariable: prec
(0) ** surface total precipitation [kg/m^2] : min=0 max=641.75
Variable: tempsfc
(0) ** surface temperature [k] : min=200.6 max=356.09
Variable: temp2m
(0) ** 2 m above ground temperature [k] : min=208.9 max=325.9
Variable: tempmax2m
(0) ** 2 m above ground maximum temperature [k] : min=0 max=0
Variable: tempmin2m
(0) ** 2 m above ground minimum temperature [k] : min=208.729 max=323.016
As it could be seen "tempmax2m" has min=0 and max=0!
Reading from thredds:
$ ncl gfsopendapplotncltalk2.ncl
Variable: prec
(0) Total precipitation (Mixed_intervals Accumulation) @ Ground or water surface (kg.m-2) : min=0 max=0
Variable: tempsfc
(0) Temperature @ Ground or water surface (K) : min=200.3 max=356.09
Variable: temphag
(0) Temperature @ Specified height level above ground (K) : min=0 max=0
Variable: tempap
(0) Apparent temperature @ Specified height level above ground (K) : min=181.281 max=328.033
Variable: tempmax
(0) Maximum temperature (Mixed_intervals Maximum) @ Specified height level above ground (K) : min=209.532 max=327
Variable: tempmin
(0) Minimum temperature (Mixed_intervals Minimum) @ Specified height level above ground (K) : min=0 max=0
As it could be seen "temphag" and "tempmin" have min=0 and max=0! Do you have any idea?May I ask whether anyone could plot different variables from these opendap servers?I'll be thankful for any help.
SincerelyEhsan
On Sunday, July 5, 2020, 08:28:05 PM GMT+4:30, Jorge Humberto Bravo Méndez <jh.herk at gmail.com> wrote:
Hello Ehsan,
you can try to plot GSF data using the nomdas opoendap server (https://nomads.ncep.noaa.gov:9090/dods/) or unidata server (https://thredds.ucar.edu/thredds/idd/forecastModels.html)
for download the gfs data you can download only your interest area in this website: (https://nomads.ncep.noaa.gov/) in grib filter, the download is much faster if you not download entire global GFS, and you can create some script for download and with a CRON job you can download periodically
El dom., 5 jul. 2020 a las 5:00, Ehsan Taghizadeh via ncl-talk (<ncl-talk at mailman.ucar.edu>) escribió:
Dear NCL group,I hope you're doing great.
I want to plot GFS data as a daily job. The problem is about downloading GFS data, because of time. So, I've tried to open GFS data from their URL (https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/), however, it seems this URL isn't an OPeNDAP server. Am I right?May I ask if there is any OPeNDAP server to plot GFS data (forecast, 0.5, and 0.25 degrees) with a minimum delay from the mentioned URL?
Thanks in advance for any help.Ehsan_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk
--
Saludos Cordiales M.Sc. Jorge H. Bravo Méndez
Investigador del Clima en eI Centro de Estudios y Pronósticos Meteorológicos
Secretaría de Protección Civil
Gobierno del Estado de Veracruz
Tel 228 1414523
Miembro de la Red de Desastres Asociados a Fenómenos Hidrometeorológicos y Climáticos (REDESClim), red temática CONACYT _______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk
--
Saludos Cordiales M.Sc. Jorge H. Bravo Méndez
Investigador del Clima en eI Centro de Estudios y Pronósticos Meteorológicos
Secretaría de Protección Civil
Gobierno del Estado de Veracruz
Tel 228 1414523
Miembro de la Red de Desastres Asociados a Fenómenos Hidrometeorológicos y Climáticos (REDESClim), red temática CONACYT
--
Saludos Cordiales M.Sc. Jorge H. Bravo Méndez
Investigador del Clima en eI Centro de Estudios y Pronósticos Meteorológicos
Secretaría de Protección Civil
Gobierno del Estado de Veracruz
Tel 228 1414523
Miembro de la Red de Desastres Asociados a Fenómenos Hidrometeorológicos y Climáticos (REDESClim), red temática CONACYT
--
Saludos Cordiales M.Sc. Jorge H. Bravo Méndez
Investigador del Clima en eI Centro de Estudios y Pronósticos Meteorológicos
Secretaría de Protección Civil
Gobierno del Estado de Veracruz
Tel 228 1414523
Miembro de la Red de Desastres Asociados a Fenómenos Hidrometeorológicos y Climáticos (REDESClim), red temática CONACYT
--
Saludos Cordiales M.Sc. Jorge H. Bravo Méndez
Investigador del Clima en eI Centro de Estudios y Pronósticos Meteorológicos
Secretaría de Protección Civil
Gobierno del Estado de Veracruz
Tel 228 1414523
Miembro de la Red de Desastres Asociados a Fenómenos Hidrometeorológicos y Climáticos (REDESClim), red temática CONACYT
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gfsurl.ncl
Type: application/octet-stream
Size: 656 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla de 2020-07-16 16-59-52.png
Type: image/png
Size: 263783 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla de 2020-10-12 10-59-51.png
Type: image/png
Size: 387412 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla de 2020-07-16 16-59-33.png
Type: image/png
Size: 266596 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0007.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla de 2020-07-16 16-59-38.png
Type: image/png
Size: 264076 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0008.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla de 2020-07-16 16-59-50.png
Type: image/png
Size: 263236 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201017/9b8dacf5/attachment-0009.png>
More information about the ncl-talk
mailing list