Simrad ar

README.md

@@ -7,5 +7,3 @@ For an overview of the WRF modeling system, along with information regarding dow
 Information regarding WRF Model citations (including a DOI) can be found here: [https://www2.mmm.ucar.edu/wrf/users/citing_wrf.html](https://www2.mmm.ucar.edu/wrf/users/citing_wrf.html).
 
 The WRF Model is open-source code in the public domain, and its use is unrestricted. The name "WRF", however, is a registered trademark of the University Corporation for Atmospheric Research. The WRF public domain notice and related information may be found here: [https://www2.mmm.ucar.edu/wrf/users/public.html](https://www2.mmm.ucar.edu/wrf/users/public.html).
-
-

Registry/registry.diags

@@ -11,7 +11,7 @@ dimspec    np      2     namelist=num_press_levels         z     num_press_level
 rconfig   integer     p_lev_diags         namelist,diags   1            0   -    "flag to process vertical interp diagnostics: 0=nope, 1=yep" "flag"
 rconfig   integer     p_lev_diags_dfi     namelist,diags   1            0   -    "when doing p_level diags and dfi, turn off diags during 'non forecast'"
 rconfig   integer     num_press_levels    namelist,diags   1            0   -    "number of pressure levels to interpolate diagnostics to"    "index"
-rconfig   real            press_levels    namelist,diags   max_plevs    0   -    "array of pressure levels to interpolate diagnostics to"     "Pa"
+rconfig   real        press_levels        namelist,diags   max_plevs    0   -    "array of pressure levels to interpolate diagnostics to"     "Pa"
 rconfig   integer     use_tot_or_hyd_p    namelist,diags   1            2   -    "1=use total pressure, 2=use hydrostatic pressure"           "flag"
 rconfig   integer     extrap_below_grnd   namelist,diags   1            1   -    "1=no extrapolation, 2=extrapolate adiabatically"            "flag"
 rconfig   real        p_lev_missing       namelist,diags   1         -999   -    "missing values below ground, no extrapolation"              "constant"

Registry/registry.em_shared_collection

@@ -30,3 +30,6 @@ include registry.solar_fields
 include registry.diags
 include registry.iau
 include registry.CMAQ
+#ifdef SIMRADAR
+include registry.simradar
+#endif

doc/README.simradAR

@@ -0,0 +1,124 @@
+simrad.AR: Implementation of a radar emulator in WRF
+
+L. Fita, V. Galligani (CIMA-IFAECI), Aregntina
+
+This effort is absolutely in shoulders of giants and is based in multiple different works:
+
+ - Bessel, F. W., (1826) On the computation of geographical longitude and latitude from geodetic measurements, Astronomische Nachrichten (Astronomical Notes), Band 4, No. 86, Spalten 241-254, Altona 1826
+ - Leinonen, J., High-level interface to T-matrix scattering calculations: architecture, capabilities and limitations, Opt. Express, vol. 22, issue 2, 1655-1660 (2014), doi: 10.1364/OE.22.001655.
+ - Mishchenko, M. I. and L. D. Travis, T-matrix computations of light scattering by large spheroidal particles, Opt. Commun., vol. 109, 16-21 (1994)
+ - Mishchenko, M. I., L. D. Travis, and A. Macke, Scattering of light by polydisperse, randomly oriented, finite circular cylinders, Appl. Opt., vol. 35, 4927-4940 (1996)
+ - Vincenty, T., (1975) Direct and inverse solutions on the ellipsoid with application of nested equations, Survey Review, Vol. 23, No. 176, pp. 88-93.
+ - Vincenty, T., (1976) Correspondence: solutions of geodesics, survey Review, Vol. 23, No. 180, p. 294.
+ - Wielaard, D. J., M. I. Mishchenko, A. Macke, and B. E. Carlson, Improved T-matrix computations for large, nonabsorbing and weakly absorbing nonspherical particles and comparison with geometrical optics approximation, Appl. Opt., vol. 36, 4305-4313 (1997)
+
+
+There is a prototype off-line working version available from CIMA's GIT repository, from which the code is based on:
+https://git.cima.fcen.uba.ar/lluis.fita/simradar/-/wikis/home
+
+This code is available from Lluis Fita's WRF fork (git@github.com:LluisFB/WRF.git) in branch 'simradAR'
+
+== Usage ==
+A new section into the WRF's namelist calles 'simradar' has to be set up
+
+=== namelist ===
+simradar&
+/
+
+integer  nradar        1     amount of surface radar"                                ""
+integer  e_dist        namelist,simradar     1              800     rh{11}    "e_dist"            "radar grid points along the radii"                      ""
+integer  e_angle       namelist,simradar     1              376     rh{11}    "e_angle"           "radar amount of horizontal angles"                      ""
+integer  e_azimuth     namelist,simradar     1               15     rh{11}    "e_azimuth"         "radar amount of azimuths"                               ""
+integer  ngrid_radar   namelist,simradar     1             1500     rh{11}    "ngrid_radar"       "amount of radar cells into a given WRF grid point"      ""
+integer  ngrid3_radar  namelist,simradar     1            25000     rh{11}    "ngrid3_radar"      "amount of vertical radar cells into a given WRF grid point"      ""
+integer  nradar_ngrid  namelist,simradar     1            30000     rh{11}    "nradar_ngrid"      "nradar*ngrid_radar"                                     ""
+integer  nradar_ngrid_coord2  namelist,simradar 1         60000     rh{11}    "nradar_ngrid_coord2" "nradar*ngrid_radar*coord2"                            ""
+integer  nradar_coord2 namelist,simradar     1               40     rh{11}    "nradar_coord2"     "nradar*coord2"                                          ""
+integer  nradar_coord4 namelist,simradar     1              160     rh{11}    "nradar_coord4"     "nradar*coord4"                                          ""
+integer  ngrid_coord2  namelist,simradar     1            60000     rh{11}    "ngrid_coord2"      "ngrid_radar*coord2"                                     ""
+integer  nradar_azimuth namelist,simradar    1              300     rh{11}    "nradar_azimuth"    "nradar*e_azimuth"                                       ""
+integer  nradar_dist   namelist,simradar     1             1600     rh{11}    "nradar_dist"       "nradar*e_dist"                                          ""
+
+integer  Nbands        namelist,simradar     1                3     rh{11}    "Nbands"            "amount of radar bands"                                  ""
+integer  Nluticebins   namelist,simradar     1              200     rh{11}    "Nluticebins"       "amount of ice bins for LUT interpolations"              "#"
+integer  Nlutqxbins    namelist,simradar     1            50000     rh{11}    "Nlutqxbins"        "amount of water species bins for LUT interpolations"    "#"
+integer  Nlutlevs      namelist,simradar     1               15     rh{11}    "Nlutlevs"          "amount of azimuths for LUT interpolations"              ""
+
+character radar_name   namelist,simradar     max_nradar    (/'-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-'/) rh{11}    "radar_name"     "name of the radar"            ""
+real     radar_clon    namelist,simradar     max_nradar    (/0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/) rh{11}    "radar_clon"           "longitude of the location of the radar"  ""
+real     radar_clat    namelist,simradar     max_nradar    (/0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/) rh{11}    "radar_clat"           "latitude of the location of the radar"  ""
+real     radar_chgt    namelist,simradar     max_nradar    (/0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/) rh{11}    "radar_chgt"           "height of the location of the radar"      ""
+character radar_band   namelist,simradar     max_nradar    (/'-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-','-'/) rh{11}    "radar_band"     "band of the radar: 'c', 's', or 'x'"            ""
+integer  radar_scan    namelist,simradar     max_nradar    (/1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1/)     rh{11}    "radar_scan" "type of scan 1:"      ""
+real     radar_distmax namelist,simradar     max_nradar    (/250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250./) rh{11}    "radar_distmax"           "maximum distance from radar"  "km"
+real     radar_hgtmax  namelist,simradar     max_nradar    (/20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20.,20./) rh{11}    "radar_hgtmax"           "maximum height from radar"  "km"
+integer  em_ray        namelist,simradar     max_nradar    (/1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1/)     rh{11}    "em_ray"     "type of propagation of em ray. 1: 4/3 Ref"                                          ""
+real     radards       namelist,simradar     max_nradar    (/250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250.,250./)    rh{11}    "radards"           "nominal horizontal resolution of radar values in m"     ""
+integer  radinterp     namelist,simradar     1                1     rh{11}    "radinterp"         "method of interpolation to the radar space (1: lineal)" ""
+real     radii_dist    namelist,simradar     1              15.     rh{11}    "radii_dist"        "radii (km) at which arc resolution is the same as minimum ddist" "km"
+integer  eloption      namelist,simradar     1                1     rh{11}    "eloption"          "Ellipsoide"                                             ""
+real     missing_value namelist,simradar     1            1.e20     rh{11}    "missing_value"     "missing value"                                          ""
+integer  doppler_scheme namelist,simradar    1                1     rh{11}    "doppler_scheme"    "method to compute the doppler velocities"               ""
+
+integer  timedbg       namelist,simradar     1              720     rh{11}    "timedbg"           "frecuency of time-debugging"                 "minute"
+real     londbg        namelist,simradar     1             -999.    rh{11}    "londbg"            "longitude of grid-point debugging (-999., not used)"  "degree"
+real     latdbg        namelist,simradar     1             -999.    rh{11}    "latdbg"            "latitude of grid-point debugging (-999., not used)"   "degree"
+
+# LUT options (maybe from netCDF file!)
+character lutfilen     namelist,simradar     1  "TMATRIX_WRF_WSM6.nc"    -    "lutfilen"          "name of the netCDF file with the T-Matrix lool-up-tables" ""
+real     q_rain_min    namelist,simradar     1              1e-7    rh{11}    "q_rain_min"        "minimum rain mixing ratio to consider"                "kgkg-1"
+real     q_rain_max    namelist,simradar     1              0.02    rh{11}    "q_rain_max"        "maximum rain mixing ratio to consider"                "kgkg-1"
+real     qn_rain_min   namelist,simradar     1              1e-8    rh{11}    "qn_rain_min"       "minimum amount of nuclii for rain to consider"        "#"
+real     qn_rain_max   namelist,simradar     1             2.6e6    rh{11}    "qn_rain_max"       "maximum amount of nuclii for rain to consider"        "#"
+real     q_snow_min    namelist,simradar     1              1e-7    rh{11}    "q_snow_min"        "minimum snow mixing ratio to consider"                "kgkg-1"
+real     q_snow_max    namelist,simradar     1             0.018    rh{11}    "q_snow_max"        "maximum snow mixing ratio to consider"                "kgkg-1"
+real     qn_snow_min   namelist,simradar     1              1e-7    rh{11}    "qn_snow_min"       "minimum amount of nuclii for snow to consider"        "#"
+real     qn_snow_max   namelist,simradar     1             2.6e6    rh{11}    "qn_snow_max"       "maximum amount of nuclii for snow to consider"        "#"
+real     q_grau_min    namelist,simradar     1              1e-7    rh{11}    "q_grau_min"        "minimum grau mixing ratio to consider"                "kgkg-1"
+real     q_grau_max    namelist,simradar     1             0.018    rh{11}    "q_grau_max"        "maximum grau mixing ratio to consider"                "kgkg-1"
+real     qn_grau_min   namelist,simradar     1              1E-7    rh{11}    "qn_grau_min"       "minimum amount of nuclii for grau to consider"        "#"
+real     qn_grau_max   namelist,simradar     1             2.6e6    rh{11}    "qn_grau_max"       "maximum amount of nuclii for grau to consider"        "#"
+real     q_ice_min     namelist,simradar     1              1e-7    rh{11}    "q_ice_min"         "minimum cloud ice mixing ratio to consider"           "kgkg-1"
+real     q_ice_max     namelist,simradar     1             0.001    rh{11}    "q_ice_max"         "maximum cloud ice mixing ratio to consider"           "kgkg-1"
+real     qn_ice_min    namelist,simradar     1              1e-8    rh{11}    "qn_ice_min"        "minimum amount of nuclii for cloud ice to consider"   "#"
+real     qn_ice_max    namelist,simradar     1             2.6e6    rh{11}    "qn_ice_max"        "maximum amount of nuclii for cloud ice to consider"   "#"
+real     q_hail_min    namelist,simradar     1            -9999.    rh{11}    "q_hail_min"        "minimum hail mixing ratio to consider"                "kgkg-1"
+real     q_hail_max    namelist,simradar     1            -9999.    rh{11}    "q_hail_max"        "maximum hail mixing ratio to consider"                "kgkg-1"
+real     qn_hail_min   namelist,simradar     1            -9999.    rh{11}    "qn_hail_min"       "minimum amount of nuclii for hail to consider"        "#"
+real     qn_hail_max   namelist,simradar     1            -9999.    rh{11}    "qn_hail_max"       "maximum amount of nuclii for hail to consider"        "#"
+real     q_clou_min    namelist,simradar     1             1e-10    rh{11}    "q_clou_min"        "minimum clou mixing ratio to consider"                "kgkg-1"
+real     q_clou_max    namelist,simradar     1             0.013    rh{11}    "q_clou_max"        "maximum clou mixing ratio to consider"                "kgkg-1"
+integer  dim_ice       namelist,simradar     1               200    rh{11}    "dim_ice"           "dimension for ice bins"                               "-"
+integer  Nqxbins       namelist,simradar     1             50000    rh{11}    "Nqxbins"           "dimension for water species bins"                     "-"
+
+
+== Installation ==
+
+1.- Get the repository from a $WORKDIR
+git clone git@github.com:LluisFB/WRF.git
+
+2.- getting submoules
+git submodule update --init --recursive
+
+3.- Switching to the simrad.AR branch
+git checkout simradAR
+
+4.- Preparing installation (remember to define WRFIO_NCD_LARGE_FILE_SUPPORT=1 and NETCDF)
+./configure
+
+5.- Editing configure.wrf to add -DSIMRADAR pre-compilation flag
+cp configure.wrf configure.simradAR.orig
+cp configure.wrf configure.simradAR.wrf
+vim configure.simradAR.wrf
+diff configure.wrf.orig configure.simradAR.wrf
+216c216
+<                        \
+---
+>                       -DSIMRADAR \
+cp configure.simradAR.wrf configure.wrf
+
+6.- Compile as usual
+./compile em_real >& compile.log
+
+
+

dyn_em/start_em.F

@@ -52,6 +52,11 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
    USE module_model_constants
    USE module_avgflx_em, ONLY : zero_avgflx
 
+#ifdef SIMRADAR
+! simrad.AR
+   USE module_simradar_tools, ONLY: lut_vars
+#endif
+
    IMPLICIT NONE
    !  Input data.
    TYPE (domain)          :: grid
@@ -135,6 +140,9 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
    REAL :: dt_s, dx_km
 
    REAL :: max_mf, max_rot_angle
+#ifdef SIMRADAR
+   TYPE(lut_vars)                                        :: lutgrid
+#endif
    CALL get_ijk_from_grid ( grid ,                              &
                            ids, ide, jds, jde, kds, kde,        &
                            ims, ime, jms, jme, kms, kme,        &
@@ -1271,6 +1279,24 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
                       ,grid%NC_CU, grid%NI_CU, grid%NR_CU, grid%NS_CU,grid%CCN_CU                    &
                       ,grid%alevsiz_cu,grid%num_months,grid%no_src_types_cu,grid%aeromcu             &
                       ,grid%aeropcu,grid%EFCG,grid%EFCS,grid%EFIG,grid%EFIS,grid%EFSG,grid%EFSS      &
+#ifdef SIMRADAR
+! simrad.AR
+                      ,grid%xlong_v, grid%xlat_u                                                     &
+                      ,config_flags%nradar                                                           &
+                      ,config_flags%e_angle, config_flags%e_azimuth, config_flags%e_dist             &
+                      ,config_flags%nradar_coord4, config_flags%nradar_ngrid_coord2                  &
+                      ,config_flags%nradar_ngrid, model_config_rec%ngrid_coord2                      &
+                      ,config_flags%londbg, config_flags%latdbg                                      &
+                      ,grid%ijraddbg                                                                 &
+                      ,grid%radars_dist, grid%radars_angle, grid%radars_azimuth                      &
+                      ,grid%ellipse_a, grid%ellipse_f, grid%ellipse_name                             &
+                      ,grid%gatehdist, grid%inradar, grid%radarnum                                   &
+                      ,grid%radarlons, grid%radarlats, grid%radardist                                & 
+                      ,grid%radarhgts, grid%radarlhgts                                               &
+                      ,grid%radarlon_bnds, grid%radarlat_bnds, grid%radarhgt_bnds                    &
+                      ,grid%ngridradar, grid%ijgridradar, grid%ijdistgridradar                       &
+                      ,lutgrid                                                                       &
+#endif
                       )
        ENDDO ! loop of tiles for phy_init
 
@@ -1297,6 +1323,68 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
 
    CALL wrf_debug ( 100 , 'start_domain_em: After call to phy_init' )
 
+#ifdef SIRADAR
+    !! Filling grid with lutgrid values
+    grid%lutelevations = grid%elevations
+    grid%lut_Zh_RAIN = lutgrid%Zh_RAIN
+    grid%lut_Zdr_RAIN = lutgrid%Zdr_RAIN
+    grid%lut_LDR_RAIN = lutgrid%LDR_RAIN
+    grid%lut_Aih_RAIN = lutgrid%Aih_RAIN
+    grid%lut_Aiv_RAIN = lutgrid%Aiv_RAIN
+    grid%lut_KDP_RAIN = lutgrid%KDP_RAIN
+
+    grid%lut_Zh_SNOW = lutgrid%Zh_SNOW
+    grid%lut_Zdr_SNOW = lutgrid%Zdr_SNOW
+    grid%lut_LDR_SNOW = lutgrid%LDR_SNOW
+    grid%lut_Aih_SNOW = lutgrid%Aih_SNOW
+    grid%lut_Aiv_SNOW = lutgrid%Aiv_SNOW
+    grid%lut_KDP_SNOW = lutgrid%KDP_SNOW
+
+    grid%lut_Zh_GRAU = lutgrid%Zh_GRAU
+    grid%lut_Zdr_GRAU = lutgrid%Zdr_GRAU
+    grid%lut_LDR_GRAU = lutgrid%LDR_GRAU
+    grid%lut_Aih_GRAU = lutgrid%Aih_GRAU
+    grid%lut_Aiv_GRAU = lutgrid%Aiv_GRAU
+    grid%lut_KDP_GRAU = lutgrid%KDP_GRAU
+
+    grid%lut_Zh_ICE = lutgrid%Zh_ICE
+    grid%lut_Zdr_ICE = lutgrid%Zdr_ICE
+    grid%lut_LDR_ICE = lutgrid%LDR_ICE
+    grid%lut_Aih_ICE = lutgrid%Aih_ICE
+    grid%lut_Aiv_ICE = lutgrid%Aiv_ICE
+    grid%lut_KDP_ICE = lutgrid%KDP_ICE
+
+    grid%lut_Zh_HAIL = lutgrid%Zh_HAIL
+    grid%lut_Zdr_HAIL = lutgrid%Zdr_HAIL
+    grid%lut_LDR_HAIL = lutgrid%LDR_HAIL
+    grid%lut_Aih_HAIL = lutgrid%Aih_HAIL
+    grid%lut_Aiv_HAIL = lutgrid%Aiv_HAIL
+    grid%lut_KDP_HAIL = lutgrid%KDP_HAIL
+
+    grid%N_tocompute = lutgrid%N_tocompute
+    grid%tocompute = lutgrid%tocompute
+    grid%N_tcmoist = lutgrid%Nmoist
+    grid%tcmoist = lutgrid%imoist
+
+    grid%q_rain_vec = lutgrid%q_rain_vec
+    grid%q_snow_vec = lutgrid%q_snow_vec
+    grid%q_grau_vec = lutgrid%q_grau_vec
+    grid%q_ice_vec = lutgrid%q_ice_vec
+    grid%q_hail_vec = lutgrid%q_rhail_vec
+
+    ! microphysics-derived Denisities
+    grid%rhoo = lutgrid%rhoo
+    grid%rho_r = lutgrid%rho_r
+    grid%n0r = lutgrid%n0r
+    grid%rho_s = lutgrid%rho_s
+    grid%n0s = lutgrid%n0s
+    grid%rho_g = lutgrid%rho_g
+    grid%n0g = lutgrid%n0g
+    grid%rho_h = lutgrid%rho_h
+    grid%n0h = lutgrid%n0h
+
+#endif
+
    IF (config_flags%do_avgflx_em .EQ. 1) THEN
       WRITE ( message , FMT = '("start_em: initializing avgflx on domain ",I3)' ) &
            & grid%id