[gradsusr] Writing a netcdf in fortran that's readable by grads

Thomas Robinson ter at hawaii.edu
Thu Jan 15 17:12:15 EST 2015


Thank you for the reply.
I added:
istatus = nf90_def_var(ncid,"lat", nf90_float, ila, ilat )
          write (6,*) trim(nf90_strerror(istatus))
i*status = nf90_put_att(ncid, ilat, 'units', 'degrees_north')*
istatus = nf90_def_var(ncid,"lon", nf90_float, ilo, ilon )
          write (6,*) trim(nf90_strerror(istatus))
*istatus = nf90_put_att(ncid, ilon, 'units', 'degrees_east')*

To give units to the latitude and longitude.  Now I have a new error:

ga-> sdfopen foo.nc
Scanning self-describing file:  foo.nc
SDF file has no discernable time coordinate -- using default values.
gadsdf: SDF file does not have any non-coordinate variables.

What does that mean?

On Thu, Jan 15, 2015 at 11:53 AM, James T. Potemra <jimp at hawaii.edu> wrote:

>  Hi Tom,
>
> I think you'll need to add attributes to your variables, specifically the
> units.  For lat and lon these
> should be "degrees_north" and "degrees_east".
>
> Jim
>
>
> On 1/15/15 11:38 AM, Thomas Robinson wrote:
>
> Aloha Grads Users,
>
>  I am bringing in some data and running calculations using fortran 90.  I
> want to output the data as a netCDF file, and I would like it to be
> readable by grads.  I have no problem creating a file called foo.nc and
> putting my lat, lon, and ref variables in the file.  The problem is if I go
> to grads and type 'sdfopen foo.nc', it says
>  ga-> sdfopen foo.nc
> Scanning self-describing file:  foo.nc
> gadsdf: SDF file has no discernable X coordinate.
>   To open this file with GrADS, use a descriptor file with an XDEF entry.
>   Documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html
>  What do I need to add to make the file self describing?
>   I'm not sure what I need to add to get this to work.  Here is what I
> have so far:
>
>  istatus = nf90_create(path = "foo.nc", cmode = nf90_noclobber, ncid =
> ncid)
>           write (6,*) trim(nf90_strerror(istatus))
>
>  istatus =nf90_def_dim(ncid,'latitude',ny,ila)
>           write (6,*) trim(nf90_strerror(istatus))
> istatus =nf90_def_dim(ncid,'longitude',nx,ilo)
>           write (6,*) trim(nf90_strerror(istatus))
> istatus =nf90_def_dim(ncid,'time',1,itt)
>           write (6,*) trim(nf90_strerror(istatus))
>
>  iids = (/ilo,ila,itt/)
>
>  istatus = nf90_def_var(ncid,"ref", nf90_float, iids, iref )
>           write (6,*) trim(nf90_strerror(istatus))
> istatus = nf90_def_var(ncid,"lat", nf90_float, ila, ilat )
>           write (6,*) trim(nf90_strerror(istatus))
> istatus = nf90_def_var(ncid,"lon", nf90_float, ilo, ilon )
>           write (6,*) trim(nf90_strerror(istatus))
> istatus = nf90_enddef(ncid)
>           write (6,*) trim(nf90_strerror(istatus))
>
>  istatus = nf90_put_var(ncid,iref,ref)
>           write (6,*) trim(nf90_strerror(istatus))
> istatus = nf90_put_var(ncid,ilat,latrad)
>           write (6,*) trim(nf90_strerror(istatus))
> istatus = nf90_put_var(ncid,ilon,lonref)
>            write (6,*) trim(nf90_strerror(istatus))
>
>  istatus = nf90_close(ncid)
>           write (6,*) trim(nf90_strerror(istatus))
>
>  --
>   Tom Robinson
> Treasurer - Graduate Student Organization
> Secretary - Student Athletic Fee Committee
> Graduate Student - Department of Atmospheric Science
>
>
> _______________________________________________
> gradsusr mailing listgradsusr at gradsusr.orghttp://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Tom Robinson
Treasurer - Graduate Student Organization
Secretary - Student Athletic Fee Committee
Graduate Student - Department of Atmospheric Science
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150115/f28d5e21/attachment-0001.html 


More information about the gradsusr mailing list