<div dir="ltr"><div>Setareh, with simple 1-D coordinate variables, you would easily do this with coordinate range subscripting. Read about that in the NCL ref manual. Here is an example that writes directly into an existing file, rather than making a new file.<br><br> f->MSEBIO_ISOP( {0,10}, {10,30} ) = 0<br></div><div><br></div><div>However, the coordinates in your particular file are 3-D, so you have to use extra code. Here is the simplest way. This is also the slowest, but if you are only working with a small amount of data, it is better to keep the code simple. This example also writes directly into the existing file. Not tested:</div><div><br></div><div> lat1 = 0</div><div> lat2 = 10</div><div> lon1 = 10</div><div> lon2 = 30</div><div><br></div><div> do t = 0, ntimes-1</div><div> do i = 0, nlats-1</div><div> do j = 0, nlons-1</div><div> if ( (XLAT(t,i,j) .ge. lat1) .and. (XLAT(t,i,j) .le. lat2) \</div><div> .and. (XLONG(t,i,j) .ge. lon1) .and. (XLONG(t,i,j) .le. lon2) ) then</div><div> f->MSEBIO_ISOP(t,i,j) = 0</div><div> ! You can add more variables here</div><div> end if</div><div> end do</div><div> end do</div><div> end do</div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 25, 2024 at 1:04 PM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com">setareh.rahimi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Dave,<div><br></div><div>So many thanks for your help. I just checked it and it worked. </div><div><br></div><div>May I ask you please to advise me on how to only convert the values of all variables to zero for a specific region and not whole the domain? (for example lat[0,10], and lon[10,30])?</div><div><br></div><div>Thanks in advance.</div><div>Kind regards,<br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 25, 2024 at 11:18 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>The problem is that Dennis's program did not copy the file global attributes, which contain part of the georeferencing. Just add this line near the top of the program, after the line "print (vNames)":</div><div><br></div><div> copy_VarAtts (f, fout)</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 25, 2024 at 11:07 AM Ehsan Taghizadeh via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><div><div dir="ltr">Dear Setareh,</div><div dir="ltr"><div><div>I suggest you read the variables from the fout file again. For example, if there is a variable named "V10" in fout, you just need to read that variable from fout, as follows:</div><div>V10_zero=fout->V10</div><div><br></div><div>This V10_zero has the same coordinates as V10 in the original file, but its values are zero.</div><div>I hope this is correct.</div></div></div><div dir="ltr"><br></div></div><div><br></div>
</div><div id="m_-6165975761838463980m_8078946574129595812m_-2853437640240134846m_5376276308783386129yahoo_quoted_6819469700">
<div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
<div>
On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:
</div>
<div><br></div>
<div><div id="m_-6165975761838463980m_8078946574129595812m_-2853437640240134846m_5376276308783386129yiv5053542009"><div><div>Dear Dennis,</div><div>I noticed an issue. When I use the attached script it converts all values to zero, but the data is no georefrenced any more. So how can I fix this please?</div><div>Thanks in advance.</div><div>Best regards, <br clear="all"><br clear="all"><div><div dir="ltr">S.Rahimi<br clear="none"><br clear="none"></div></div></div><div><br></div><div id="m_-6165975761838463980m_8078946574129595812m_-2853437640240134846m_5376276308783386129yiv5053542009yqt39142"><div><div><div dir="ltr">On Sun, Jan 21, 2024 at 18:44 Dennis Shea <<a rel="nofollow noopener noreferrer" shape="rect" href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Attached is a script that does what you requested.</div><div><br clear="none"></div></div><br clear="none"><div></div><div><div dir="ltr">On Sat, Jan 20, 2024 at 6:40 AM Setareh Rahimi via ncl-talk <<a rel="nofollow noopener noreferrer" shape="rect" href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br clear="none"></div></div><div><blockquote style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"></blockquote></div><div><blockquote style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br clear="all"><div>Dear NCL users,</div><div><br clear="none"></div><div>I have a NetCDF file (attached) that contains some variables. I would like to change the values of all variables to zero for all grid points. So how can I do this using NCL, please?</div><div><br clear="none"></div><div>This file is an input file for the WRF-Chem model, and I do not want to change the format of the file, only want to change the values to zero.</div><div><br clear="none"></div><div>Please kindly advise me in this regard.</div><div><br clear="none"></div><div>Thanks in advance.</div><div>Kind regards,</div><span>-- </span><br clear="none"><div dir="ltr">S.Rahimi</div></div></blockquote></div></blockquote></div></div></div></div></div></div></div></div></div></blockquote></div></div></blockquote></div></div></div>
</blockquote></div></div>