search

anikfal / postwrf

Visualization and postprocessing of the WRF and ERA5 data
MIT License
54 stars 11 forks source link

Contour shading #53

Closed SownalChand closed 3 weeks ago

SownalChand commented 3 weeks ago

Hello sir,

I have try to draw the precipitation by using RAINC+RAINNC.

I also tried to reinstall the postwrf. But still not working.

Please provide details on why is this happening

Kind regards Sownalc

anikfal commented 3 weeks ago

Hi. Please provide more information, including the errors and namelist.wrf. I'll check them out.

SownalChand commented 3 weeks ago

================== WRF/ERA5 DATA EXTRACTION (ON SURFACE) ====================

WRF_Extract_ON_OFF = 0 ERA5_Extract_ON_OFF = 0

Variable_names = sp Location_names = Anzali, Berlin, Cairo, Houston, Shanghai Location_latitudes = 37.46, 52.52, 30.04, 29.76, 31.23 Location_longitudes = 49.48, 13.40, 31.23, -95.37, 121.47

+Vertical_profile_plot_on_off = 0 +Vprofile_X_axis_decimals = 3

============================= WRF CONTOUR PLOT ==============================

CONTOUR_ON_OFF = 1

!!! Variable_1 by line contours !!! 1st_Variable_on_off = 0 1st_Variable_name = T2 1st_Variable_pressure_level = 0 1st_Variable_intervals = 5 1st_Variable_line_color = red 1st_Variable_line_thickness = 0.5 1st_Variable_label_size = 0.5

!!! Variable_2 by (broken) line contours !!! 2nd_Variable_on_off = 0 2nd_Variable_name = pvo 2nd_Variable_pressure_level = 700 2nd_Variable_intervals = 0 2nd_Variable_line_color = red 2nd_Variable_line_thickness = 0.5 2nd_Variable_label_size = 0.4 2nd_Variable_brokenline_on_off = 1

!!! Variable_3 by shaded contours !!! 3rd_Variable_on_off = 1 3rd_Variable_name = RAINC+RAINNC 3rd_Variable_pressure_level = 0,

Wind_Vectors_on_off = 1 Wind_Vectors_pressure_level = 850 Wind_Vectors_speed (m/s) = 5 Wind_Vectors_density = 0.5 Wind_Vectors_color = black Wind_Vectors_thickness = 0.5 Wind_Barbs_on_off = 1

Location_point_mark_on_off = 0 Location_labels_on_off = 1 Location_labels_size = 0.2 Location_labels_color = black Shapefile_on-off = 0 Shapefile_path = /home/anikfal/shapefiles/anazali_swamp.shp HighRes_coastline_on_off = 0 HighRes_coastline_data_path = /home/anikfal/PostWRF/highres

=========================== WRF CROSS SECTION PLOT ==========================

CROSSSECTION_ON-OFF = 0

CrossSection_path_on_off = 1 Topographic_map_on_off = 1

!!! 3D Variable_1 by line contours !!! 1st_var_on_off = 0 1st_var_name = (ua^2 + va^2)^0.5 1st_var_intervals = 10 1st_var_line_color = black 1st_var_line_thickness = 0.5 1st_var_label_size = 0.5

!!! 3D Variable_2 by (broken) line contours !!! 2nd_var_on_off = 1 2nd_var_name = tc 2nd_var_intervals = 20 2nd_var_line_color = red 2nd_var_line_thickness = 0.5 2nd_var_label_size = 0.5 2nd_var_brokenline_on_off = 1

!!! 3D Variable_3 by shaded contours !!! 3rd_var_on_off = 1 3rd_var_name = QCLOUD

======================== WRF STATISTICAL DIAGRAMS ===========================

STATISTICAL_DIAGRAMS_ON_OFF = 0

Timeseries_ON_OFF = 0 Timeseries_Line_on_off = 1 Timeseries_Marker_on_off = 1

Scatterplot_ON_OFF = 1

Taylor_diagram_ON_OFF = 0 Taylor_labels_size = 0.5 Taylor_markers_size = 0.5

============================== WRF DOMAIN PLOT ==============================

Domain_On-Off = 0

namelist.wps_path (schematic) = /home/sownal/Build_WRF/postwrf-master/vut_shp/vutvut_admbndl_admALL_spc_itos_20180824.shp background_color (schematic) = forestgreen lables_On-Off (schematic) = 0

+Domain_on_off (topography) = 1 +number_of_domains (topography) = 2 +box_color = white

=========================== WRF WIND ROSE DIAGRAM ===========================

windrose_On-Off = 0 wind_location_name = shiraz wind_location_latitude = 29.5 wind_location_longitude = 52.5

============================ WRF SKEWT DIAGRAM ==============================

skewT_On-Off = 0 skewt_location_name = shiraz skewt_location_latitude = 29.5 skewt_location_longitude = 52.5

