<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi! <div class=""><br class=""></div><div class="">Welcome to NCL! This task might be quite a steep learning curve, but it is possible. I have not worked with your exact kind of data, but others might chime in with more specific steps and examples.</div><div class=""><br class=""></div><div class="">To separate your data into severe and very severe, you could load the original data twice and use an if-condition and do-loops (page 2 <a href="https://www.ncl.ucar.edu/Document/Reference_Cards/NCL_scripting_language_reference_card_legal.pdf" class="">here</a>) based on the wind speed number to check every data point and delete data points that do not meet the criteria. Then save those two variables with your original time/latitude/longitude coordinates in a new file. To plot them both in one <a href="https://www.ncl.ucar.edu/Applications/vector.shtml" class="">map</a> you can use an <a href="https://www.ncl.ucar.edu/Applications/overlay.shtml" class="">overlay plot</a>, which means you plot one map first and then another map on top of it with the lines/arrows in a different color. Finally, add a legend to define the two colors.</div><div class=""><br class=""></div><div class="">Below and attached are a very basic NCL script which I used to create a spatial subset of global daily precipitation data. It might give you a general idea of the structure of an NCL script. First you load the data, then you manipulate it in the way you need to, then you save it in a new file, and in this case it finally plots the data (however, not using an overlay plot).</div><div class=""><br class=""></div><div class="">I hope this will help to get you started.</div><div class=""><br class=""></div><div class="">Best of luck,</div><div class="">Toni</div><div class=""><br class=""><div class=""><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div class=""><font face="Calibri" class=""><br class=""></font></div><div class=""><font face="Calibri" class=""><b class="">Toni Klemm, Ph.D.<br class=""></b>Postdoctoral Research Associate</font></div><div class="">Department of Ecosystem Science and Management</div><div class="">College of Agriculture and Life Sciences</div><div class="">Texas A&M University, College Station, TX</div><div class=""><font face="Calibri" class="">Contributor to the <a href="http://www.eccforum.org" class="">Early Career Climate Forum</a><br class=""><a href="http://www.toni-klemm.de" class="">www.toni-klemm.de</a> | <a href="http://twitter.com/toniklemm" class="">@toniklemm</a></font></div><div class=""><font face="Calibri" class=""><br class=""></font></div><div class=""><font face="Calibri" class=""><br class=""></font></div></div></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">; ---------------------------------------------------</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">; ---------------------------------------------------</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">; ***** GFDL FLOR-B01 daily precipitation data ******</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">; ---------------------------------------------------</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">begin</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> do year = 1980,2013<span class="Apple-tab-span" style="white-space: pre;"> </span>; do-loop through the record years, 1980-2013</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> do m = 1,12<span class="Apple-tab-span" style="white-space: pre;"> </span> ; do-loop through the 12 calendar months</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> do r = 1,12<span class="Apple-tab-span" style="white-space: pre;"> </span>; do-loop through the 12 model runs</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> <span class="Apple-tab-span" style="white-space: pre;"> </span>if (m.lt.10) then</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><span class="Apple-tab-span" style="white-space: pre;"> </span> m_zero = "0" + m<span class="Apple-tab-span" style="white-space: pre;"> </span>; turn "1"..."9" into "01"..."09" to match filename structure (months have 2 digits)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><span class="Apple-tab-span" style="white-space: pre;"> </span>end if</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><span class="Apple-tab-span" style="white-space: pre;"> </span>if(m.gt.9) then</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><span class="Apple-tab-span" style="white-space: pre;"> </span> m_zero = m<span class="Apple-tab-span" style="white-space: pre;"> </span>; leave "10..."12" unchanged</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><span class="Apple-tab-span" style="white-space: pre;"> </span>end if</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("**** START: " + year + m_zero + ", run " + r + " ****")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; READ IN MULTIPLE FILES </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> nmme_files = systemfunc("ls pr_day_GFDL*-v3.1-" + m_zero + year + "_r" + r + "i1p1" + "*.nc")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> nmme_f = addfiles(nmme_files,"r")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ListSetType(nmme_f,"join")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 1/8 Read in multiple files")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; READ IN A SINGLE FILE</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; year = 1980</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; m = 3</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; r = 1</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; nmme_f = addfile("pr_day_GFDL-FLORB01_FLORB01-P1-ECDA-v3.1-031980_r1i1p1_19800301-19810228.nc","r")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lat = nmme_f[0]->lat</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lon = nmme_f[0]->lon </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("Finished: 1/8 Read in a single file")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; PRINT THE DIMENSIONS OF THE READ-IN FILE</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("the dimensions of latitude are:")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(dimsizes(lat))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("the dimensions of longitude are:")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(dimsizes(lon))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; printMinMax(lon,True) <span class="Apple-tab-span" style="white-space: pre;"> </span>; prints the minimum and maximum longitude (0.3125 ... 359.6875)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(lat) <span class="Apple-tab-span" style="white-space: pre;"> </span>; prints all 360 latitude entries</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 2/8 Print the dimensions of the read-in file")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; READ IN THE VARIABLE (precipitation)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip = nmme_f[:]->pr*3600*24<span class="Apple-tab-span" style="white-space: pre;"> </span> <span class="Apple-tab-span" style="white-space: pre;"> </span>; extract precipitation into new variable, output is mm/day</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("variable dimensions")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(dimsizes(precip))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; num_precip = num(precip)<span class="Apple-tab-span" style="white-space: pre;"> </span>; finds out the number of total values in all 12 model runs</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; max_precip = max(precip)<span class="Apple-tab-span" style="white-space: pre;"> </span>; finds the smallest value in any of the 12 model runs</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("The total number of values in the 12 1980 GFDL model runs are:")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(num_precip)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("Largest number in the GFDL 1980 files:")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(max_precip)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; since we used join command above, the dimensions of precip should be (number files,ndays,nlat,nlon)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 3/8 Read in the variable (precipitation)")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; CREATE US-SUBSET OF THE GLOBAL DATA</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; let's choose latitudes of </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> latind = (/20,51/)<span class="Apple-tab-span" style="white-space: pre;"> </span>; latind is a variable</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; let's choose longitudes of </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lonind = (/230,306/)<span class="Apple-tab-span" style="white-space: pre;"> </span>; lonind is a variable</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lti = ind_nearest_coord(latind,lat,0)<span class="Apple-tab-span" style="white-space: pre;"> </span>; lti, lni = variables, ind_nearest_coord = pre-defined functions to determine the indices of locations closest to the coordinate array, needed later to define precip_subdomain_4D</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lni = ind_nearest_coord(lonind,lon,0)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lat_subdomain = lat(lti(0):lti(1))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> lon_subdomain = lon(lni(0):lni(1))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_3D = precip(:,lti(0):lti(1),lni(0):lni(1))<span class="Apple-tab-span" style="white-space: pre;"> </span>; define the subdomain based on latind and lonind</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(lat_subdomain_4D) <span class="Apple-tab-span" style="white-space: pre;"> </span>; displays the latitudes of the chosen area</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(lon_subdomain_4D) <span class="Apple-tab-span" style="white-space: pre;"> </span>; displays the longitudes of the chosen area</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(dimsizes(precip_subdomain_3D)) <span class="Apple-tab-span" style="white-space: pre;"> </span>; prints all lats and lons of the subset </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> day = nmme_f[:]->time</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ensemble_member = ispan(1,12,1)<span class="Apple-tab-span" style="white-space: pre;"> </span>; there are 12 ensemble members</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 4/8 Create US-subset of the global data")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; WRITE THE SUBSET INTO NEW netCDF FILE</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> x = (/year,100/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> y = product(x)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> z = (/y,m/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> init = sum(z)<span class="Apple-tab-span" style="white-space: pre;"> </span>; filename element in YYYYMM format to indicate initiation month</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> netCDF1 = True<span class="Apple-tab-span" style="white-space: pre;"> </span>; Output format is NetCDF</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> if (netCDF1) then</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> diro = "./subset/nc/"<span class="Apple-tab-span" style="white-space: pre;"> </span>; output in subdirectory "subset" within input directory</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filename = "GFDL_US_pr_daily_init_" + init + "run_" + r +".nc"<span class="Apple-tab-span" style="white-space: pre;"> </span>; output filename for subset using initialization year and month to differentiate</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filo = filename</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> end if</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> setfileoption("nc","Format","LargeFile")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> if (netCDF1) then</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> system("/bin/rm -f "+diro+filo)<span class="Apple-tab-span" style="white-space: pre;"> </span>; removes old files with the same name if they are present</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf = addfile(diro+filo,"c")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> setfileoption(ncdf,"DefineMode",True)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; create attributes</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fAtt = True</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fAtt@title = "Data: Daily GFDL FLOR-B01, run: " + r + " i1p1, for the Contiguous U.S., initialized " + init</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fAtt@source = "Vecchi et al. 2014"</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fAtt@Conventions = "None"</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fAtt@creation_date = systemfunc("date")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> fileattdef(ncdf,fAtt)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> dimNames = (/"year","ensemble_member","days","lat","lon"/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> dimSizes = (/-1,12,dimsizes(day),dimsizes(lat_subdomain),dimsizes(lon_subdomain)/)<span class="Apple-tab-span" style="white-space: pre;"> </span>; r = model run (1-12), was "12" before</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> dimUnlim = (/True,False,False,False,False/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filedimdef(ncdf,dimNames,dimSizes,dimUnlim)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"year",typeof(year),(/"year"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"lat",typeof(lat_subdomain),(/"lat"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"lon",typeof(lon_subdomain),(/"lon"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print("Variable day: ")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; Print(day)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"days",typeof(day),(/"days"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"ensemble_member",typeof(ensemble_member),(/"ensemble_member"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; get attributes of our variable</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> filevardef(ncdf,"precipitation",typeof(precip_subdomain_3D),(/"days","lat","lon"/))</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; now write all the variables to the file</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->year = (/year/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->ensemble_member = (/ensemble_member/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->days = (/day/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->lat = (/lat_subdomain/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->lon = (/lon_subdomain/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ncdf->precipitation = (/precip_subdomain_3D/)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> end if</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 5/8 Write the subset into new netCDF file")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; REDUCE 3D ARRAY TO PLOTABLE 2D ARRAY</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(dimsizes(precip_subdomain_3D)) <span class="Apple-tab-span" style="white-space: pre;"> </span>; prints all lats and lons of the subset </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; printVarSummary(precip_subdomain_3D)<span class="Apple-tab-span" style="white-space: pre;"> </span>; print summary of subdomain precip: var name, type, byte size, # of values, dimensions (x,y,z,temp)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D = dim_avg_n(precip_subdomain_3D, 0)<span class="Apple-tab-span" style="white-space: pre;"> </span>; averaging all daily precipitation values into one annual average for every lat and lon, creating 2D array, needed to plot single map</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; printVarSummary(precip_subdomain_3D)<span class="Apple-tab-span" style="white-space: pre;"> </span>; print size (byte, values) and dimensions (x, y, z)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(precip_subdomain_3D)<span class="Apple-tab-span" style="white-space: pre;"> </span>; prints *ALL* precip values for the subset, about 5 million lines total!!!</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; print(precip_subdomain_3D(:,:,:))<span class="Apple-tab-span" style="white-space: pre;"> </span>; array subscripting for all subset lats and lons for 1980, all 12 model runs</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 6/8 Reduce 3D array to plottable 2D array")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; CREATE MAP</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; Define boundaries for map</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D!0 = "lat_subdomain"<span class="Apple-tab-span" style="white-space: pre;"> </span>; latitude information come from variable lat_subdomain</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D&lat_subdomain = lat_subdomain</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D!1 = "lon_subdomain"<span class="Apple-tab-span" style="white-space: pre;"> </span>; longitude information come from variable lon_subdomain</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D&lon_subdomain = lon_subdomain</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D@long_name = "Precipitation"<span class="Apple-tab-span" style="white-space: pre;"> </span>; map subtitle, printed on the top left</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> precip_subdomain_2D@units = "mm/day"<span class="Apple-tab-span" style="white-space: pre;"> </span>; map units, printed on the top right</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> diro = "./subset/png/"<span class="Apple-tab-span" style="white-space: pre;"> </span>; rename output directory to save map in png folder</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> wks = gsn_open_wks("png", diro + "GFDL_US_pr_daily_init_" + init + "run_" + r)<span class="Apple-tab-span" style="white-space: pre;"> </span>; map name and file type</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> cmap = read_colormap_file("BlueDarkRed18")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res = True<span class="Apple-tab-span" style="white-space: pre;"> </span>; plot mods desired</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@tiMainString = "Run " + r + ", init. " + m + "/" + year + " - GFDL-Hindcast, daily"<span class="Apple-tab-span" style="white-space: pre;"> </span>; main title</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@cnFillOn = True<span class="Apple-tab-span" style="white-space: pre;"> </span>; turn on color fill, works with gsn_csm_map() function 3 lines down</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@gsnAddCyclic = False</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@cnLinesOn = False<span class="Apple-tab-span" style="white-space: pre;"> </span>; contour lines of the map fill on or off </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpProjection = "Mercator"<span class="Apple-tab-span" style="white-space: pre;"> </span>; Mercator map projection</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@cnFillPalette = cmap(::-1,:)<span class="Apple-tab-span" style="white-space: pre;"> </span>; reverse color map</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpLimitMode = "LatLon"<span class="Apple-tab-span" style="white-space: pre;"> </span>; define plotted area by lats and lons, namely those of the selection (20-50N, 230-300E)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpMaxLatF = 50</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpMinLatF = 20</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; res@mpCenterLatF = 35<span class="Apple-tab-span" style="white-space: pre;"> </span>; define center latitude, only use with res@mpPrjection is conical equidistant</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpMaxLonF = 305</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpMinLonF = 230</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; res@mpCenterLonF = 265<span class="Apple-tab-span" style="white-space: pre;"> </span>; define center longitude, only use with res@mpPrjection is conical equidistant</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpOutlineBoundarySets = "GeophysicalAndUSStates"<span class="Apple-tab-span" style="white-space: pre;"> </span>; lines and boundaries of USA and US States</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpGeophysicalLineThicknessF = 2<span class="Apple-tab-span" style="white-space: pre;"> </span>; thickness of the USA lines</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpUSStateLineThicknessF = 2<span class="Apple-tab-span" style="white-space: pre;"> </span>; thickness of the state lines</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> res@mpSpecifiedFillColors = (/0,100/)<span class="Apple-tab-span" style="white-space: pre;"> </span>; fill with background color, works with gsn_map() 2 lines down</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 7/8 Create map")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; PLOT MAP</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; ***********************************************</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; plot = gsn_csm_contour_map_ce(wks,precip_subdomain_2D, res)<span class="Apple-tab-span" style="white-space: pre;"> </span>; conical equidistant projection, if used take out res@mpProjection line </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> plot = gsn_csm_contour_map(wks,precip_subdomain_2D, res)<span class="Apple-tab-span" style="white-space: pre;"> </span>; universal map projection, define in res@mpProjection</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("Finished: 8/8 Plot map")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> print("**** ALL DONE: " + init + ", run " + r + " ****")</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> delete(precip)<span class="Apple-tab-span" style="white-space: pre;"> </span>; delete before the next loop or program will fail</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> delete(precip_subdomain_3D)<span class="Apple-tab-span" style="white-space: pre;"> </span>; delete before the next loop or program will fail</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> ; delete(precip_subdomain_2D)<span class="Apple-tab-span" style="white-space: pre;"> </span>; delete before the next loop or program will fail</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> delete(day)<span class="Apple-tab-span" style="white-space: pre;"> </span>; delete before the next loop or program will fail</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> </div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> end do<span class="Apple-tab-span" style="white-space: pre;"> </span>; end the loop through the 12 model runs (1-12)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> end do<span class="Apple-tab-span" style="white-space: pre;"> </span>; end the loop through the 12 calendar months (1-12)</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"> end do<span class="Apple-tab-span" style="white-space: pre;"> </span>; end the loop through the years of 1980-2013</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">end</div><div dir="auto" class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><br class=""></div></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline"></div><div class=""><br class=""><blockquote type="cite" class=""><div class="">On Jan 29, 2019, at 11:50 AM, prithvi raj <<a href="mailto:rajjallu@gmail.com" class="">rajjallu@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hello all,<br class=""> I am rather new to NCL. I have cyclone track data from 1945-2018 as ascii files. In those files the date, lat, lon, wind speed at that particular lat and lon is mentioned. I want to segregate the them based on the intensity of cyclones into severe and very severe and plot them as bars in a time series. Can it be done using NCL? <br class="">_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class=""><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br class=""></div></div></blockquote></div></div><br class=""><br class=""><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div></div></div></div></div></div>
</div>
<br class=""></body></html>