[gradsusr] displaying data for negative longitudes

Adams, Jennifer M. (GSFC-610.2)[ADNET SYSTEMS INC] jennifer.m.adams at nasa.gov
Wed Jun 5 13:42:48 EDT 2019


Hi, Alan -- 
Defined variables never get the 'wrap-the-globe' flag set to 'true', so you have to extend the domain before you define your variables, then restrict it when you are ready to display. 

set lon -180 360   
myvar = myexpr
set lon -170 190
display myvar

--Jennifer


-- 
Jennifer Miletta Adams
ADNET Systems, Inc.
NASA/GSFC, Code 610.2
Building 32, Room S159
(301) 614-6070
 

On 6/5/19, 1:12 PM, "gradsusr on behalf of Alan Robock" <gradsusr-bounces at gradsusr.org on behalf of robock at envsci.rutgers.edu> wrote:

    Hi,
    
    When I
    
    set lon -170 190
    
    I have no problem drawing global maps of data without separating the 
    continents.  But when I ran a script to calculate averages from the 
    data, the resulting data will only display at longitudes 0 to 360, and 
    the data for other longitudes, if I set them, are blank.  When I did the 
    calculation without running a script, the resulting averages can be 
    plotted at any longitude.
    
    So, is there a setting I can use to plot data at any longitude?
    
    Or is there a setting I need to tell the script to produce output that 
    can be displayed at any longitude?
    
    Here are some details.  Here is the script:
    
    i = 1
    'fiymn = 0'
    while ( i <= 42 )
      'fiymne'i' = ( ( ave(ta.'i', time = 1dec1986,time=1feb1987)+ 
    ave(ta.'i', time = 1dec1987,time=1feb1988)+ ave(ta.'i', time = 
    1dec1988,time=1feb1989)+ ave(ta.'i', time = 1dec1989,time=1feb1990)+ 
    ave(ta.'i', time = 1dec1990,time=1feb1991 ))/5)'
    'fiymn = fiymn + fiymne'i
    i = i + 1
    endwhile
    'fiymn = fiymn/42'
    
    When I plot
    
    d fiymne2 (or any number)
    
    or
    
    d fiymn
    
    or
    
    d fiymne2-fiymn
    
    the data only show up between lon 0 and lon 360.
    
    If I do, for example,
    
    d ta.2
    
    I can plot it for any lon.
    
    Here is an example ddf
    
    DSET 
    /glade/collections/cdg/data/cesmLE/CESM-CAM5-BGC-LE/atm/proc/tseries/monthly/TREFHT/b.e11.B20TRC5CNBDRD.f09_g16.002.cam.h0.TREFHT.192001-200512.nc
    DTYPE netcdf
    TITLE historical ensemble member 2 for LE
    UNDEF -99999e8
    OPTIONS 365_day_calendar
    XDEF 288 linear 0.0 1.25
    YDEF 192 linear -90 0.9375
    ZDEF 30 levels 3.64346569404006, 7.59481964632869, 14.3566322512925,  
    24.6122200042009, 38.2682997733355, 54.5954797416925, 72.0124505460262, 
    87.8212302923203, 103.317126631737, 121.547240763903, 142.994038760662, 
    168.225079774857, 197.908086702228, 232.828618958592, 273.910816758871, 
    322.241902351379, 379.100903868675,
      445.992574095726, 524.687174707651, 609.778694808483, 
    691.389430314302, 763.404481112957, 820.858368650079, 859.53476652503, 
    887.020248919725, 912.644546944648, 936.198398470879, 957.485479535535, 
    976.325407391414, 992.556095123291
    TDEF 1032 linear 1jan1920 1mo
    VARS 1
    TREFHT=>TA 0 t,y,x K Surface air temp
    ENDVARS
    
    Thanks.
    
    -- 
    Alan
    _________________________________________________________________________
    Alan Robock, Distinguished Professor
       Associate Editor, Reviews of Geophysics
    Department of Environmental Sciences             Phone: +1-848-932-5751
    Rutgers University                    E-mail: robock at envsci.rutgers.edu
    14 College Farm Road            https://urldefense.proofpoint.com/v2/url?u=http-3A__people.envsci.rutgers.edu_robock&d=DwIGaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=M2mFZsxo0ZoztoK2Lyoypu1kKkBVKyEoM9CXGWO42zk&m=01pdYNWJh4oo5Zb85hF88gGxQJnQ_ESN3Kh-d-XvhqU&s=kKXbgYaTNZoCyUIAAZ4AX7LxhNSMnf-Ln2Y1rrsU_xg&e=
    New Brunswick, NJ 08901-8551  USA      ☮  https://urldefense.proofpoint.com/v2/url?u=http-3A__twitter.com_AlanRobock&d=DwIGaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=M2mFZsxo0ZoztoK2Lyoypu1kKkBVKyEoM9CXGWO42zk&m=01pdYNWJh4oo5Zb85hF88gGxQJnQ_ESN3Kh-d-XvhqU&s=1sialPClpN38BVR5il5xQF7pg3xbiam6c_5AtjL9fuc&e=
    
    
    _______________________________________________
    gradsusr mailing list
    gradsusr at gradsusr.org
    https://urldefense.proofpoint.com/v2/url?u=http-3A__gradsusr.org_mailman_listinfo_gradsusr&d=DwIGaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=M2mFZsxo0ZoztoK2Lyoypu1kKkBVKyEoM9CXGWO42zk&m=01pdYNWJh4oo5Zb85hF88gGxQJnQ_ESN3Kh-d-XvhqU&s=2RWDRqKOnwVcM68SwkL2P1Sa-dlfeQJ4KLhNT8ZXkFE&e=
    




More information about the gradsusr mailing list