writeusrmsg.c

/** FILE NAME: writeusrmsg.c
________________________________________________________________________
   2017 U.S. Government

   Licensed under the U.S. Army Research Laboratory CC0 1.0 Universal (CC0_1.0) Public Domain Dedication 
   and the Contributor License Agreement (CLA) (together known as the "License"); you may not use this 
   file except in compliance with the License.  You may obtain a copy of the License at 

       https://github.com/USArmyResearchLab.

   Unless required by applicable law or agreed to in writing, software
   distributed under the License is distributed on an "AS IS" BASIS,
   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
   See the License for the specific language governing permissions and
   limitations under the License. 
_________________________________________________________________________

This function prints height based layer type output to the appropriate file.  It writes an ouput sounding for 
a user defined height or pressure level profile.  Input data are taken from the University of Wyoming web site, 
NOAA archive, or other source of upper air files. If wind speed is in m/s, "uncomment" the appropriate lines 
near the end. 
*/

#include "convert.h"

int writeusrmsg(struct sound *user, char *outpath)
{

   FILE *fuserout;
   int i, size;
   float zone, midheight;
   char outfile[121], infilename[71];
   
  /* Get output file name. ********************/

    strcpy(outfile, outpath);
    sscanf(user->site.filename, "%29s", infilename);
    /*printf("infilename: %s\n\n", infilename);*/ 
    strcat(outfile, infilename);
    strcat(outfile, "_USRMSG");

  /* Open output file.  ************************************/
     
     if(!(fuserout = fopen(outfile,"w")))
      {
        fprintf(stderr,"Unable to open user data file.\n");
        exit(1);
      }

   /* Print output to files. ************************************/
  
     size = user->nht+1; /*** loop indices ***/

   /*   USER DEFINED layer output output ***********************/

    fprintf(fuserout,"USER DEFINED LAYER OUTPUT\n\n");

    fprintf(fuserout,"Date: %s   Time: %s  Latitude: %9.4f     Longitude: %9.4f\n", 
                        user->site.date, user->site.time, user->site.lat, user->site.lon);
    fprintf(fuserout,"Elevation: %7.2f     Ceiling: %7.1f     Visibility: %7.1f\n\n", 
                        user->site.elev, user->site.ceil, user->site.vis);  

    fprintf(fuserout,"Layer   Height  Wind Direction  Wind Speed   Virt Temp  Pressure  Temperature\n");
    fprintf(fuserout,"          (m)      (degrees)      (kn)          (K)      (hPa)        (K)\n\n");

    for(i=0;i<size;i++)
      {
         if(user->level[i].hgt != ERROR && i > 0)       /*Mean value for midpoint.*/
            midheight = (user->level[i].hgt + user->level[i-1].hgt) * 0.5;
         else if (user->level[i].hgt != ERROR)
            midheight = user->level[0].hgt;

         /*if(user->level[i].spd != ERROR)
           {*/
             /*user->level[i].spd *= NM;*/ /*U Wyoming data already in knots. No change needed.*/
         /*}*/

         fprintf(fuserout,"%3d   %8.0f    %7.0f      %8.1f      %8.1f  %8.1f   %8.1f\n",
                 i, midheight, user->level[i].dir, user->level[i].spd,
                 user->level[i].vtmp, user->level[i].prs, user->level[i].tmp);
      }

    printf("\n User defined message printed to USRMSG output.\n");

/**********End of output statements**********************/

   fclose(fuserout);

   return;
}