search

anikfal / postwrf

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

ERA5 and WRF simulation #11

Closed sstikle closed 2 years ago

sstikle commented 2 years ago

Hi, I want to analyse the simulated WRF output with ERA5. Could you please tell the simple method. I am new to WRF and intial WRF output share my previous colleague.

Thanks Suvarna

anikfal commented 2 years ago

Hi,

Please look at the PostWRF documentation. It has several examples which can help you to easily visualize and post-process the WRF (and ERA5) outputs: https://postwrf.readthedocs.io/en/latest/

In case of further problems, you can discuss it here.

sstikle commented 2 years ago

Ok. Thanks. I have already started to read it.

sstikle commented 2 years ago

when I run ./postwrf.sh Error massage is as follows

More than one task is on Select only one task or section in namelist.wrf and run again

anikfal commented 2 years ago

The message is clear. You can activate only one task-section at a time. For example, if both WRF_Extract_ON_OFF = 1 and CONTOUR_ON_OFF = 1, you will face this error message

If you want to extract data, only WRF_Extract_ON_OFF = 1, and turn off other sections.

sstikle commented 2 years ago

ok, Thanks

sstikle commented 2 years ago

Its run successfully. I have choose Bilinear. Not able to see the output ascii file.

anikfal commented 2 years ago

outputs will be stored in new directories Their names are started with output

sstikle commented 2 years ago

Output folder is saved but it is empty. I will try it again and get back to you.

anikfal commented 2 years ago

You can send me any warning or error messages you get. Is your data WRF or ERA5? For the WRF data you shouldn't face such problems. But for ERA5, you may have chosen some variables that are not inside the file.

sstikle commented 2 years ago

It is for WRF Now not able to run upto interpolation , although work path is correct

================ PostWRF (Version 1.1) ===================== Run './postwrf.sh -h' to get help
In case of problems, please contact <ah.nikfal....com> 

  EXTRACT VARIABLE             EXTRACT VARIABLE
  No WRF output files in the current directory.
  You can link or copy one or more files to the current directory.
anikfal commented 2 years ago

Where is your WRF output files? Copy or link them into the PostWRF directory, and run postwrf.sh again.

Please pay attention to the messages. It says:

No WRF output files in the current directory. You can link or copy one or more files to the current directory.

So it is actually mentioning how to deal with this problem.

sstikle commented 2 years ago

I already linked it. my data is in same directory which I linked but still this massage is coming.

anikfal commented 2 years ago

What are the names of your WRF output files?

You shouldn't change their names. They must have their generic names, starting with wrfout_d01 or wrfout_d02, ...

Starting of the names must be always wrfout_d0. You can change the rest of the names however.

sstikle commented 2 years ago

WRF outfiles starting as wrfout_d01

anikfal commented 2 years ago

Are you sure about your WRF file? Did you run it correctly? Maybe your desired variable is not simulated and stored in the WRF output.

Can you send your namelist.wrf or copy the contents of WRF/ERA5 DATA EXTRACTION (ON SURFACE)?

sstikle commented 2 years ago

I am sure about my WRF data. I repeat the entire process again but still problem exist. Surprisingly yesterday it run on same path but as I informed you outfolder was empty. In your namelist file I have made following changes. WRF_Extract_ON_OFF = 1 ERA5_Extract_ON_OFF = 0

Variable_names = sp Location_names = Indore Location_latitudes = 22.7196 Location_longitudes = 75.8577

+Vertical_profile_plot_on_off = 0 +Vprofile_X_axis_decimals = 3

sstikle commented 2 years ago

Do you have any video of doc?

anikfal commented 2 years ago

Variable_names = sp

What is sp? I have not seen this variable before. If this variable is inside your wrfoutput, you should see it by the command below:

./postwrf.sh -f wrfout_d01_2022-01-15_all

List of variables inside wrfout_d01_2022-01-15_all:

