SCALE-RM
Classes | Macros | Functions
read_toshiba_mpr.h File Reference
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes

struct  mppawr_header
 

Macros

#define RDIM   800 /* max number of Range bins */
 
#define AZDIM   320 /* max number of AZ angles */
 
#define ELDIM   120 /* max number of EL angles */
 
#define NFNAME   128 /* number of characters in file names */
 
#define DMISS   -327.68 /* missing data in output (& input data offset) */
 
#define DNOISE   -327.00 /* noise level data in output */
 
#define RAD   (3.141592653589/180.0)
 
#define LEN_BUFRDHD   512
 
#define LEN_BUFPBHD   10368
 
#define LEN_BUFELHD   432
 
#define XYRESO   0.25 /* horizontal resolution in km */
 
#define ZRESO   0.25 /* vertical resolution in km */
 
#define CXDIM   481 /* X-dim of CAPPI output (include center) */
 
#define CYDIM   481 /* Y-dim of CAPPI output (include center) */
 
#define CZDIM   65 /* Z-dim of CAPPI output (include surface) */
 
#define CXKM0   -60.0 /* start X coords from the radar (in km) */
 
#define CYKM0   -60.0 /* start Y coords from the radar (in km) */
 
#define CZKM0   0.0 /* start Z coords from the radar (in km) */
 
#define PXDIM   481 /* X-dim of PPI output (include center) */
 
#define PYDIM   481 /* Y-dim of PPI output (include center) */
 
#define PELDIM   20 /* EL-dim of output EL number (< ELDIM) */
 
#define PXKM0   -60.0 /* start X coords from the radar (in km) */
 
#define PYKM0   -60.0 /* start Y coords from the radar (in km) */
 
#define RXDIM   241 /* X-dim of RHI output (include center) */
 
#define RZDIM   65 /* Z-dim of RHI output (include surface) */
 
#define RAZDIM   300 /* AZ-dim of RHI output (< AZDIM) */
 
#define RZKM0   0.0 /* start Z coords from the radar (in km) */
 

Functions

int16_t char2int16 (void *input)
 
uint16_t char2uint16 (void *input)
 
int32_t char2int32 (void *input)
 
int read_toshiba_mpr (char *in_file, int opt_verbose, mppawr_header *hd, float az[ELDIM][AZDIM], float el[ELDIM][AZDIM], float rtdat[ELDIM][AZDIM][RDIM])
 
int decode_toshiba_mpr (size_t bufsize, unsigned char *buf, int opt_verbose, mppawr_header *hd, float az[ELDIM][AZDIM], float el[ELDIM][AZDIM], float rtdat[ELDIM][AZDIM][RDIM])
 
int ggrads_ctl (char *c_name, char *g_name, int nx, int ny, int nz, int nt, double xmin, double ymin, double zmin, double tmin, double dx, double dy, double dz, double dt, double no_data, int ihour, int imin, int iday, int imonth, int iyear, char *c_var, char *c_memo)
 

Macro Definition Documentation

◆ RDIM

#define RDIM   800 /* max number of Range bins */

Definition at line 5 of file read_toshiba_mpr.h.

◆ AZDIM

#define AZDIM   320 /* max number of AZ angles */

Definition at line 6 of file read_toshiba_mpr.h.

◆ ELDIM

#define ELDIM   120 /* max number of EL angles */

Definition at line 7 of file read_toshiba_mpr.h.

◆ NFNAME

#define NFNAME   128 /* number of characters in file names */

Definition at line 8 of file read_toshiba_mpr.h.

◆ DMISS

#define DMISS   -327.68 /* missing data in output (& input data offset) */

Definition at line 9 of file read_toshiba_mpr.h.

◆ DNOISE

#define DNOISE   -327.00 /* noise level data in output */

Definition at line 10 of file read_toshiba_mpr.h.

◆ RAD

#define RAD   (3.141592653589/180.0)

Definition at line 12 of file read_toshiba_mpr.h.

◆ LEN_BUFRDHD

#define LEN_BUFRDHD   512

Definition at line 15 of file read_toshiba_mpr.h.

◆ LEN_BUFPBHD

#define LEN_BUFPBHD   10368

Definition at line 16 of file read_toshiba_mpr.h.

◆ LEN_BUFELHD

#define LEN_BUFELHD   432

Definition at line 17 of file read_toshiba_mpr.h.

◆ XYRESO

#define XYRESO   0.25 /* horizontal resolution in km */

Definition at line 20 of file read_toshiba_mpr.h.

◆ ZRESO

#define ZRESO   0.25 /* vertical resolution in km */

Definition at line 21 of file read_toshiba_mpr.h.

◆ CXDIM

#define CXDIM   481 /* X-dim of CAPPI output (include center) */

