[ncl-talk] How to convert NCL script to a NetCDF file?

Dennis Shea shea at ucar.edu
Sun Oct 30 17:47:03 MDT 2022


"The [R] software can not recognize the lat/lon. So I need to convert my
NCL script to a NetCDF file in an easily readable way for the R."
====


ncl-talk does not have much experience with R.

The description "The [R] software can not recognize the lat/lon."

Perhaps, the fact the the XLAT and XLONG are three dimensional  (Time,
south_north, west_east)  causes the problem?
I do not think ncl-talk is going to learn R's reading of netCDF files.
***That is your responsibility. ***
======================
NOT TESTED
====================


              load
"$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"


       inNC = "wef_chem.nc"         ; original file
       f = addfile(inNC, "r")


       Times    = f->Times                ; char Times(Time, DateStrLen) ;
    ;
    ; https://www.ncl.ucar.edu/Document/Functions/WRF_contributed/wrf_times_c.shtml
    ;
       *time      = wrf_times_c*( Times, 3 )   ; yyyymmddhh   (integer)
       time!0    = "time"                    ; name the dimension

       time&time = time

       TAUAER2   = f->TAUAER2
       TAUAER2!0 = "time"
       TAUAER2&time = time

       TAUAER3  = f->TAUAER3
       TAUAER3!0 = "time"
       TAUAER3&time = time

       xlat_2d  = f->XLAT(0,:,:)             ; 3D->2D
       xlong_2d = f->XLONG(0,:,:)

       printVarSummaary(time)

       printVarSummaary(xlat_2d)
       printVarSummaary(xlong_2d)

       printVarSummaary(TAUAER2)
       printVarSummaary(TAUAER3)

;-----------------------------

       newNC = "setarah.nc"

       *system*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm
-f + newNC)     ; remove any pre-existing file
       ncOUT = *addfile*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>(newNC
,"c")     ; open OUTPUT netCDF file

    ;===================================================================
    ; create global attributes of the file (optional)
    ;===================================================================
       fAtt               = True            ; assign file attributes
       fAtt at title         = "REWRITE of " + inNC
       fAtt at source_file   =  inNC
     ;;fAtt at Conventions   = "None"
       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>(
ncOUT, fAtt )            ; add new attributes

       ncOUT->time  = time

       ncOUT->XLAT  = xlat_2d

       ncOUT->XLONG = xlong_2d

       ncOUT->TAUAER2 = TAUAER2
       ncOUT->TAUAER3 = TAUAER3


On Sun, Oct 30, 2022 at 11:31 AM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all NCL users,
> Using R software, I want to plot aerosol optical depth (AOD) from the
> WRF-chem model. Once import the model output in the R, the software can not
> recognize the lat/lon. So I need to convert my NCL script to a NetCDF file
> in an easily readable way for the R. How can I convert my NCL script (which
> includes converted time, lat/lon, and calculate AOD) to a new  NetCDF file,
> please? The script and the model output have been attached.
>  wef_chem.nc.gz
> <https://drive.google.com/file/d/1K7lEFyBtr3B88PHF6QWhMRB9yc_ORn3e/view?usp=drive_web>
>
> --
> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221030/774c90ce/attachment.htm>


More information about the ncl-talk mailing list