Times, AOD_OUT, AOD2D_OUT, ATOP2D_OUT, ICN_DIAG, NC_DIAG, CLDFRA2, UST_T, ROUGH_COR, SMOIS_COR, DRYDEPVEL, E_VASH1, E_VASH2, DRY_DEP_LEN, EBIO_ISO, EBIO_API, LAI_VEGMASK, dvel_o3, PM2_5_DRY, PM10, DMS_0, PHOTR2, PHOTR4, SNU, SAC, so2, vash_1, vash_2, vash_3, vash_4, vash_5, vash_6, vash_7, vash_8, vash_9, vash_10, PV, GAMN2O5, CN2O5, KN2O5, YCLNO2, XLAT, XLONG, LU_INDEX, ZNU, ZNW, ZS, DZS, VAR_SSO, BATHYMETRY_FLAG, U, V, W, PH, PHB, T, THM, HFX_FORCE, LH_FORCE, TSK_FORCE, HFX_FORCE_TEND, LH_FORCE_TEND, TSK_FORCE_TEND, MU, MUB, NEST_POS, P, PB, FNM, FNP, RDNW, RDN, DNW, DN, CFN, CFN1, THIS_IS_AN_IDEAL_RUN, P_HYD, Q2, T2, TH2, PSFC, U10, V10, RDX, RDY, AREA2D, DX2D, RESM, ZETATOP, CF1, CF2, CF3, ITIMESTEP, XTIME, QVAPOR, QCLOUD, QRAIN, SHDMAX, SHDMIN, SNOALB, TSLB, SMOIS, SH2O, SEAICE, XICEM, SFROFF, UDROFF, IVGTYP, ISLTYP, VEGFRA, GRDFLX, ACGRDFLX, ACSNOM, SNOW, SNOWH, CANWAT, SSTSK, WATER_DEPTH, COSZEN, LAI, DTAUX3D, DTAUY3D, DUSFCG, DVSFCG, VAR, CON, OA1, OA2, OA3, OA4, OL1, OL2, OL3, OL4, MAPFAC_M, MAPFAC_U, MAPFAC_V, MAPFAC_MX, MAPFAC_MY, MAPFAC_UX, MAPFAC_UY, MAPFAC_VX, MF_VX_INV, MAPFAC_VY, F, E, SINALPHA, COSALPHA, HGT, TSK, P_TOP, GOT_VAR_SSO, T00, P00, TLP, TISO, TLP_STRAT, P_STRAT, MAX_MSFTX, MAX_MSFTY, RAINC, RAINSH, RAINNC, SNOWNC, GRAUPELNC, HAILNC, CLDFRA, SWDOWN, GLW, SWNORM, ACSWUPT, ACSWUPTC, ACSWDNT, ACSWDNTC, ACSWUPB, ACSWUPBC, ACSWDNB, ACSWDNBC, ACLWUPT, ACLWUPTC, ACLWDNT, ACLWDNTC, ACLWUPB, ACLWUPBC, ACLWDNB, ACLWDNBC, SWUPT, SWUPTC, SWDNT, SWDNTC, SWUPB, SWUPBC, SWDNB, SWDNBC, LWUPT, LWUPTC, LWDNT, LWDNTC, LWUPB, LWUPBC, LWDNB, LWDNBC, OLR, XLAT_U, XLONG_U, XLAT_V, XLONG_V, ALBEDO, CLAT, ALBBCK, EMISS, NOAHRES, TMN, XLAND, UST, PBLH, HFX, QFX, LH, ACHFX, ACLHF, SNOWC, SR, SAVE_TOPO_FROM_REAL, ISEEDARR_SPPT, ISEEDARR_SKEBS, ISEEDARR_RAND_PERTURB, ISEEDARRAY_SPP_CONV, ISEEDARRAY_SPP_PBL, ISEEDARRAY_SPP_LSM, ISNOW, TV, TG, CANICE, CANLIQ, EAH, TAH, CM, CH, FWET, SNEQVO, ALBOLD, QSNOWXY, QRAINXY, WSLAKE, ZWT, WA, WT, TSNO, ZSNSO, SNICE, SNLIQ, LFMASS, RTMASS, STMASS, WOOD, STBLCP, FASTCP, XSAI, TAUSS, T2V, T2B, Q2V, Q2B, TRAD, NEE, GPP, NPP, FVEG, QIN, RUNSF, RUNSB, ECAN, EDIR, ETRAN, FSA, FIRA, APAR, PSN, SAV, SAG, RSSUN, RSSHA, BGAP, WGAP, TGV, TGB, CHV, CHB, SHG, SHC, SHB, EVG, EVB, GHV, GHB, IRG, IRC, IRB, TR, EVC, CHLEAF, CHUC, CHV2, CHB2, CHSTAR, SMCWTD, RECH, QRFS, QSPRINGS, QSLAT, ACINTS, ACINTR, ACDRIPR, ACTHROR, ACEVAC, ACDEWC, FORCTLSM, FORCQLSM, FORCPLSM, FORCZLSM, FORCWLSM, ACRAINLSM, ACRUNSB, ACRUNSF, ACECAN, ACETRAN, ACEDIR, ACQLAT, ACQRF, ACETLSM, ACSNOWLSM, ACSUBC, ACFROC, ACFRZC, ACMELTC, ACSNBOT, ACSNMELT, ACPONDING, ACSNSUB, ACSNFRO, ACRAINSNOW, ACDRIPS, ACTHROS, ACSAGB, ACIRB, ACSHB, ACEVB, ACGHB, ACPAHB, ACSAGV, ACIRG, ACSHG, ACEVG, ACGHV, ACPAHG, ACSAV, ACIRC, ACSHC, ACEVC, ACTR, ACPAHV, ACSWDNLSM, ACSWUPLSM, ACLWDNLSM, ACLWUPLSM, ACSHFLSM, ACLHFLSM, ACGHFLSM, ACPAHLSM, ACCANHS, SOILENERGY, SNOWENERGY, ACEFLXB, GRAIN, GDD, CROPCAT, PGS, QTDRAIN, IRNUMSI, IRNUMMI, IRNUMFI, IRSIVOL, IRMIVOL, IRFIVOL, IRELOSS, IRRSPLH, C1H, C2H, C1F, C2F, C3H, C4H, C3F, C4F, PCB, PC, LANDMASK, LAKEMASK, SST, SST_INPUT

