[ncl-talk] Save output as a NetCDF or tiff file

[Setareh Rahimi]

Thank you so much for your help. How can I assign lat/lon to the AOD,
please?
I used the below lines in the script but received the following error:
" fatal:(XLAT)
is not a named dimension in variable (a2)".



       system("/bin/rm -f AOD550_2D.nc")    ; remove any pre-existing file

       ncdf = addfile("AOD550_2D.nc" ,"c")  ; open output netCDF file


    ;===================================================================

    ; create global attributes of the file (optional)

    ;===================================================================

       fAtt               = True            ; assign file attributes

       fAtt at title         = "Derived AOD"

       fAtt at creation_date = systemfunc ("date")


       fileattdef( ncdf, fAtt )            ; copy file attributes


    ;===================================================================

    ; make time an UNLIMITED dimension; recommended  for most applications

    ;===================================================================

       filedimdef(ncdf,"time",-1,True)


    ;===================================================================

    ; output variables directly; NCL will call appropriate functions

    ; to write the meta data associated with each variable

    ;===================================================================

       ncdf->AOD550_2D  = AOD550_2D


      ncdf->a2  = a2(0,:,:,:)

      AOD550_2D!0    = "lat"

      AOD550_2D!1    = "lon"

             AOD550_2D&lat  = a2&XLAT
===========================================================================================
here is printVarSummary of "a2"

Variable: a2

Type: float

Total Size: 144434880 bytes

            36108720 values

Number of Dimensions: 4

Dimensions and sizes: [Time | 24] x [bottom_top | 30] x [south_north | 229]
x [west_east | 219]

Coordinates:

            Time: [   0..  23]

Number Of Attributes: 6

  coordinates : XLONG XLAT XTIME

  stagger :

  units : ?

  description : bin 2 layer optical thickness

  MemoryOrder : XYZ

  FieldType : 104



On Fri, Nov 10, 2023 at 8:17 PM Dennis Shea <shea at ucar.edu> wrote:

> netCDF operator  'ncks' to extract a variable from a netCDF file and copy
> it to another
>
> %> ncks -v xxx Model.nc XXX.nc
>
> where 'xxx' is the name of the variable on the model file
> ===============================
>     AOD
>
>        *system* <https://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm -f AOD.nc")    ; remove any pre-existing file
>        ncdf = *addfile* <https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>("AOD.nc" ,"c")  ; open output netCDF file
>
>     ;===================================================================
>     ; create global attributes of the file (optional)
>     ;===================================================================
>        fAtt               = True            ; assign file attributes
>        fAtt at title         = "Derived AOD"
>        fAtt at creation_date = *systemfunc* <https://www.ncl.ucar.edu/Document/Functions/Built-in/systemfunc.shtml> ("date")
>        *fileattdef* <https://www.ncl.ucar.edu/Document/Functions/Built-in/fileattdef.shtml>( ncdf, fAtt )            ; copy file attributes
>
>     ;===================================================================
>     ; make time an UNLIMITED dimension; recommended  for most applications
>     ;===================================================================
>        *filedimdef* <https://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml>(ncdf,"time",-1,True)
>
>     ;===================================================================
>     ; output variables directly; NCL will call appropriate functions
>     ; to write the meta data associated with each variable
>     ;===================================================================
>        ncdf->AOD  = AOD
>
> ================================
> See
> *https://www.ncl.ucar.edu/Applications/method_1.shtml
> <https://www.ncl.ucar.edu/Applications/method_1.shtml>*
>
> On Fri, Nov 10, 2023 at 1:02 AM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear all,
>> I have calculated Aerosol Optical Depth (AOD) from WRF-Chem output using
>> "TAUAER2" and "TAUAER3" variables (script attached). Instead of using the
>> created AOD to plot it, I would like to save the AOD as a NetCDF or tiff
>> file. So how can I do this, please?
>> I have another question. How can I extract one variable from the model
>> output and save it as a NetCDF or tiff file?
>> Please help me with these issues.
>> Kind regards,
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231111/aa2e6186/attachment.htm>