bufr2tac_synop.c File Reference

file with the code to parse a sequence as a SYNOP FM-12, FM-13 and FM-14 More...

#include "bufr2tac.h"

Include dependency graph for bufr2tac_synop.c:

Go to the source code of this file.

Functions
int	synop_YYYYMMDDHHmm_to_YYGG (struct synop_chunks *syn)
	Sets YYGG from YYYYMMDDHHmm extended group. More...
char *	guess_WMO_region_synop (struct synop_chunks *syn)
	Try to find WMO region if it is not already set and WMO Block and number index are known. More...
int	parse_subset_as_synop (struct metreport *m, struct bufr2tac_subset_state *s, struct bufr_subset_sequence_data *sq, char *err)
	parses a subset sequence as an Land fixed SYNOP FM-12, SHIP FM-13 or SYNOP-mobil FM-14 report More...

Detailed Description

file with the code to parse a sequence as a SYNOP FM-12, FM-13 and FM-14

Definition in file bufr2tac_synop.c.

Function Documentation

guess_WMO_region_synop()

char * guess_WMO_region_synop

(

struct synop_chunks *

syn

)

Try to find WMO region if it is not already set and WMO Block and number index are known.

Parameters

syn	pointer to the struct synop_chunks with all known data for a synop

It returns a pointer to the string with WMO region

Definition at line 93 of file bufr2tac_synop.c.

{
  if ( syn->s0.A1[0] )
    {
      return syn->s0.A1;
    }
 
  if ( syn->s0.II[0] == 0  || syn->s0.iii[0] == 0 )
    {
      return syn->s0.A1;
    }
 
  return guess_WMO_region ( syn->s0.A1, syn->s0.Reg, syn->s0.II, syn->s0.iii );
}

References synop_sec0::A1, guess_WMO_region(), synop_sec0::II, synop_sec0::iii, synop_sec0::Reg, and synop_chunks::s0.

Here is the call graph for this function:

parse_subset_as_synop()

int parse_subset_as_synop	(	struct metreport *	m,
		struct bufr2tac_subset_state *	s,
		struct bufr_subset_sequence_data *	sq,
		char *	err
	)

parses a subset sequence as an Land fixed SYNOP FM-12, SHIP FM-13 or SYNOP-mobil FM-14 report

Parameters

m	pointer to a struct metreport where set some results
s	pointer to a struct bufr2tac_subset_state
sq	pointer to a struct bufr_subset_sequence_data with the parsed sequence on input
err	string with detected errors, if any

It return 0 if all is OK. Otherwise returns 1 and it also fills the err string

Definition at line 119 of file bufr2tac_synop.c.