Instead of wrfout_d01_2022-01-15_all, put your own wrfoutput and run it. Do you see sp among the variable?

You can also select diagnostic variables (not inside the WRF output) in the the list bellow:

./postwrf.sh -d

List of choosable diagnostic variables 1) avo ( Absolute Vorticity [10-5 s-1]) 2) eth ( Equivalent Potential Temperature [K] ) 3) cape_surface ( 2D fields of CAPE [J/kg] ) 4) cin_surface ( 2D fields of CAPE [J/kg] ) 5) cape_3d ( 2D fields of CAPE [J/kg] ) 6) cin_3d ( 2D fields of CAPE [J/kg] ) 7) ctt ( Cloud Top Temperature [degC] ) 8) dbz ( Reflectivity [dBZ] ) 9) mdbz ( Maximum reflectivity [dBZ] ) 10) geopt ( Full model geopotential [m2 s-2] ) 11) geopotential ( Full model geopotential [m2 s-2] ) 12) helicity ( Storm Relative Helicity [m-2/s-2] ) 13) omg ( Omega [C] ) 14) tk ( Temperature [K] ) 15) p/pres ( Full model pressure [Pa] ) 16) pressure ( Full model pressure [hPa] ) 17) pvo ( Potential Vorticity [PVU] ) 18) pw ( Precipitable Water ) 19) rh2 ( 2m Relative Humidity [%] ) 20) rh ( Relative Humidity [%] ) 21) slp ( Sea Level Pressure [hpa] ) 22) ter ( Model terrain height [m] ) 23) td2 ( 2m dew point temperature [C] ) 24) td ( Dew point temperature [C] ) 25) tc ( Temperature [C] ) 26) th/theta ( Potential temperature [K] ) 27) tv ( Virtual temperature [K] ) 28) twb ( Wet bulb temperature [K] ) 29) updraft_helicity ( Updraft helicity [m-2/s-2] ) 30) ua ( U component of wind on mass points ) 31) va ( V component of wind on mass points ) 32) wa ( W component of wind on mass points ) 33) z/height ( Full Model Height [m] ) 34) lfc ( Level of Free Convection ) 35) lcl ( Lifted Condensation Level ) 36) dust_total ( Total Dust Concentration [Microgram/Kg_dry_air] ) 37) dust_pm10 ( Dust Concentration equal to PM10 [Microgram/Kg_dry_air] ) 38) dust_pm2.5 ( Dust Concentration equal to PM2.5 [Microgram/Kg_dry_air] ) 39) wind_s ( NCL6.5: Wind Speed [m/s] ) 40) wind_d ( NCL6.5: Wind Direction [degrees from north] )

sstikle commented 2 years ago

Sorry by mistake I shared wrong. My variable is tr18_24. Whether it will run or not. Correct information

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

WRF_Extract_ON_OFF = 1 ERA5_Extract_ON_OFF = 0

Variable_names = tr18_24 Location_names = Indore Location_latitudes = 22.7196 Location_longitudes = 75.8577

+Vertical_profile_plot_on_off = 0 +Vprofile_X_axis_decimals = 3

anikfal commented 2 years ago

What is tr18_24 exactly? I want to know if tr18_24 is inside your WRF output?

Run the command below (put your WRF file-name instead of wrfout). What do you see after running the command? ./postwrf.sh -f wrfout | grep -i tr18_24

sstikle commented 2 years ago

run the command, no massage.

anikfal commented 2 years ago

Show the whole command you ran. "No message" is not possible, if the WRF file and name is okay.

sstikle commented 2 years ago

./postwrf.sh -f wrfout_d01_2019-01-18_00:00:00 | grep -i tr18_24. I checked wrfout put bu ncview, variable is their.

anikfal commented 2 years ago

Are you working on Ubuntu?