========================== WRF CONVERT TO GEOTIFF ===========================

Geotiff_ON_OFF = 0 Variable_name_to_Geotiff = rh Geotiff_pressure_level = 500 tiff_resolution (km) = 10

========================== RTTOV-WRF INPUT/OUTPUT ===========================

RTTOV_On-Off = 0

WRF2RTTOV_profiles_OnOff = 0 aerosol_profile_OnOff = 0

RTTOV_OUTPUT_OnOff = 1 rttov_output_prefix = outputexample Brightness_temperature = 1 Reflectance = 0 Radiance = 0 Surface_emissivity = 0

============================= ERA5 CONTOUR PLOT =============================

ERA5_ON_OFF = 0

!!! Variable_1 by line contours !!! 1st_ERA5_Var_on_off = 1 1st_ERA5_Var_name = z 1st_ERA5_Var_pressure_level = 500 1st_ERA5_Var_intervals = 5 1st_ERA5_Var_line_color = blue 1st_ERA5_Var_line_thickness = 0.5 1st_ERA5_Var_label_size = 0.3

!!! Variable_2 by (broken) line contours !!! 2nd_ERA5_Var_on_off = 0 2nd_ERA5_Var_name = t 2nd_ERA5_Var_pressure_level = 500 2nd_ERA5_Var_intervals = 0 2nd_ERA5_Var_line_color = red 2nd_ERA5_Var_line_thickness = 0.5 2nd_ERA5_Var_label_size = 0.3 2nd_ERA5_Var_brokenline_on_off = 0

!!! Variable_3 by shaded contours !!! 3rd_ERA5_Var_on_off = 1 3rd_ERA5_Var_name = t 3rd_ERA5_Var_pressure_level = 0

Wind_ERA5_on_off = 1 Wind_ERA5_pressure_level = 500 Wind_ERA5_speed (m/s) = 7 Wind_ERA5_density = 0.5 Wind_ERA5_color = black Wind_ERA5_thickness = 0.5 wind_barbs_on_off = 1

ERA5_point_mark_on_off = 0 ERA5_labels_on_off = 0 ERA5_labels_size = 0.3 ERA5_labels_color = black ERA5_Shapefile_on_off = 0 ERA5_shapefile_path = /home/anikfal/shapefiles/anzali_swamp.shp

ERA5_SubDomain_on_off = 0 Min_lat = 30 Max_lat = 50 Min_lon = 55 Max_lon = 65

============================= GENERAL SETTINGS ==============================

spin-up_time (hours) = 0 averaging_on_off = 0 average_time (hours) = 12 panel_ploton_off = 0 panel_plotrows_number = 2 panel_plot__columns_number = 3 smooth_maps_on_off = 1 borderline_on_off = 1 map_borderline_color = black map_borderline_thickness = 0.5 map_filled_ocean_on_off = 1 map_gridlines_on_off = 0 plot_titles_on_off = 1

SownalChand commented 3 weeks ago

wrfout_d01_2023-10-20_000000 https://drive.google.com/file/d/1H8Rib3-Dd57SMWHN4I9KRGDfMsKI2LlK/view?usp=drive_web Dear sir,

Please find attached the wrfout file.

Thank you for the assistance.

Kind regards Sownalc

On Sun, Oct 13, 2024 at 11:33 PM Amirhossein Nikfal < @.***> wrote:

Hi. Please provide more information, including the errors and namelist.wrf. I'll check them out.

— Reply to this email directly, view it on GitHub https://github.com/anikfal/postwrf/issues/53#issuecomment-2408941962, or unsubscribe https://github.com/notifications/unsubscribe-auth/AVN7SISPCKV6ZJRIYU4SEG3Z3JK65AVCNFSM6AAAAABP3LCPC2VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDIMBYHE2DCOJWGI . You are receiving this because you authored the thread.Message ID: @.***>

SownalChand commented 3 weeks ago

sir,

I have emailed you the namelist.wrf and wrfout file.

thank you and kind regards sownalc

anikfal commented 3 weeks ago

This is a sample output image from your file and namelist: image

They seem to be okay. You can also provide the errors raised after you run postwrf.

SownalChand commented 3 weeks ago

Sir,

./pos twrf.sh

======================== PostWRF ============================ Run './postwrf.sh -h' to get help Reference any problems on GitHub <github.com/anikfal>

   CONTOUR PLOT                CONTOUR PLOT

There are multiple wrf-files in the current directory:

