how to use different initial soil moisture as the input

Questions/topics related to the first-guess input meteorological data, and Vtables
Post Reply
Iamwrfuser
Posts: 1
Joined: Mon Jul 01, 2019 8:20 pm

how to use different initial soil moisture as the input

Post by Iamwrfuser » Mon Jul 01, 2019 8:33 pm

Hi All,
I want to change the initial condition with only with different soil moisture, but I have no clue how to do it?
Could you give me a guide?
Thank you!
Mandy

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Tue Jul 02, 2019 10:48 pm

Hi Mandy,
Do you have some soil moisture data that you want to use, in addition to your other first-guess data, or are you wanting to manually modify data points to new values for soil moisture?
NCAR/MMM

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Mon Sep 23, 2019 3:37 pm

Hi,
I had some problems like that. May I ask whether you solved yours?
I want to use SMOPS data (including just soil moisture) as initializing. I wrote its Vtable and run (ungribe.exe and metgrid.exe) correctly (kwerner helped me about it: posted as SMOPS Vtable).
After running ./real.exe I found not suitable output for SMOIS in wrfinput files.
Some points I should mention: GFS data linked per 6hr (00,06,12,18) but SMOPS data linked daily (only 00).
So, for GFS and ./ungrib.exe (interval_seconds = 21600 in namelist.wps), but for SMOPS and ./ungrib.exe (interval_seconds = 86400 in namelist.wps). While ./real.exe (interval_seconds = 21600, in namelist.input).
However, I've test for GFS daily linked (only 00 hr and interval_seconds = 86400), but there is still bad difference between with and without SMOPS initialization.
I've attached 4 files for:
1- Only GFS 6rh (00,06,12,18)
2- GFS 6 hr (00,06,12,18) and SMOPS 24 hr (00)
3- Only GFS 24 hr (00)
4- GFS 24 hr (00) and SMOPS 24 hr (00)

I hope I could describe my issue clearly.
I'll be thankful for any help.

Sincerely
Ehsan
Attachments
4-SMOIS_GFS24hr_SMOPS24hr.png
3-SMOIS_GFS24hr.png
2-SMOIS_GFS6hr_SMOPS24hr.pdf
(21.69 KiB) Downloaded 14 times
1-SMOIS_GFS6hr.pdf
(33.54 KiB) Downloaded 16 times

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Tue Sep 24, 2019 10:30 am

About my previous post, I should mention SMOPS data are daily. So I think using daily data as input for 00 hr, may resulted these bad SMOIS values, may be not. However, SM in met_em* files doesn't show very significant differences.
Also, GFS first layer soil moisture data are for 0-10 cm and SMOPS soil moisture data are for 0-5cm!
Any solution will be appreciated.

Sincerely
Ehsan

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Thu Sep 26, 2019 8:37 pm

Hi Ehsan,
Are there missing points where you need GFS data to be used for soil moisture?
Can you attach the namelist.wps file that you are using, and if possible, can you attach 2 (consecutive times) of your SMOPS data files? Thanks.
NCAR/MMM

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Fri Sep 27, 2019 11:06 am

Dear kwerner,
Thank you so much for your attention.
1- I don't think there is missing point.
2- I've attached namelist.wps, namelist.input, Vtable.SMOPS, guidance of SMOPS data (gribtab_SMOPS3.0.dat), GFSp25 and SMOPS daily data.

I did following instruction:

$ ln -s ungrib/Variable_Tables/Vtable.GFS Vtable
$ ./link_grib.csh /home/x_292026/98_6_15/data_GFSp25/gfs.0p25.201704*
$ vi namelist.wps
interval_seconds = 21600,
prefix = 'GFS'
$ ./ungrib.exe

$ ./link_grib.csh /home/x_292026/98_6_15/data_SMOPS/NPR_SMOPS_CMAP_D2017041*.gr2
$ ln -sf ungrib/Variable_Tables/Vtable.SMOPS Vtable
$ vi namelist.wps
interval_seconds = 86400,
prefix = 'SMOPS'
fg_name = 'GFS', 'SMOPS'
$ ./ungrib.exe

$ vi namelist.wps
interval_seconds = 21600,
$ vi METGRID.TBL.ARW
name=SM000010
change missing_value=-1.E30
to
missing_value=0.
$ ./metgrid.exe

$ cd ../WRF/run/
$ ln -sf ../../WPS/met_em.d0* .
$ ./real.exe
$ ./wrf.exe

