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

Southwest pacific no map_proj works Namelist.wps

tolux

New member
I am trying to run a wrf model for TC in the southwest pacific in this geolocation (-5 to -55 lon) and (145 E to 170 W lon). Because this longitude goes through the 180, it is given problem to define, when try to define it the longitude goes reverse and run from 170 W to 140 E making my interest area empty. Please someone help with suitable map projection and how to define the ref longitude and others to capture my interested area.

to add to this,
when I set the domain using mercator and extend the longitude to cross 180, geogrid.exe run and when I try to visualize the geo_*.nc using python. for example If I just contour(landmask) I will see the domain I specify, however if I contourf(Xlat,xlong,landmask) it shows the other part and leave area of interest empty. Then when I try to run metgrid the error comes up as Pres at 200100 is missing (i=682 j=1) something like that.

Thank you

Tolux
 
Last edited:
when I tried to extend the longitude beyond the 180 degree, this happened. I have try using lat-lon map_proj as well. got similar problem
 

Attachments

  • 1727665909655.png
    1727665909655.png
    43.9 KB · Views: 5
  • 1727665992548.png
    1727665992548.png
    325.4 KB · Views: 3
  • 1727666043515.png
    1727666043515.png
    303.9 KB · Views: 4
for me this problem looks complicated, the maps will reverse if the set number of longitude above 165, it because it will cross international date line and will back to west,

actually i never try this, but you can check this site Pacific Northwest Mesoscale Model Numerical Forecast Information , this is for northwest pacific but i hope u can get any clue about lon lat setting from this site, the namelist.wps detail attach on this site to, you can check it, i hope this site can give you some clue to solve your problem
 
can i know the dataset that you use for meteorological data while run ungrib? like gdas or gfs or era5 or other data? also your before do metgrid, can i know your ungrib process? like the output of ungrib (FILE (domain) (date) ), because looks like the model cant interpolate that data due to missing value on some level
 
@tolux

As a test, I set up my namelist.wps to match yours, and there are no issues with the domain. I used the util/plotgrids_new.ncl script to plot the domain, and it looks as expected (see the attachment). I also ran all wps executables and didn't run into any errors. The metgrid error you're seeing is likely more related to the input data you've obtained to use for your case. If you are only grabbing a subset of data, often, due to the map projections, some data can be missing - especially along the edges. It is recommended to obtain the full global dataset to avoid such issues.
 

Attachments

  • wps_show_dom.pdf
    38.8 KB · Views: 1
Top