<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Sh Zhou,<div><br></div><div>On first glance at your program, it looks like you are replacing part of your array with the value of 5000. If you have not specifically set the contour levels in GrADS, then it is choosing the levels for you. In your map labeled "former", the contour interval is 200, while in the "after-changed", the interval is 500 (to accommodate the new value of 5000) . You can choose your contour levels with "set clevs" to show your original contours.</div><div><br></div><div>Colleen</div><div><br></div><div><div>On May 9, 2009, at 11:34 PM, sh zhou wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><br clear="all">Hi All,<br><br>Did someone come across a problem like this: when wrote data back , the display changed unnormal in Grads? because many "contour lines" disappeared , Can you tell me why? Please help me, thank you.<br> <br>In accessory , there are some files of my data and pictures. Please have a look, the code I used as following :<br><br>###<br> program read_DOMAIN_INFO<br> implicit none<br> integer jx,iy,i,j<br>! real yy,xx,z<br> parameter (iy=110,jx=140) <br> real ht(jx,iy)<br> open(10,file='<a href="http://DOMAIN.INFO">DOMAIN.INFO</a>',form='unformatted'<br> $ ,recl=jx*iy*4, access='direct',convert='big_endian') ! replace 4 to 1 if IFC, I use PGI<br> <br> read(10,rec=2) ht <br><br> do i=65,75<br><br> ht(i,30) = 5000.<br> ht(i,55)=5000.<br> enddo<br><br> write(10,rec=2) ht<br><br> print *, 'end of program'<br> close(10)<br> <br> end<br>###<br><br>Sincerely.<br><br><br><br><br> <span><former.jpg></span><span><after-changed.jpg></span><span><DOMAIN.CTL></span><span><DOMAIN.INFO></span></blockquote></div><br></body></html>