Page 1 of 1

ERROR 7 IN AREDAM - AREA-MAP ARRAY OVERFLOW

PostPosted: Thu Apr 04, 2019 4:38 am
by domagojp
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

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

PostPosted: Fri Apr 12, 2019 2:01 am
by dcvz
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.

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

PostPosted: Mon May 06, 2019 10:50 am
by domagojp
Playing with niam did not help but increasing cint did the trick. Thanks!