a little help with reading binary: coordinates from Erik Noble on 2008-07-15 (ncl-talk 2008 archive)

From: Erik Noble <enoble_at_nyahnyahspammersnyahnyah>
Date: Tue, 15 Jul 2008 20:23:07 -0400

Hi. I have a flat binary file that contains 16 records of gridded
precipitation data.
Today was my first time to use NCL to read a binary file.

Although I was successful at reading the data, I am unsure if I am
creating the lat and lon coordinates correctly. When I do a
printVarSummary on both variables in the file, I get:

Variable: cpc_combined_MW
Type: float
Total Size: 2764800 bytes
            691200 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 480] x [lon | 1440]
Coordinates:
            lat: [ 0.. 0]
            lon: [ 0.. 0]
Number Of Attributes: 2
  units : mm/hr
  long_name : merged microwave precipitation

Should I be getting lat: [ 0.. 0] and lon: [ 0.. 0]?

-Erik

;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
     setfileoption ("bin", "ReadByteOrder", "BigEndian")
   do i = 0, 14 ; CMORPH binary file contains 16 records. NCL starts
at counting at 0
     lat = fbindirread ("./20060901_3hr-025deg_cpc+comb",0, 480, "double")
     lon = fbindirread ("./20060901_3hr-025deg_cpc+comb",1, 1440, "double")
     cpc_combined_MW = fbindirread
("./20060901_3hr-025deg_cpc+comb",i,(/480,1440/),"float")
     CMORPH = fbindirread
("./20060901_3hr-025deg_cpc+comb",i+1,(/480,1440/),"float")

    ;=============================
    ; create lat and long coordinate variables
    ;============================
     lon!0 = "lon"
     lon_at_long_name = "lon"
     lon_at_units = "degrees-east"
     lon&lon = lon

     lat!0 = "lat"
     lat_at_long_name = "lat"
     lat_at_units = "degrees_north"
     lat&lat = lat
    ;=============================
    ; name dimensions of variables
    ;============================
     cpc_combined_MW!0 = "lat"
     cpc_combined_MW!1 = "lon"
     cpc_combined_MW&lat = lat
     cpc_combined_MW&lon = lon
     cpc_combined_MW_at_long_name = "merged microwave precipitation"
     cpc_combined_MW_at_units = "mm/hr"
     print(i)
     printVarSummary(cpc_combined_MW)

     CMORPH!0 = "lat"
     CMORPH!1 = "lon"
     CMORPH&lat = lat
     CMORPH&lon = lon
     CMORPH_at_long_name = "CMORPH"
     CMORPH_at_units = "mm/hr"
     print(i+1)
     printVarSummary(CMORPH)

end do
end

ENOBLE_at_ATHENA: /usr/people/enoble> ncl bin3.ncl
Copyright (C) 1995-2007 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 4.3.1
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 0

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 1

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 2

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 3

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 4

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 5

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 6

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 7

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 8

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 9

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 10

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 11

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 12

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 13

Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 14

Variable: CMORPH
Type: float
Total Size: 2764800 bytes
            691200 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 480] x [lon | 1440]
Coordinates:
            lat: [ 0.. 0]
            lon: [ 0.. 0]
Number Of Attributes: 2
  units : mm/hr
  long_name : CMORPH
ENOBLE_at_ATHENA: /usr/people/enoble>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 15 2008 - 18:23:07 MDT

This archive was generated by hypermail 2.2.0 : Fri Jul 18 2008 - 08:51:52 MDT