<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><meta http-equiv="Content-Type" content="text/html charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Alan, <div class="">Unfortunately, I get a series of fatal errors if I don’t define them apriori. </div><div class=""><br class=""></div><div class=""><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (speed) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 32 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (latvals) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 33 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (lonvals) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 34 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (speed) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 36 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (coord) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 37 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (coord) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 38 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Variable (latvals) is undefined</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 42 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (wnout) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 43 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Undefined identifier: (wnout) is undefined, can't continue</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 44 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Variable (var0) is undefined</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 45 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255); min-height: 16px;" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class=""></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Variable (var1) is undefined</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 46 in file varcalc.ncl</span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class=""></span></div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Monaco; background-color: rgb(255, 255, 255);" class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class=""></span></div><div class=""><blockquote type="cite" class=""><div class="">On Oct 14, 2017, at 12:44 PM, Alan Brammer <<a href="mailto:abrammer@albany.edu" class="">abrammer@albany.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="gmail_default" style="font-family:verdana,sans-serif"><font color="#0061ff" style="font-family:arial,sans-serif;font-size:12px" class="">> I don’t know what my dimensions of </font><font color="#ff2600" style="font-family:arial,sans-serif;font-size:12px" class="">wnout</font><font color="#0061ff" style="font-family:arial,sans-serif;font-size:12px" class=""> will be. Will it be a </font><font color="#ff2600" style="font-family:arial,sans-serif;font-size:12px" class="">251x251x73</font><font color="#0061ff" style="font-family:arial,sans-serif;font-size:12px" class=""> array? </font><br class=""></div><div class="gmail_default" style="font-family:verdana,sans-serif"><font color="#0061ff" style="font-family:arial,sans-serif;font-size:12px" class=""><br class=""></font></div>It will<div class="gmail_default" style="font-family:verdana,sans-serif;display:inline">​ be a list, there won't be array dimensions.  You do not need to define wnout before hand. </div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline">likewise, you do not need to define ​var0 or var1 before hand.  </div></div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline"><br class=""></div></div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline">e.g.</div></div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline"><br class=""></div></div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline"><br class=""></div></div><div class=""><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline"><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> a = addfile(filename,"r")</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> u10 = a->u10</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> v10 = a->v10</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> slp = a->slp</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;; Constants ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> R         = 6371 </span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> d2r       = atan(1.0)/45; Equivalent of pi/180 - 0.01745329</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> r2d       = 1.0/d2r</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> toknots   = 1.94384</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> pi        = 3.14159</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class="">;; this section wasn't really needed</div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> slpsub    = slp(time,0,:,:) ; X is now a 2D array</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> speed     = wind_speed(u10(time,0,:,:),v10<wbr class="">(time,0,:,:)) * toknots</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> latvals   = u10&lat(0:800)</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lonvals   = u10&lon(0:809)</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> coord        = center_find(speed(:,:),slpsub)<wbr class=""> </span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lat1         = coord[0]</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lon1         = coord[1]</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;; Find all variables here and store into allvars ;;;;;;;</span></div><span class="m_5598048521643050020gmail-im" style="font-family:Monaco;font-size:12px"><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> wnout        = wn01(latvals,lonvals,lat1,lon1<wbr class="">,speed)   </font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> var0         = wnout[0] </font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> var1         = wnout[1]  </font></span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal;min-height:16px" class="">printVarSummary(var0)</div><div style="margin:0px;line-height:normal;min-height:16px" class="">printVarSummary(var1)</div><div style="margin:0px;line-height:normal;min-height:16px" class=""><br class=""></div></span><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal;" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;; Write to a netcdf file ;;;;;;;;;;;;;;;;;;</span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class=""><br class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span></div><div style="font-family: Monaco; font-size: 12px; margin: 0px; line-height: normal; min-height: 16px;" class="">;;;;   I do not think the lines that were here will work. </div></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Oct 14, 2017 at 10:26 AM, Prashanth Bhalachandran <span dir="ltr" class=""><<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank" class="">prashanth.bhalachandran@<wbr class="">gmail.com</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class=""><div style="word-wrap:break-word" class=""><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class="">Dear Marston, </font></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class="">Thank you for your response. I will ensure that the missing values are the default ones from NCL. I hope my attaching the code will give you a little more context. </font></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class=""><br class=""></font></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class="">My concern is this : If I am returning two arrays var0  and var01 (please see the red highlighted portions for array dimensions), how will I receive it in the main wrapper script? </font></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class=""><br class=""></font></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#3a88fe" class=""><span style="font-size:12px" class="">That is, if my var0 dimension is : </span></font><span style="color:rgb(255,38,0);font-family:Monaco;font-size:12px" class="">(/251/) </span><span style="color:rgb(0,97,255);font-size:12px" class="">and my var01 dimension is </span><span style="font-size:12px" class=""><font color="#ff2600" class="">(/251,73/)</font><font color="#0061ff" class="">, and I return the variable from the function using </font></span><span style="font-size:12px" class=""><font color="#ff2600" class="">return([/var0,var01/])</font></span><span style="color:rgb(0,97,255);font-size:12px" class="">, when I collect it in the wrapper script using </span><span style="color:rgb(0,97,255);font-size:12px" class=""> </span></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#ff2600" class=""><span style="font-size:12px" class="">wnout</span><span style="font-size:12px" class="">        </span><span style="font-size:12px" class="">= wn01(latvals,lonvals,lat1,lon1<wbr class="">,speed), </span></font></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><span style="font-size:12px" class=""><font color="#0061ff" class=""><br class=""></font></span></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><span style="font-size:12px" class=""><font color="#0061ff" class="">I don’t know what my dimensions of </font><font color="#ff2600" class="">wnout</font><font color="#0061ff" class=""> will be. Will it be a </font><font color="#ff2600" class="">251x251x73</font><font color="#0061ff" class=""> array? I do know that I want to extract it as: </font></span></div><span class=""><div style="font-size:12px;margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#0061ff" class=""> </font><font color="#ff2600" class="">var0         = wnout[0] </font></span></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#ff2600" class=""><span style="font-size:12px" class=""> var1</span><span style="font-size:12px" class=""> </span><span style="font-size:12px" class="">       </span><span style="font-size:12px" class=""> </span><span style="font-size:12px" class="">= wnout[1]</span><span style="font-size:12px" class=""> </span><span style="font-size:12px" class=""> </span></font><span style="font-size:12px" class=""><font color="#ff2600" class=""> </font><font color="#0061ff" class=""> </font></span></div><div style="margin:0px;line-height:normal;background-color:rgb(255,255,255)" class=""><font color="#ff2600" face="Monaco" class=""><span style="font-size:12px" class=""> </span></font></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><br class=""></div></span><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><b class=""><u class="">CODE</u></b><font color="#ff2600" style="font-family:Monaco" class=""> </font></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr class="">pts/csm/contributed.ncl"</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> load "/scratch/lustreD/s/sbhalach/D<wbr class="">ATA/func_center.ncl"</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> load "/scratch/lustreD/s/sbhalach/D<wbr class="">ATA/func_rtheta.ncl"</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> a = addfile(filename,"r")</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> u10 = a->u10</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> v10 = a->v10</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> slp = a->slp</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;; Constants ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> R         = 6371 </span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> d2r       = atan(1.0)/45; Equivalent of pi/180 - 0.01745329</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> r2d       = 1.0/d2r</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> toknots   = 1.94384</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> pi        = 3.14159</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> latvals   = new((/801/),float,1e+30)</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lonvals   = new((/810/),float,1e+30)</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> speed     = new((/801,810/),float,1e+30)  </span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> <font color="#ff2600" class="">var01     = new((/251,73/),float,1e+30)</font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> var0      = new((/251/),float,1e+30)</font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> </span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr class="">;</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> slpsub       = slp(time,0,:,:) ; X is now a 2D array</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> speed(:,:)   = wind_speed(u10(time,0,:,:),v10<wbr class="">(time,0,:,:)) * toknots</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> latvals(:)   = u10&lat(0:800)</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lonvals(:)   = u10&lon(0:809)</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> coord        = center_find(speed(:,:),slpsub)<wbr class=""> </span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lat1         = coord[0]</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> lon1         = coord[1]</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;; Find all variables here and store into allvars ;;;;;;;</span></div><span class=""><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> wnout        = wn01(latvals,lonvals,lat1,lon1<wbr class="">,speed)   </font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> var0         = wnout[0] </font></span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""> var1         = wnout[1]  </font></span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div></span><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;;;;;;;;;;;;;;;;;;;;;; Write to a netcdf file ;;;;;;;;;;;;;;;;;;</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> system("rm -f $f1_$<a href="http://f2.nc/" target="_blank" class="">f2.nc</a>")</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> ncdf = addfile("$f1_$<a href="http://f2.nc/" target="_blank" class="">f2.nc</a>","c")</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> ncdf->allvars = allvars</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> ncdf->yymmdd  = $f1</span></div><div style="margin:0px;line-height:normal" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> ncdf->time    = $f2</span></div><div style="margin:0px;line-height:normal;min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> </span><br class="m_5598048521643050020m_-8056821693561989853webkit-block-placeholder"></div><div class=""><span style="font-variant-ligatures:no-common-ligatures" class="">------------------------------<wbr class="">------------------------------<wbr class="">-------</span></div></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div style="margin:0px;font-size:12px;line-height:normal;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><b class=""><u class="">FUNCTION </u></b></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><br class=""></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">load "$NCARG_ROOT/lib/ncarg/nclscri<wbr class="">pts/csm/contributed.ncl"</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">load "/scratch/lustreD/s/sbhalach/D<wbr class="">ATA/func_center.ncl"</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><font color="#ff2600" class="">undef</font>("wn01")</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> function wn01(latvals:float,lonvals:flo<wbr class="">at,lat1:float,lon1:float,var:<wbr class="">float)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> <font color="#ff2600" class="">local</font> R,d2r,r2d,toknots,pi,delr,Xi,Y<wbr class="">i,xo,yo,nearxo,nearyo,ovar,var<wbr class="">01,var0,c1,b1,orad,otha,arr,<wbr class="">a1,a2,a3,a4,a5,a6,a7,a8,a9,r,<wbr class="">weights </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> begin </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> R         = 6371 </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> d2r       = atan(1.0)/45; Equivalent of pi/180 or 0.01745329</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> r2d       = 1.0/d2r</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> toknots   = 1.94384</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> pi        = 3.14159 </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> delr      = 12</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> Xi        = new((/810/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> Yi        = new((/801/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> ovar      = new((/251,73/),float,1e+30)  </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> var01     = new((/251,73/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> var0      = new((/251/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> c1        = new((/251/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> b1        = new((/251/),float,1e+30)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> orad      = ispan(0,1000,delr)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> otha      = fspan(0,6.283185,73)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> arr       = new((/3,3/),float,1e+30)  </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> Xi(:) = R *cos(lat1*d2r)*(lonvals(:)-lon<wbr class="">1) * d2r</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> Yi(:) = R *(latvals(:)-lat1) * d2r</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> </span><br class="m_5598048521643050020m_-8056821693561989853webkit-block-placeholder"></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> do radius=0,250,1 </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  do theta=0,72,1   </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    xo  = (radius*delr)*cos((theta*5)*d2<wbr class="">r)   ; reference points </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    yo  = (radius*delr)*sin((theta*5)*d2<wbr class="">r)   ; reference points</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  </span><br class="m_5598048521643050020m_-8056821693561989853webkit-block-placeholder"></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    nearxo = ind_nearest_coord (xo, Xi(:), 0)  ;Find the index of the bounds</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    nearyo = ind_nearest_coord (yo, Yi(:), 0) </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    if(nearxo .eq. 0 .or. nearxo .eq. 809 .or. nearyo .eq. 0 .or. nearyo .eq. 800) then</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">    ovar(radius,theta) = var(nearyo,nearxo)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  else</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a1 = var(nearyo+1,nearxo+1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a2 = var(nearyo+1,nearxo)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a3 = var(nearyo-1,nearxo-1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a4 = var(nearyo,nearxo+1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a5 = var(nearyo,nearxo)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a6 = var(nearyo,nearxo-1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a7 = var(nearyo-1,nearxo+1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a8 = var(nearyo-1,nearxo)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   a9 = var(nearyo-1,nearxo-1)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   arr = (/ (/a1,a2,a3/), (/a4,a5,a6/), (/a7,a8,a9/)/)   ; Written as nD arrays for visual clarity of neighbors</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   weights = (/ (/0.0625,0.0625,0.0625/), (/0.0625,0.5,0.0625/), (/0.0625,0.0625,0.0625/)/) </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   ovar(radius,theta) =  dim_avg_wgt(ndtooned(arr),ndto<wbr class="">oned(weights),1)  </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">;  Give 50% weight to the nearest location and distribute the other 50% equally amongst the neighbours  </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">   end if</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  end do ; Loop over every radius and theta </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> end do    ; End of r-theta loop</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  do r=0,250,1 </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">     <font color="#ff2600" class="">var0(r)</font>    =  1/(2*pi)*simpeq(ovar(r,:),5.0*<wbr class="">d2r)                    ; Dimension </span>(/251/)</div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">     c1(r)      =  1/(2*pi)*simpeq(ovar(r,:)*cos(<wbr class="">otha),5.0*d2r)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">     b1(r)      =  1/(2*pi)*simpeq(ovar(r,:)*sin(<wbr class="">otha),5.0*d2r)</span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">     <font color="#ff2600" class="">var01(r,:)</font> =  var0(r) + (c1(r)*cos(otha)) + (b1(r)*sin(otha))<span class="m_5598048521643050020m_-8056821693561989853Apple-tab-span" style="white-space:pre-wrap">      </span> ; Dimension </span>(/251,73/)</div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">  end do    </span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> <font color="#ff2600" class="">return([/var0,var01/])</font></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><br class=""></span></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255)" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""> end </span></div><div class=""><div class="m_5598048521643050020h5"><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><br class=""></div><div style="margin:0px;font-size:12px;line-height:normal;font-family:Monaco;background-color:rgb(255,255,255);min-height:16px" class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class="">On Oct 14, 2017, at 12:42 AM, Marston Johnston <<a href="mailto:shejo284@gmail.com" target="_blank" class="">shejo284@gmail.com</a>> wrote:</div><br class="m_5598048521643050020m_-8056821693561989853Apple-interchange-newline"><div class=""><div class="m_5598048521643050020m_-8056821693561989853WordSection1" style="font-family:PalatinoLinotype-Roman;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)"><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">There is an error in your code:        var0 = new((/251/),float,1e+30<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">As a preference, I always find it best to use the default FillValues when creating arrays. Read up on missing values/fillvalues on the NCL wedpage.<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class=""><u class=""></u> <u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">There is nothing wrong with returning a list from a function, even if the dimensions are different. Why would the dimensions matter?<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">This doesn’t seem to be your problem. You have not provided enough information about your code and variable information (printVarSumary) to give a definite answer.<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">It seems you need to do some more debugging of your code yourself.<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class=""><u class=""></u> <u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">When defining functions, it is good practice to set the define the intrinsic variables as “local” as well as using the “undef” function.<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span><u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class="">Using the same name for the variables inside and outside the function, if you are indeed using a function that is properly designed, can mask a bug in your code.<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span><u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class=""><u class=""></u> <u class=""></u></span></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">Marston S. Ward, PhD<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">Department of Earth Sciences<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">University of Gothenburg, Sweden<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">Email: <a href="mailto:marston.johnston@gu.se" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" class=""><span style="color:rgb(0,112,192)" class="">marston.johnston@gu.se</span></a><u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">SkypeID: marston.johnston <u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">Phone: <a href="tel:+46%2031%20786%2049%2001" value="+46317864901" target="_blank" class="">+46-31-7864901</a> <u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">Only the fruitful thing is true!<u class=""></u><u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span style="font-size:12pt;font-family:'Times New Roman',serif;color:rgb(0,112,192)" class="">~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<u class=""></u><u class=""></u></span></div></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class=""><u class=""></u> <u class=""></u></span></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><span class=""><u class=""></u> <u class=""></u></span></div><div style="border-style:solid none none;border-top-width:1pt;border-top-color:rgb(181,196,223);padding:3pt 0cm 0cm" class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><b class=""><span style="font-size:12pt" class="">From:<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span></span></b><span style="font-size:12pt" class="">ncl-talk <<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank" class="">ncl-talk-bounces@ucar.edu</a>> on behalf of Prashanth Bhalachandran <<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank" class="">prashanth.bhalachandran@gmail<wbr class="">.com</a>><br class=""><b class="">Date:<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span></b>Saturday, 14 October 2017 at 03:57<br class=""><b class="">To:<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span></b><<a href="mailto:ncl-talk@ucar.edu" target="_blank" class="">ncl-talk@ucar.edu</a>><br class=""><b class="">Subject:<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span></b>Re: [ncl-talk] Function return<u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><u class=""></u> <u class=""></u></div></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">Hello, <u class=""></u><u class=""></u></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">I have a quick question regarding returning of arrays from functions. <u class=""></u><u class=""></u></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class=""><u class=""></u> <u class=""></u></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">In my code, I have two arrays var0 and var01 that I want to return from my function. Their dimensionalities are as follows: <u class=""></u><u class=""></u></div></div><div class=""><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""><br class=""><br class=""><u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class="">var0      = new((/251/),float,1e+30<u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class="">var01     = new((/251,73/),float,1e+30)<u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""> <span class="m_5598048521643050020m_-8056821693561989853apple-tab-span">   <span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span></span><u class=""></u><u class=""></u></span></div></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""><u class=""></u> <u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class="">At present, I am trying to return this using the command : <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""><u class=""></u> <u class=""></u></span></div></div><div class=""><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(77,47,45)" class="">return([/var0,var01/]) ; Note that the dimensions of var0 and var1 are different. <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(77,47,45)" class=""><br class=""><br class=""><u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(77,47,45)" class="">As a result, in my wrapper script, when I receive the variables, it is becoming a confusing task. <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(77,47,45)" class=""><br class=""><br class=""><u class=""></u><u class=""></u></span></div></div><div class=""><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""> wnout        = wn01(latvals,lonvals,lat1,lon1<wbr class="">,speed)   <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""> var0         = wnout[0] <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""> var1         = wnout[1]  <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""><br class=""><br class=""><u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""> The above return procedure is giving out an error since I don’t know how to define the dimensions of wnout. Is there any better way to return two arrays of varying dimensions in NCL? <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class=""><br class=""><br class=""><u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class="">Thank you, <u class=""></u><u class=""></u></span></div></div><div class=""><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif;background-color:rgb(237,234,212)" class=""><span style="font-size:12pt;font-family:Courier,serif;color:rgb(184,45,24)" class="">Prashanth  <u class=""></u><u class=""></u></span></div></div></div></div><div style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">______________________________<wbr class="">_________________ ncl-talk mailing list<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span><a href="mailto:ncl-talk@ucar.edu" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" class="">ncl-talk@ucar.edu</a><span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span>List instructions, subscriber options, unsubscribe:<span class="m_5598048521643050020m_-8056821693561989853Apple-converted-space"> </span><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" class="">http://mailman.uc<wbr class="">ar.edu/mailman/listinfo/ncl-ta<wbr class="">lk</a><u class=""></u><u class=""></u></div></div></div></blockquote></div><br class=""></div></div></div></div><br class="">______________________________<wbr class="">_________________<br class="">
ncl-talk mailing list<br class="">
<a href="mailto:ncl-talk@ucar.edu" target="_blank" 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" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/mailma<wbr class="">n/listinfo/ncl-talk</a><br class="">
<br class=""></blockquote></div><br class=""></div></div>
</div></blockquote></div><br class=""></div></div></body></html>