<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<!--[if lte mso 15 || CheckWebRef]-->
<div id="OwaReferenceAttachments" contenteditable="false">
<table style="padding-bottom: 13px; border-width: 0px; border-style: none;">
<tbody>
<tr valign="top">
<td>
<table style="border-width: 0px 0px 1px 0px; border-color:#C7C7C7; border-style: none none dotted none;">
<tbody>
<tr valign="top">
<td style="padding-bottom:7px;">
<table align="left" style="padding-right: 28px; border-width: 0px; background-color: rgb(255, 255, 255); border-spacing: 0px">
<tbody>
<tr valign="top">
<td style="padding: 0px;">
<div id="OwaReferenceAttachmentDescription" style="padding-left: 3px; font-size: 14px; font-family: 'Segoe UI', 'Segoe WP', 'Segoe UI WPC', Tahoma, Arial, sans-serif; color: rgb(102, 102, 102);">
Liu, Xiaoli has shared OneDrive for Business files with you. To view them, click the links below.
</div>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
<tr valign="top">
<td><a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=043cc05571c414a5db1454ae19a28cc36&authkey=AYoZc-0X9kH7gkzGJO7GBMM&expiration=2017-09-29T15%3a03%3a51.000Z" target="_blank">
<table align="left" style="padding-right: 28px; padding-bottom:10px; border-width: 0px; height:20px; background-color: rgb(255, 255, 255); border-spacing: 0px">
<tbody>
<tr valign="top">
<td style="padding: 0px;">
<div style="background-color: rgb(255, 255, 255); height: 20px; width: 20px; max-height: 20px;">
<a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=043cc05571c414a5db1454ae19a28cc36&authkey=AYoZc-0X9kH7gkzGJO7GBMM&expiration=2017-09-29T15%3a03%3a51.000Z" target="_blank"><img width="20" style="border:0px;" src="https://r1.res.office365.com/owa/prem/images/dc-png_20.png"></a></div>
</td>
<td>
<div id="OwaReferenceAttachmentFileName2" style="padding: 0px 0px 0px 5px; font-size: 14px; font-family: 'Segoe UI', 'Segoe WP', 'Segoe UI WPC', Tahoma, Arial, sans-serif; color: rgb(0, 114, 198);">
<a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=043cc05571c414a5db1454ae19a28cc36&authkey=AYoZc-0X9kH7gkzGJO7GBMM&expiration=2017-09-29T15%3a03%3a51.000Z" target="_blank" style="text-decoration: none; margin: 0px; font-size: 14px; font-family: 'Segoe UI', 'Segoe WP', 'Segoe UI WPC', Tahoma, Arial, sans-serif; color: rgb(0, 114, 198);">Orbit
View.png</a></div>
</td>
<td style="display:none;visibility:hidden;" width="0" height="0"><img width="0" height="0" style="visibility:hidden;border:0px;display:none" title="Orbit View.png" src="cid:6feea805-81b8-411b-aa66-2ae4d4b40d6f"></td>
</tr>
</tbody>
</table>
</a><a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=1a7f340c75074400a8b8d8c29d359140d&authkey=AaHal-OUvrqFxwfnh5iAIpE&expiration=2017-09-29T15%3a03%3a56.000Z" target="_blank">
<table align="left" style="padding-right: 28px; padding-bottom:10px; border-width: 0px; height:20px; background-color: rgb(255, 255, 255); border-spacing: 0px">
<tbody>
<tr valign="top">
<td style="padding: 0px;">
<div style="background-color: rgb(255, 255, 255); height: 20px; width: 20px; max-height: 20px;">
<a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=1a7f340c75074400a8b8d8c29d359140d&authkey=AaHal-OUvrqFxwfnh5iAIpE&expiration=2017-09-29T15%3a03%3a56.000Z" target="_blank"><img width="20" style="border:0px;" src="https://r1.res.office365.com/owa/prem/images/dc-generic_20.png"></a></div>
</td>
<td>
<div id="OwaReferenceAttachmentFileName2" style="padding: 0px 0px 0px 5px; font-size: 14px; font-family: 'Segoe UI', 'Segoe WP', 'Segoe UI WPC', Tahoma, Arial, sans-serif; color: rgb(0, 114, 198);">
<a href="https://sooners-my.sharepoint.com/personal/xiaoli_liu_ou_edu/_layouts/15/guestaccess.aspx?docid=1a7f340c75074400a8b8d8c29d359140d&authkey=AaHal-OUvrqFxwfnh5iAIpE&expiration=2017-09-29T15%3a03%3a56.000Z" target="_blank" style="text-decoration: none; margin: 0px; font-size: 14px; font-family: 'Segoe UI', 'Segoe WP', 'Segoe UI WPC', Tahoma, Arial, sans-serif; color: rgb(0, 114, 198);">file.out</a></div>
</td>
<td style="display:none;visibility:hidden;" width="0" height="0"><img width="0" height="0" style="visibility:hidden;border:0px;display:none" title="file.out" src="cid:8a52ad5e-886e-4fe8-b595-5f69ff11aefc"></td>
</tr>
</tbody>
</table>
</a></td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</div>
<div id="OwaReferenceAttachmentsEnd" style="display:none;visibility:hidden;"></div>
<!--[endif]-->
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" dir="ltr">
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" dir="ltr">
<p>Hello,</p>
<p>I tried to print the max and min value, and found them all zero. But when I use the Orbit Viewer software, I found there shoul be some values not equal zero. So I don't know why. I sent the data file "<span>2A.GPM.DPR.V7-20170308.20150627-S000754-E014026.007536.V05A</span>",
the NCL scripy"gpm.ncl" via ftp.</p>
And I aosl attached the printscreen of Orbit viewer and the output file of ncl_filedump command.</div>
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" dir="ltr">
Thanks!</div>
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" dir="ltr">
Xiaoli<br>
<br>
<div style="color: rgb(0, 0, 0);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> ncl-talk <ncl-talk-bounces@ucar.edu> on behalf of ncl-talk-request@ucar.edu <ncl-talk-request@ucar.edu><br>
<b>Sent:</b> Wednesday, August 30, 2017 4:39 PM<br>
<b>To:</b> ncl-talk@ucar.edu<br>
<b>Subject:</b> ncl-talk Digest, Vol 165, Issue 49</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Send ncl-talk mailing list submissions to<br>
ncl-talk@ucar.edu<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" id="LPlnk827561" previewremoved="true">
http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
<div id="LPBorder_GT_15041048177300.7009210482165977" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_15041048177280.4822220164087656" role="presentation" cellspacing="0" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);">
<tbody>
<tr valign="top" style="border-spacing: 0px;">
<td id="TextCell_15041048177290.5510832957142069" colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;">
<div id="LPRemovePreviewContainer_15041048177290.8928577921665406"></div>
<div id="LPTitle_15041048177290.9119805796271108" style="top: 0px; color: rgb(0, 120, 215); font-weight: normal; font-size: 21px; font-family: wf_segoe-ui_light, "Segoe UI Light", "Segoe WP Light", "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; line-height: 21px;">
<a id="LPUrlAnchor_15041048177290.6906814983226757" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank" style="text-decoration: none;">ncl-talk Info Page - University Corporation for ...</a></div>
<div id="LPMetadata_15041048177300.30595630142900165" style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 14px;">
mailman.ucar.edu</div>
<div id="LPDescription_15041048177300.13891612456844538" style="display: block; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;">
This email list is for NCL users to exchange information, ask questions, and report bugs on NCL or other related software if appropriate. It is also for the ...</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
or, via email, send a message with subject or body 'help' to<br>
ncl-talk-request@ucar.edu<br>
<br>
You can reach the person managing the list at<br>
ncl-talk-owner@ucar.edu<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: fatal:["Execute.c":8640]:Execute: Error occurred at or<br>
near line 897 in file<br>
$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl (Dennis Shea)<br>
2. Fw: Can't got value of variables other than Lat and Lon<br>
(Liu, Xiaoli)<br>
3. Re: problem with adding file (Mary Haley)<br>
4. Re: Fw: Can't got value of variables other than Lat and Lon<br>
(Mary Haley)<br>
5. Remove SST data for specific area (Dinuka Silva)<br>
6. Re: problem with adding file (Setareh Rahimi)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 29 Aug 2017 14:57:53 -0600<br>
From: Dennis Shea <shea@ucar.edu><br>
To: Guido Cioni <guidocioni@gmail.com><br>
Cc: ncl-talk <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] fatal:["Execute.c":8640]:Execute: Error<br>
occurred at or near line 897 in file<br>
$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl<br>
Message-ID:<br>
<CAOF1d_4aUYmTyaL+DHvAxhaMpVfhvHwM0bzUQBNNQbGeNG4sxw@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
I fixed an error and added writing a netCDF file.<br>
<br>
<a href="http://www.ncl.ucar.edu/Document/Functions/Crop/refevt_hargreaves_fao56.shtml" id="LPlnk168526" previewremoved="true">http://www.ncl.ucar.edu/Document/Functions/Crop/refevt_hargreaves_fao56.shtml</a>
<div id="LPBorder_GT_15041048177640.6989375302528256" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_15041048177620.4166410927812778" role="presentation" cellspacing="0" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);">
<tbody>
<tr valign="top" style="border-spacing: 0px;">
<td id="TextCell_15041048177630.1424655387912861" colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;">
<div id="LPRemovePreviewContainer_15041048177630.5378948963196042"></div>
<div id="LPTitle_15041048177630.9122508154070359" style="top: 0px; color: rgb(0, 120, 215); font-weight: normal; font-size: 21px; font-family: wf_segoe-ui_light, "Segoe UI Light", "Segoe WP Light", "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; line-height: 21px;">
<a id="LPUrlAnchor_15041048177630.055519843389450596" href="http://www.ncl.ucar.edu/Document/Functions/Crop/refevt_hargreaves_fao56.shtml" target="_blank" style="text-decoration: none;">refevt_hargreaves_fao56 - CISL's NCAR Command Language (NCL)</a></div>
<div id="LPMetadata_15041048177630.2845136189906168" style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 14px;">
www.ncl.ucar.edu</div>
<div id="LPDescription_15041048177640.015519475576827757" style="display: block; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;">
NCL Home > Documentation > Functions > Meteorology, Crop refevt_hargreaves_fao56. Use the Hargreaves ETo equation to derive reference evapotranspiration as described ...</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
Example 2<br>
<br>
On Tue, Aug 29, 2017 at 8:37 AM, Dennis Shea <shea@ucar.edu> wrote:<br>
<br>
> THX for sending the files. I have added Example 2 for<br>
><br>
> <a href="http://www.ncl.ucar.edu/Document/Functions/Crop/" id="LPlnk851481" previewremoved="true">
http://www.ncl.ucar.edu/Document/Functions/Crop/</a>
<div id="LPBorder_GT_15041049393560.2314944841096187" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_15041049393560.025477141462297848" role="presentation" cellspacing="0" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);">
<tbody>
<tr valign="top" style="border-spacing: 0px;">
<td id="TextCell_15041049393560.576358635688127" colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;">
<div id="LPRemovePreviewContainer_15041049393560.22118093675577555"></div>
<div id="LPTitle_15041049393560.550755152357469" style="top: 0px; color: rgb(0, 120, 215); font-weight: normal; font-size: 21px; font-family: wf_segoe-ui_light, "Segoe UI Light", "Segoe WP Light", "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; line-height: 21px;">
<a id="LPUrlAnchor_15041049393560.15126809833468724" href="http://www.ncl.ucar.edu/Document/Functions/Crop/" target="_blank" style="text-decoration: none;">crop.ncl functions</a></div>
<div id="LPMetadata_15041049393560.35157321941459263" style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 14px;">
www.ncl.ucar.edu</div>
<div id="LPDescription_15041049393560.04649399816237576" style="display: block; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;">
Description of crop functions in crop.ncl</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
> refevt_hargreaves_fao56.shtml<br>
><br>
> ----<br>
> Note: Given the size of the files (each > 4GB), this calculation takes<br>
> some time (about 5 minutes)<br>
><br>
> radext = radext_fao56(jday, lat2d, radunit)<br>
> ; <=== 186 seconds<br>
><br>
> refevt = refevt_hargreaves_fao56( tmin, tmax, radext, (/0,0,0/) ) ; <===<br>
> 113 seconds<br>
><br>
> =============<br>
> A tip when testing scripts on large files, it may be best to implement a<br>
> 'TEST' option:<br>
><br>
> TEST = True<br>
> ...<br>
> if (TEST) then<br>
> tStrt= 0<br>
> tLast= 100<br>
> time = fmin->time(tStrt:tLast) ; small subset<br>
> else<br>
> time = fmin->time ; entire file<br>
> end if<br>
><br>
> ...<br>
><br>
> ;---Read daily min & max<br>
> if (TEST) then<br>
> tmin = fmin->tasmin(tStrt:tLast,:,:)<br>
> tmax = fmax->tasmax(tStrt:tLast,:,:)<br>
> else<br>
> tmin = fmin->tasmin<br>
> tmax = fmax->tasmax<br>
> end if<br>
><br>
> ====<br>
><br>
> This would provide musch faster turn around for debugging.<br>
><br>
> Good Luck<br>
><br>
><br>
><br>
><br>
> On Wed, Aug 23, 2017 at 10:30 AM, Dennis Shea <shea@ucar.edu> wrote:<br>
><br>
>> Guido's response is correct.<br>
>><br>
>> Still, a 'trick' when encountering a malloc error is to use a small<br>
>> subset of the data<br>
>><br>
>> Rather than:<br>
>><br>
>> tasmax = a->tasmax<br>
>><br>
>> tasmin = b->tasmin<br>
>><br>
>> Pr_annual_total = c->pr<br>
>><br>
>> lon = a->longitude<br>
>><br>
>> lat = a->latitude<br>
>> time = a->time<br>
>><br>
>> Use<br>
>> NT = 100 ; arbitrary<br>
>><br>
>> tasmax = a->tasmax(0:NT,:,:)<br>
>><br>
>> tasmin = b->tasmin(0:NT,:,:)<br>
>><br>
>> Pr_annual_total = c->pr(0:NT,:,:)<br>
>><br>
>> lon = a->longitude<br>
>><br>
>> lat = a->latitude<br>
>><br>
>><br>
>> Note: rather than reading the 'time' variable directly<br>
>> time = a->time(0:NT)<br>
>> I suggest using the 'time' coordinate associated with the variable<br>
>> time = tasmax&time<br>
>><br>
>><br>
>> ====<br>
>><br>
>><br>
>> On Wed, Aug 23, 2017 at 5:15 AM, Guido Cioni <guidocioni@gmail.com><br>
>> wrote:<br>
>><br>
>>> You have to increase the maximum memory (ws maximum size) used by ncl to<br>
>>> plot. See here (<a href=""></a>https://www.ncl.ucar.edu/Document/Graphics/hlures.shtml<br>
>>> ).<br>
>>><br>
>>> Additionally, depending on the machine, you may need to increase the<br>
>>> maximum memory available for ncl at execution time (like ulimit in Linux).<br>
>>> Cheers<br>
>>><br>
>>> Il 23 ago 2017 1:09 PM, "Be?ta Szab?-Tak?cs" <szabo.b@czechglobe.cz> ha<br>
>>> scritto:<br>
>>><br>
>>>> Dear NCL Users,<br>
>>>><br>
>>>><br>
>>>><br>
>>>> I would like to calculate Aridity index base on Hargreaves ETo<br>
>>>> equation. I use daily maximum and minimum temperature data from 1971 to<br>
>>>> 2000 and yearly sum precipitation (30 years). The dimension of latitude is<br>
>>>> (201x464) and the day of year is (10957) vector for radext_fao56 function.<br>
>>>> When I run the scipt I get the following error message:<br>
>>>><br>
>>>> fatal:NclMalloc Failed:[errno=12]<br>
>>>><br>
>>>> fatal:New: could not create new array:[errno=12]<br>
>>>><br>
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 897<br>
>>>> in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl<br>
>>>><br>
>>>><br>
>>>><br>
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 162<br>
>>>> in file lambert_PET.ncl<br>
>>>><br>
>>>><br>
>>>><br>
>>>> ; Main code.<br>
>>>><br>
>>>> ;----------------------------------------------------------------------<br>
>>>><br>
>>>> begin<br>
>>>><br>
>>>> ;************************************************<br>
>>>><br>
>>>> ; read in netCDF file<br>
>>>><br>
>>>> ;************************************************<br>
>>>><br>
>>>> a = addfile("E-OBS_tasmax.nc","r")<br>
>>>><br>
>>>> b = addfile("E-OBS_tasmin.nc","r")<br>
>>>><br>
>>>> c = addfile("E-OBS_1971-2000_ys.nc","r")<br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>> tasmax = a->tasmax<br>
>>>><br>
>>>> tasmin = b->tasmin<br>
>>>><br>
>>>> Pr_annual_total = c->pr<br>
>>>><br>
>>>> lon = a->longitude<br>
>>>><br>
>>>> lat = a->latitude<br>
>>>><br>
>>>> time = a->time<br>
>>>><br>
>>>> minlat = 30 ; min lat to mask<br>
>>>><br>
>>>> maxlat = 75 ; max lat to mask<br>
>>>><br>
>>>> minlon = -20 ; min lon to mask<br>
>>>><br>
>>>> maxlon = 40 ; max lon to mask<br>
>>>><br>
>>>><br>
>>>><br>
>>>> nlat = dimsizes(lat)<br>
>>>><br>
>>>> mlon = dimsizes(lon)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> lat_new=conform_dims( (/nlat,mlon/), lat, 0)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> time2 = calendar_decode2(time,-5)<br>
>>>><br>
>>>> doy = day_of_year(time2(:,0),time2(:,1),time2(:,2))<br>
>>>><br>
>>>> print(doy)<br>
>>>><br>
>>>> radext = radext_fao56(doy, lat_new, 0)<br>
>>>><br>
>>>> radext = where(ismissing(radext), 0, radext)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> evtH_0 = refevt_hargreaves_fao56( tmin, tmax, radext, (/0,0,0/) )<br>
>>>><br>
>>>><br>
>>>><br>
>>>> print(evtH_0)<br>
>>>><br>
>>>> opt = True<br>
>>>><br>
>>>> opt@nval_crit = 12<br>
>>>><br>
>>>><br>
>>>><br>
>>>> PET_m = calculate_monthly_values (evtH_0, "sum", 0,opt)<br>
>>>><br>
>>>> print(PET_m)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> PET_annual_total = month_to_annual(PET_m, 0)<br>
>>>><br>
>>>> ; Pr_annual_total = month_to_annual(Pr_m, 0)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> AI = Pr_annual_total/PET_annual_total<br>
>>>><br>
>>>> AI2 = dim_avg_n(AI,0)<br>
>>>><br>
>>>> IM = new(dimsizes(AI2),typeof(AI2))<br>
>>>><br>
>>>><br>
>>>><br>
>>>> IM = where(AI2 .lt. 0.03,1.0,IM)<br>
>>>><br>
>>>> IM = where(AI2 .ge. 0.03 .and. AI2 .lt. 0.2,2.0,IM)<br>
>>>><br>
>>>> IM = where(AI2 .ge. 0.2 .and. AI2 .lt. 0.5,3.0,IM)<br>
>>>><br>
>>>> IM = where(AI2 .ge. 0.5 .and. AI2 .lt. 0.65,4.0,IM)<br>
>>>><br>
>>>> IM = where(AI2 .ge. 0.65,5.0,IM)<br>
>>>><br>
>>>><br>
>>>><br>
>>>> The error refer to the radext_fao56 in my script which is attached as<br>
>>>> well. Could someone suggest a solution?<br>
>>>><br>
>>>><br>
>>>><br>
>>>> Thank you for your help in advance!<br>
>>>><br>
>>>><br>
>>>><br>
>>>> Best regards,<br>
>>>><br>
>>>> Beata<br>
>>>><br>
>>>> _______________________________________________<br>
>>>> ncl-talk mailing list<br>
>>>> ncl-talk@ucar.edu<br>
>>>> List instructions, subscriber options, unsubscribe:<br>
>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>>>><br>
>>>><br>
>>> _______________________________________________<br>
>>> ncl-talk mailing list<br>
>>> ncl-talk@ucar.edu<br>
>>> List instructions, subscriber options, unsubscribe:<br>
>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>>><br>
>>><br>
>><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d0d92219/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d0d92219/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 29 Aug 2017 21:38:25 +0000<br>
From: "Liu, Xiaoli" <xiaoli.liu@ou.edu><br>
To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: [ncl-talk] Fw: Can't got value of variables other than Lat<br>
and Lon<br>
Message-ID:<br>
<DM5PR03MB31776AD4934CF7B8F7B5DDE7FE9F0@DM5PR03MB3177.namprd03.prod.outlook.com><br>
<br>
Content-Type: text/plain; charset="iso-2022-jp"<br>
<br>
<br>
<br>
<br>
________________________________<br>
From: Liu, Xiaoli <xiaoli.liu@ou.edu><br>
Sent: Wednesday, August 30, 2017 6:18 AM<br>
To: ncl-talk@ucar.edu<br>
Subject: Can't got value of variables other than Lat and Lon<br>
<br>
<br>
Hello,<br>
<br>
I want use the data of GPM 2A.GPM.DPR.V7 with NCL. I found that I can get the value of Lat and lon correctlly. But as to the other variables, the value got is very unusal. When I use orbit view, I found there are normal values for these variables.<br>
<br>
Here is my script?<br>
<br>
dir="/scratch/liu4119/gpm/data/"<br>
FILES=systemfunc(" ls -1 "+dir+"*.HDF5")<br>
numFILES = dimsizes(FILES)<br>
; print("numFILES = " + numFILES)<br>
; print(FILES)<br>
; print (" ")<br>
do i=32,32<br>
h5_file = addfile(FILES(i),"r")<br>
; g = h5_file=>/NS<br>
; print(g)<br>
; g1=g=>/SLV<br>
; print(g1)<br>
; list_filevars(g)<br>
data_raw = h5_file->/NS/SLV/precipRateESurface<br>
print(data_raw(5396,:))<br>
data = data_raw<br>
data@long_name = "precipRate"<br>
data@_FillValue = -9999.9<br>
data@units = data_raw@Units<br>
longitude=h5_file->/NS/Longitude<br>
printVarSummary(longitude)<br>
latitude=h5_file->/NS/Latitude<br>
print(latitude(5396,:))<br>
print(longitude(5396,:))<br>
printVarSummary(latitude)<br>
<br>
[cid:c1ec0f44-513a-44d6-964c-a1000a7b32ce]<br>
[cid:891609e8-6d88-4b50-ace4-52853a685e3f]<br>
[cid:ab77d561-deb8-4ef8-82d6-06293786eb5d]<br>
Thanks!<br>
<br>
Xiaoli<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: data_raw.png<br>
Type: image/png<br>
Size: 25095 bytes<br>
Desc: data_raw.png<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0003.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0003.png</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: latitude.png<br>
Type: image/png<br>
Size: 33981 bytes<br>
Desc: latitude.png<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0004.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0004.png</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: longitude.png<br>
Type: image/png<br>
Size: 32193 bytes<br>
Desc: longitude.png<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0005.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/fd739654/attachment-0005.png</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 29 Aug 2017 15:54:03 -0600<br>
From: Mary Haley <haley@ucar.edu><br>
To: Setareh Rahimi <setareh.rahimi@gmail.com><br>
Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] problem with adding file<br>
Message-ID:<br>
<CACNN_CJAeJ9g4sCNizTh4x=EHMTwS7PEmZ6XvvHinX-U1EkUAg@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
As another test, please see if you can look at the file using<br>
"ncl_filedump" from the UNIX command line. But, first make sure you have<br>
the file correctly spelled:<br>
<br>
ls precip.mon.total.v6.nc<br>
<br>
then:<br>
<br>
ncl_filedump precip.mon.total.v6.nc<br>
<br>
?This is the same thing as the 2-line script that Dennis asked you to try.<br>
If this works, then you can also verify at the same time if "precip" is a<br>
variable in that file.<br>
<br>
If you get any errors from the ncl_filedump command, respond back to<br>
ncl-talk with all of the error messages included.<br>
<br>
--Mary<br>
<br>
<br>
<br>
On Tue, Aug 29, 2017 at 11:45 AM, Dennis Shea <shea@ucar.edu> wrote:<br>
<br>
> Hello ... Please respond to ncl-talk and not directly to me.<br>
><br>
> Your script has:<br>
><br>
> *f * = addfile("precip.mon.total.v6.nc","r")<br>
> prc =* f*->precip({tStrt:tLast},:,:) ; whole grid<br>
><br>
> The error message is:<br>
><br>
> fatal: file (*f*) isn't defined<br>
><br>
> I am not sure what the problem would be. Someone else will have to respond.<br>
><br>
> ---<br>
> What does the following yoeld? Likely the same error.<br>
><br>
> *f * = addfile("precip.mon.total.v6.nc","r")<br>
> print(*f*)<br>
><br>
> =====================================================<br>
> NOTE:<br>
> You do *not* have to load the following libraries from NCL 6.2.0 onward.<br>
> So for 6.3.0, they are *automatically loaded.*<br>
><br>
> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
><br>
> From 6.4.0, the following are automatically loaded:<br>
><br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/bootstrap.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/extval.ncl"<br>
> "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br>
><br>
><br>
> On Tue, Aug 29, 2017 at 9:27 AM, Setareh Rahimi <setareh.rahimi@gmail.com><br>
> wrote:<br>
><br>
>> dear Dennis,<br>
>><br>
>> Thank you for your advice, but both data file and script are in the same<br>
>> directory.<br>
>><br>
>> Regards,<br>
>><br>
>> On Tue 29 Aug 2017 at 19:12, Dennis Shea <shea@ucar.edu> wrote:<br>
>><br>
>>> By default, NCL always looks in the 'current directory'.<br>
>>> My guess is that the data file is not in the same directory as the NCL<br>
>>> script.<br>
>>><br>
>>> dirp = "/some/path/to/file/"<br>
>>> filp = "precip.mon.total.v6.nc"<br>
>>><br>
>>> pthp = dirp + filp<br>
>>> f = addfile(pthp, "r")<br>
>>><br>
>>> or<br>
>>><br>
>>> f = addfile(dirp + filp, "r")<br>
>>><br>
>>><br>
>>> On Tue, Aug 29, 2017 at 8:26 AM, Setareh Rahimi <<br>
>>> setareh.rahimi@gmail.com> wrote:<br>
>>><br>
>>>> Dear all,<br>
>>>><br>
>>>> Once I put my NCL command line in NCL environment, it returns me an<br>
>>>> error.<br>
>>>><br>
>>>> Please have a look at the attached files, which contain a screenshot<br>
>>>> of the error and my script.<br>
>>>><br>
>>>> Please help me as soon as possible.<br>
>>>><br>
>>>> Best wishes,<br>
>>>> --<br>
>>>> S.Rahimi<br>
>>>><br>
>>>><br>
>>>> _______________________________________________<br>
>>>> ncl-talk mailing list<br>
>>>> ncl-talk@ucar.edu<br>
>>>> List instructions, subscriber options, unsubscribe:<br>
>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>>>><br>
>>>><br>
>>> --<br>
>> S.Rahimi<br>
>><br>
>><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> ncl-talk@ucar.edu<br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/262ac813/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/262ac813/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 29 Aug 2017 16:02:50 -0600<br>
From: Mary Haley <haley@ucar.edu><br>
To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] Fw: Can't got value of variables other than<br>
Lat and Lon<br>
Message-ID:<br>
<CACNN_CKZgYJGMZLw9iupKzofhmy8SmKH014DgqNvDe_YUgD2-g@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Xiaoli,<br>
<br>
Just because you have a bunch of values equal to 0.0, it doesn't mean your<br>
data is incorrect.<br>
? You should ?<br>
also look at the min/max of your<br>
?whole ?<br>
data<br>
?variable ?<br>
with:<br>
<br>
printMinMax(data_raw,0)<br>
<br>
Also, to be on the safe side, rather than hard-coding data@_FillValue to<br>
-9999.9, you should use the attribute that is attached to data:<br>
<br>
data@_FillValue = data@CodeMissingValue<br>
?data_raw = data?<br>
<br>
Please do this before the "printMinMax" call.<br>
<br>
?If you still think something is wron?g with your data after looking at the<br>
printMinMax output, then respond back to this message and cc<br>
ncl-talk@ucar.edu. It helps if you can provide more information about why<br>
you think the data is incorrect.<br>
<br>
Also, if you continue to have problems, then it also helps if you can<br>
provide your data file. You can put it on our anonymous ftp:<br>
<br>
<a href="http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP">http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP</a><br>
<br>
?Thanks,<br>
?<br>
<br>
--Mary<br>
<br>
On Tue, Aug 29, 2017 at 3:38 PM, Liu, Xiaoli <xiaoli.liu@ou.edu> wrote:<br>
<br>
><br>
><br>
><br>
> ------------------------------<br>
> *From:* Liu, Xiaoli <xiaoli.liu@ou.edu><br>
> *Sent:* Wednesday, August 30, 2017 6:18 AM<br>
> *To:* ncl-talk@ucar.edu<br>
> *Subject:* Can't got value of variables other than Lat and Lon<br>
><br>
><br>
> Hello,<br>
><br>
> I want use the data of GPM 2A.GPM.DPR.V7 with NCL. I found that I can get<br>
> the value of Lat and lon correctlly. But as to the other variables, the<br>
> value got is very unusal. When I use orbit view, I found there are<br>
> normal values for these variables.<br>
><br>
> Here is my script?<br>
><br>
> dir="/scratch/liu4119/gpm/data/"<br>
> FILES=systemfunc(" ls -1 "+dir+"*.HDF5")<br>
> numFILES = dimsizes(FILES)<br>
> ; print("numFILES = " + numFILES)<br>
> ; print(FILES)<br>
> ; print (" ")<br>
> do i=32,32<br>
> h5_file = addfile(FILES(i),"r")<br>
> ; g = h5_file=>/NS<br>
> ; print(g)<br>
> ; g1=g=>/SLV<br>
> ; print(g1)<br>
> ; list_filevars(g)<br>
> data_raw = h5_file->/NS/SLV/precipRateESurface<br>
> print(data_raw(5396,:))<br>
> data = data_raw<br>
> data@long_name = "precipRate"<br>
> data@_FillValue = -9999.9<br>
> data@units = data_raw@Units<br>
> longitude=h5_file->/NS/Longitude<br>
> printVarSummary(longitude)<br>
> latitude=h5_file->/NS/Latitude<br>
> print(latitude(5396,:))<br>
> print(longitude(5396,:))<br>
> printVarSummary(latitude)<br>
><br>
><br>
><br>
><br>
> Thanks!<br>
><br>
> Xiaoli<br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> ncl-talk@ucar.edu<br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: data_raw.png<br>
Type: image/png<br>
Size: 25095 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0003.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0003.png</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: longitude.png<br>
Type: image/png<br>
Size: 32193 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0004.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0004.png</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: latitude.png<br>
Type: image/png<br>
Size: 33981 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0005.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170829/d43fa209/attachment-0005.png</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Wed, 30 Aug 2017 14:30:12 +0900<br>
From: Dinuka Silva <dinu.kera@gmail.com><br>
To: ncl-talk@ucar.edu<br>
Subject: [ncl-talk] Remove SST data for specific area<br>
Message-ID:<br>
<CAOjtyXkBfSPhsUstxzsMVxtibOvDfR1y0yqtdTGc_0jXXpVhng@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi all,<br>
<br>
I want to remove SST data of one box ({03N - 12N}, {90E- 140E}) and then<br>
find the correlation between remaining SST data and rainfall index. I want<br>
to know how can I remove the SST data of that box. SST data has 3<br>
dimensions. I used mask function but I got error message.<br>
<br>
x = fin2->sst(:,:,:)<br>
<br>
x@lat2d = fin2->latitude<br>
x@lon2d = fin2->longitude<br>
<br>
<br>
latMin = 03<br>
latMax = 12<br>
lonMin = 90<br>
lonMax = 140<br>
<br>
x = mask (x, conform(x, (x@lat2d.ge.latMin .and.x@lat2d.le.latMax<br>
.and.x@lon2d.ge.lonMin .and.x@lon2d.le.lonMax), (/1,2/)), True)<br>
printVarSummary(x)<br>
<br>
<br>
The error message : Dimension size, for dimension number 0, of operands<br>
does not match, can't continue<br>
<br>
<br>
Thanks for your help.<br>
Dinuka<br>
<br>
<br>
<<a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon">https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon</a>><br>
Virus-free.<br>
<a href="http://www.avast.com">www.avast.com</a><br>
<<a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link">https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link</a>><br>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/53da195e/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/53da195e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Wed, 30 Aug 2017 11:09:06 +0430<br>
From: Setareh Rahimi <setareh.rahimi@gmail.com><br>
Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] problem with adding file<br>
Message-ID:<br>
<CAN+DFgzHevxBcc6XZWEu_4kJ3ArpHsTc4tCVSViKEApUoP_ATg@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear all,<br>
<br>
Please have a look at the attached file, which is a screenshot of<br>
ncl_filedump command .<br>
<br>
what is your idea?<br>
<br>
Best wishes,<br>
<br>
On Wed, Aug 30, 2017 at 2:24 AM, Mary Haley <haley@ucar.edu> wrote:<br>
<br>
> As another test, please see if you can look at the file using<br>
> "ncl_filedump" from the UNIX command line. But, first make sure you have<br>
> the file correctly spelled:<br>
><br>
> ls precip.mon.total.v6.nc<br>
><br>
> then:<br>
><br>
> ncl_filedump precip.mon.total.v6.nc<br>
><br>
> ?This is the same thing as the 2-line script that Dennis asked you to try.<br>
> If this works, then you can also verify at the same time if "precip" is a<br>
> variable in that file.<br>
><br>
> If you get any errors from the ncl_filedump command, respond back to<br>
> ncl-talk with all of the error messages included.<br>
><br>
> --Mary<br>
><br>
><br>
><br>
> On Tue, Aug 29, 2017 at 11:45 AM, Dennis Shea <shea@ucar.edu> wrote:<br>
><br>
>> Hello ... Please respond to ncl-talk and not directly to me.<br>
>><br>
>> Your script has:<br>
>><br>
>> *f * = addfile("precip.mon.total.v6.nc","r")<br>
>> prc =* f*->precip({tStrt:tLast},:,:) ; whole grid<br>
>><br>
>> The error message is:<br>
>><br>
>> fatal: file (*f*) isn't defined<br>
>><br>
>> I am not sure what the problem would be. Someone else will have to<br>
>> respond.<br>
>><br>
>> ---<br>
>> What does the following yoeld? Likely the same error.<br>
>><br>
>> *f * = addfile("precip.mon.total.v6.nc","r")<br>
>> print(*f*)<br>
>><br>
>> =====================================================<br>
>> NOTE:<br>
>> You do *not* have to load the following libraries from NCL 6.2.0 onward.<br>
>> So for 6.3.0, they are *automatically loaded.*<br>
>><br>
>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
>><br>
>> From 6.4.0, the following are automatically loaded:<br>
>><br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/bootstrap.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/extval.ncl"<br>
>> "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br>
>><br>
>><br>
>> On Tue, Aug 29, 2017 at 9:27 AM, Setareh Rahimi <setareh.rahimi@gmail.com<br>
>> > wrote:<br>
>><br>
>>> dear Dennis,<br>
>>><br>
>>> Thank you for your advice, but both data file and script are in the<br>
>>> same directory.<br>
>>><br>
>>> Regards,<br>
>>><br>
>>> On Tue 29 Aug 2017 at 19:12, Dennis Shea <shea@ucar.edu> wrote:<br>
>>><br>
>>>> By default, NCL always looks in the 'current directory'.<br>
>>>> My guess is that the data file is not in the same directory as the NCL<br>
>>>> script.<br>
>>>><br>
>>>> dirp = "/some/path/to/file/"<br>
>>>> filp = "precip.mon.total.v6.nc"<br>
>>>><br>
>>>> pthp = dirp + filp<br>
>>>> f = addfile(pthp, "r")<br>
>>>><br>
>>>> or<br>
>>>><br>
>>>> f = addfile(dirp + filp, "r")<br>
>>>><br>
>>>><br>
>>>> On Tue, Aug 29, 2017 at 8:26 AM, Setareh Rahimi <<br>
>>>> setareh.rahimi@gmail.com> wrote:<br>
>>>><br>
>>>>> Dear all,<br>
>>>>><br>
>>>>> Once I put my NCL command line in NCL environment, it returns me an<br>
>>>>> error.<br>
>>>>><br>
>>>>> Please have a look at the attached files, which contain a screenshot<br>
>>>>> of the error and my script.<br>
>>>>><br>
>>>>> Please help me as soon as possible.<br>
>>>>><br>
>>>>> Best wishes,<br>
>>>>> --<br>
>>>>> S.Rahimi<br>
>>>>><br>
>>>>><br>
>>>>> _______________________________________________<br>
>>>>> ncl-talk mailing list<br>
>>>>> ncl-talk@ucar.edu<br>
>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>>>>><br>
>>>>><br>
>>>> --<br>
>>> S.Rahimi<br>
>>><br>
>>><br>
>><br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> ncl-talk@ucar.edu<br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>><br>
>><br>
><br>
<br>
<br>
-- <br>
S.Rahimi<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/9e4394ae/attachment.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/9e4394ae/attachment.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: Screenshot from 2017-08-30 11-00-59.png<br>
Type: image/png<br>
Size: 69507 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/9e4394ae/attachment.png">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/9e4394ae/attachment.png</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
ncl-talk@ucar.edu<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 165, Issue 49<br>
*****************************************<br>
</div>
</span></font></div>
</div>
</div>
</body>
</html>