Definition at line 23 of file read_toshiba_mpr.h.

◆ CYDIM

#define CYDIM   481 /* Y-dim of CAPPI output (include center) */

Definition at line 24 of file read_toshiba_mpr.h.

◆ CZDIM

#define CZDIM   65 /* Z-dim of CAPPI output (include surface) */

Definition at line 25 of file read_toshiba_mpr.h.

◆ CXKM0

#define CXKM0   -60.0 /* start X coords from the radar (in km) */

Definition at line 26 of file read_toshiba_mpr.h.

◆ CYKM0

#define CYKM0   -60.0 /* start Y coords from the radar (in km) */

Definition at line 27 of file read_toshiba_mpr.h.

◆ CZKM0

#define CZKM0   0.0 /* start Z coords from the radar (in km) */

Definition at line 28 of file read_toshiba_mpr.h.

◆ PXDIM

#define PXDIM   481 /* X-dim of PPI output (include center) */

Definition at line 30 of file read_toshiba_mpr.h.

◆ PYDIM

#define PYDIM   481 /* Y-dim of PPI output (include center) */

Definition at line 31 of file read_toshiba_mpr.h.

◆ PELDIM

#define PELDIM   20 /* EL-dim of output EL number (< ELDIM) */

Definition at line 32 of file read_toshiba_mpr.h.

◆ PXKM0

#define PXKM0   -60.0 /* start X coords from the radar (in km) */

Definition at line 33 of file read_toshiba_mpr.h.

◆ PYKM0

#define PYKM0   -60.0 /* start Y coords from the radar (in km) */

Definition at line 34 of file read_toshiba_mpr.h.

◆ RXDIM

#define RXDIM   241 /* X-dim of RHI output (include center) */

Definition at line 36 of file read_toshiba_mpr.h.

◆ RZDIM

#define RZDIM   65 /* Z-dim of RHI output (include surface) */

Definition at line 37 of file read_toshiba_mpr.h.

◆ RAZDIM

#define RAZDIM   300 /* AZ-dim of RHI output (< AZDIM) */

Definition at line 38 of file read_toshiba_mpr.h.

◆ RZKM0

#define RZKM0   0.0 /* start Z coords from the radar (in km) */

Definition at line 39 of file read_toshiba_mpr.h.

Function Documentation

◆ char2int16()

int16_t char2int16 ( void *  input)

Definition at line 22 of file read_toshiba.c.

23 {
24  int16_t output;
25 
26 #if __BYTE_ORDER == __LITTLE_ENDIAN
27  output = *((int16_t*)input);
28 #elif __BYTE_ORDER == __BIG_ENDIAN
29  output = ((int16_t)(((unsigned char*)input)[1]) << 8) | (int16_t)(((unsigned char*)input)[0]);
30 #else
31  output = ((char*)input)[1] * 0x100 + ((unsigned char*)input)[0];
32 #endif
33  return output;
34 }

◆ char2uint16()

uint16_t char2uint16 ( void *  input)

Definition at line 36 of file read_toshiba.c.

37 {
38  uint16_t output;
39 
40 #if __BYTE_ORDER == __LITTLE_ENDIAN
41  output = *((uint16_t*)input);
42 #elif __BYTE_ORDER == __BIG_ENDIAN
43  output = ((uint16_t)(((unsigned char*)input)[1]) << 8) | (uint16_t)(((unsigned char*)input)[0]);
44 #else
45  output = ((unsigned char*)input)[1] * 0x100 + ((unsigned char*)input)[0];
46 #endif
47  return output;
48 }

◆ char2int32()

int32_t char2int32 ( void *  input)

Definition at line 50 of file read_toshiba.c.

51 {
52  int32_t output;
53 #if __BYTE_ORDER == __LITTLE_ENDIAN
54  output = *((int32_t*)input);
55 #elif __BYTE_ORDER == __BIG_ENDIAN
56  output = ((int32_t)(((unsigned char*)input)[3]) << 24) | ((int32_t)(((unsigned char*)input)[2]) << 16) | ((int32_t)(((unsigned char*)input)[1]) << 8) | (int32_t)(((unsigned char*)input)[0]);
57 #else
58  output = ((char*)input)[3] * 0x1000000 + ((unsigned char*)input)[2] * 0x10000 + ((unsigned char*)input)[1] * 0x100 + ((unsigned char*)input)[0];
59 #endif
60  return output;
61 }

◆ read_toshiba_mpr()

int read_toshiba_mpr ( char *  in_file,
int  opt_verbose,
mppawr_header *  hd,
float  az[ELDIM][AZDIM],
float  el[ELDIM][AZDIM],
float  rtdat[ELDIM][AZDIM][RDIM] 
)

