ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

The RIP4 graphics package.

ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

Postby domagojp » Thu Apr 04, 2019 4:38 am

Hi,

I am trying to plot PBLH in feet by modifying ripdp_wrfarw.f to read in PBLH and convert and output in feet:

c
c Get PBL height (PBLH), and convert to feet
c
nf_status = nf_inq_varid (ncid, 'PBLH', varid)
if ( nf_status .ne. 0 ) then ! we have PBLH?
nf_status = nf_inq_varid (ncid, 'PBLH', varid)
endif
call handle_err(777.,nf_status)
iprocvarid(varid)=1
nf_status = nf_get_vara_real (ncid, varid, nf_tstart2,
& nf_tcount2, nf_tarr2)
call handle_err(778.,nf_status)
do j=1,mjx-1
do i=1,miy-1
pblh(i,j)=nf_tarr2(j,i)*3.28
enddo
enddo
plchun='ft'
vardesc='PBL height, ft'
call writefile_rdp(pblh,'pblh ',2,1,vardesc,plchun,
& fname,iendf1,ihrip,rhrip,chrip,iexpanded,
& iexpandedout,ioffexp,joffexp,miy,mjx,mkzh_out)

Ripdp runs fine and creates the *.pblh files but when I try plotting them with rip I get this error:
ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

which indicates problem in arinam function. I tried increasing niam and ncs all the way to 3000000 in driver.f but still get the same error.

Any help appreciated.

Thanks,
Domagoj
domagojp
 
Posts: 8
Joined: Wed Mar 10, 2010 6:30 pm

Re: ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

Postby dcvz » Fri Apr 12, 2019 2:01 am

The first question to ask is can you plot the PBLH field in meters using the standard ripdp? Can you post your plspec (.in) file?

The problem is that your PBLH field is too complicated for aredam to handle. But, there are ways to get around it. First, try making niam even larger. The latest version of rip has niam equal to 10 million and ncs equal to 1 million. You can keep making these larger until you run out of memory or it takes forever to finish. A better option is to make the field simpler to plot. 1) Increase your cint. Maybe try cint=1000. If that works, then reduce to a better value. 2) Run the PBLH through one of rip's smoothers. e.g. smth=104. 3) subset the domain using xwin, ywin.

As a last resort there is a method described in arex02.f that you can try. (Run 'ncargex -areas' to get the source). But, that would require non-trivial source code mods to rip.
dcvz
 
Posts: 183
Joined: Tue Apr 15, 2008 12:02 am

Re: ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

Postby domagojp » Mon May 06, 2019 10:50 am

Playing with niam did not help but increasing cint did the trick. Thanks!
domagojp
 
Posts: 8
Joined: Wed Mar 10, 2010 6:30 pm


Return to RIP4

Who is online

Users browsing this forum: No registered users and 2 guests