int main() { int appid, wid, pid; int rlist; char *labels[5]; NhlColorIndex colors[5]; NhlMarkerIndex item_ind[5]; float mkthik[5]; char const *wks_type = "x11"; /* * Initialize data values */ labels[0] = "Marker_0"; labels[1] = "Marker_1"; labels[2] = "Marker_2"; labels[3] = "Marker_3"; labels[4] = "Marker_4"; colors[0] = 3; colors[1] = 5; colors[2] = 7; colors[3] = 9; colors[4] = 11; mkthik[0] = 2.; mkthik[1] = 3.; mkthik[2] = 4.; mkthik[3] = 5.; mkthik[4] = 6.; item_ind[0] = 2; item_ind[1] = 3; item_ind[2] = 4; item_ind[3] = 5; item_ind[4] = 6; /* * Initialize the high level utility library */ NhlInitialize(); /* * Create an application context. Set the app dir to the current directory * so the application looks for a resource file in the working directory. * In this example the resource file supplies the plot title only. */ rlist = NhlRLCreate(NhlSETRL); NhlRLClear(rlist); NhlRLSetString(rlist,NhlNappDefaultParent,"True"); NhlRLSetString(rlist,NhlNappUsrDir,"./"); NhlCreate(&appid,"lg02",NhlappClass,NhlDEFAULT_APP,rlist); if (!strcmp(wks_type,"ncgm") || !strcmp(wks_type,"NCGM")) { /* * Create a meta file workstation. */ NhlRLClear(rlist); NhlRLSetString(rlist,NhlNwkMetaName,"./lg02c.ncgm"); NhlCreate(&wid,"lg02Work",NhlncgmWorkstationClass,NhlDEFAULT_APP, rlist); } else if (!strcmp(wks_type,"x11") || !strcmp(wks_type,"X11")) { /* * Create an X Workstation. */ NhlRLClear(rlist); NhlRLSetInteger(rlist,NhlNwkPause,True); NhlCreate(&wid,"lg02Work",NhlcairoWindowWorkstationClass,NhlDEFAULT_APP, rlist); } else if (!strcmp(wks_type,"oldps") || !strcmp(wks_type,"OLDPS")) { /* * Create an older-style PostScript workstation. */ NhlRLClear(rlist); NhlRLSetString(rlist,NhlNwkPSFileName,"./lg02c.ps"); NhlCreate(&wid,"lg02Work",NhlpsWorkstationClass,NhlDEFAULT_APP, rlist); } else if (!strcmp(wks_type,"oldpdf") || !strcmp(wks_type,"OLDPDF")) { /* * Create an older-style PDF workstation. */ NhlRLClear(rlist); NhlRLSetString(rlist,NhlNwkPDFFileName,"./lg02c.pdf"); NhlCreate(&wid,"lg02Work",NhlpdfWorkstationClass,NhlDEFAULT_APP, rlist); } else if (!strcmp(wks_type,"pdf") || !strcmp(wks_type,"PDF") || !strcmp(wks_type,"ps") || !strcmp(wks_type,"PS")) { /* * Create a cairo PS/PDF workstation. */ NhlRLClear(rlist); NhlRLSetString(rlist,NhlNwkFileName,"./lg02c"); NhlRLSetString(rlist,NhlNwkFormat,(char*)wks_type); NhlCreate(&wid,"lg02Work",NhlcairoDocumentWorkstationClass,NhlDEFAULT_APP, rlist); } else if (!strcmp(wks_type,"png") || !strcmp(wks_type,"PNG")) { /* * Create a cairo PNG workstation. */ NhlRLClear(rlist); NhlRLSetString(rlist,NhlNwkFileName,"./lg02c"); NhlRLSetString(rlist,NhlNwkFormat,(char*)wks_type); NhlCreate(&wid,"lg02Work",NhlcairoImageWorkstationClass,NhlDEFAULT_APP, rlist); } /* * Specify the viewport extent of the object. */ NhlRLClear(rlist); NhlRLSetFloat(rlist,NhlNvpXF,0.); NhlRLSetFloat(rlist,NhlNvpYF,1.); NhlRLSetFloat(rlist,NhlNvpWidthF,1.); NhlRLSetFloat(rlist,NhlNvpHeightF,1.); /* * Specify the type of markers for the legend. */ NhlRLSetInteger(rlist,NhlNlgItemCount,5); NhlRLSetInteger(rlist,NhlNlgMonoItemType,True); NhlRLSetInteger(rlist,NhlNlgItemType,NhlMARKERS); NhlRLSetStringArray(rlist,NhlNlgLabelStrings,labels,5); NhlRLSetIntegerArray(rlist,NhlNlgMarkerColors,colors,5); NhlRLSetIntegerArray(rlist,NhlNlgMarkerIndexes,item_ind,5); NhlRLSetString(rlist,NhlNlgMonoMarkerThickness,"False"); NhlRLSetFloatArray(rlist,NhlNlgMarkerThicknesses,mkthik,5); NhlRLSetFloat(rlist,NhlNlgMarkerSizeF,.05); NhlCreate(&pid,"Legend",NhllegendClass,wid,rlist); NhlDraw(pid); NhlFrame(wid); NhlDestroy(pid); NhlDestroy(wid); NhlDestroy(appid); NhlClose(); exit(0); }
void QNCLContourLine::plot(QGRIB2Message *grb2Msg, QString outputPath, QString outputFileName) { /*float xbvalues1[7] = {0.,60.,120.,180.,240.,300.,360.}; char *xblabels1[7] = {"0","60E","120E","180","120W","60W","0"}; float xbvalues2[11] = {0.,3.,6.,9.,12.,15.,18.,21.,24.,27.,30.}; float ylvalues[7] = {-90., -60.,-30.,0.,30.,60.,90.}; char *xblabels2[11] = {"0","3","6","9","12","15","18","21","24","27","30"}; char *yllabels[7] = {"90S","60S","30S","EQ","30N","60N","90N"};*/ /*rules setting*/ float xbValues[7] = {0}; char *xbLabels[7]; float ybValues[7] = {0}; char *ybLabels[7]; int xLabelCount = 0, yLabelCount = 0; double * dData; float * fData; //float cmap[255][3] = {0}; int ni, nj; int mi, mj; int xStart, xEnd; int yStart, yEnd; int appid, jan, cn, mp, tx, srlist; ng_size_t length[2]; float cmap[256][3]; //float * cmap; size_t nSize; ng_size_t icount[2]; int k = 0; int ps2;//1, ps2, png1; char setting[256]= {0}; /*x label setting*/ memset(xbLabels, 0, sizeof(xbLabels)); memset(ybLabels, 0, sizeof(ybLabels)); xLabelCount = width/100 + 1; xLabelCount = (xLabelCount < 2)?2:((xLabelCount>7)?7:xLabelCount); for (int i = 0; i < xLabelCount; i ++) { xbValues[i] = leftLongitude + (rightLongitude - leftLongitude) * i / (xLabelCount - 1); xbLabels[i] = (char *)calloc(10, 1); if ((xbValues[i] == 0) || (xbValues[i] == 180) || (xbValues[i]) == 360) sprintf(xbLabels[i], "%d", ((int)xbValues[i])%360); else if ((xbValues[i] > 0) && (xbValues[i] < 180)) { if ((xbValues[i]-(int)xbValues[i]) == 0) sprintf(xbLabels[i], "%dE", ((int)xbValues[i])%360); else sprintf(xbLabels[i], "%3.4fE", xbValues[i]); } else if ((xbValues[i] > 180) && (xbValues[i] < 360)) { if ((xbValues[i]-(int)xbValues[i]) == 0) sprintf(xbLabels[i], "%dW", (360 - (int)xbValues[i])); else sprintf(xbLabels[i], "%3.4fW", 360. - xbValues[i]); } } yLabelCount = height/100+1; yLabelCount = (yLabelCount < 2)?2:((yLabelCount>7)?7:yLabelCount); for (int i = 0; i < yLabelCount; i++) { ybValues[i] = topLatitude - (topLatitude - bottomLatitude) * i / (yLabelCount - 1); ybLabels[i] = (char *)calloc(10, 1); if (ybValues[i] == 0) sprintf(ybLabels[i], "%d", ((int)ybValues[i])%360); else if ((ybValues[i] > 0) && (ybValues[i] <= 90)) { if ((ybValues[i]-(int)ybValues[i]) == 0) sprintf(ybLabels[i], "%dN", ((int)ybValues[i])%360); else sprintf(ybLabels[i], "%3.4fN", ybValues[i]); } else if ((ybValues[i] < 0) && (ybValues[i] >= -90)) { if ((ybValues[i]-(int)ybValues[i]) == 0) sprintf(ybLabels[i], "%dW", ((int)ybValues[i])%360); else sprintf(ybLabels[i], "%3.4fW", ybValues[i]); } } /*get data from grib2 message between in the contour setting*/ grib_handle * hgrib = grb2Msg->getGRIBHandle(); QString workName; workName = outputFileName; workName = workName.left(workName.indexOf('.')); ni = grb2Msg->getGrid()->xp(); nj = grb2Msg->getGrid()->yp(); nSize = ni * nj; dData = (double *)calloc(ni*nj, sizeof(double)); grib_get_double_array(hgrib, "values", dData, &nSize); mi = ((rightLongitude - leftLongitude) * 1000000 / grb2Msg->getGrid()->dx()) + 1; mj = ((topLatitude - bottomLatitude) * 1000000 / grb2Msg->getGrid()->dy()) + 1; if (mi > ni) mi = ni; if (mj > nj) mj = nj; fData = (float *)calloc(mi * mj, sizeof(float)); xStart = (int)((leftLongitude * 1000000 - grb2Msg->getGrid()->lon1())/grb2Msg->getGrid()->dx()); if ((xStart * grb2Msg->getGrid()->dx()) > (leftLongitude * 1000000)) xStart --; if (xStart < 0) xStart = 0; xEnd = (int)((rightLongitude * 1000000 - grb2Msg->getGrid()->lon1())/grb2Msg->getGrid()->dx()) + 1; if ((xEnd * grb2Msg->getGrid()->dx()) < (rightLongitude * 1000000)) xEnd++; if (xEnd > ni) xEnd = ni; yStart = (int)((grb2Msg->getGrid()->lat1() - topLatitude * 1000000)/grb2Msg->getGrid()->dy()); if ((grb2Msg->getGrid()->lat1() - yStart * grb2Msg->getGrid()->dy()) < (topLatitude * 1000000)) yStart--; if (yStart < 0) yStart = 0; yEnd = (int)((grb2Msg->getGrid()->lat1() - bottomLatitude * 1000000)/grb2Msg->getGrid()->dy()) + 1; if ((grb2Msg->getGrid()->lat1() - yEnd * grb2Msg->getGrid()->dy()) > (bottomLatitude * 1000000)) yEnd++; if (yEnd > nj) yEnd = nj; for (int jj = yStart; jj < yEnd; jj++) { for (int ii = xStart; ii < xEnd; ii ++) { fData[k++] = dData[jj * ni + ii]; } } NhlInitialize(); srlist = NhlRLCreate(NhlSETRL); NhlRLClear(srlist); NhlRLSetString(srlist, "appDefaultParent", "True"); memset(setting, 0, sizeof(setting)); strncpy(setting, outputPath.toStdString().c_str(), sizeof(setting)-1); NhlRLSetString(srlist, "appUsrDir", setting); NhlCreate(&appid, grb2Msg->getName().toStdString().c_str(), NhlappClass, 0, srlist); NhlRLClear(srlist); QString outputpsFileName = outputPath + "/" + outputFileName; memset(setting, 0, sizeof(setting)); strncpy(setting, outputpsFileName.toStdString().c_str(), sizeof(setting)-1); NhlRLSetString(srlist, NhlNwkFileName, setting); /*NhlRLSetString(srlist, NhlNwkFormat, "ps");*/ NhlRLSetString(srlist, NhlNwkFormat, "png"); /*NhlRLSetString(srlist, NhlNwkOrientation, "landscape"); NhlRLSetString(srlist, NhlNwkOrientation, "portrait"); NhlRLSetInteger(srlist, NhlNwkDeviceLowerX, 0); NhlRLSetInteger(srlist, NhlNwkDeviceLowerY, 60); NhlRLSetInteger(srlist, NhlNwkDeviceUpperX, width); NhlRLSetInteger(srlist, NhlNwkDeviceUpperY, height); */ float fcolor[3] = {0}; float bcolor[3] = {0}; int r, g, b; foregroundColor.getRgb(&r, &g, &b); fcolor[0] = ((float)r)/255.; fcolor[1] = ((float)g)/255.; fcolor[2] = ((float)b)/255.; length[0] = 3; NhlRLSetFloatArray(srlist, NhlNwkForegroundColor, fcolor, length[0]); backgroundColor.getRgb(&r, &g, &b); bcolor[0] = ((float)r)/255.; bcolor[1] = ((float)g)/255.; bcolor[2] = ((float)b)/255.; length[0] = 3; NhlRLSetFloatArray(srlist, NhlNwkBackgroundColor, bcolor, length[0]); length[0] = loadColormap(colormap, &cmap[0][0], 256); length[1] = 3; if (length[0] > 0) { NhlRLSetMDFloatArray(srlist,NhlNwkColorMap,&cmap[0][0], 2, length); } else { memset(setting, 0, sizeof(setting)); strncpy(setting, colormap.toStdString().c_str(), sizeof(setting)-1); NhlRLSetString(srlist, NhlNwkColorMap, setting); } /*NhlCreate(&ps2, workName.toStdString().c_str(), NhlpsWorkstationClass, 0, srlist);*/ NhlCreate(&ps2, workName.toStdString().c_str(), NhlcairoImageWorkstationClass, 0, srlist); icount[1] = xEnd - xStart; icount[0] = yEnd - yStart; /**/ NhlRLClear(srlist); NhlRLSetMDFloatArray(srlist, NhlNsfDataArray, fData, 2, icount); NhlRLSetFloat(srlist, NhlNsfXCStartV, leftLongitude); NhlRLSetFloat(srlist, NhlNsfXCEndV, rightLongitude); NhlRLSetFloat(srlist, NhlNsfYCStartV, topLatitude); NhlRLSetFloat(srlist, NhlNsfYCEndV, bottomLatitude); NhlCreate(&jan, "sf", NhlscalarFieldClass, appid, srlist); //NhlRLSetString(srlist, NhlNwkOrientation, "landscape"); /*NhlRLSetString(srlist, NhlNwkOrientation, "portrait");*/ NhlRLSetString(srlist, NhlNcnLineDrawOrder, "predraw"); /** Create a ContourPlot object.*/ NhlRLClear(srlist); NhlRLSetInteger(srlist, NhlNcnScalarFieldData, jan); NhlRLSetFloat(srlist, NhlNvpXF, .1); NhlRLSetFloat(srlist, NhlNvpYF, .75); NhlRLSetFloat(srlist, NhlNvpWidthF, 0.8); NhlRLSetFloat(srlist, NhlNvpHeightF, 0.4); NhlRLSetString(srlist, NhlNcnInfoLabelOn, "False"); NhlRLSetString(srlist, NhlNcnHighLabelsOn, "False"); NhlRLSetString(srlist, NhlNcnLowLabelsOn, "False"); //NhlRLSetIntegerArray(srlist, NhlNcnLineColors, linecolors, 15); NhlRLSetString(srlist, NhlNcnLineDrawOrder, "predraw"); NhlRLSetString(srlist, NhlNcnFillDrawOrder, "predraw"); NhlRLSetString(srlist, NhlNcnLabelDrawOrder, "predraw"); memset(setting, 0, sizeof(setting)); strncpy(setting, levelSelectionMode.toStdString().c_str(), sizeof(setting)-1); NhlRLSetString(srlist, NhlNcnLevelSelectionMode, setting); NhlRLSetInteger(srlist, NhlNcnMaxLevelCount, maxLevelCount); if (levelSelectionMode == "ManualLevels") { NhlRLSetFloat(srlist, NhlNcnMinLevelValF, minLevel); NhlRLSetFloat(srlist, NhlNcnMaxLevelValF, maxLevel); NhlRLSetFloat(srlist, NhlNcnLevelSpacingF, levelSpacing); } else if (levelSelectionMode == "ExplicitLevels") { length[0] = maxLevelCount; NhlRLSetFloatArray(srlist, NhlNcnLevelSpacingF, lineLevels, length[0]); } else if (levelSelectionMode == "EqualSpacedLevels") { NhlRLSetFloat(srlist, NhlNcnLevelSpacingF, levelSpacing); } if (plotMode == 0) { NhlRLSetString(srlist, NhlNcnFillOn, "False"); NhlRLSetString(srlist, NhlNcnLinesOn, "True"); NhlRLSetString(srlist, NhlNcnLineLabelsOn, "True"); NhlRLSetInteger(srlist, NhlNcnLineLabelInterval, 2); NhlRLSetString(srlist, NhlNcnLineLabelPlacementMode, "computed"); if (monoColor == false) { NhlRLSetString(srlist, NhlNcnMonoLineColor, "False"); NhlRLSetIntegerArray(srlist, NhlNcnLineColors, lineColors, maxLevelCount); } else { NhlRLSetString(srlist, NhlNcnMonoLineColor, "True"); /*NhlRLSetString(srlist, NhlNcnLineColor, "Black");*/ NhlRLSetInteger(srlist, NhlNcnLineColor, lineColor); } } else { NhlRLSetString(srlist, NhlNcnFillOn, "True"); NhlRLSetString(srlist, NhlNcnLinesOn, "False"); NhlRLSetString(srlist, NhlNcnFillMode, "AreaFill"); NhlRLSetString(srlist, NhlNcnLineLabelsOn, "False"); length[0] = maxLevelCount; NhlRLSetIntegerArray(srlist, NhlNcnFillColors, lineColors, length[0]); } NhlRLSetString(srlist, NhlNtmXBMode, "EXPLICIT"); NhlRLSetFloatArray(srlist, NhlNtmXBValues, xbValues, xLabelCount); NhlRLSetStringArray(srlist, NhlNtmXBLabels, &xbLabels[0], xLabelCount); NhlRLSetString(srlist, NhlNtmYLMode, "EXPLICIT"); NhlRLSetFloatArray(srlist, NhlNtmYLValues, ybValues, yLabelCount); NhlRLSetStringArray(srlist, NhlNtmYLLabels, &ybLabels[0], yLabelCount); NhlRLSetString(srlist, NhlNtmXTLabelsOn, "False"); NhlRLSetString(srlist, NhlNtmYRLabelsOn, "False"); NhlRLSetFloat(srlist, NhlNtmXBLabelFontHeightF, .010); NhlRLSetFloat(srlist, NhlNtmYLLabelFontHeightF, .010); NhlRLSetFloat(srlist, NhlNtmXBMajorOutwardLengthF, .006); NhlRLSetFloat(srlist, NhlNtmXBMajorLengthF, .006); NhlRLSetFloat(srlist, NhlNtmXTMajorLengthF, 0.); NhlRLSetFloat(srlist, NhlNtmXTMajorOutwardLengthF, 0.); NhlRLSetFloat(srlist, NhlNtmYLMajorOutwardLengthF, .006); NhlRLSetFloat(srlist, NhlNtmYLMajorLengthF, .006); NhlRLSetString(srlist, NhlNtmXBMinorOn, "False"); NhlRLSetString(srlist, NhlNtmXTMinorOn, "False"); NhlRLSetString(srlist, NhlNtmYLMinorOn, "False"); NhlRLSetString(srlist, NhlNtmYRMinorOn, "False"); if (monoThickness == true) { NhlRLSetString(srlist, NhlNcnMonoLineThickness, "True"); NhlRLSetFloat(srlist, NhlNcnLineThicknessF, thickness); } else { NhlRLSetString(srlist, NhlNcnMonoLineThickness, "False"); NhlRLSetFloatArray(srlist, NhlNcnLineThicknesses, lineThicknesses, maxLevelCount); } NhlCreate(&cn, "cn", NhlcontourPlotClass, ps2, srlist); /* Create a MapPlot object. */ NhlRLClear(srlist); NhlRLSetFloat(srlist, NhlNvpXF, .1); NhlRLSetFloat(srlist, NhlNvpYF, .75); NhlRLSetFloat(srlist, NhlNvpWidthF, 0.8); NhlRLSetFloat(srlist, NhlNvpHeightF, 0.4); /*NhlRLSetString(srlist, NhlNmpFillOn, "True");*/ NhlRLSetString(srlist, NhlNmpFillOn, "False"); NhlRLSetString(srlist, NhlNmpLabelsOn, "False"); /*NhlRLSetString(srlist, NhlNmpGeophysicalLineColor, "Black"); NhlRLSetString(srlist, NhlNmpDefaultFillColor, "DarkSalmon"); NhlRLSetString(srlist, NhlNmpLandFillColor, "DarkSalmon"); NhlRLSetString(srlist, NhlNmpOceanFillColor, "Blue"); NhlRLSetString(srlist, NhlNmpInlandWaterFillColor, "Blue");*/ NhlRLSetString(srlist, NhlNmpOutlineOn, "True"); NhlRLSetString(srlist, NhlNgsLineColor, "Black"); /*NhlRLSetString(srlist, NhlNmpOutlineOn, "True"); NhlRLSetString(srlist, NhlNmpOutlineBoundarySets, "Geophysical");*/ /*NhlRLSetString(srlist, NhlNmpOutlineDrawOrder, "predraw");*/ NhlRLSetString(srlist, NhlNmpOutlineBoundarySets, "Geophysical"); NhlRLSetFloat(srlist, NhlNmpGeophysicalLineThicknessF, 1); NhlRLSetString(srlist, NhlNmpGeophysicalLineColor, "Black"); /*NhlRLSetString(srlist, NhlNmpPerimDrawOrder, "predraw");*/ memset(setting, 0, sizeof(setting)); strncpy(setting, mapProject.toStdString().c_str(), sizeof(setting)-1); if (mapProject == "Polar Stereographic") { NhlRLSetString(srlist, NhlNmpProjection, "Stereographic"); if (polarPosition == "NH") { NhlRLSetFloat(srlist, NhlNmpMinLatF, 0); NhlRLSetFloat(srlist, NhlNmpMaxLatF, 90); NhlRLSetFloat(srlist, NhlNmpMinLonF, leftLongitude); NhlRLSetFloat(srlist, NhlNmpMaxLonF, rightLongitude); NhlRLSetFloat(srlist, NhlNmpCenterLonF, (leftLongitude+rightLongitude)/2); NhlRLSetFloat(srlist, NhlNmpCenterLatF, 90); } else { NhlRLSetFloat(srlist, NhlNmpMinLatF, -90); NhlRLSetFloat(srlist, NhlNmpMaxLatF, 0); NhlRLSetFloat(srlist, NhlNmpMinLonF, leftLongitude); NhlRLSetFloat(srlist, NhlNmpMaxLonF, rightLongitude); NhlRLSetFloat(srlist, NhlNmpCenterLonF, (leftLongitude+rightLongitude)/2); NhlRLSetFloat(srlist, NhlNmpCenterLatF, -90); } } else { NhlRLSetString(srlist, NhlNmpProjection, setting); NhlRLSetFloat(srlist, NhlNmpMinLatF, bottomLatitude); NhlRLSetFloat(srlist, NhlNmpMaxLatF, topLatitude); NhlRLSetFloat(srlist, NhlNmpMinLonF, leftLongitude); NhlRLSetFloat(srlist, NhlNmpMaxLonF, rightLongitude); NhlRLSetFloat(srlist, NhlNmpCenterLonF, (leftLongitude+rightLongitude)/2); NhlRLSetFloat(srlist, NhlNmpCenterLatF, (topLatitude + bottomLatitude)/2); } NhlRLSetString(srlist, NhlNmpGridAndLimbOn, "False"); NhlRLSetString(srlist, NhlNmpLimitMode, "latlon"); NhlCreate(&mp, "mp", NhlmapPlotClass, ps2, srlist); /* Create a TextItem object. */ QString title = grb2Msg->getParameter()->nclName(); QString level; level.setNum(grb2Msg->getLevel()->firstLevel()); switch (grb2Msg->getLevel()->firstLevelType()) { case 100: case 101: level += "Pa"; break; case 0: case 106: level += "m"; break; } title += " " + level + " " + grb2Msg->getTimeLevel()->getName() ; NhlRLClear(srlist); NhlRLSetFloat(srlist, NhlNtxPosXF, 0.5); NhlRLSetFloat(srlist, NhlNtxPosYF, 0.8); NhlRLSetString(srlist, NhlNtxJust, "CENTERCENTER"); memset(setting, 0, sizeof(setting)); strncpy(setting, title.toStdString().c_str(), sizeof(setting)-1); NhlRLSetString(srlist, NhlNtxString, setting); NhlRLSetFloat(srlist, NhlNtxFontHeightF, .030); NhlRLSetInteger(srlist, NhlNtxFont, 25); NhlCreate(&tx, "tx", NhltextItemClass, ps2, srlist); NhlDraw(mp); NhlDraw(cn); NhlDraw(tx); NhlFrame(ps2); NhlDestroy(ps2); NhlDestroy(appid); for (int idx = 0; idx < xLabelCount; idx ++) { free(xbLabels[idx]); } for (int idx = 0; idx < yLabelCount; idx ++) { free(ybLabels[idx]); } free(dData); free(fData); }