Innermost (horizontal) gridlength 22 m: divide by zero

Questions and other topics related to UEMS 15.
Post Reply
j1d1w1
Posts: 34
Joined: Wed Mar 13, 2013 2:02 am

Innermost (horizontal) gridlength 22 m: divide by zero

Post by j1d1w1 » Fri Jan 05, 2018 7:03 pm

I'm hoping to use UEMS to diagnose the near-ground wind field over a domain of about a square kilometer with a horizontal resolution of the order of 20 m; this innermost domain (no. 3) is nested within a domain 2 that (in turn, in the usual way) is nested within domain 1, where the outermost domain (domain 1) covers about 20 by 20 km with gridlength 3 x (3 x 20) = 180 m. I'm using the ARW core, and the lat-lon grid. The goal is to simply spin up the fields over 6 hours onto the fine mesh, using an analysis or reanalysis field for initial and boundary conditions. Eventually, i.e. if I can get this to "go", I will provide high resolution (lidar-derived) topography for domain 3.

In my experience dwiz has proven worse than useless to create such a domain, or even to view such a domain that one had created via the command line.

But sidestepping that point, I have used ems_domain.pl to create my domain, adjusted the namelist.wpl file to reflect what I want, then localized successfully. But when I run (ems_run), I immediately get the error: "Illegal division by zero at .../NLphysics.pm line 4412" and that line reads

Code: Select all

my@opts=();
Any suggestions would be appreciated. I'm not trying LES, and have the YSU boundary-layer scheme enabled. No Cu parameterization. Thanks, John

Opsasy
Posts: 3
Joined: Fri Jun 16, 2017 10:21 am

Re: Innermost (horizontal) gridlength 22 m: divide by zero

Post by Opsasy » Sat Jan 06, 2018 7:46 pm

Try to change time step for 60 (seconds). It worked in my case. If I used time step equal to more than 60 (e. g. 120) or Auto_R or Auto_S the I had the same error :)

Opsasy
Posts: 3
Joined: Fri Jun 16, 2017 10:21 am

Re: Innermost (horizontal) gridlength 22 m: divide by zero

Post by Opsasy » Sat Jan 06, 2018 8:45 pm

You should also change radt value for more or equal 1.

j1d1w1
Posts: 34
Joined: Wed Mar 13, 2013 2:02 am

Re: Innermost (horizontal) gridlength 22 m: divide by zero

Post by j1d1w1 » Mon Jan 08, 2018 3:06 am

Thanks -- I'll play with that. I had set TIME_STEP = AUTO_S and UEMS had chosen "a large timestep of 0.5 seconds". Maybe I'll prescribe an even smaller value. John

Post Reply

Who is online

Users browsing this forum: No registered users and 7 guests