See addfile doc
http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml
Use the "w" option.
I suggest making a copy of the original file ... "just in case"
f = addfile("met_em.d01.1998-03-20_12.nc", "w")
si = f->SEAICE ; SEAICE(Time, south_north, west_east)
lon = f->XLONG_M
lat = f->XLAT_M
si_1d = ndtooned( si )
lon_1d = ndtooned( lon )
lat_1d = ndtooned( lat )
if (lon_1d.ge.-134 .and. lon_1d.le.169) then
si_new = ....
end if
si_new = onedtond(si_1d, dimsizes(si) )
f->SEAICE = (/ si_new /)
Brandon Fisel wrote:
> Hi NCL'ers,
>
> I'm attempting to read in all the variables from a met_em.d01 file and
> write all of them to a new met_em.d01 file with the exception of a
> variable, SEAICE, which will need to be modified before being written to
> the new file.
>
> The variable, SEAICE, contains data between 0.0 - 1.0. The constraints
> are that the variable will not need to be changed between longitudes
> (-135 - -179) and (170 - 179). The negative longitudes are West and
> positive longitudes are East. The remainder longitudes (-134 - 0 & 0 -
> 169) will need to be modified for values between 0-0.49 and 0.50-1.
>
> Attached are the contents of the met_em.d01 file.
>
> -------------------------------------------------------------
> Brandon Fisel
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 21 20:51:02 2009
This archive was generated by hypermail 2.1.8 : Fri Oct 23 2009 - 11:36:30 MDT