Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

how to use different initial soil moisture as the input

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

Iamwrfuser

New member
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
 
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?
 
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

  • 1-SMOIS_GFS6hr.pdf
    33.5 KB · Views: 56
  • 2-SMOIS_GFS6hr_SMOPS24hr.pdf
    21.7 KB · Views: 112
  • 3-SMOIS_GFS24hr.png
    3-SMOIS_GFS24hr.png
    153.6 KB · Views: 4,019
  • 4-SMOIS_GFS24hr_SMOPS24hr.png
    4-SMOIS_GFS24hr_SMOPS24hr.png
    106 KB · Views: 4,014
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
 
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.
 
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

  • namelist.wps
    1.3 KB · Views: 56
  • namelist.input
    4.4 KB · Views: 53
  • GFSp25.tar
    87.5 MB · Views: 56
  • SMOPS.tar
    38.1 MB · Views: 111
  • gribtab_SMOPS3.0.dat
    3.3 KB · Views: 60
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.
 
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
    n1_SM000010_GFS6hr_SMOPS6hr.png
    27.9 KB · Views: 3,980
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.
 
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

  • namelist.wps
    1.3 KB · Views: 976
  • ungrib.log
    408.8 KB · Views: 57
Hi Ehsan,
In the namelist you sent, in the &metgrid section, I only see that you are using GFS data:
Code:
&metgrid
 fg_name = 'GFS',

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

How are you actually using both datasets?
 
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

  • namelist.wps
    1.3 KB · Views: 58
  • ncview.SM000010.png
    ncview.SM000010.png
    23.8 KB · Views: 3,954
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

  • gfs_and_smops.png
    gfs_and_smops.png
    373.7 KB · Views: 3,946
  • only_gfs.png
    only_gfs.png
    269.1 KB · Views: 3,946
  • namelist.wps
    754 bytes · Views: 21
Top