1) wrfout_d01_2023-10-20_00:00:00 2) wrfout_d01_2023-10-21_00:00:00 3) wrfout_d01_2023-10-22_00:00:00 4) wrfout_d01_2023-10-23_00:00:00 5) wrfout_d01_2023-10-24_00:00:00 6) wrfout_d01_2023-10-25_00:00:00 7) wrfout_d01_2024-10-06_00:00:00 8) wrfout_d01_2024-10-07_00:00:00 9) wrfout_d01_2024-10-08_00:00:00 10) wrfout_d01_2024-10-09_00:00:00 11) wrfout_d02_2023-10-20_00:00:00 12) wrfout_d02_2023-10-21_00:00:00 13) wrfout_d02_2023-10-22_00:00:00 14) wrfout_d02_2023-10-23_00:00:00 15) wrfout_d02_2023-10-24_00:00:00 16) wrfout_d02_2023-10-25_00:00:00 17) wrfout_d02_2024-10-06_00:00:00 18) wrfout_d02_2024-10-07_00:00:00 19) wrfout_d02_2024-10-08_00:00:00 20) wrfout_d02_2024-10-09_00:00:00

Enter the number of the favored file: 1

wrfout_d01_2023-10-20_00:00:00 has been selected

Plotting contour maps ...

Specify Time_Step(s) Between The Images (Default=1): Time_Step Has Been Set To 1

Specify The Method Of Drawing Contours (1 or 2) for RAINC+RAINNC: 1) Automatic 2) Manual

? 2

Contour Lines Min Value (press Enter for the best value): Got no input. Best Min value will be set automatically.

Contour Lines Max Value (press Enter for the best value): Got no input. Best Max value will be set automatically.

Select a color pattern from the list below (a number from 1 to 10) for RAINC+RAINNC:: 1) Rainbow-start_from_blue 2) Rainbow-start_from_red 3) Rainbow-start_from_white 4) Blue..Red 5) Blue 6) Red 7) Green 8) White..Yellow..Orange..Red 9) White..Black 10) Black..White

? 5

Format Of The Images: 1, 2, 3, or 4? 1) x11 3) png 2) pdf 4) animated_gif

? 3

Name of the Output File? (Press Enter for the default name):OOOout Please wait. NCL is processing ...

Please wait. NCL is processing ...

Contour 3: Name of the polynomial RAINC+RAINNC ? Rain Unit-scale of the polynomial RAINC+RAINNC ? mm (0) Working on time: 2023-10-20_00:00:00 warning:ContourPlotInitialize: scalar field is constant; no contour lines will appear; use cnConstFEnableFill to enable fill warning:ContourPlotSetValues: Data values out of range of levels set by MANUALLEVELS mode (0) Working on time: 2023-10-20_01:00:00 (0) Working on time: 2023-10-20_02:00:00 (0) Working on time: 2023-10-20_03:00:00 (0) Working on time: 2023-10-20_04:00:00 (0) Working on time: 2023-10-20_05:00:00 (0) Working on time: 2023-10-20_06:00:00 (0) Working on time: 2023-10-20_07:00:00 (0) Working on time: 2023-10-20_08:00:00 (0) Working on time: 2023-10-20_09:00:00 (0) Working on time: 2023-10-20_10:00:00 (0) Working on time: 2023-10-20_11:00:00 (0) Working on time: 2023-10-20_12:00:00 (0) Working on time: 2023-10-20_13:00:00 (0) Working on time: 2023-10-20_14:00:00 (0) Working on time: 2023-10-20_15:00:00 (0) Working on time: 2023-10-20_16:00:00 (0) Working on time: 2023-10-20_17:00:00 (0) Working on time: 2023-10-20_18:00:00 (0) Working on time: 2023-10-20_19:00:00 (0) Working on time: 2023-10-20_20:00:00 (0) Working on time: 2023-10-20_21:00:00 (0) Working on time: 2023-10-20_22:00:00 (0) Working on time: 2023-10-20_23:00:00 (0) Working on time: 2023-10-21_00:00:00 (base) sownal@LAPTOP-DSEPU3MN:~/Build_WRF/postwrf-master

Oout-wrfout_d01 000001

Sir attached is the output and png image file.

please can you check why its saying constant value for rain ?

thanking you in advance Sownalc

SownalChand commented 3 weeks ago

Oout-wrfout_d01 000003 Oout-wrfout_d01 000002

anikfal commented 3 weeks ago

In the line Specify The Method Of Drawing Contours (1 or 2) for RAINC+RAINNC:, choose the first option (Automatic).

SownalChand commented 3 weeks ago

Sir,

I will check and confirm you soon as i am back to my desk.

Kind regards sownalc

SownalChand commented 3 weeks ago

Dear sir,

========================  PostWRF  ============================
              Run './postwrf.sh -h' to get help
     Reference any problems on GitHub <github.com/anikfal>
===============================================================

---------------------------------------------------------
       CONTOUR PLOT                CONTOUR PLOT
---------------------------------------------------------

There are multiple wrf-files in the current directory:

   1) wrfout_d01_2023-10-20_00:00:00
   2) wrfout_d01_2023-10-21_00:00:00
   3) wrfout_d01_2023-10-22_00:00:00
   4) wrfout_d01_2023-10-23_00:00:00
   5) wrfout_d01_2023-10-24_00:00:00
   6) wrfout_d01_2023-10-25_00:00:00
   7) wrfout_d01_2024-10-06_00:00:00
   8) wrfout_d01_2024-10-07_00:00:00
   9) wrfout_d01_2024-10-08_00:00:00
   10) wrfout_d01_2024-10-09_00:00:00
   11) wrfout_d02_2023-10-20_00:00:00
   12) wrfout_d02_2023-10-21_00:00:00
   13) wrfout_d02_2023-10-22_00:00:00
   14) wrfout_d02_2023-10-23_00:00:00
   15) wrfout_d02_2023-10-24_00:00:00
   16) wrfout_d02_2023-10-25_00:00:00
   17) wrfout_d02_2024-10-06_00:00:00
   18) wrfout_d02_2024-10-07_00:00:00
   19) wrfout_d02_2024-10-08_00:00:00
   20) wrfout_d02_2024-10-09_00:00:00

Enter the number of the favored file: 6

wrfout_d01_2023-10-25_00:00:00 has been selected

Plotting  contour maps ...

Specify Time_Step(s) Between The Images (Default=1):
Time_Step Has Been Set To 1

 Specify The Method Of Drawing Contours (1 or 2) for RAINC+RAINNC:
1) Automatic
2) Manual
#? 1

Contour Lines Will Be Set Automatically ...

Select a color pattern from the list below (a number from 1 to 10) for RAINC+RAINNC::
1) Rainbow-start_from_blue       5) Blue                         9) White..Black
2) Rainbow-start_from_red        6) Red                         10) Black..White
3) Rainbow-start_from_white      7) Green
4) Blue..Red                     8) White..Yellow..Orange..Red
#? 5

 Format Of The Images: 1, 2, 3, or 4?
1) x11
2) pdf
3) png
4) animated_gif
#? 1
Please wait. NCL is processing ...

Please wait. NCL is processing ...

Contour 3: Name of the polynomial RAINC+RAINNC ?
Rain
Unit-scale of the polynomial RAINC+RAINNC ?
mm
(0)     Working on time: 2023-10-25_00:00:00
warning:ContourPlotInitialize: scalar field is constant; no contour lines will appear; use cnConstFEnableFill to enable fill
(0)     Working on time: 2023-10-25_01:00:00
(0)     Working on time: 2023-10-25_02:00:00
(0)     Working on time: 2023-10-25_03:00:00
(0)     Working on time: 2023-10-25_04:00:00
(0)     Working on time: 2023-10-25_05:00:00
(0)     Working on time: 2023-10-25_06:00:00

Still having the same issue and this time the warning messge is seen as above:

Sir any assistance will be greatly appreciated

Kind regards sownalc

anikfal commented 3 weeks ago

warning:ContourPlotInitialize: scalar field is constant; no contour lines will appear; use cnConstFEnableFill to enable fill

This warning message is not a serious issue. It simply means that, in the first time-slot, there is no rain values available to plot.

After you selected x11 as the Format of the images, you need to click on each image, so the next image will be plotted.

Next time try png or pdf as the Format of the images.

SownalChand commented 3 weeks ago

Dear sir,

The images are simply like below: Rain-wrfout_d01 000001 Rain-wrfout_d01 000002 Rain-wrfout_d01 000003

It doesnt appear also when i try to use png or pdf.

I am still confused as how this could happen?

hoping for your expertise in this area to assist me.

Kind regards Sownalc

SownalChand commented 3 weeks ago

sir i think its due to that i am running Ubuntu using wsl on windows 11.

I am closing the issue shortly

Kind regards sownalc

SownalChand commented 3 weeks ago

Thanks sir for all the assistance

Kind regards sownalc

anikfal commented 3 weeks ago

If you are serious in atmospheric modeling and its associated tasks, Windows WSL can never be a good option. You can at least install a virtual Linux OS on your Windows, using VMware or similar software.

SownalChand commented 3 weeks ago

Thanks, I'll check it out.

On Tue, Oct 15, 2024, 22:45 Amirhossein Nikfal @.***> wrote:

If you are serious in atmospheric modeling and its associated tasks, Windows WSL can never be a good option. You can at least install a virtual Linux OS on your Windows, using VMware or similar software.

— Reply to this email directly, view it on GitHub https://github.com/anikfal/postwrf/issues/53#issuecomment-2413540173, or unsubscribe https://github.com/notifications/unsubscribe-auth/AVN7SITS53BTU5CKMZU3HILZ3TW4TAVCNFSM6AAAAABP3LCPC2VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDIMJTGU2DAMJXGM . You are receiving this because you modified the open/close state.Message ID: @.***>