{
  size_t is;
  char aux[32];
  struct synop_chunks *syn;
 
  // An auxiliar pointer
  syn = &m->synop;
 
  // clean data
  bufr2tac_clean_synop_chunks ( syn );
 
  // reject if still not coded type
  if ( strcmp ( s->type_report,"AAXX" ) && strcmp ( s->type_report,"BBXX" ) && strcmp ( s->type_report,"OOXX" ) )
    {
      sprintf ( err,"bufr2tac: parse_subset_as_synop(): '%s' reports still not decoded in this software", s->type_report );
      return 1;
    }
 
  syn->s0.MiMi[0] = s->type_report[0];
  syn->s0.MiMi[1] = s->type_report[1];
  syn->s0.MjMj[0] = s->type_report[2];
  syn->s0.MjMj[1] = s->type_report[3];
 
  strcpy ( m->type, s->type_report );
 
  /**** First pass, sequential analysis *****/
  for ( is = 0; is < sq->nd; is++ )
    {
      // check if is a significance qualifier
      if ( sq->sequence[is].desc.x == 8 )
        {
          s->i = is;
          s->a = &sq->sequence[is];
          s->ival = ( int ) sq->sequence[is].val;
          s->val = sq->sequence[is].val;
          syn_parse_x08 ( syn, s );
        }
      if ( s->isq )  // case of a significance qualifier
        {
          continue;
        }
 
      s->i = is;
      s->ival = ( int ) ( sq->sequence[is].val );
      s->val = sq->sequence[is].val;
      s->a = &sq->sequence[is];
      if ( is > 0 )
        {
          s->a1 = &sq->sequence[is - 1];
        }
 
      switch ( sq->sequence[is].desc.x )
        {
        case 1: //localization descriptors
          syn_parse_x01 ( syn, s );
          break;
 
        case 2: //Type of station descriptors
          syn_parse_x02 ( syn, s );
          break;
 
        case 4: //Date and time descriptors
          syn_parse_x04 ( syn, s );
          break;
 
        case 5: // Horizontal position. Latitude
          syn_parse_x05 ( syn, s );
          break;
 
        case 6: // Horizontal position. Longitude
          syn_parse_x06 ( syn, s );
          break;
 
        case 7: // Vertical position
          syn_parse_x07 ( syn, s );
          break;
 
        case 10: // Air Pressure descriptors
          syn_parse_x10 ( syn, s );
          break;
 
        case 11: // wind  data
          syn_parse_x11 ( syn, s );
          break;
 
        case 12: //Temperature descriptors
          syn_parse_x12 ( syn, s );
          break;
 
        case 13: // Humidity and precipitation data
          syn_parse_x13 ( syn, s );
          break;
 
        case 14: // Radiation
          syn_parse_x14 ( syn, s );
          break;
 
        case 20: // Cloud data
          syn_parse_x20 ( syn, s );
          break;
 
        case 22: // Oceanographic data
          syn_parse_x22 ( syn, s );
          break;
 
        case 31: // Replicators
          syn_parse_x31 ( syn, s );
          break;
 
 
        default:
          break;
        }
    }
 
  /* Check about needed descriptors */
  if ( ( ( s->mask & SUBSET_MASK_HAVE_LATITUDE ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_LONGITUDE ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_YEAR ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_MONTH ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_DAY ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_HOUR ) == 0 ) ||
       ( ( s->mask & SUBSET_MASK_HAVE_MINUTE ) == 0 )
     )
    {
      sprintf ( err,"bufr2tac: %s(): lack of mandatory descriptor in sequence", __func__ );
      return 1;
    }
 
  // copy WIGOS ID 
  memcpy(&m->g.wid, &m->synop.wid, sizeof (struct wigos_id) );
  
  if ( strcmp ( s->type_report, "AAXX" ) == 0 && ( syn->s0.II[0] == 0 || syn->s0.iii[0] == 0 ) )
    {
      if ( s->mask & SUBSET_MASK_HAVE_WIGOS_ID &&
           ( syn->wid.issuer != 0 || syn->wid.issue != 0 ) &&
           syn->wid.local_id[0] )
        {
          strcpy ( syn->s0.II, "00" );
          strcpy ( syn->s0.iii, "000" );
        }
      else
        {
          sprintf ( err,"bufr2tac: %s(): lack of mandatory index station for AAXX report", __func__ );
          return 1;
        }
    }
 
 
  /****** Second pass. Global results and consistence analysis ************/
  // set ir
  if ( syn->s1.RRR[0] == 0 && syn->s3.RRR[0] == 0 && syn->s3.RRRR24[0] == 0 )
    {
      // no prec data case
      syn->s1.ir[0] = '4';
    }
  else if ( syn->s1.RRR[0] == 0 && syn->s3.RRR[0] != 0 )
    {
      syn->s1.ir[0] = '2';
    }
  else if ( syn->s1.RRR[0] != 0 && syn->s3.RRR[0] == 0 )
    {
      syn->s1.ir[0] = '1';
    }
  else
    {
      syn->s1.ir[0] = '0';
    }
 
  // Supress 7wwW1W2 if no significant weather
  if ( ( syn->s1.ww[0] || syn->s1.W1[0] || syn->s1.W2[0] ) &&
       ( syn->s1.ww[0] == 0 || syn->s1.ww[0] == '/' || ( strcmp ( syn->s1.ww,"04" ) < 0 ) ) &&
       ( syn->s1.W1[0] == 0 || syn->s1.W1[0] == '/' || ( strcmp ( syn->s1.W1,"3" ) < 0 ) ) &&
       ( syn->s1.W2[0] == 0 || syn->s1.W2[0] == '/' || ( strcmp ( syn->s1.W2,"3" ) < 0 ) ) )
    {
      syn->s1.ww[0] = 0;
      syn->s1.W1[0] = 0;
      syn->s1.W2[0] = 0;
      s->mask |= ( SUBSET_MASK_HAVE_NO_SIGNIFICANT_WW | SUBSET_MASK_HAVE_NO_SIGNIFICANT_W1 | SUBSET_MASK_HAVE_NO_SIGNIFICANT_W2 );
      if ( syn->s1.ix[0] == '1' )
        {
          syn->s1.ix[0] = '2';
        }
    }
 
  // adjust iw
  if ( syn->s0.iw[0] == '/' && syn->s1.ff[0] != '/' )
    {
      syn->s0.iw[0] = '1';
    }
 
  // adjust ix
  if ( syn->s1.ix[0] == '/' /*&& syn->s1.ww[0] == 0*/ )
    {
      if ( ( s->mask & SUBSET_MASK_HAVE_NO_SIGNIFICANT_WW ) &&
           ( s->mask & SUBSET_MASK_HAVE_NO_SIGNIFICANT_W1 ) &&
           ( s->mask & SUBSET_MASK_HAVE_NO_SIGNIFICANT_W2 )
         )
        {
          if ( s->mask & SUBSET_MASK_HAVE_TYPE_STATION )
            {
              if ( s->type == 1 || s->type == 2 )
                {
                  strcpy ( syn->s1.ix,"2" );
                }
              else if ( s->type == 0 )
                {
                  strcpy ( syn->s1.ix,"5" );
                }
            }
          else
            {
              strcpy ( syn->s1.ix,"6" );  //NOTE: here we assume an automatic station without W data
            }
        }
      else if ( s->mask & SUBSET_MASK_HAVE_TYPE_STATION )
        {
          if ( s->type == 1 )
            {
              strcpy ( syn->s1.ix,"3" );
            }
          else if ( s->type == 0 )
            {
              strcpy ( syn->s1.ix,"6" );
            }
        }
      else
        {
          strcpy ( syn->s1.ix,"6" );  //NOTE: here we assume an automatic station without W data
        }
    }
 
  // adjust radiation
  if ( ( syn->mask & SYNOP_SEC3 ) && ( syn->s3.SS[0] == '\0' ) &&
       ( syn->s3.j5[0][0] || syn->s3.j5[1][0] ||
         syn->s3.j5[2][0] || syn->s3.j5[3][0] ||
         syn->s3.j5[4][0] || syn->s3.j5[5][0] ||
         syn->s3.j5[6][0] ) )
    {
      sprintf ( syn->s3.SS, "//" );
    }
 
  if ( ( syn->mask & SYNOP_SEC3 ) && ( syn->s3.SSS[0] == '\0' ) &&
       ( syn->s3.j524[0][0] || syn->s3.j524[1][0] ||
         syn->s3.j524[2][0] || syn->s3.j524[3][0] ||
         syn->s3.j524[4][0] || syn->s3.j524[5][0] ||
         syn->s3.j524[6][0] ) )
    {
      sprintf ( syn->s3.SSS, "///" );
    }
 
  /****** Final Adjust ***********/
 
  // fix YYGG according with YYYYMMDDHHmm
  if ( syn->e.mm[0] == 0 )
    sprintf ( syn->e.mm, "00" );
  synop_YYYYMMDDHHmm_to_YYGG ( syn );
  if ( strcmp ( syn->e.mm,"00" ) )
    {
      strcpy ( syn->s1.GG, syn->e.HH );
      strcpy ( syn->s1.gg, syn->e.mm );
    }
  syn->mask |= SYNOP_EXT;
 
  // Marc as a synop from bufr
  syn->mask |= SYNOP_BUFR;
 
  // Fill some metreport fields
  if ( strlen ( syn->s0.II ) )
    {
      strcpy ( m->g.index, syn->s0.II );
      strcat ( m->g.index, syn->s0.iii );
    }
  else if ( strlen ( syn->s0.D_D ) )
    {
      strcpy ( m->g.index, syn->s0.D_D );
    }
  else if ( strlen ( syn->s0.IIIII ) )
    {
      strcpy ( m->g.index, syn->s0.IIIII );
    }
 
  if ( s->mask & SUBSET_MASK_HAVE_LATITUDE )
    {
      m->g.lat = s->lat;
    }
  if ( s->mask & SUBSET_MASK_HAVE_LONGITUDE )
    {
      m->g.lon = s->lon;
    }
  if ( s->mask & SUBSET_MASK_HAVE_ALTITUDE )
    {
      m->g.alt = s->alt;
    }
  if ( s->mask & SUBSET_MASK_HAVE_NAME )
    {
      strcpy ( m->g.name, s->name );
    }
  if ( s->mask & SUBSET_MASK_HAVE_COUNTRY )
    {
      strcpy ( m->g.country, s->country );
    }
 
  sprintf ( aux,"%s%s%s%s%s", syn->e.YYYY, syn->e.MM, syn->e.DD, syn->e.HH, syn->e.mm );
  YYYYMMDDHHmm_to_met_datetime ( &m->t, aux );
 
  if ( check_date_from_future ( m ) )
    {
      return 1;  // Bad date/time . Is a report from future!
    }
 
  // If finally we arrive here, It succeded
  return 0;
}

References bufr2tac_subset_state::a, bufr2tac_subset_state::a1, bufr2tac_subset_state::alt, met_geo::alt, bufr2tac_clean_synop_chunks(), check_date_from_future(), bufr2tac_subset_state::country, met_geo::country, synop_sec0::D_D, report_date_ext::DD, synop_chunks::e, synop_sec1::ff, metreport::g, synop_sec1::GG, synop_sec1::gg, report_date_ext::HH, bufr2tac_subset_state::i, synop_sec0::II, synop_sec0::iii, synop_sec0::IIIII, met_geo::index, synop_sec1::ir, bufr2tac_subset_state::isq, wigos_id::issue, wigos_id::issuer, bufr2tac_subset_state::ival, synop_sec0::iw, synop_sec1::ix, synop_sec3::j5, synop_sec3::j524, bufr2tac_subset_state::lat, met_geo::lat, wigos_id::local_id, bufr2tac_subset_state::lon, met_geo::lon, bufr2tac_subset_state::mask, synop_chunks::mask, synop_sec0::MiMi, synop_sec0::MjMj, report_date_ext::MM, report_date_ext::mm, bufr2tac_subset_state::name, met_geo::name, synop_sec1::RRR, synop_sec3::RRR, synop_sec3::RRRR24, synop_chunks::s0, synop_chunks::s1, synop_chunks::s3, synop_sec3::SS, synop_sec3::SSS, SUBSET_MASK_HAVE_ALTITUDE, SUBSET_MASK_HAVE_COUNTRY, SUBSET_MASK_HAVE_DAY, SUBSET_MASK_HAVE_HOUR, SUBSET_MASK_HAVE_LATITUDE, SUBSET_MASK_HAVE_LONGITUDE, SUBSET_MASK_HAVE_MINUTE, SUBSET_MASK_HAVE_MONTH, SUBSET_MASK_HAVE_NAME, SUBSET_MASK_HAVE_NO_SIGNIFICANT_W1, SUBSET_MASK_HAVE_NO_SIGNIFICANT_W2, SUBSET_MASK_HAVE_NO_SIGNIFICANT_WW, SUBSET_MASK_HAVE_TYPE_STATION, SUBSET_MASK_HAVE_WIGOS_ID, SUBSET_MASK_HAVE_YEAR, syn_parse_x01(), syn_parse_x02(), syn_parse_x04(), syn_parse_x05(), syn_parse_x06(), syn_parse_x07(), syn_parse_x08(), syn_parse_x10(), syn_parse_x11(), syn_parse_x12(), syn_parse_x13(), syn_parse_x14(), syn_parse_x20(), syn_parse_x22(), syn_parse_x31(), metreport::synop, SYNOP_BUFR, SYNOP_EXT, SYNOP_SEC3, synop_YYYYMMDDHHmm_to_YYGG(), metreport::t, bufr2tac_subset_state::type, metreport::type, bufr2tac_subset_state::type_report, bufr2tac_subset_state::val, synop_sec1::W1, synop_sec1::W2, met_geo::wid, synop_chunks::wid, synop_sec1::ww, report_date_ext::YYYY, and YYYYMMDDHHmm_to_met_datetime().

