[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