Definition at line 109 of file read_toshiba_mpr.c.

115 {
116  const size_t bufsize = AZDIM * (ELDIM * (RDIM * 2 + LEN_BUFELHD) + LEN_BUFPBHD) + LEN_BUFRDHD; // max size
117  size_t bsize; // actual data size
118  int ierr;
119  unsigned char *buf;
120  FILE *fp;
121  char *is_gzip;
122 
123  buf = malloc(bufsize);
124  if(buf == NULL){
125  printf("failed to allocate memory in read_toshiba");
126  return -99;
127  }
128  if((fp = fopen(in_file, "r")) == NULL){
129  is_gzip = "true";
130  if((fp = fopen(strcat(in_file,".gz"), "r")) == NULL){
131  //printf("file not found: %s\n", in_file);
132  return -9;
133  }
134  }
135  bsize = fread(buf, 1, bufsize, fp);
136  if(bsize == 0){
137  printf("file size is 0: %s\n", in_file);
138  return -9;
139  }
140 
141  if(is_gzip == "true") bsize = ungzip_toshiba_mpr(bufsize, bsize, buf);
142 
143  ierr = decode_toshiba_mpr(bsize, buf, opt_verbose, hd, az, el, rtdat);
144  if(ierr != 0) return ierr;
145 
146  free(buf);
147 
148  return 0;
149 }

References AZDIM, ELDIM, LEN_BUFELHD, LEN_BUFPBHD, LEN_BUFRDHD, and RDIM.

Referenced by scale_letkf::letkf_add_inflation_setup().

Here is the caller graph for this function:

◆ decode_toshiba_mpr()

int decode_toshiba_mpr ( size_t  bufsize,
unsigned char *  buf,
int  opt_verbose,
mppawr_header *  hd,
float  az[ELDIM][AZDIM],
float  el[ELDIM][AZDIM],
float  rtdat[ELDIM][AZDIM][RDIM] 
)

Definition at line 152 of file read_toshiba_mpr.c.

