[gradsusr] displaying data for negative longitudes

Alan Robock robock at envsci.rutgers.edu
Wed Jun 5 12:57:53 EDT 2019


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            http://people.envsci.rutgers.edu/robock
New Brunswick, NJ 08901-8551  USA      ☮  http://twitter.com/AlanRobock




More information about the gradsusr mailing list