vgalligani / SimRad.AR

Herramienta capaz de simular las variables observadas por los radar polarimétrico del Sistema Nacional de Radares Meteorológicos (SINERAME) de Argentina, a partir de las salidas de modelos numéricos en escala de nubes, particularmente del Weather Research and Forecasting model (WRF) el cual es utilizado operativamente en el Servicio Meteorológico Nacional. El simulador está escrito en Python y está disponible para su uso libremente.
0 stars 0 forks source link

error when use "mp == 10" #1

Open zzwei1 opened 3 years ago

zzwei1 commented 3 years ago

Hi, nice work! I'm doing some research about simulating polarimetric weather radar variables from WRF. I'm very interesting in your repo SimRad.AR, a python version radar simulator. I also found CRSIM, however, I'm not familar with fortran, so I wonna use a python version. However, after generating lookup tables for a double moment microphysics scheme, I wonna run wrf2radar.py. When I set "mp == 10", I found there's no ineterpolation functions ( line 258 in wrf2radar.py, elif (mp == 10): display('get from wrf2radar_old')display('get from wrf2radar_old')).

So where could I find wrf2radar_old ? Is it another file ? Is there more information and reference about this repo? thanks in advance!

vgalligani commented 3 years ago

CONGREMET_GALLIGANi.pdf https://drive.google.com/file/d/1EX96aNHffNeNu0U8fQl5MU8AvZZdG5PS/view?usp=drive_web Hi ! How are you? Thanks for contacting me. What kind of work are you planning on doing? The current distribution does not include the Morrison scheme, but I have previously worked with it in the past so it shouldn't pose too much of a problem. The only problem being that I am very busy at this very moment, but I would definitely put on my list of things to do in the near future. Would that be ok? I am attaching a poster I presented once (sorry it's in spanish) but understanding the figures should be straightforward. In the poster I compare my radar fw model to CRSIM. One important thing is that SimRad.AR also needs some information regarding the radar you are trying to simulate. Did you see that? Let me know, Victoria

On Tue, Apr 13, 2021 at 7:16 AM zzwei1 @.***> wrote:

Hi, nice work! I'm doing some research about simulating polarimetric weather radar variables from WRF. I'm very interesting in your repo SimRad.AR, a python version radar simulator. I also found CRSIM, however, I'm not familar with fortran, so I wonna use a python version. However, after generating lookup tables for a double moment microphysics scheme, I wonna run wrf2radar.py. When I set "mp == 10", I found there's no ineterpolation functions ( line 258 in wrf2radar.py, elif (mp == 10): display('get from wrf2radar_old')display('get from wrf2radar_old')).

So where could I find wrf2radar_old ? Is it another file ? Is there more information and reference about this repo? thanks in advance!

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/vgalligani/SimRad.AR/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/AEOEKN3MO66AQTOIM37SJVDTIQKZBANCNFSM423BLIHA .

zzwei1 commented 3 years ago

CONGREMET_GALLIGANi.pdf https://drive.google.com/file/d/1EX96aNHffNeNu0U8fQl5MU8AvZZdG5PS/view?usp=drive_web Hi ! How are you? Thanks for contacting me. What kind of work are you planning on doing? The current distribution does not include the Morrison scheme, but I have previously worked with it in the past so it shouldn't pose too much of a problem. The only problem being that I am very busy at this very moment, but I would definitely put on my list of things to do in the near future. Would that be ok? I am attaching a poster I presented once (sorry it's in spanish) but understanding the figures should be straightforward. In the poster I compare my radar fw model to CRSIM. One important thing is that SimRad.AR also needs some information regarding the radar you are trying to simulate. Did you see that? Let me know, Victoria On Tue, Apr 13, 2021 at 7:16 AM zzwei1 @.***> wrote: Hi, nice work! I'm doing some research about simulating polarimetric weather radar variables from WRF. I'm very interesting in your repo SimRad.AR, a python version radar simulator. I also found CRSIM, however, I'm not familar with fortran, so I wonna use a python version. However, after generating lookup tables for a double moment microphysics scheme, I wonna run wrf2radar.py. When I set "mp == 10", I found there's no ineterpolation functions ( line 258 in wrf2radar.py, elif (mp == 10): display('get from wrf2radar_old')display('get from wrf2radar_old')). So where could I find wrf2radar_old ? Is it another file ? Is there more information and reference about this repo? thanks in advance! — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#1>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AEOEKN3MO66AQTOIM37SJVDTIQKZBANCNFSM423BLIHA .

CONGREMET_GALLIGANi.pdf https://drive.google.com/file/d/1EX96aNHffNeNu0U8fQl5MU8AvZZdG5PS/view?usp=drive_web Hi ! How are you? Thanks for contacting me. What kind of work are you planning on doing? The current distribution does not include the Morrison scheme, but I have previously worked with it in the past so it shouldn't pose too much of a problem. The only problem being that I am very busy at this very moment, but I would definitely put on my list of things to do in the near future. Would that be ok? I am attaching a poster I presented once (sorry it's in spanish) but understanding the figures should be straightforward. In the poster I compare my radar fw model to CRSIM. One important thing is that SimRad.AR also needs some information regarding the radar you are trying to simulate. Did you see that? Let me know, Victoria On Tue, Apr 13, 2021 at 7:16 AM zzwei1 @.***> wrote: Hi, nice work! I'm doing some research about simulating polarimetric weather radar variables from WRF. I'm very interesting in your repo SimRad.AR, a python version radar simulator. I also found CRSIM, however, I'm not familar with fortran, so I wonna use a python version. However, after generating lookup tables for a double moment microphysics scheme, I wonna run wrf2radar.py. When I set "mp == 10", I found there's no ineterpolation functions ( line 258 in wrf2radar.py, elif (mp == 10): display('get from wrf2radar_old')display('get from wrf2radar_old')). So where could I find wrf2radar_old ? Is it another file ? Is there more information and reference about this repo? thanks in advance! — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#1>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AEOEKN3MO66AQTOIM37SJVDTIQKZBANCNFSM423BLIHA .

Thanks for your reply!

I'm doing some work about simulating airship-borne polarimetric weather radar variables from WRF outpus. It would be great if you can help me in your free time !

I notice you use a lookuptable method for calculating polarimetric weather radar variables. The current method I can think of resolving the error I mesioned above is that, I directly use WRF outputs to calculate weather radar variables instead of using a lookuptable method. Of course, it takes more computaion and time than a lookuptable method.

Yes, I have noticed that SimRad.AR needs some information regarding the radar that I'm trying to simulate. And I have some confusions about it. Maybe you can help me in your free time.

There are differences between the geometric relationship between the ground_based and the airborne weather radar, beacuse airborne weather radar always have a “look-down” geometry. I think I can modify it based on your method. This involves recalculating the range and elevation angle using different formulas.

However, the first confusion is, how to interpolate WRF variables to a radar coordinate. My WRF variables have a resolution about 2km in horizontal resolution and a vertical resoltion about 500m. However, the radar system I woona simulate has a radial resolution of 375m. The spatial scale are not match. What's you WRF resolution and radar radial resolution, maybe I can take as an reference.

The second is, what does the "stacked_column_mass" (in wrf2radar.py, line 354) mean, when I set "apply_mask == 1" ? By the way, using mask in your code is such an efficient method of reducing comptation! I was inspired by it!

The third is, what does the "mode" (in wrf2radar.py, line 495) mean ? Actually, I have not understood the code about attenea pattern weighting, maybe I can ask you more about it after I understanding it.

Thanks in advance!

Hope to keep communication!