What do you see by running the command below? ncl_filedump -v tr18_24 wrfout_d01_2019-01-18_00:00:00.nc

sstikle commented 2 years ago

0.3134E-19 0.1302E-20 0.3929E-21 0.3757E-20 0.7027E-20 0.1201E-18 0.6118E-19 0.7452E-21 0.1043E-20 0.8516E-20 0.6330E-19 0.6330E-1

values like this

anikfal commented 2 years ago

Go to the directory of wrfout_d01_2019-01-18_00:00:00:

Run the commands below:

ncl
a = addfile("wrfout_d01_2019-01-18_00:00:00", "r")
printVarSummary(a->tr18_24)

What do you see now?

sstikle commented 2 years ago

all details of variable like Variable: tr18_24 (file variable) Type: float Total Size: 4145568 bytes 1036392 values Number of Dimensions: 4 Dimensions and sizes: [Time | 1] x [bottom_top | 42] x [south_north | 124] x [west_east | 199] Coordinates: Number Of Attributes: 6 etc..........................

sstikle commented 2 years ago

it means wrfout is correct. then why it is showing that No WRF output files in the current directory. You can link or copy one or more files to the current directory.

anikfal commented 2 years ago

If wrfout_d01_2019-01-18_00:00:00 is your WRF file, and if you can display it by ncview, then your file is okay and ./postwrf.sh should recognize your file and work on it.

There is one point:

Dimensions and sizes: [Time | 1] x [bottom_top | 42] x [south_north | 124] x [west_east | 199]

What kind of physical variable is tr18_24? What is the name of it? This variable doesn't have any time series (because of [Time | 1]).

sstikle commented 2 years ago

it is emission variables so we used different code for different 23 pollutants. tr18_24 this is total. how to resolve the issue?

anikfal commented 2 years ago

Can you extract T2 from your wrfoutput?

WRF_Extract_ON_OFF = 1
ERA5_Extract_ON_OFF = 0

Variable_names = T2
Location_names = Indore
Location_latitudes = 22.7196
Location_longitudes = 75.8577
sstikle commented 2 years ago

will try

On Wed, Nov 9, 2022 at 3:37 PM Amirhossein Nikfal @.***> wrote:

Can you extract T2 from you wrfoutput?

WRF_Extract_ON_OFF = 1 ERA5_Extract_ON_OFF = 0

Variable_names = T2 Location_names = Indore Location_latitudes = 22.7196 Location_longitudes = 75.8577

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

-- Dr. Suvarna Tikle Postdoctoral researcher Environmental Modeling Division Max Planck Institute for Meteorology Bundesstr. 53 D-20146 Hamburg @.*** +91 9822082084

sstikle commented 2 years ago

same massage ================= PostWRF (Version 1.1) ===================== Run './postwrf.sh -h' to get help In case of problems, please contact @.***>

 EXTRACT VARIABLE             EXTRACT VARIABLE
 No WRF output files in the current directory.
 You can link or copy one or more files to the current directory.

On Wed, Nov 9, 2022 at 3:38 PM Dr.Suvarna Tikle @.***> wrote:

will try

On Wed, Nov 9, 2022 at 3:37 PM Amirhossein Nikfal < @.***> wrote:

Can you extract T2 from you wrfoutput?

WRF_Extract_ON_OFF = 1 ERA5_Extract_ON_OFF = 0

Variable_names = T2 Location_names = Indore Location_latitudes = 22.7196 Location_longitudes = 75.8577

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

-- Dr. Suvarna Tikle Postdoctoral researcher Environmental Modeling Division Max Planck Institute for Meteorology Bundesstr. 53 D-20146 Hamburg @.*** +91 9822082084

-- Dr. Suvarna Tikle Postdoctoral researcher Environmental Modeling Division Max Planck Institute for Meteorology Bundesstr. 53 D-20146 Hamburg @.*** +91 9822082084

anikfal commented 2 years ago

Can you send your WRF file to me?

sstikle commented 2 years ago

Sorry but I can not because of our Maxplank policy.

anikfal commented 2 years ago

Where is your WRF output files? Copy or link them into the PostWRF directory, and run postwrf.sh again.

So in our meeting we realized this was the problem.

sstikle commented 2 years ago

Yes. Thanks for resolving the issue.

On Thu, Nov 10, 2022 at 10:19 AM Amirhossein Nikfal < @.***> wrote:

Where is your WRF output files? Copy or link them into the PostWRF directory, and run postwrf.sh again.

So in our meeting we realized this was the problem.

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

-- Dr. Suvarna Tikle Postdoctoral researcher Environmental Modeling Division Max Planck Institute for Meteorology Bundesstr. 53 D-20146 Hamburg @.*** +91 9822082084