[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20010925: 20010924: Gempak Omega plotting problem



Hi Mike.
I have to work from home for the time being....so the 56K dialup is
really a pinch when trying to load your annimation applet.

Anyhow, let me know how the images are tomorrow. I can always
launch the browser and go sift the deluge of resumes were
getting for job openings. Thanks for resending your original
question. I lost sight of it in my work box with all the
status reports, committee meeting, hiring panel, and other operating
system holy war messages I have pending at the moment.

Steve




>From: Mike Voss <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200109251754.f8PHs3120565

>Steve,
>
>The irony kills me, there has been data to fill the plots every day since I do
> n't know when, and of course the day you look at it nothing is there....arrrg
> h!
>
>I will try to set the SCALE and see if that fixes the problem. If not I'll sen
> d an email in the AM (when I know the plots are there) so you can check out w
> hat we're talking about. Much thanks,
>Mike
>
>At 06:09 AM 9/25/2001 -0600, you wrote:
>
>>Mike,
>>
>>I couldn't see much in the urls. Most frames do not have any color fill.
>>The CONUS had none.
>>
>>One thing I notice is that you are specifying all your FINT & FLINE
>>values. But, you are using a SCALE=999.
>>
>>When scale is 9999, it allows GEMPAK to choose the appropriate scaling based 
> on the
>>values in the data. What you might be seeing is the scale value changing betw
> een frames.
>>
>>If you are doing a loop and expecting certain FINT values, then you should se
> t SCALE appropriately
>>so that you have the same value throughout you loops.
>>
>>Steve Chiswell
>>Unidata User Support
>>
>>
>>
>>>From: Mike Voss <address@hidden>
>>>Organization: UCAR/Unidata
>>>Keywords: 200109242314.f8ONED117319
>>
>>>Hello,
>>>I sent this message last Monday (9/17), but maybe it got lost in the shuffle
>  o
>>> r never made it. Anyway, we're still having the same problem so I'll try ag
> ai
>>> in. Thanks for any help you can provide!!
>>>
>>>--------
>>>Hi Steve,
>>>We add an omega grid to the MRF201 data daily with the gdomeg program.  We d
> is
>>> play this grid in a polar hemispheric output:
>>>http://www.met.sjsu.edu/weather/models/mrfp/omeg700-ani.html
>>>and a USA (CONUS) output: 
>>>http://www.met.sjsu.edu/weather/models/conusmrf/omeg700-ani.html
>>>
>>>The only differences in the .nts files which create these images are with GA
> RE
>>> A and PROJ parameters (see below).  The polar output plots fine. The CONUS 
> ou
>>> tput files, however, have a strange color scheme that develops. This does n
> ot
>>>  happen consistently as you can see by looking at today's output, but does 
> te
>>> nd to occur in the later forecast periods, i.e., days 7 through 10. In comp
> ar
>>> ing the same forecast hour image in
>>>the two different outputs we have found that the max, min, and zero values
>>>are in the same place geographically in both images.  The difference is
>>>that the CONUS output (when having trouble) seems to use all available color
> s 
>>> to contour from zero to the max or min value. In doing this the incremental
>  s
>>> pacing is too small and therefore shows false values when compared to the c
> ol
>>> orbar.  Do you have some insight into what may be causing our problem?  Bel
> ow
>>>  I have a sample .nts file from each script and a copy of the script itself
>  u
>>> sed to generate the images:
>>>
>>>.nts from the polar plot:
>>>GDFILE  $METDAT/grids/2001010900_mrf201.gem                                 
>   
>>>                                                           
>>>GDATTIM first                                                               
>   
>>>                                                           
>>>GLEVEL  700                                                                 
>   
>>>                                                           
>>>GVCORD  PRES                                                                
>   
>>>                                                           
>>>PANEL   0                                                                   
>   
>>>                                                           
>>>SKIP    0/3;3                                                               
>   
>>>                                                           
>>>SCALE   999                                                                 
>   
>>>                                                           
>>>GDPFUN  mul(-1,omeg)!quo(hght,10)   !kntv(wnd)                              
>   
>>>                                                           
>>>TYPE    f!c                         !b                                      
>   
>>>                                                           
>>>CONTUR  0                                                                   
>   
>>>                                                           
>>>CINT    3                                                                   
>   
>>>                                                           
>>>LINE    1/1/1                                                               
>   
>>>                                                           
>>>FINT    -6;-5;-4;-3;-2;-1;0;1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19 
>   
>>>                                                           
>>>FLINE   21;22;23;24;25;26;27;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20   
>   
>>>                                                           
>>>HILO                                                                        
>   
>>>                                                           
>>>HLSYM   2;.5//21;1//hw                                                      
>   
>>>                                                           
>>>CLRBAR  1/v/ll/.01;.25/.6;.03                                               
>   
>>>                                                           
>>>WIND    bk1/.4                                                              
>   
>>>                                                           
>>>REFVEC                                                                      
>   
>>>                                                           
>>>TITLE   1/-3/@ Omega(-ub/s) & Hght(dm)!1/-2/ Valid ~                        
>   
>>>                                                           
>>>TEXT    1                                                                   
>   
>>>                                                           
>>>CLEAR   y                                                                   
>   
>>>                                                           
>>>GAREA   10;-90;10;-90                                                       
>   
>>>                                                           
>>>PROJ    str/90;-90;0                                                        
>   
>>>                                                           
>>>MAP     32                                                                  
>   
>>>                                                           
>>>LATLON  32/1/1/30/30;90                                                     
>   
>>>                                                           
>>>------------------------------------
>>>
>>>And .nts for the conus plot:
>>>GDFILE  $METDAT/grids/2001010900_mrf201.gem                                 
>   
>>>                                                           
>>>GDATTIM first                                                               
>   
>>>                                                           
>>>GLEVEL  700                                                                 
>   
>>>                                                           
>>>GVCORD  PRES                                                                
>   
>>>                                                           
>>>PANEL   0                                                                   
>   
>>>                                                           
>>>SKIP    0/3;3                                                               
>   
>>>                                                           
>>>SCALE   999                                                                 
>   
>>>                                                           
>>>GDPFUN  mul(-1,omeg)!quo(hght,10)   !kntv(wnd)                              
>   
>>>                                                           
>>>TYPE    f!c                         !b                                      
>   
>>>                                                           
>>>CONTUR  0                                                                   
>   
>>>                                                           
>>>CINT    3                                                                   
>   
>>>                                                           
>>>LINE    1/1/2                                                               
>   
>>>                                                           
>>>FINT    -6;-5;-4;-3;-2;-1;0;1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19 
>   
>>>                                                           
>>>FLINE   21;22;23;24;25;26;27;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20   
>   
>>>                                                           
>>>HILO                                                                        
>   
>>>                                                           
>>>HLSYM   2;.5//21;1//hw                                                      
>   
>>>                                                           
>>>CLRBAR  1/v/ll/.01;.25/.6;.03                                               
>   
>>>                                                           
>>>WIND    bk1/.4                                                              
>   
>>>                                                           
>>>REFVEC                                                                      
>   
>>>                                                           
>>>TITLE   1/-3/@ Omega(-ub/s) & Hght(dm)!1/-2/ Valid ~                        
>   
>>>                                                           
>>>TEXT    1                                                                   
>   
>>>                                                           
>>>CLEAR   y                                                                   
>   
>>>                                                           
>>>GAREA   10;-140;55;-55                                                      
>   
>>>                                                           
>>>PROJ    lcc//6;2;2;1                                                        
>   
>>>                                                  
>>>MAP     32                                                                  
>   
>>>                                                           
>>>LATLON  32/1/1/30/30;90                                                     
>   
>>>                                                           
>>>---------------------------------------------------
>>>
>>>And the script (this is identical for both plots except for the directory in
> fo
>>> ):
>>>
>>>#!/bin/csh -f
>>>#This script uses GEMPAK to make gif files of virt. velocity and wind @700mb
> s.
>>>
>>> setenv DISPLAY rossby:1
>>> source /usr/local/gempak/gempak5.6/Gemenviron
>>> set CURRENT_MRF=`ls -t $METDAT/grids/*_mrf201.gem | head -1`
>>> cd /export/home/baggett/gempak/mrfp
>>>
>>>#if (-e omeg700f00.gif ) then
>>>     rm -f omeg700*.gif*
>>>#endif
>>>
>>>rm gemglb.nts
>>>set a=00 b=06 c=12 d=18 e=24 f=30 g=36 h=42 i=48 j=54 k=60      
>>>set l=72 m=84 n=96 o=108 p=120 q=132 r=144 s=156 t=168 u=180
>>>set v=192 w=204 x=216 y=228 z=240
>>>foreach time($a $e $i $l $n $p $r $t $v $x $z)
>>>#set a=00 b=06
>>>#foreach time($a $b)
>>>set thing = f$time
>>>echo $thing
>>>    gdplot2<< EOF
>>> restore omeg700.nts
>>>    gdfile=$CURRENT_MRF
>>>    gdattim=$thing
>>>    device=gf|omeg700$thing.gif|800;540
>>>    clear=y
>>>    text=1/2/1/hw
>>>    title=1/-3/@ Omega(-ub/s) & Hght(dm)!1/-2/ Valid ~
>>>    r
>>>
>>>    e
>>>
>>>EOF
>>>
>>>rm gemglb.nts
>>>
>>>gpend
>>>
>>>end
>>>-------------------------------------------------
>>>
>>>thanks for your help!!
>>>
>>>Mike
>>>
>>>
>>>
>>>--------------------------
>>>Mike Voss                                
>>>Department of Meteorology               
>>>San Jose State University                        
>>>One Washington Square                      
>>>San Jose, CA 95192-0104   
>>>             
>>>408.924.5204 voice
>>>408.924.5191 fax   
>>>
>>
>>**************************************************************************** 
>>Unidata User Support                                    UCAR Unidata Program 
>>(303)497-8644                                                  P.O. Box 3000 
>>address@hidden                                   Boulder, CO 80307 
>>---------------------------------------------------------------------------- 
>>Unidata WWW Service                        http://www.unidata.ucar.edu/      
>>**************************************************************************** 
> < 
>
>--------------------------
>Mike Voss                                
>Department of Meteorology               
>San Jose State University                        
>One Washington Square                      
>San Jose, CA 95192-0104   
>             
>408.924.5204 voice
>408.924.5191 fax   
>



pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy