create_rsel2_file.c

Go to the documentation of this file.

#include <stdio.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <X11/Xos.h>
#include <ctype.h>
#include "mm2rfs.h"
 
#include <Xm/Text.h>
#include <Xm/TextF.h>
 
extern struct time_cntl tc;
extern struct fgroup_info fg[20];
extern fgroup_num;
extern Widget text_edit;
extern Widget text_output;
extern int dur_toggle_item_set;
extern int snow_toggle_item_set;
 
int create_rsel2_file(int remove)
{
FILE *fp;
char *text;
char fg_last[9];
char line[80];
char tmp_file[80];
int i,ii,j,k,kk,last_per;
int y,m,d,h;
int vy,vm,vd,vh;
int nhz;
int rsel;
int len,len_actual;
int hz2rsel=-1000;
float fhz;
struct fmaz_struct fmaz[100];
 
  last_per=3+(4*(dur_toggle_item_set-1));
 
  /* load structure with fmaz informaiton */
  i=fmaz_load(&fmaz[0]);
 
 
  /* identify snow level offset */
  if(snow_toggle_item_set < 1 || snow_toggle_item_set >= 7)
    hz2rsel=-1000;
  else
    hz2rsel=1000-500*snow_toggle_item_set;
 
  /* make sure that there are postable fmazs */
  ii=0;
  if(i)
  {
    for(j=0;j<i;j++)
      if(strlen(fmaz[j].fmazid))
        ii=1;
  }
        
  if(i && ii)
  {
    /* set time control */
    adj_time(-6,tc.year,tc.month,tc.day,tc.hour,&y,&m,&d,&h);
    if(tc.period==1) h=10;
    if(tc.period==2) h=16;
    if(tc.period==3) h=22;
    if(tc.period==4) h=4;
    /* mods are valid for the next 3 computational periods */
    /* mods are valid for the current period only */
    adj_time(0,y,m,d,h,&vy,&vm,&vd,&vh);
    y=y%100;
    vy=vy%100;
 
    /* create tmp filename */
    sprintf(tmp_file,"rsel_fg.%d%d%d.%d%d%d",
            tc.year,tc.month,tc.day,tc.hour,tc.minute,tc.second);
 
    if((fp = fopen (tmp_file, "w")))
    {
      strcpy(fg_last,"\0");
      for(j=0;j<i;j++)
      {
        if( strlen(fmaz[j].fmazid) && 
            (!strcmp(fmaz[j].pos,"O") || !strcmp(fmaz[j].pos,"U")) 
          )
        {
          if(strcmp(fg_last,fmaz[j].fgroupid))
          {
            if(j)
              fprintf(fp,"$\n");
            for(k=0;k<fgroup_num;k++)
              if(!strcmp(fg[k].id,fmaz[j].fgroupid))
                fprintf(fp,"$%s\n.TSCHNG "
                           "%02d%02d%02d%02dPST %02d%02d%02d%02dPST\n",
                            fg[k].desc,m,d,y,h,vm,vd,vy,vh);
          }
          strcpy(fg_last,fmaz[j].fgroupid);
  
          sprintf(line,"%s        ",fmaz[j].segid);
          *(line+8)='\0';
          fprintf(fp,"%s ",line);
  
          sprintf(line,"%s        ",fmaz[j].fmazid);
          *(line+8)='\0';
          fprintf(fp,"%s ",line);
  
/*
          sprintf(line,"%s    ",rsel[j].ts_dt);
*/
      
          sprintf(line,"RSEL    ");
          *(line+4)='\0';
          fprintf(fp,"%s  6 ",line);
  
          for(k=0;k<=last_per;k++)
          {
            if(k<tc.period-1)
              strcpy(line,"     ");
            else
            {
              
              fhz=0.;nhz=0;
              for(kk=k;kk<=k+1;kk++)
              {
                if(strcmp(fmaz[j].per[kk].qual,"M"))
                {
                  fhz=fhz+fmaz[j].per[kk].amt;
                  nhz++;
                }
              }
              if(nhz)
              {
                fhz=fhz/(float)nhz;
                rsel=100*((int)(10*(fhz+.0499)))+hz2rsel;
                if(rsel<0) rsel=0;
                sprintf(line,"%5d",rsel);
              }
              else
                sprintf(line,"-9999");
            }
            fprintf(fp," %s",line);
 
            if(!((k+1)%4) && k != last_per)
              fprintf(fp,"  &\n                          ");
          }
          fprintf(fp,"    SNOW-17\n");
/*
          fprintf(fp,"    %s\n",rsel[j].oper);
*/
        }
      }
 
      fclose(fp);
      if(remove)
      {
        if((fp = fopen (tmp_file, "r")))
        { 
          len=10000;
          if (!(text = XtMalloc ((unsigned)(len+1)))) /* +1 for NULL */
          {
            sprintf (line, "%s: XtMalloc(%ld) failed", len, tmp_file);
            if(text_output)
              XmTextFieldSetString (text_output, line);
            return(0);
          }
          else
          {
            len_actual=fread (text, sizeof (char), len, fp);
            text[len_actual] = 0; /* NULL-terminate */
            if(text_edit);
              XmTextSetString (text_edit, text);
            sprintf (line, "Click 'Save as *' to store as '* File'");
            if(text_output);
              XmTextFieldSetString (text_output, line);
          }
          /* close and rm file */
          fclose(fp);
          sprintf(line,"rm %s",tmp_file);
          system(line);
        }
        else
        {
          sprintf(line,"Couldn't open temp file %s",tmp_file);
          if(text_output)
            XmTextFieldSetString (text_output, line);
          return(0);
        }
      }
    }
    else
    {
      sprintf(line,"Couldn't open temp file %s",tmp_file);
      if(text_output)
        XmTextFieldSetString (text_output, line);
      return(0);
    }
  }
  else
  {
    sprintf(line,"No Freezing Level forecasts in table fmaqtrly to read.");
    if(text_output)
      XmTextFieldSetString (text_output, line);
    return(0);
  }
   
  return(1);
}