* +-======-+ * Copyright (c) 2003-2018 United States Government as represented by * the Admistrator of the National Aeronautics and Space Administration. * All Rights Reserved. * * THIS OPEN SOURCE AGREEMENT ("AGREEMENT") DEFINES THE RIGHTS OF USE, * REPRODUCTION, DISTRIBUTION, MODIFICATION AND REDISTRIBUTION OF CERTAIN * COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS * REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY"). * THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN * INTENDED THIRD-PARTY BENEFICIARY OF ALL SUBSEQUENT DISTRIBUTIONS OR * REDISTRIBUTIONS OF THE SUBJECT SOFTWARE. ANYONE WHO USES, REPRODUCES, * DISTRIBUTES, MODIFIES OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED * HEREIN, OR ANY PART THEREOF, IS, BY THAT ACTION, ACCEPTING IN FULL THE * RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT. * * Government Agency: National Aeronautics and Space Administration * Government Agency Original Software Designation: GSC-15354-1 * Government Agency Original Software Title: GEOS-5 GCM Modeling Software * User Registration Requested. Please Visit http://opensource.gsfc.nasa.gov * Government Agency Point of Contact for Original Software: * Dale Hithon, SRA Assistant, (301) 286-2691 * * +-======-+ function g (args) string = subwrd(args,1) cint = subwrd(args,2) minval = subwrd(args,3) 'run getenv GEOSUTIL' geosutil = result 'getinfo lat' lat = result 'getinfo lon' lon = result 'getinfo time' time = result 'getinfo zpos' zpos = result 'getinfo file' curfile = result 'open 'geosutil'/plots/grads_util/lwmask1440721.tabl' 'getinfo numfiles' newfile = result 'set dfile 'newfile 'set t 1' 'set z 1' 'define mask = lwmask' 'close 'newfile 'set dfile 'curfile 'set t 'time 'set z 'zpos 'set grads off' 'set gxout shaded' 'set ccols 92 0' 'set clevs 0.5' 'd mask' cval = cint if( cint = '' & minval = '' ) 'shades 'string' 0' cval = result endif if( cint != '' & minval = '' ) 'shades 'string' 0 -cint 'cint endif if( cint = '' & minval != '' ) 'shades 'string' 0 -minval 'minval cval = result endif if( cint != '' & minval != '' ) 'shades 'string' 0 -cint 'cint' -minval 'minval endif 'define dummy = regrid2( 'string', .25, .25, bs_p1, 'lon', 'lat' )' 'd maskout(dummy,abs(dummy)-'cval')' say 'd maskout(dummy,abs(dummy)-'cval')' 'cbarn -snum 0.8'