20020619: 20020618: nested gfunc

Hmmm... 73x144 doesn't sound very big at all. I was thinking more like
greater than 100,000 grid points.

Strange,

Steve

>From: Chris Hennon <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200206191555.g5JFt2614031

>Steve -
>
>First I tried to sub vecr(urel,vrel) in for wnd and got the same problem.
>Since I have quite a few scripts that are having this issue I decided to
>up the memory in GEMPRM.PRM (I doubled it) and recompile.  Now all the
>scripts seem to be performing normally.  I am working with
>global reanalysis data so I did have quite a big grid (73 x 144).  Thanks
>for your help.
>
>Chris
>
>
>On Tue, 18 Jun 2002, Unidata Support wrote:
>
>> 
>> Chris,
>> 
>> Sounds like you are having trouble composing WND from the components).
>> Not sure why without seeing what your grid navigation and size is since 
>> the function you show works for me using current ETA data.
>> 
>> Could be a problem with the stack size. Are these big grids (kxky)?
>> If so, it could be that you don't have enough memory to load in
>> all the grids (remember that you probably have a conversion from
>> dwpc=sub(dwpk,273.15) and vecr(urel,vrel), before the mixr() computation
>> and the resulting output grid.
>> 
>> Presumably you have UREL and VREL at 925mb. What happens if you
>> try using sdiv(mixr(dwpc,pres),vecr(urel,vrel))? 
>> 
>> If the result works with smaller grids, then you may have to
>> do the computation in smaller pieces using gddiag.
>> 
>> Similarly, your tmpc calculation is also undergoing a tmpc=tmpk-273.15.
>> 
>> Otherwise, we'd have to recompile with a bigger heap defined for
>> LLMDGG in $GEMPAK/include/GEMPRM.PRM
>> 
>> 
>> Steve Chiswell
>> 
>> 
>> >From: Chris Hennon <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200206182105.g5IL5L617501
>> 
>> >Steve -
>> >
>> >Can you tell me why these arguments:
>> >
>> >gddiag << EOF
>> > GDFILE   = /home/chennon/data/grib/mdivrepo_1998.grd
>> > GDOUTF   = /home/chennon/data/grib/mdiv_1998.grd
>> > GFUNC    = sdiv(mixr(dwpc,pres),wnd) 
>> > GDATTIM  = 980501/0000 
>> > GLEVEL   = 925
>> > GVCORD   = PRES
>> > GRDNAM   = MDIV 
>> > GPACK    =  
>> >r
>> >
>> >e
>> >EOF
>> >
>> >produce this error:
>> >
>> >[DG -7]  Input grid WND ^980501/0000 @925 %PRES in SDIV cannot be found.
>> >
>> >I also have this problem when I try to add three scalar fields, for
>> >example:
>> >
>> >add(add(tmpc+1,tmpc+2),tmpc+3)
>> >
>> >All grids that I need to compute moisture convergence exist in the
>> >mdivrepo grid file.  It isn't liking the nesting for some reason - I have
>> >used this technique successfully in the past.  I'm using version
>> >5.6.e.1.  Thanks.
>> >
>> >Chris
>> >
>> >================================================
>> >| Chris Hennon             Ohio State University   |
>> >| Tropical Meteorology      address@hidden   |
>> >|                                              |
>> >| Dept of Geography   Office: 1155 Derby Hall  |
>> >| 1036 Derby Hall     Phone : (614) 292-2704   |
>> >| Columbus, OH 43210  Fax   : (614) 292-6213   |
>> >================================================
>> >
>> 
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8643                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>> 
>