Code for Evaluation of cloud microphysics schemes using a satellite radar simulator over the Oklahoma and Korea regions

Intro The zip file contains the source codes for the WRF simulations, and the plot figures used in our manuscript "Evaluation of cloud microphysics schemes using a satellite radar simulator over the Oklahoma and Korea regions" The Joint Simulator is available on the Joint simulator group's wiki (https://sites.google.com/site/jointsimulator/home) Contents of the ZIP file WRF The source codes for the WRF simulations are as follows. module_ra_effectiveradius.F namelist.wps.kor, namelist.wps.usa namelist.input Forecast_period_kor.xlsx, Forecast_period_usa.xlsx The file namelist.wps.domain are the namelist used for WPS in each domain, namelist.input is the namelist used for running WRF for both domains, and the xlsx files specifiy the simulation period for each domain. Python scirpt The scripts for pre-processing HDF5 are as follows. preku.py preka.py q10.py plot_cfads_obs.py prewrf.py plot_cfads_wrf.py prejsim.py plot_cfads_jsim.py diff.py calcr.py OBS The preku.py, and preka.py preprocess the 2ADPR data for each band by counting number of pixels where the "precipRateNearSurface" variable exceeds 10 mm h-1. The q10.py extracts the "zFactorFinal" values from the files in which the number of pixels exceeding 10 mm h-1 ranks in the top 10% across all cases The plot_cfads_obs.py plot the Ku-band, and Ka-band Contoured Frequency by Altitude Diagrams (CFADs) using results of q10.py WRF The prewrf.py processes the wrfout data by extracting the "REFL_10CM" profiles into a csv file when the "RAINNCV" variable exceeds 10 mm h-1. The plot_cfads_wrf.py plot the S-Band CFADs using results of prewrf.py JSIM The prejsim.py processes the JSIM output data by extracting the "dbz_13", and "dbz_35" profiles into a csv file when the "RAINNCV" variable exceeds 10 mm h-1. The plot_cfads_jsim.py plot the simultated Ku-band, and Ka-band Contoured Frequency by Altitude Diagrams (CFADs) using results of prejsim.py The diff.py plot the difference between Ku-band and Ka-band. The calcr.py calculates the correlation (Pearson R) between the obs and the model.