! +-======-+ ! Copyright (c) 2003-2007 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 cloud (args) expid = subwrd(args,1) season = subwrd(args,2) output = subwrd(args,3) expdsc = subwrd(args,4) 'set vpage off' 'set parea off' 'set grid off' 'set mproj scaled' 'set frame on' 'set xlopts 1 3 .11' 'set ylopts 1 3 .11' 'rgbset' 'set dfile 1' 'set t 1' 'define mod = cldlo'season'mod' 'getinfo lon' lon = result 'set dfile 2' 'define obs = cldlow'season'isccp*0.01' 'set vpage 0 8.5 0.0 11' 'set parea 1.5 7.0 7.70 10.50' 'set grads off' 'set clevs .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6 .65 .7 .75 .8 .85 .9 .95 1.0 290 300' 'set clevs .03 .07 .10 .13 .17 .2 .23 .27 .3 .33 .37 .4 .43 .47 .5 .53 .57 .6 .63 .67 .7 .73' 'set ccols 59 58 57 56 55 49 47 45 44 37 36 34 33 32 31 21 22 23 24 25 26 27 28' 'd mod' 'set parea 0 8.5 7.0 11' 'cbarn -vert' 'set parea off' 'set vpage 0 8.5 0.0 11' 'set parea 1.5 7.0 4.30 7.10' 'set grads off' 'set clevs .03 .07 .10 .13 .17 .2 .23 .27 .3 .33 .37 .4 .43 .47 .5 .53 .57 .6 .63 .67 .7 .73' 'set ccols 59 58 57 56 55 49 47 45 44 37 36 34 33 32 31 21 22 23 24 25 26 27 28' 'd obs' 'set parea off' 'set vpage 0 8.5 0.0 11' 'set parea 1.5 7.0 0.90 3.70' 'set grads off' 'define obsg = regrid2( obs,1,1,bs_p1,'lon',-90)' 'define modg = regrid2( mod,1,1,bs_p1,'lon',-90)' 'set clevs -.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9' 'set clevs -.45 -.4 -.35 -.3 -.25 -.2 -.15 -.1 -.05 .05 .1 .15 .2 .25 .3 .35 .4 .45' 'set ccols 55 49 47 45 44 36 34 33 32 0 21 22 23 24 25 26 27 28 69' 'd maskout( modg-obsg,abs(obsg) )' 'cbarn' 'stats maskout(modg,abs(obsg))' avgmod = subwrd(result,1) stdmod = subwrd(result,2) 'stats maskout(obsg,abs(obsg))' avgobs = subwrd(result,1) stdobs = subwrd(result,2) 'stats maskout(modg-obsg,abs(obsg))' avgdif = subwrd(result,1) stddif = subwrd(result,2) 'set vpage off' 'set string 1 c 6' 'set strsiz .11' 'draw string 4.25 10.9 'expdsc 'draw string 4.25 10.6 Low-Level Clouds (fraction) 'season 'draw string 4.25 7.2 Observations (ISCCP) 'season 'draw string 4.25 3.80 Model - Observations' 'set string 1 c 4' 'set strsiz .08' 'draw string 0.52 10 Mean: 'avgmod 'draw string 0.52 9.8 Std: 'stdmod 'draw string 0.52 6.6 Mean: 'avgobs 'draw string 0.52 6.4 Std: 'stdobs 'draw string 0.52 3.2 Mean: 'avgdif 'draw string 0.52 3.0 Std: 'stddif 'myprint -name 'output'/cldlo.'season 'set mproj latlon'