Please let me know if I should prepare anything else.
I'll be thankful if I have your idea.

Sincerely
Ehsan
Attachments
gribtab_SMOPS3.0.dat
(3.3 KiB) Downloaded 16 times
SMOPS.tar
(38.06 MiB) Downloaded 15 times
GFSp25.tar
(87.49 MiB) Downloaded 16 times
namelist.input
(4.37 KiB) Downloaded 16 times
namelist.wps
(1.29 KiB) Downloaded 14 times

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Sat Sep 28, 2019 11:45 am

I had forgotten to upload Vtable.SMOPS. I've attached it in this post.
Attachments
Vtable.SMOPS.tar
(10 KiB) Downloaded 15 times

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Wed Oct 02, 2019 5:53 pm

Hi,
Even though your SMOPS data aren't available as often as the GFS data, in your namelist.wps, you should still be able to set interval_seconds = 21600 when ungribbing your SMOPS data, and the data will be interpolated to 6 hour intervals. This should give you a GFS and SMOPS file every 6 hours, and then when you run metgrid, it will use the interpolated SMOPS data at each time, instead of filling in GFS data during the gaps, which is very likely causing your problem. Give that a try and see if it makes any difference.
NCAR/MMM

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Fri Oct 04, 2019 3:18 pm

Dear kwerner,
Thank you for your nice reply.
I tried with interval_seconds = 21600. It seems it could reduce the problem of SMIOS in wrfinput files. However, about SM000010 in met_em files, for 06, 12 and 18 hr, which there isn't SMOPS data, it doesn't show a good result. Attached file is SM000010 for 12hr.

Sincerely
Ehsan
Attachments
n1_SM000010_GFS6hr_SMOPS6hr.png

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Wed Oct 09, 2019 4:12 pm

Hi,
Can you attach the namelist.wps file you are using to run metgrid, along with the ungrib.log file for your GFS run? Thanks.
NCAR/MMM

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Sun Oct 13, 2019 7:46 pm

Dear kwerner,
Pardon me to check your reply with a delay.
I've attached namelist.wps and ungrib.log with for run GFS. I hope these files be correct. However, the date in these files may differ with previous files.
Please let me know if I should prepare any other files.

Sincerely
Ehsan
Attachments
ungrib.log
(408.8 KiB) Downloaded 13 times
namelist.wps
(1.29 KiB) Downloaded 14 times

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Fri Oct 18, 2019 12:39 pm

Hi Ehsan,
In the namelist you sent, in the &metgrid section, I only see that you are using GFS data:

Code: Select all

&metgrid
 fg_name = 'GFS',
If you are using both GFS and SMOPS data, then that line should read something like:

Code: Select all

&metgrid
 fg_name = 'GFS','SMOPS'
How are you actually using both datasets?
NCAR/MMM

ehsantaghizadeh
Posts: 14
Joined: Sun Oct 21, 2018 8:42 am

Re: how to use different initial soil moisture as the input

Post by ehsantaghizadeh » Thu Nov 07, 2019 4:31 pm

Dear kwerner,
Pardon me to back too late. Actually, I left that job for a few weeks. I was in a course.
You are right. I had uploaded wrong namelist.wps. I've attached the correct one here.
First, I've used
&ungrib
prefix = 'GFS',
and run ungrib.exe
then I've used:
&ungrib
prefix = 'SMOPS',

&metgrid
fg_name = 'GFS', 'SMOPS'

As, you said before, I've used
interval_seconds = 21600
for both "GFS" and "SMOPS"
However, I've attached png file for SM000010 in met_em.d02.2017-04-13_00:00:00.nc
May I ask if you test it. Did you get right picture?

Sincerely
Ehsan
Attachments
ncview.SM000010.png
namelist.wps
(1.3 KiB) Downloaded 15 times

kwerner
Posts: 1859
Joined: Wed Feb 14, 2018 9:21 pm

Re: how to use different initial soil moisture as the input

Post by kwerner » Fri Nov 15, 2019 10:03 pm

Hi,
I apologize for the delay. I did run it using the SMOPS data you previously sent, along with GFS data. I didn't get the same image you are seeing. I'm attaching screenshots of my images (one with only GFS, and the other is with GFS+SMOPS). I'm also attaching the namelist.wps file I used.
Attachments
namelist.wps
(754 Bytes) Downloaded 6 times
only_gfs.png
gfs_and_smops.png
NCAR/MMM

Post Reply

Return to “Meteorological Input Data”