external/io_grib1/WGRIB/ec_ext.c

   1 #include <stdio.h>
   2 #include <stdlib.h>
   3 #include <stddef.h>
   4 #include <string.h>
   5
   6 #include "pds4.h"
   7 #include "grib.h"
   8
   9 /*
  10  * EC_ext       v1.0 wesley ebisuzaki
  11  *
  12  * prints something readable from the EC stream parameter
  13  *
  14  * prefix and suffix are only printed if EC_ext has text
  15  */
  16
  17 void EC_ext(unsigned char *pds, char *prefix, char *suffix) {
  18
  19     /* int i;
  20     printf("\n");
  21     for (i=0; i < PDS_LEN(pds); i++) {
  22       printf("%x ",pds[i]);
  23     }
  24     */
  25     /*
  26      10/03/2000: R.Rudsar : subroutine changed.
  27                  Tests for EcType and extra test for EcStream 1035
  28     */
  29     if (PDS_Center(pds) == ECMWF && PDS_LEN(pds) >= 43) {
  30
  31         switch(PDS_EcType(pds)) {
  32             case 10:
  33                 printf("%sControl forecast%s", prefix, suffix);
  34                 break;
  35             case 11:
  36                 printf("%sPerturbed forecasts%s", prefix, suffix);
  37                 break;
  38             case 14:
  39                 printf("%sCluster means%s", prefix, suffix);
  40                 break;
  41             case 15:
  42                 printf("%sCluster std. dev.%s", prefix, suffix);
  43                 break;
  44             case 16:
  45                 printf("%sForecast probability%s", prefix, suffix);
  46                 break;
  47             case 17:
  48                 printf("%sEnsemble means%s", prefix, suffix);
  49                 break;
  50             case 18:
  51                 printf("%sEnsemble std. dev.%s", prefix, suffix);
  52                 break;
  53             default:
  54                 printf("%sECMWF type?%s", prefix, suffix);
  55                 break;
  56         }
  57     }
  58     if (PDS_Center(pds) == ECMWF && PDS_LEN(pds) >= 45) {
  59
  60         switch(PDS_EcStream(pds)) {
  61             case 1035:
  62                 printf("%sensemble forecasts%s", prefix, suffix);
  63                 break;
  64             case 1043:
  65                 printf("%smon mean%s", prefix, suffix);
  66                 break;
  67             case 1070:
  68                 printf("%smon (co)var%s", prefix, suffix);
  69                 break;
  70             case 1071:
  71                 printf("%smon mean from daily%s", prefix, suffix);
  72                 break;
  73             default:
  74                 printf("%sECMWF stream?%s", prefix, suffix);
  75                 break;
  76         }
  77     }
  78 }