<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000099'><br><br><hr id="zwchr"><br><div style="color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><style>p { margin: 0; }</style><div style="font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000099"><div>Dear Sir/Madam,</div><div><br></div><div>I wan't to replace default landcover class cropland to evergreenbroadleaf in lis_input.d01.nc over paticular region.</div><div><br></div><div>Landcover variable is 3-D</div><div>-----------------------------------------------------------------------------------------------------------------<br></div><div>Type: float<br>Total Size: 53763120 bytes<br> 13440780 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [sfctypes | 20] x [north_south | 801] x [east_west | 839]<br>Coordinates: <br>Number Of Attributes: 9<br> standard_name : MODIS-IGBP (NCEP-modified) landcover map<br> units : <br> scale_factor : 1<br> add_offset : 0<br> missing_value : -9999<br> vmin : 0<br> vmax : 0<br> num_bins : 20<br> _FillValue : -9999<br>----------------------------------------------------------------------------------------------------------------</div><div><font color="#000000">I created a new file nc having only landcover variable like this</font></div><div>---------------------------------------------------------<br></div><div>f=addfile("lis_input.d01.nc","r")<br>a=f->LANDCOVER(:,:,:)<br>a1=new((/20,801,839/),float)<br>a1=a(:,:,:)<br>do i=341,541 <br>do j=166,491<br>if(a(12,i,j) .eq. 1) then<br>a1(2,i,j)=1<br>a1(12,i,j)=0<br>end if<br>end do<br>end do<br><br>m=addfile("new.nc","c")<br>m->t=a1<br>---------------------------------------------------------------------------</div><div><br></div><div><font color="#000000">Now in order to modify my original lis_input.d01.nc I did</font></div><div><br></div><div>f1=addfile("lis_input.d01.nc","w")<br>f2=addfile("new.nc","r")<br>f1->LANDCOVER(:,:,:)=(/f2->t(:,:,:)/)<br></div><div><br></div><div><font color="#000000">The error is</font></div><div><font color="#000000"> <br></font></div><div>copyright (C) 1995-2015 - All Rights Reserved<br> University Corporation for Atmospheric Research<br> NCAR Command Language Version 6.3.0<br> The use of this software is governed by a License Agreement.<br> See http://www.ncl.ucar.edu/ for more details.<br>Segmentation fault (core dumped)</div><div><br></div><div><font color="#000000">ncl -x command gives</font></div><div>f1=addfile("lis_input.d01.nc","w")<br>+ f2=addfile("new.nc","r")<br>+ f1->LANDCOVER(:,:,:)=(/f2->t(:,:,:)/)<br>Segmentation fault (core dumped)<br><br></div><div><br></div><div><font color="#000000">Why this segmentation fault error is coming and there is approximately 200 GB free space .Kindly help me to solve this issue.</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">If there is any other way to replace landcover class, please explain me. </font></div><div><br></div><div><font color="#000000">Thanks</font>,</div><div><font color="#000000">G.Pavani</font><br></div></div></div><br></div></body></html>