158 {
159  unsigned char *bufrdhd; /* RAW DATA hdr */
160  unsigned char *bufpbhd; /* POLAR BLOCK hdr */
161  unsigned char *bufelhd; /* EL hdr */
162  unsigned char *bufdat; /* 2 byte x 800 range (max) */
163  int i, j, k, iss, i0, i1, i2, i3;
164  int aznum, elnum, rnum, sp_rnum, lp_rnum;
165  int r_byte, irb, mesh_size, data_size;
166  int ideg, ideg_tc, inum, ibunshi, ibunbo;
167  int s1_yr, s1_mn, s1_dy, s1_hr, s1_mi, s1_sc;
168  int s2_yr, s2_mn, s2_dy, s2_hr, s2_mi, s2_sc;
169  int s3_yr, s3_mn, s3_dy, s3_hr, s3_mi, s3_sc;
170  int el_num_sec, cpi_num, sec_num, rnum1;
171  int hit_num, spa_num, lpa_num, spv_num, lpv_num;
172  int pw_hsp, pw_hlp, pw_vsp, pw_vlp, prf;
173  float mesh_lsb, mesh_offset;
174  float start_az_sec, end_az_sec, start_el_sec, end_el_sec;
175  float start_az_cpi, end_az_cpi, start_el_cpi, end_el_cpi;
176  float lower_el_angle, upper_el_angle;
177  float tx_freq_h, tx_freq_v, tx_power_h, tx_power_v;
178  float tx_plen_hs, tx_plen_hl, tx_plen_vs, tx_plen_vl;
179  float tx_antgain_h, tx_antgain_v, rx_antgain_h, rx_antgain_v;
180  float hori_beamwidth_h, hori_beamwidth_v, vert_beamwidth_h, vert_beamwidth_v;
181  float start_obs_az, start_obs_el, end_obs_az, end_obs_el;
182  // added by Otsuka
183  size_t buf_offset;
184  const float ideg2deg = 360.0 / 65536.0;
185  const float ideg_tc2deg = 180.0 / 32768.0;
186  int tmp_range_res;
187 
188  // init param
189  hd->range_res = -1;
190 
191  /*----- open RAW DATA file -----*/
192  buf_offset = 0;
193 
194  /*----- read RAW DATA header (512 byte)-----*/
195  if(buf_offset == bufsize) {
196  goto READEND;
197  } else if(buf_offset + LEN_BUFRDHD > bufsize) {
198  printf("# abnormal record length (RDhd)\n");
199  goto READERR;
200  }
201  bufrdhd = buf + buf_offset;
202  buf_offset += LEN_BUFRDHD;
203 
204  strncpy(hd->data_name, &bufrdhd[16], 32);
205  hd->s_yr = 1000 * (bufrdhd[48] >> 4) + 100 * (bufrdhd[48] & 15) + 10 * (bufrdhd[49] >> 4) + (bufrdhd[49] & 15);
206  hd->s_mn = 10 * (bufrdhd[50] >> 4) + (bufrdhd[50] & 15);
207  hd->s_dy = 10 * (bufrdhd[51] >> 4) + (bufrdhd[51] & 15);
208  hd->s_hr = 10 * (bufrdhd[52] >> 4) + (bufrdhd[52] & 15);
209  hd->s_mi = 10 * (bufrdhd[53] >> 4) + (bufrdhd[53] & 15);
210  hd->s_sc = 10 * (bufrdhd[54] >> 4) + (bufrdhd[54] & 15);
211  hd->e_yr = 1000 * (bufrdhd[56] >> 4) + 100 * (bufrdhd[56] & 15) + 10 * (bufrdhd[57] >> 4) + (bufrdhd[57] & 15);
212  hd->e_mn = 10 * (bufrdhd[58] >> 4) + (bufrdhd[58] & 15);
213  hd->e_dy = 10 * (bufrdhd[59] >> 4) + (bufrdhd[59] & 15);
214  hd->e_hr = 10 * (bufrdhd[60] >> 4) + (bufrdhd[60] & 15);
215  hd->e_mi = 10 * (bufrdhd[61] >> 4) + (bufrdhd[61] & 15);
216  hd->e_sc = 10 * (bufrdhd[62] >> 4) + (bufrdhd[62] & 15);
217  /*--- Scan-start(s1), Sequence-start(s2), Mode-start(s3) time-- */
218  s1_yr = 1000 * (bufrdhd[64] >> 4) + 100 * (bufrdhd[64] & 15) + 10 * (bufrdhd[65] >> 4) + (bufrdhd[65] & 15);
219  s1_mn = 10 * (bufrdhd[66] >> 4) + (bufrdhd[66] & 15);
220  s1_dy = 10 * (bufrdhd[67] >> 4) + (bufrdhd[67] & 15);
221  s1_hr = 10 * (bufrdhd[68] >> 4) + (bufrdhd[68] & 15);
222  s1_mi = 10 * (bufrdhd[69] >> 4) + (bufrdhd[69] & 15);
223  s1_sc = 10 * (bufrdhd[70] >> 4) + (bufrdhd[70] & 15);
224  s2_yr = 1000 * (bufrdhd[72] >> 4) + 100 * (bufrdhd[72] & 15) + 10 * (bufrdhd[73] >> 4) + (bufrdhd[73] & 15);
225  s2_mn = 10 * (bufrdhd[74] >> 4) + (bufrdhd[74] & 15);
226  s2_dy = 10 * (bufrdhd[75] >> 4) + (bufrdhd[75] & 15);
227  s2_hr = 10 * (bufrdhd[76] >> 4) + (bufrdhd[76] & 15);
228  s2_mi = 10 * (bufrdhd[77] >> 4) + (bufrdhd[77] & 15);
229  s2_sc = 10 * (bufrdhd[78] >> 4) + (bufrdhd[78] & 15);
230  s3_yr = 1000 * (bufrdhd[80] >> 4) + 100 * (bufrdhd[80] & 15) + 10 * (bufrdhd[81] >> 4) + (bufrdhd[81] & 15);
231  s3_mn = 10 * (bufrdhd[82] >> 4) + (bufrdhd[82] & 15);
232  s3_dy = 10 * (bufrdhd[83] >> 4) + (bufrdhd[83] & 15);
233  s3_hr = 10 * (bufrdhd[84] >> 4) + (bufrdhd[84] & 15);
234  s3_mi = 10 * (bufrdhd[85] >> 4) + (bufrdhd[85] & 15);
235  s3_sc = 10 * (bufrdhd[86] >> 4) + (bufrdhd[86] & 15);
236  strncpy(hd->site_name, &bufrdhd[136], 32);
237  hd->latitude = 0.000001 * char2int32(bufrdhd + 168 + 0);
238  hd->longitude = 0.000001 * char2int32(bufrdhd + 172 + 0);
239  hd->altitude = 0.01 * char2int32(bufrdhd + 180 + 0);
240  strncpy(hd->sq_name, &bufrdhd[216], 16);
241 
242  hd->el_num = bufrdhd[236 + 1];
243  hd->ray_num = char2uint16(bufrdhd + 238); //unsigned
244  hd->start_az = char2uint16(bufrdhd + 240) * ideg2deg; //unsigned
245  hd->start_el = char2int16(bufrdhd + 242) * ideg_tc2deg;
246  hd->end_az = char2int16(bufrdhd + 244) * ideg2deg; //unsigned
247  hd->end_el = char2int16(bufrdhd + 246) * ideg_tc2deg;
248  mesh_size = char2int16(bufrdhd + 384) * 8; /* mesh_size stored in byte? */
249  ibunshi = char2int32(bufrdhd + 388);
250  ibunbo = char2int32(bufrdhd + 392);
251  mesh_lsb = (float)ibunshi / (float)ibunbo;
252  ibunshi = char2int32(bufrdhd + 396);
253  ibunbo = char2int32(bufrdhd + 400);
254  mesh_offset = (float)ibunshi / (float)ibunbo;
255  hd->range_num = RDIM; /* <== shoud be fixed num, becasue of rtdat's base dim */
256 
257  hd->mesh_offset = mesh_offset; //save offset
258 
259  /*--- print RAW DATA header ---*/
260  if (opt_verbose >= 1){
261  //printf("## input file: %s\n", in_file);
262  printf("## data name: %s\n", hd->data_name);
263  printf("## scan collection start/end time: ");
264  printf("%d/%d/%d, %d:%d:%d - %d/%d/%d, %d:%d:%d\n",
265  hd->s_yr, hd->s_mn, hd->s_dy, hd->s_hr, hd->s_mi, hd->s_sc,
266  hd->e_yr, hd->e_mn, hd->e_dy, hd->e_hr, hd->e_mi, hd->e_sc);
267  printf("## scan start time (s1): %d/%d/%d, %d:%d:%d\n",
268  s1_yr, s1_mn, s1_dy, s1_hr, s1_mi, s1_sc);
269  printf("## sequence start time (s2): %d/%d/%d, %d:%d:%d\n",
270  s2_yr, s2_mn, s2_dy, s2_hr, s2_mi, s2_sc);
271  printf("## mode start time (s3): %d/%d/%d, %d:%d:%d\n",
272  s3_yr, s3_mn, s3_dy, s3_hr, s3_mi, s3_sc);
273  printf("## site name: %s", hd->site_name);
274  printf(" lat=%7.4f, lon=%8.4f, alt=%6.2f\n",
275  hd->latitude, hd->longitude, hd->altitude);
276  printf("## sequence name: %s", hd->sq_name);
277  printf(" el_num=%d, ray_num=%d\n", hd->el_num, hd->ray_num);
278  /*---
279  printf("## start/end AZ angle=%5.1f-%5.1f ",hd->start_az, hd->end_az);
280  printf(" start/end EL angle=%5.2f-%5.2f\n",hd->start_el, hd->end_el);
281  ---*/
282  printf("## mesh size=%d, LSB=%7.3f, offset=%7.3f\n",
283  mesh_size, mesh_lsb, mesh_offset);
284  }
285 
286  if(mesh_size != 16){ /* 2-byte data */
287  printf("# abnormal mesh size: %d bit\n", mesh_size);
288  printf("# this program (read_toshiba_mpr) is not supported !\n");
289  goto READERR;
290  }
291 
292  /*----- RAY data (PPI) loop -----*/
293  aznum = hd->ray_num;
294  if(aznum > AZDIM){
295  printf("# abnormal AZ dimension: %d (> %d)\n", aznum, AZDIM);
296  goto READERR;
297  }
298  for(j = 0; j < aznum; j++) { /*= AZ loop uses j =*/
299  /*----- read POLAR BLOCK (RAY) header (10368 byte)-----*/
300  if(buf_offset == bufsize) {
301  goto READEND;
302  } else if(buf_offset + LEN_BUFPBHD > bufsize) {
303  printf("# abnormal record length (PB header)\n");
304  goto READERR;
305  }
306  bufpbhd = buf + buf_offset;
307  buf_offset += LEN_BUFPBHD;
308 
309  /*--- get PB header info for print ---*/
310  if (opt_verbose >= 2){
311  el_num_sec = bufpbhd[9]; //unsigned
312  cpi_num = char2uint16(bufpbhd + 12); //unsigned
313  sec_num = char2uint16(bufpbhd + 14); //unsigned
314  start_az_sec = char2uint16(bufpbhd + 16) * ideg2deg; //unsigned
315  start_el_sec = char2int16(bufpbhd + 18) * ideg_tc2deg;
316  end_az_sec = char2uint16(bufpbhd + 20) * ideg2deg; //unsigned
317  end_el_sec = char2int16(bufpbhd + 22) * ideg_tc2deg;
318  lower_el_angle = char2int16(bufpbhd + 26) * ideg_tc2deg;
319  upper_el_angle = char2int16(bufpbhd + 30) * ideg_tc2deg;
320  if(lower_el_angle > 180.0) lower_el_angle -= 360.0;
321 
322  /*--- print PB header ---*/
323  printf("## Polar Block (RAY) header: j=%d\n", j);
324  printf(" el_num_sec=%d, cpi_num=%d, sec_num=%d | rnum=%d\n",
325  el_num_sec, cpi_num, sec_num, rnum);
326  printf(" start/end AZ in sector=%6.2f-%6.2f ",
327  start_az_sec, end_az_sec);
328  printf(" start/end EL of pedestal=%6.2f-%6.2f ",
329  start_el_sec, end_el_sec);
330  printf(" lower/upper EL angle=%6.2f-%6.2f\n",
331  lower_el_angle, upper_el_angle);
332  }
333 
334  /*----- EL data (RHI) loop -----*/
335  elnum = hd->el_num;
336  if(elnum > ELDIM){
337  printf("# abnormal EL dimension: %d (> %d)\n", elnum, ELDIM);
338  goto READERR;
339  }
340  for(k = 0; k < elnum; k++) { /*= EL loop uses k =*/
341 
342  /*----- read EL header (432 byte)-----*/
343  if(buf_offset == bufsize) {
344  goto READEND;
345  } else if(buf_offset + LEN_BUFELHD > bufsize) {
346  printf("# abnormal record length (EL header)\n");
347  goto READERR;
348  }
349  bufelhd = buf + buf_offset;
350  buf_offset += LEN_BUFELHD;
351 
352  /*** if (opt_verbose >= 3 && (j==0 || j==1) && (k==3 || k==4)){ ***/
353  if (opt_verbose >= 3){
354 
355  /*--- get EL header info for print ---*/
356  data_size = char2uint16(bufelhd + 0); //unsigned
357  hit_num = char2uint16(bufelhd + 8); //unsigned
358  cpi_num = char2uint16(bufelhd + 10); //unsigned
359  spa_num = char2uint16(bufelhd + 12); /*short pulse align range num*/ //unsigned
360  lpa_num = char2uint16(bufelhd + 14); /*long pulse align range num*/ //unsigned
361  spv_num = char2uint16(bufelhd + 16); /*short pulse valid range num*/ //unsigned
362  lpv_num = char2uint16(bufelhd + 18); /*long pulse valid range num*/ //unsigned
363  start_az_cpi = char2uint16(bufelhd + 128) * ideg2deg; //unsigned
364  start_el_cpi = char2int16(bufelhd + 130) * ideg_tc2deg;
365  end_az_cpi = char2uint16(bufelhd + 132) * ideg2deg; //unsigned
366  end_el_cpi = char2int16(bufelhd + 134) * ideg_tc2deg;
367  /*short pulse: 0=0.5us, 1=1.0us, 2=2.0us, 3=4.0us*/
368  /*long pulse: 0=24us, 1=32us, 2=48us, 3=72us*/
369  pw_hsp = bufelhd[144]; /*pulse width of H-pol short pulse*/
370  pw_hlp = bufelhd[145]; /*pulse width of H-pol long pulse*/
371  pw_vsp = bufelhd[146]; /*pulse width of V-pol short pulse*/
372  pw_vlp = bufelhd[147]; /*pulse width of V-pol long pulse*/
373  prf = char2uint16(bufelhd + 148); //unsigned
374  tx_freq_h = 0.0001 * char2int32(bufelhd + 176);
375  tx_freq_v = 0.0001 * char2int32(bufelhd + 180);
376  tx_power_h = 0.0001 * char2int32(bufelhd + 184);
377  tx_power_v = 0.0001 * char2int32(bufelhd + 188);
378  tx_plen_hs = 0.0001 * char2int32(bufelhd + 192);
379  tx_plen_hl = 0.0001 * char2int32(bufelhd + 196);
380  tx_plen_vs = 0.0001 * char2int32(bufelhd + 200);
381  tx_plen_vl = 0.0001 * char2int32(bufelhd + 204);
382  tx_antgain_h = 0.0001 * char2int32(bufelhd + 208);
383  tx_antgain_v = 0.0001 * char2int32(bufelhd + 212);
384  rx_antgain_h = 0.0001 * char2int32(bufelhd + 216);
385  rx_antgain_v = 0.0001 * char2int32(bufelhd + 220);
386  hori_beamwidth_h = 0.0001 * char2int32(bufelhd + 224);
387  hori_beamwidth_v = 0.0001 * char2int32(bufelhd + 228);
388  vert_beamwidth_h = 0.0001 * char2int32(bufelhd + 232);
389  vert_beamwidth_v = 0.0001 * char2int32(bufelhd + 236);
390 
391  /*--- print EL header ---*/
392  printf("## EL header (%d byte): j=%d k=%d | ", data_size, j, k);
393  printf(" hit_num=%d, cpi_num=%d, prf=%d,",
394  hit_num, cpi_num, prf);
395  printf(" spv/a=%d/%d, lpv/a=%d/%d\n",
396  spv_num, spa_num, lpv_num, lpa_num);
397 
398  if (opt_verbose >= 4){
399  printf(" start/end AZ in cpi=%6.2f-%6.2f\n",
400  start_az_sec, end_az_sec);
401  printf(" start/end EL in cpi=%6.2f-%6.2f\n",
402  start_el_sec, end_el_sec);
403  printf(" < pw_short pulse: 0=0.5us, 1=1.0us, 2=2.0us, 3=4.0us >\n");
404  printf(" < pw_long pulse: 0=24us, 1=32us, 2=48us, 3=72us >\n");
405  printf(" pw_hsp=%d, pw_hlp=%d, pw_vsp=%d, pw_vlp=%d\n",
406  pw_hsp, pw_hlp, pw_vsp, pw_vlp);
407  printf(" tx_freq_h=%.1f, tx_freq_v=%.1f, ",
408  tx_freq_h, tx_freq_v);
409  printf("tx_power_h=%.1f, tx_power_v=%.1f\n",
410  tx_power_h, tx_power_v);
411  printf(" tx_plen_hs=%.1f, tx_plen_hl=%.1f, ",
412  tx_plen_hs, tx_plen_hl);
413  printf("tx_plen_vs=%.1f, tx_plen_vl=%.1f\n",
414  tx_plen_vs, tx_plen_vl);
415  printf(" tx_antgain_h=%.2f, tx_antgain_v=%.2f, ",
416  tx_antgain_h, tx_antgain_v);
417  printf("rx_antgain_h=%.2f, rx_antgain_v=%.2f\n",
418  rx_antgain_h, rx_antgain_v);
419  printf(" hori_beamwidth_h=%.2f, hori_beamwidth_v=%.2f, ",
420  hori_beamwidth_h, hori_beamwidth_v);
421  printf("vert_beamwidth_h=%.2f, vert_beamwidth_v=%.2f\n",
422  vert_beamwidth_h, vert_beamwidth_v);
423  }
424  } /* end of if (opt_verbose >= 3) */
425 
426  /* range resolution (RAY) */
427  tmp_range_res = char2uint16(bufelhd + 372); //unsigned
428  if((hd->range_res > 0) && (hd->range_res != tmp_range_res)){
429  printf("!!! causion: non-uniform range resolution is not supported !!! %d, %d\n", tmp_range_res, hd->range_res);
430  goto READERR;
431  }
432  hd->range_res = tmp_range_res; //overwrite
433 
434  /*--- get range num and angle data---*/
435  sp_rnum = char2uint16(bufelhd + 376); //unsigned
436  lp_rnum = char2uint16(bufelhd + 378); //unsigned
437  rnum = sp_rnum + lp_rnum;
438  /* hd->range_num = rnum; */
439  start_obs_az = char2uint16(bufelhd + 380) * ideg2deg; //unsigned
440  start_obs_el = char2int16(bufelhd + 382) * ideg_tc2deg;
441  end_obs_az = char2uint16(bufelhd + 384) * ideg2deg; //unsigned
442  end_obs_el = char2int16(bufelhd + 386) * ideg_tc2deg;
443 
444  if(end_obs_az < start_obs_az){
445  end_obs_az += 360.0;
446  }
447  if(start_obs_el > 180.0){
448  start_obs_el -= 360.0;
449  }
450 
451  /*--- store AZ and EL ---*/
452  az[k][j] = 0.5 * (start_obs_az + end_obs_az);
453  el[k][j] = 0.5 * (start_obs_el + end_obs_el);
454 
455  if (opt_verbose >= 3){
456  printf("=== short/long rnum=%d %d (%d)",sp_rnum, lp_rnum, rnum);
457  printf(" start/end AZ=%6.2f-%6.2f (%6.2f) ",start_obs_az, end_obs_az, az[k][j]);
458  printf(" start/end EL=%5.2f-%5.2f (%5.2f)\n",start_obs_el, end_obs_el, el[k][j]);
459  }
460 
461  /*----- read data block -----*/
462  irb = rnum * (mesh_size / 8);
463  bufdat = buf + buf_offset;
464  buf_offset += irb;
465 
466  /*--- store raw data ---*/
467  for(i = 0; i < rnum; i++)
468  rtdat[k][j][i] = mesh_offset + mesh_lsb * (float)(char2uint16(bufdat + 2 * i));
469 
470  } /* end of k-loop */
471  } /* end of j-loop */
472 
473  goto READEND;
474 
475  READERR:
476  return -8;
477 
478  READEND:
479  return 0;
480 }

