Reverted some mistakes andersFeatureBranch master jenkins-wrwp-43
authorAnders Henja <anders@henjab.se>
Thu, 18 Oct 2018 18:32:16 +0000 (20:32 +0200)
committerAnders Henja <anders@henjab.se>
Thu, 18 Oct 2018 18:32:16 +0000 (20:32 +0200)
lib/wrwp.c
lib/wrwp.h

index 6fe0c7a..f0d621e 100644 (file)
@@ -458,7 +458,7 @@ VerticalProfile_t* Wrwp_generate(Wrwp_t* self, PolarVolume_t* inobj, const char*
        become the HGHT array */
     centerOfLayer = iz + (self->dz / 2.0);
     firstInit = 0;
-    fprintf(stderr, "0\n");
+
     // loop over scans
     for (is = 0; is < nscans; is++) {
       char* malfuncString = NULL;
@@ -470,7 +470,6 @@ VerticalProfile_t* Wrwp_generate(Wrwp_t* self, PolarVolume_t* inobj, const char*
       
       if (elangleForThisScan * RAD2DEG >= self->emin) { /* We only do the calculation for scans with elangle >= the minimum one */
         RaveDateTime_t* startDTofThisScan = WrwpInternal_getStartDateTimeFromScan(scan);
-       fprintf(stderr, "1: REFCNT = %d\n", RAVE_OBJECT_REFCNT(startDTofThisScan));
         RaveDateTime_t* endDTofThisScan = WrwpInternal_getEndDateTimeFromScan(scan);        
         RaveAttribute_t* malfuncattr = PolarScan_getAttribute(scan, "how/malfunc");
         if (malfuncattr != NULL) {
@@ -484,16 +483,12 @@ VerticalProfile_t* Wrwp_generate(Wrwp_t* self, PolarVolume_t* inobj, const char*
             firstStartDT = RAVE_OBJECT_COPY(startDTofThisScan);
             lastEndDT = RAVE_OBJECT_COPY(endDTofThisScan);
             firstInit = 1;
-           fprintf(stderr, "2: REFCNT = %d\n", RAVE_OBJECT_REFCNT(startDTofThisScan));
-           fprintf(stderr, "3: REFCNT = %d\n", RAVE_OBJECT_REFCNT(firstStartDT));
           }
           if (firstInit == 1 && countAcceptedScans > 1 && iz==0) {
             if (RaveDateTime_compare(startDTofThisScan, firstStartDT) < 0) {
               /* if start datetime of this scan is before the first saved start datetime, save this one instead */
               RAVE_OBJECT_RELEASE(firstStartDT);
               firstStartDT = RAVE_OBJECT_COPY(startDTofThisScan);
-             fprintf(stderr, "4: REFCNT = %d\n", RAVE_OBJECT_REFCNT(startDTofThisScan));
-             fprintf(stderr, "5: REFCNT = %d\n", RAVE_OBJECT_REFCNT(firstStartDT));
             }
             if (RaveDateTime_compare(endDTofThisScan, lastEndDT) > 0) {
               /* If end datetime of this scan is after the last saved end datetime, save this one instead */
@@ -525,8 +520,6 @@ VerticalProfile_t* Wrwp_generate(Wrwp_t* self, PolarVolume_t* inobj, const char*
                     (val != nodata) &&
                     (val != undetect) &&
                     (abs(offset + gain * val) >= self->vmin)) {
-                  fprintf(stderr, "NV=%d, NOR=%d\n",nv, NOR);
-                  if (nv > NOR*NOC-2) fprintf(stderr, "WARNING WARNING, nv=%d\n", nv);
                   *(v+nv) = offset+gain*val;
                   *(az+nv) = 360./nrays*ir*DEG2RAD;
                   *(A+nv*NOC) = sin(*(az+nv));
@@ -567,17 +560,12 @@ VerticalProfile_t* Wrwp_generate(Wrwp_t* self, PolarVolume_t* inobj, const char*
             RAVE_OBJECT_RELEASE(dbz);         
           }        
         }
-        fprintf(stderr, "6: REFCNT = %d\n", RAVE_OBJECT_REFCNT(startDTofThisScan));
-
         RAVE_OBJECT_RELEASE(startDTofThisScan);
-       fprintf(stderr, "6.01\n");
         RAVE_OBJECT_RELEASE(endDTofThisScan);
-       fprintf(stderr, "6.02\n");
       }
       RAVE_OBJECT_RELEASE(scan); 
-      fprintf(stderr, "6.03\n");
     }
-    fprintf(stderr, "6.1\n");
+
     if (countAcceptedScans == 0) { /* Emergency exit if no accepted scans were found */
       result = NULL;               /* If this is the case, we don't bother with checking */
       RAVE_FREE(A);                /* the same thing for the other atmospheric layers, */
@@ -775,8 +763,6 @@ done:
   RAVE_OBJECT_RELEASE(hght_field);
   RAVE_OBJECT_RELEASE(uwnd_field);
   RAVE_OBJECT_RELEASE(vwnd_field);
-  fprintf(stderr, "RELEASING FIRST START DT\n");
-  fprintf(stderr, "7: REFCNT = %d\n", RAVE_OBJECT_REFCNT(firstStartDT));
   RAVE_OBJECT_RELEASE(firstStartDT);
   RAVE_OBJECT_RELEASE(lastEndDT);
   RaveList_freeAndDestroy(&wantedFields);
index 57cb235..a3f9cbc 100644 (file)
@@ -56,7 +56,7 @@ along with baltrad-wrwp.  If not, see <http://www.gnu.org/licenses/>.
 
 #define DEG2RAD     DEG_TO_RAD      /* Degrees to radians. From PROJ.4 */
 #define RAD2DEG     RAD_TO_DEG      /* Radians to degrees. From PROJ.4 */
-#define NOR         30000           /* Number of rows in matrix A used in the computation */
+#define NOR         20000           /* Number of rows in matrix A used in the computation */
 #define NOC         3               /* Number of columns in matrix A used in the computation */
 #define NRHS        1               /* Number of right-hand sides; that is, the number of columns in matrix B used in the computation */
 #define LDA         NOC             /* Leading dimension of the array specified for a */
@@ -64,8 +64,7 @@ along with baltrad-wrwp.  If not, see <http://www.gnu.org/licenses/>.
 #define DMIN        4000            /* Minimum distance for deriving a profile [m] */
 #define DMAX        40000           /* Maximum distance for deriving a profile [m] */
 #define NMIN        36              /* Minimum sample size */
-#define EMIN        0.5             /* Minimum elevation angle [deg] */
-//#define EMAX        45.0            /* Maximum elevation angle [deg] */
+#define EMIN        2.5             /* Minimum elevation angle [deg] */
 #define VMIN        2.0             /* Radial velocity threshold [m/s] */
 #define DZ          200             /* Height interval for deriving a profile [m] */
 #define HMAX        12000           /* Maximum height of the profile [m] */
@@ -211,20 +210,6 @@ void Wrwp_setEMIN(Wrwp_t* self, double emin);
 double Wrwp_getEMIN(Wrwp_t* self);
 
 /**
- * Sets maximum elevation angle [deg]
- * @param[in] self - self
- * @param[in] emax - maximum elevation angle [deg]
- */
-//void Wrwp_setEMAX(Wrwp_t* self, double emax);
-
-/**
- * Returns maximum elevation angle [deg]
- * @param[in] self - self
- * @return Maximum elevation angle [deg] (default EMAX)
- */
-//double Wrwp_getEMAX(Wrwp_t* self);
-
-/**
  * Sets radial velocity threshold [m/s]
  * @param[in] self - self
  * @param[in] vmin - radial velocity threshold [m/s]