Re: NCL webpage problem

Dear Muhammad,

You didn't include any code, so I'm not sure exactly what you need to do with the lat/lon arrays on the file. It sounds like you want to plot them, which I already showed you what you needed to do below.

When you set up your resource list for plotting, you will need to do this:

   sst@lon2d = f->lon
   sst@lat2d = f->lat

Where "f" needs to be replaced with whatever variable you used when you called "addfile" and "sst" needs to be replaced with whatever your SST variable is called.

You might also need to set:

> res@gsnAddCyclic = False

If this doesn't help, then please include your script and a clear description of what is not working, or what your exact error messages are, if any. Even better, if you can include your data so we can run the script, that helps. But, please try what I recommended before you send me any scripts or data.

Thanks,

--Mary

On Nov 20, 2012, at 7:06 PM, Muhammad Afzaal wrote:

> Dear Mary!
>
> Thank you very much for reply.
> The file has lat lon coordinates in the form lat(j,i), lon(j,i).
> Following is the header information of data file. Please advise in the light of this information.
>
> Thanking you in anticipation!
>
> P.S: problem of NCL page loading is solved. Thanks again.
>
> ==========================================
>
> dimensions:
> time = 120 // unlimited
> j = 292
> i = 362
> bnds = 2
> vertices = 4
> variables:
> double time ( time )
> bounds : time_bnds
> units : days since 1850-01-01 00:00:00
> calendar : gregorian
> axis : T
> long_name : time
> standard_name : time
>
> double time_bnds ( time, bnds )
>
> integer j ( j )
> units : 1
> long_name : cell index along second dimension
>
> integer i ( i )
> units : 1
> long_name : cell index along first dimension
>
> float lat ( j, i )
> standard_name : latitude
> long_name : latitude coordinate
> units : degrees_north
> bounds : lat_vertices
>
> float lon ( j, i )
> standard_name : longitude
> long_name : longitude coordinate
> units : degrees_east
> bounds : lon_vertices
>
> float lat_vertices ( j, i, vertices )
> units : degrees_north
>
> float lon_vertices ( j, i, vertices )
> units : degrees_east
>
> float tos ( time, j, i )
> standard_name : sea_surface_temperature
> long_name : Sea Surface Temperature
> comment : "this may differ from ""surface temperature"" in regions of sea ice."
> units : K
> original_name : tos
> original_units : degC
> history : 2011-10-06T23:34:59Z altered by CMOR: Converted units from 'degC' to 'K'.
> cell_methods : time: mean
> cell_measures : area: areacello
> missing_value : 1e+20
> _FillValue : 1e+20
>
> ====================================================
>
> MUHAMMAD AFZAAL KARORI
> Institute of Atmospheric Physics
> Chinese Academy of Sciences
> Beijing, China
>
>
> On Tue, Nov 20, 2012 at 11:37 PM, Mary Haley <haley@ucar.edu> wrote:
>
> On Nov 18, 2012, at 5:05 AM, Muhammad Afzaal wrote:
>
> > Dear NCL team!
> >
> > 1. Thanks for your efforts for always providing a prompt and valueable help on NCL issues. However, It has been noticed that the webpage http://www.ncl.ucar.edu/ takes very very long time to open, while the internet speed ok. is it a problem at my end or at yours? Please see into the matter.
>
> Hi Muhammad,
>
> I haven't had any recent reports of people dealing with slow times on the NCL website. We had some problems last week with power, and our websites were behaving erratically. We hope to have this fixed by now.
>
> >
> > 2. I have SST data with dimensions (time, j, i) instead of (time, lat, lon). How can I assign lat lon to i j. I tried lat2d lon2d function, but the data is not placed properly in the display. Any help will be highly appreciated.
>
> It's possible the file that contains your SST data also contains lat/lon information. If you look at your SST variable:
>
> printVarSummary(sst)
>
> Does it have a "coordinates" attribute? If so, does the coordinates attribute say something like "lat lon" or "latitude longitude"? If it does, then this likely means that you have variables "lat" and "lon" (or "latitude" "longitude") on the file, and you can read these just like you did sst:
>
> lat = f->lat
> lon = f->lon
>
> Once you have these, then you can use them for plotting:
>
> sst@lat2d = lat
> sst@lon2d = lon
>
> You may also need to set (for plotting):
>
> res@gsnAddCyclic = False
>
> if your data is not global.
>
> --Mary
>
> >
> > Thanks a lot!
> >
> > MUHAMMAD AFZAAL KARORI
> > Institute of Atmospheric Physics
> > Chinese Academy of Sciences
> > Beijing, China
> > _______________________________________________
> > 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 21 10:59:26 2012