References scale_file::i, mod_atmos_vars::j, and scale_tracer::k.

◆ ggrads_ctl()

int ggrads_ctl ( char *  c_name,
char *  g_name,
int  nx,
int  ny,
int  nz,
int  nt,
double  xmin,
double  ymin,
double  zmin,
double  tmin,
double  dx,
double  dy,
double  dz,
double  dt,
double  no_data,
int  ihour,
int  imin,
int  iday,
int  imonth,
int  iyear,
char *  c_var,
char *  c_memo 
)
mppawr_header::start_el
float start_el
Definition: read_toshiba_mpr.h:48
mppawr_header::e_mn
int e_mn
Definition: read_toshiba_mpr.h:44
mppawr_header::s_yr
int s_yr
Definition: read_toshiba_mpr.h:43
LEN_BUFRDHD
#define LEN_BUFRDHD
Definition: read_toshiba_mpr.h:15
mppawr_header::altitude
double altitude
Definition: read_toshiba_mpr.h:47
mppawr_header::ray_num
int ray_num
Definition: read_toshiba_mpr.h:45
mppawr_header::range_num
int range_num
Definition: read_toshiba_mpr.h:45
RDIM
#define RDIM
Definition: read_toshiba.h:6
mppawr_header::s_dy
int s_dy
Definition: read_toshiba_mpr.h:43
mppawr_header::latitude
double latitude
Definition: read_toshiba_mpr.h:47
char2uint16
uint16_t char2uint16(void *input)
Definition: read_toshiba.c:36
mppawr_header::longitude
double longitude
Definition: read_toshiba_mpr.h:47
mppawr_header::el_num
int el_num
Definition: read_toshiba_mpr.h:45
decode_toshiba_mpr
int decode_toshiba_mpr(size_t bufsize, unsigned char *buf, int opt_verbose, mppawr_header *hd, float az[ELDIM][AZDIM], float el[ELDIM][AZDIM], float rtdat[ELDIM][AZDIM][RDIM])
Definition: read_toshiba_mpr.c:152
char2int16
int16_t char2int16(void *input)
Definition: read_toshiba.c:22
mppawr_header::start_az
float start_az
Definition: read_toshiba_mpr.h:48
mppawr_header::sq_name
char sq_name[16]
Definition: read_toshiba_mpr.h:42
float
typedef float(real32_t)
char2int32
int32_t char2int32(void *input)
Definition: read_toshiba.c:50
mod_atmos_vars::j
real(rp), allocatable, target, public j
Definition: mod_atmos_vars.F90:141
mppawr_header::data_name
char data_name[32]
Definition: read_toshiba_mpr.h:42
mppawr_header::range_res
int range_res
Definition: read_toshiba_mpr.h:46
mppawr_header::e_hr
int e_hr
Definition: read_toshiba_mpr.h:44
mppawr_header::e_sc
int e_sc
Definition: read_toshiba_mpr.h:44
scale_tracer::k
real(rp), public k
Definition: scale_tracer.F90:45
scale_file::i
logical, public i
Definition: scale_file.F90:196
mppawr_header::s_sc
int s_sc
Definition: read_toshiba_mpr.h:43
mppawr_header::s_hr
int s_hr
Definition: read_toshiba_mpr.h:43
mppawr_header::end_el
float end_el
Definition: read_toshiba_mpr.h:48
ungzip_toshiba_mpr
size_t ungzip_toshiba_mpr(size_t outbufsize, size_t bufsize, unsigned char *buf)
Definition: read_toshiba_mpr.c:72
mppawr_header::mesh_offset
float mesh_offset
Definition: read_toshiba_mpr.h:49
AZDIM
#define AZDIM
Definition: read_toshiba.h:7
mppawr_header::e_yr
int e_yr
Definition: read_toshiba_mpr.h:44
mppawr_header::s_mn
int s_mn
Definition: read_toshiba_mpr.h:43
mppawr_header::e_dy
int e_dy
Definition: read_toshiba_mpr.h:44
LEN_BUFELHD
#define LEN_BUFELHD
Definition: read_toshiba_mpr.h:17
mppawr_header::s_mi
int s_mi
Definition: read_toshiba_mpr.h:43
mppawr_header::site_name
char site_name[32]
Definition: read_toshiba_mpr.h:42
mppawr_header::end_az
float end_az
Definition: read_toshiba_mpr.h:48
LEN_BUFPBHD
#define LEN_BUFPBHD
Definition: read_toshiba_mpr.h:16
mppawr_header::e_mi
int e_mi
Definition: read_toshiba_mpr.h:44
ELDIM
#define ELDIM
Definition: read_toshiba.h:8
Generated by   doxygen 1.8.17