search

gustfrontar / LETKF_WRF

Letkf-wrf source code and scripts
GNU General Public License v3.0
3 stars 2 forks source link

How to use the scripts to convert WRF outputs to radar coordinates? #7

Open zzwei1 opened 4 years ago

zzwei1 commented 4 years ago

Dear gustfrontar, I'm a freshman in radar meteorology and Linux, and now I want to convert the WRF outputs(some .nc files) to radar coordinates, it's something like interpolating the WRF outputs to radar beam center along the beam.I found the scripts in 'LETKF_WRF/wrf/obs/radar/wrf_to_radar/src' may meet my requirements.But I don't know how to use your scripts(like xxx.sh), I tried to study to use the gfortran to compile the scripts but get some errors. Could you tell me how to use your scripts to achieve my goal? Thank you in advance! Best regards, zzwei

gustfrontar commented 4 years ago

Dear Zzwei, Sorry for my late response. Currently, the code allows you to do that (for reflectivity and Doppler velocity), but you need at least one volume of radar data in cfradial format (a netcdf file which some metadata of the geometry of the data). If you have something like this, then the code can be used to interpolate WRF outputs to this radar geometry. You can convert radar data in many formats to the cfradial format using the python module pyart. Let me know if you have data in cfradial or if your data can be converted in cfradial and I can explain to you how to use this code.

All the best Juan

Juan Ruiz Centro de Investigaciones del Mar y la Atmosfera (CIMA/CONICET-UBA) Dpto. de Ciencias de la Atmosfera y los Oceanos, FCEN, Universidad de Buenos Aires Instituto Franco-Argentino para el Estudio del Clima y sus Impactos (UMI IFAECI/CNRS-CONICET-UBA)

Intendente Guiraldes 2160 - Ciudad Universitaria Pabellon II - 2do. Piso (C1428EGA) Ciudad Autonoma de Buenos Aires - Argentina Phone: (54)(11) 5285-8453 Email: jruiz@cima.fcen.uba.ar

On 26.02.2020 05:09, zzwei1 wrote:

Dear gustfrontar, I'm a freshman in radar meteorology and Linux, and now I want to convert the WRF outputs(some .nc files) to radar coordinates, it's something like interpolating the WRF outputs to radar beam center along the beam.I found the scripts in 'LETKF_WRF/wrf/obs/radar/wrf_to_radar/src' may meet my requirements.But I don't know how to use your scripts(like xxx.sh), I tried to study to use the gfortran to compile the scripts but get some errors. Could you tell me how to use your scripts to achieve my goal? Thank you in advance! Best regards, zzwei

-- You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub [1], or unsubscribe [2]. [ { "@context": "http://schema.org", "@type": "EmailMessage", "potentialAction": { "@type": "ViewAction", "target": "https://github.com/gustfrontar/LETKF_WRF/issues/7?email_source=notifications\u0026email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA", "url": "https://github.com/gustfrontar/LETKF_WRF/issues/7?email_source=notifications\u0026email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA", "name": "View Issue" }, "description": "View this Issue on GitHub", "publisher": { "@type": "Organization", "name": "GitHub", "url": "https://github.com" } } ]

Links:

[1] https://github.com/gustfrontar/LETKF_WRF/issues/7?email_source=notifications&email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA [2] https://github.com/notifications/unsubscribe-auth/AEBADYX63SQAWGMRGKL45ULREYPSTANCNFSM4K37IRVA

zzwei1 commented 4 years ago

Dear Zzwei, Sorry for my late response. Currently, the code allows you to do that (for reflectivity and Doppler velocity), but you need at least one volume of radar data in cfradial format (a netcdf file which some metadata of the geometry of the data). If you have something like this, then the code can be used to interpolate WRF outputs to this radar geometry. You can convert radar data in many formats to the cfradial format using the python module pyart. Let me know if you have data in cfradial or if your data can be converted in cfradial and I can explain to you how to use this code. All the best Juan … --- Juan Ruiz Centro de Investigaciones del Mar y la Atmosfera (CIMA/CONICET-UBA) Dpto. de Ciencias de la Atmosfera y los Oceanos, FCEN, Universidad de Buenos Aires Instituto Franco-Argentino para el Estudio del Clima y sus Impactos (UMI IFAECI/CNRS-CONICET-UBA) Intendente Guiraldes 2160 - Ciudad Universitaria Pabellon II - 2do. Piso (C1428EGA) Ciudad Autonoma de Buenos Aires - Argentina Phone: (54)(11) 5285-8453 Email: jruiz@cima.fcen.uba.ar On 26.02.2020 05:09, zzwei1 wrote: Dear gustfrontar, I'm a freshman in radar meteorology and Linux, and now I want to convert the WRF outputs(some .nc files) to radar coordinates, it's something like interpolating the WRF outputs to radar beam center along the beam.I found the scripts in 'LETKF_WRF/wrf/obs/radar/wrf_to_radar/src' may meet my requirements.But I don't know how to use your scripts(like xxx.sh), I tried to study to use the gfortran to compile the scripts but get some errors. Could you tell me how to use your scripts to achieve my goal? Thank you in advance! Best regards, zzwei -- You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub [1], or unsubscribe [2]. [ { @.": "http://schema.org", @.": "EmailMessage", "potentialAction": { @.": "ViewAction", "target": "#7?email_source=notifications\u0026email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA", "url": "#7?email_source=notifications\u0026email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA", "name": "View Issue" }, "description": "View this Issue on GitHub", "publisher": { @.": "Organization", "name": "GitHub", "url": "https://github.com" } } ] Links: ------ [1] #7?email_source=notifications&email_token=AEBADYWSOWW6ARSKLYT4ZQLREYPSTA5CNFSM4K37IRVKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4IQK6BRA [2] https://github.com/notifications/unsubscribe-auth/AEBADYX """ Input radar filename and use appropriate PyART functions to read file and return radar object63SQAWGMRGKL45ULREYPSTANCNFSM4K37IRVA

Thanks for your reply! I'm sorry to say that I don't have cfradial data for now. I wonna simulate the observations of a 'fake' radar. I have read the scripts in /LETKF_WRF-master/wrf/obs/radar/scripts, and I found that there's some defination of the 'fake radar' in H_run_model_to_radar_20140122_control.sh. Could I use the 'fake radar defination(e.g., I define a radar at a specified latitude, longitude and height)' instead of the cfradial data of a real radar? If I have to use cfradial data ,I think may be I could use open source radar data (e.g. wsr88d), to get the cfradial data.

Thanks in advance!

Best regards, Zzwei