Re: Modify variable by longitude

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