Referenced by parse_subset_sequence().

Here is the call graph for this function:

Here is the caller graph for this function:

synop_YYYYMMDDHHmm_to_YYGG()

int synop_YYYYMMDDHHmm_to_YYGG

(

struct synop_chunks *

syn

)

Sets YYGG from YYYYMMDDHHmm extended group.

Parameters

syn	pointer to the target struct synop_chunks

Definition at line 31 of file bufr2tac_synop.c.

{
  char aux[20];
  char tz[256], *c;
  time_t t;
  struct tm tim;
 
  if ( strlen ( syn->e.YYYY ) &&
       strlen ( syn->e.MM ) &&
       strlen ( syn->e.DD ) &&
       strlen ( syn->e.HH ) &&
       strlen ( syn->e.mm ) )
    {
      sprintf ( aux,"%s%s%s%s%s", syn->e.YYYY, syn->e.MM, syn->e.DD, syn->e.HH, syn->e.mm );
    }
 
  if ( strlen ( aux ) != 12 )
    {
      return 1;
    }
 
  // Get current TZ
  tz[0] = '\0';
  c = getenv ( "TZ" );
  if ( c != NULL && c[0] && strlen ( c ) < 256 )
    {
      strcpy ( tz, c );
    }
 
  // Set TZ to UTC
  setenv ( "TZ", "UTC", 1 );
  tzset();
 
  memset ( &tim, 0, sizeof ( struct tm ) );
  strptime ( aux, "%Y%m%d%H%M", &tim );
 
  t = mktime ( &tim ) + 1799 ; // rounding trick
  gmtime_r ( &t, &tim );
  sprintf ( syn->s0.YY, "%02d", tim.tm_mday );
  sprintf ( syn->s0.GG, "%02d", tim.tm_hour );
 
  // Revert TZ changes
  if ( tz[0] )
    {
      setenv ( "TZ", tz, 1 );
      tzset();
    }
  else
    {
      unsetenv ( "TZ" );
    }
 
  return 0;
}

References report_date_ext::DD, synop_chunks::e, synop_sec0::GG, report_date_ext::HH, report_date_ext::MM, report_date_ext::mm, synop_chunks::s0, synop_sec0::YY, and report_date_ext::YYYY.

Referenced by parse_subset_as_synop().

Here is the caller graph for this function: