TOSSMTimeValue* LoadTOSSMTimeValue(TMover *theOwner, short unitsIfKnownInAdvance)
{
char path[256],shortFileName[256];
char tempStr[256];
Point where = CenteredDialogUpLeft(M38d);
WorldPoint p;
OSType typeList[] = { 'NULL', 'NULL', 'NULL', 'NULL' };
MySFReply reply;
OSErr err = 0;
#if TARGET_API_MAC_CARBON
mysfpgetfile(&where, "", -1, typeList,
(MyDlgHookUPP)0, &reply, M38d, MakeModalFilterUPP(STDFilter));
if (!reply.good) return 0;
strcpy(path, reply.fullPath);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#else
sfpgetfile(&where, "",
(FileFilterUPP)0,
-1, typeList,
(DlgHookUPP)0,
&reply, M38d,
(ModalFilterUPP)MakeUPP((ProcPtr)STDFilter, uppModalFilterProcInfo));
if (!reply.good) return nil; // user canceled
my_p2cstr(reply.fName);
#ifdef MAC
GetFullPath(reply.vRefNum, 0, (char *)reply.fName, path);
strcpy(shortFileName,(char*) reply.fName);
#else
strcpy(path, reply.fName);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#endif
#endif
// return CreateTOSSMTimeValue(theOwner,path,shortFileName,kUndefined); // ask user for units
return CreateTOSSMTimeValue(theOwner,path,shortFileName,unitsIfKnownInAdvance); // ask user for units
}
ADCPTimeValue* ADCPMover_c::AddADCP(OSErr *err)
{ // might send in path for first adcp
*err = 0;
char tempStr[128], shortFileName[64], givenPath[256], givenFileName[64], s[256], fileName[64], path[256];
short unitsIfKnownInAdvance = 0;
Boolean askForFile = true;
ADCPTimeValue *timeValObj = 0;
Point where = CenteredDialogUpLeft(M38c);;
OSType typeList[] = { 'NULL', 'NULL', 'NULL', 'NULL' };
MySFReply reply;
//if(askForFile || !givenPath || !givenFileName)
{
#if TARGET_API_MAC_CARBON
mysfpgetfile(&where, "", -1, typeList,
(MyDlgHookUPP)0, &reply, M38c, MakeModalFilterUPP(STDFilter));
if (!reply.good) return 0;
strcpy(path, reply.fullPath);
#else
sfpgetfile(&where, "",
(FileFilterUPP)0,
-1, typeList,
(DlgHookUPP)0,
&reply, M38c,
(ModalFilterUPP)MakeUPP((ProcPtr)STDFilter, uppModalFilterProcInfo));
if (!reply.good) return 0;
my_p2cstr(reply.fName);
#ifdef MAC
GetFullPath(reply.vRefNum, 0, (char *)reply.fName, path);
#else
strcpy(path, reply.fName);
#endif
#endif
strcpy (s, path);
SplitPathFile (s, fileName);
}
/*else
{ // don't ask user, we were provided with the path
strcpy(path,givenPath);
strcpy(fileName,givenFileName);
}*/
if (IsADCPFile(path))
{
timeValObj = new ADCPTimeValue(dynamic_cast<ADCPMover *>(this));
//timeDep = new ADCPTimeValue(this);
if (!timeValObj)
{ TechError("TextRead()", "new ADCPTimeValue()", 0); return nil; }
*err = timeValObj->InitTimeFunc();
if(*err) {delete timeValObj; timeValObj = nil; return nil;}
*err = timeValObj->ReadTimeValues (path, M19REALREAL, unitsIfKnownInAdvance);
if(*err) { delete timeValObj; timeValObj = nil; return nil;}
timeValObj->SetTimeFileName(fileName);
//return timeValObj;
//AddTimeDep(timeValObj,0);
}
// code goes here, add code for OSSMHeightFiles, need scale factor to calculate derivative
else
{
sprintf(tempStr,"File %s is not a recognizable ADCP time file.",shortFileName);
printError(tempStr);
}
return timeValObj;
}
OSErr AddMapsDialog2()
{
char path[256], nameStr [256], shortFileName[256], tempStr[256];
OSErr err = noErr;
long n;
Point where = CenteredDialogUpLeft(M38b);
TVectorMap *vMap;
TOSSMMap *oMap;
Map3D *gMap;
GridMap_c *gridMap = 0;
OSType typeList[] = { 'NULL', 'NULL', 'NULL', 'NULL' };
MySFReply reply;
WorldRect theRect = emptyWorldRect;
short gridType;
Boolean isESI = false;
float arrowDepth = 0;
#if TARGET_API_MAC_CARBON
mysfpgetfile(&where, "", -1, typeList,
(MyDlgHookUPP)0, &reply, M38b, MakeModalFilterUPP(STDFilter));
if (!reply.good) return USERCANCEL;
strcpy(path, reply.fullPath);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#else
sfpgetfile(&where, "",
(FileFilterUPP)0,
-1, typeList,
(DlgHookUPP)0,
&reply, M38b,
(ModalFilterUPP)MakeUPP((ProcPtr)STDFilter, uppModalFilterProcInfo));
if (!reply.good) return USERCANCEL;
my_p2cstr(reply.fName);
#ifdef MAC
GetFullPath(reply.vRefNum, 0, (char *)reply.fName, path);
strcpy(shortFileName,(char*) reply.fName);
#else
strcpy(path, reply.fName);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#endif
#endif
if (IsVectorMap (path, &isESI))
{
if (isESI) {printNote("File is ESI segments not map polygons");err=-1; return err;}
strcpy (nameStr, "Vector Map: ");
strcat (nameStr, shortFileName);
vMap = (TVectorMap*) new TVectorMap (nameStr, theRect);
if (!vMap)
{ TechError("AddMapsDialog()", "new TVectorMap()", 0); return -1; }
if (err = vMap -> InitMap(path)) { delete vMap; return err; }
if (err = model->AddMap(vMap, 0))
{
vMap -> Dispose ();
delete vMap;
err = -1;
}
}
else if (IsGridMap (path))
{
strcpy (nameStr, "Grid Map: ");
strcat (nameStr, shortFileName);
oMap = new TOSSMMap(nameStr, voidWorldRect);
if (!oMap)
{ TechError("AddMapsDialog()", "new TOSSMMap()", 0); return -1; }
if (err = ((TOSSMMap*) oMap) -> InitMap(path)) { delete oMap; return err; }
if (err = model->AddMap(oMap, 0))
{
oMap -> Dispose ();
delete oMap;
err = -1;
}
}
else if (IsPtCurFile (path))
{
TMap *newMap = 0;
TCurrentMover *newMover = CreateAndInitCurrentsMover (model->uMap,false,path,"ptcurfile",&newMap); // already have path
if (newMover)
{
PtCurMover *ptCurMover = dynamic_cast<PtCurMover*>(newMover);
err = ptCurMover -> SettingsDialog();
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (!err) err = AddMoverToMap(newMap, timeFileChanged, newMover);
//if(!err) err = ((PtCurMap*)newMap)->MakeBitmaps();
if (!err) newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
else
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
}
}
}
}
else if (IsNetCDFFile (path, &gridType))
{
TMap *newMap = 0;
char s[256],fileName[256], outPath[256];
WorldRect bounds = theWorld;
strcpy(s,path);
SplitPathFile (s, fileName);
strcat (nameStr, fileName);
gMap = CreateAndInitMap3D (path,bounds); // need to fix bounds afterwards
if (gMap && gridType!=REGULAR && gridType!=REGULAR_SWAFS)
{
FLOATH depthPtsH = 0;
WORLDPOINTFH vertexPtsH=0;
long numRows=0, numCols=0, numNodes=0, numTri=0, numBoundaryPts=0;
char errmsg[256];
gridMap = new GridMap_c();
#if TARGET_API_MAC_CARBON
err = ConvertTraditionalPathToUnixPath((const char *) path, outPath, kMaxNameLen) ;
//if (!err) strcpy(path,outPath);
#else
strcpy(outPath,path);
#endif
if (gridType == CURVILINEAR)
{
DOUBLEH maskH = 0;
//err = gMap->GetPointsAndMask(path,&maskH,&vertexPtsH,&depthPtsH,&numRows, &numCols); //Text read
//if (!err) err = gMap->SetUpCurvilinearGrid(maskH,numRows,numCols,vertexPtsH,depthPtsH,errmsg); //Reorder points
err = gridMap->GetPointsAndMask(outPath,&maskH,&vertexPtsH,&depthPtsH,&numRows, &numCols); //Text read
if (!err) err = gridMap->SetUpCurvilinearGrid(maskH,numRows,numCols,vertexPtsH,depthPtsH,errmsg); //Reorder points
if(maskH) {DisposeHandle((Handle)maskH); maskH = 0;}
}
else if (gridType == TRIANGULAR)
{ // check if topology is included
LONGPTR bndry_indices=0, bndry_nums=0, bndry_type=0, tri_verts=0, tri_neighbors=0;
//err = gMap->GetPointsAndBoundary(path,&vertexPtsH,&depthPtsH, &numNodes, &bndry_indices, &bndry_nums, &bndry_type, &numBoundaryPts, &tri_verts, &tri_neighbors, &numTri); //Text read
err = gridMap->GetPointsAndBoundary(outPath,&vertexPtsH,&depthPtsH, &numNodes, &bndry_indices, &bndry_nums, &bndry_type, &numBoundaryPts, &tri_verts, &tri_neighbors, &numTri); //Text read
if (!err)
{ // separate points and boundary
if (numTri == 0)
//err = gMap->SetUpTriangleGrid(numNodes,numTri,vertexPtsH,depthPtsH, bndry_indices, bndry_nums, bndry_type, numBoundaryPts); //Reorder points
err = gridMap->SetUpTriangleGrid(numNodes,numTri,vertexPtsH,depthPtsH, bndry_indices, bndry_nums, bndry_type, numBoundaryPts); //Reorder points
else
//err = gMap->SetUpTriangleGrid2(numNodes,numTri,vertexPtsH,depthPtsH, bndry_indices, bndry_nums, bndry_type, numBoundaryPts,tri_verts, tri_neighbors); //Reorder points
err = gridMap->SetUpTriangleGrid2(numNodes,numTri,vertexPtsH,depthPtsH, bndry_indices, bndry_nums, bndry_type, numBoundaryPts,tri_verts, tri_neighbors); //Reorder points
}
if (bndry_indices) delete [] bndry_indices;
if (bndry_nums) delete [] bndry_nums;
if (bndry_type) delete [] bndry_type;
if (tri_verts) delete [] tri_verts;
if (tri_neighbors) delete [] tri_neighbors;
}
if(vertexPtsH) {DisposeHandle((Handle)vertexPtsH); vertexPtsH = 0;}
if(depthPtsH) {DisposeHandle((Handle)depthPtsH); depthPtsH = 0;}
if( !err)
{
gMap->SetBoundarySegs(gridMap->GetBoundarySegs());
gMap->SetWaterBoundaries(gridMap->GetWaterBoundaries());
gMap->SetBoundaryPoints(gridMap->GetBoundaryPoints());
gMap->SetGrid(gridMap->GetGrid());
gMap->SetMapBounds(gridMap->GetMapBounds());
err = gMap->MakeBitmaps();
}
if (!err) err = model->AddMap(gMap, 0);
if (err)
{
gMap -> Dispose ();
delete gMap;
err = -1;
}
}
else
{
err = true;
sprintf(tempStr,"File %s is a current file and should be input as a universal mover.",shortFileName);
printNote(tempStr);
}
}
else if (IsTriCurFile (path))
{
TMap *newMap = 0;
TCurrentMover *newMover = CreateAndInitCurrentsMover (model->uMap,false,path,"TriCurFile",&newMap); // already have path
if (newMover)
{
TriCurMover *triCurMover = dynamic_cast<TriCurMover *>(newMover);
err = triCurMover -> SettingsDialog();
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
err = AddMoverToMap(newMap, timeFileChanged, newMover);
if(err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
/*err = ((PtCurMap*)newMap)->MakeBitmaps();
if(err)
{
newMap->Dispose(); delete newMap; newMap = 0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}*/
newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
}
}
}
else
{
err = true;
sprintf(tempStr,"File %s is a current file and should be input as a universal mover.",shortFileName);
printNote(tempStr);
}
}
else if (IsCATS3DFile (path)) // for any CATS?
{
char s[256],fileName[256];
WorldRect bounds = theWorld;
strcpy(s,path);
SplitPathFile (s, fileName);
strcat (nameStr, fileName);
gMap = CreateAndInitMap3D (path,bounds); // need to fix bounds afterwards
if (gMap)
{
gridMap = new GridMap_c();
if (!gridMap) {err = gMap->ReadCATSMap(path);/*err = -1;*/}
else
//err = gMap->ReadCATSMap(path);
err = gridMap->ReadCATSMap(path);
if(!err)
{
if (gridMap)
{
gMap->SetBoundarySegs(gridMap->GetBoundarySegs());
gMap->SetWaterBoundaries(gridMap->GetWaterBoundaries());
gMap->SetBoundaryPoints(gridMap->GetBoundaryPoints());
gMap->SetGrid(gridMap->GetGrid());
gMap->SetMapBounds(gridMap->GetMapBounds());
}
err = gMap->MakeBitmaps();
}
if (err = model->AddMap(gMap, 0))
{
gMap -> Dispose ();
delete gMap;
err = -1;
}
}
}
else
{
WorldRect bounds = theWorld;
gMap = CreateAndInitMap3D (path,bounds); // need to fix bounds afterwards
if (gMap)
{
gridMap = new GridMap_c();
if (!gridMap) {err = gMap -> ReadTopology(path);/*err = -1;*/}
else
//err = gMap -> ReadTopology(path);
err = gridMap -> ReadTopology(path);
if(!err)
{
if (gridMap)
{
gMap->SetBoundarySegs(gridMap->GetBoundarySegs());
gMap->SetWaterBoundaries(gridMap->GetWaterBoundaries());
gMap->SetBoundaryPoints(gridMap->GetBoundaryPoints());
gMap->SetGrid(gridMap->GetGrid());
gMap->SetMapBounds(gridMap->GetMapBounds());
}
err = gMap->MakeBitmaps();
}
if (err = model->AddMap(gMap, 0))
{
gMap -> Dispose ();
delete gMap;
err = -1;
}
}
//err = true;
if (err)
{sprintf(tempStr,"File %s is not a recognizable map file.",shortFileName);
printError(tempStr);}
}
if (!err)
{
model->NewDirtNotification();
}
return err;
}
OSErr AddMapsDialog()
{
char path[256], nameStr [256], shortFileName[256], tempStr[256];
OSErr err = noErr;
long n;
Point where = CenteredDialogUpLeft(M38b);
TVectorMap *vMap;
TOSSMMap *oMap;
OSType typeList[] = { 'NULL', 'NULL', 'NULL', 'NULL' };
MySFReply reply;
WorldRect theRect = emptyWorldRect;
short gridType;
Boolean isESI = false;
float arrowDepth = 0;
#if TARGET_API_MAC_CARBON
mysfpgetfile(&where, "", -1, typeList,
(MyDlgHookUPP)0, &reply, M38b, MakeModalFilterUPP(STDFilter));
if (!reply.good) return USERCANCEL;
strcpy(path, reply.fullPath);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#else
sfpgetfile(&where, "",
(FileFilterUPP)0,
-1, typeList,
(DlgHookUPP)0,
&reply, M38b,
(ModalFilterUPP)MakeUPP((ProcPtr)STDFilter, uppModalFilterProcInfo));
if (!reply.good) return USERCANCEL;
my_p2cstr(reply.fName);
#ifdef MAC
GetFullPath(reply.vRefNum, 0, (char *)reply.fName, path);
strcpy(shortFileName,(char*) reply.fName);
#else
strcpy(path, reply.fName);
strcpy(tempStr,path);
SplitPathFile(tempStr,shortFileName);
#endif
#endif
if (IsVectorMap (path, &isESI))
{
if (isESI) {printNote("File is ESI segments not map polygons");err=-1; return err;}
strcpy (nameStr, "Vector Map: ");
strcat (nameStr, shortFileName);
vMap = (TVectorMap*) new TVectorMap (nameStr, theRect);
if (!vMap)
{ TechError("AddMapsDialog()", "new TVectorMap()", 0); return -1; }
if (err = vMap -> InitMap(path)) { delete vMap; return err; }
if (err = model->AddMap(vMap, 0))
{
vMap -> Dispose ();
delete vMap;
err = -1;
}
}
else if (IsGridMap (path))
{
strcpy (nameStr, "Grid Map: ");
strcat (nameStr, shortFileName);
oMap = new TOSSMMap(nameStr, voidWorldRect);
if (!oMap)
{ TechError("AddMapsDialog()", "new TOSSMMap()", 0); return -1; }
if (err = ((TOSSMMap*) oMap) -> InitMap(path)) { delete oMap; return err; }
if (err = model->AddMap(oMap, 0))
{
oMap -> Dispose ();
delete oMap;
err = -1;
}
}
else if (IsPtCurFile (path))
{
TMap *newMap = 0;
TCurrentMover *newMover = CreateAndInitCurrentsMover (model->uMap,false,path,"ptcurfile",&newMap); // already have path
if (newMover)
{
PtCurMover *ptCurMover = dynamic_cast<PtCurMover*>(newMover);
err = ptCurMover -> SettingsDialog();
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (!err) err = AddMoverToMap(newMap, timeFileChanged, newMover);
//if(!err) err = ((PtCurMap*)newMap)->MakeBitmaps();
if (!err) newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
else
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
}
}
}
}
else if (IsNetCDFFile (path, &gridType))
{
TMap *newMap = 0;
TCurrentMover *newMover = CreateAndInitCurrentsMover (model->uMap,false,path,"NetCDFfile",&newMap); // already have path
if (newMover && gridType!=REGULAR && gridType!=REGULAR_SWAFS) // should probably get rid of the SWAFS grid type and generalize a regular navy
{
NetCDFMover *netCDFMover = dynamic_cast<NetCDFMover *>(newMover);
err = netCDFMover -> SettingsDialog();
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (!err) err = AddMoverToMap(newMap, timeFileChanged, newMover);
// if (!err) err = ((PtCurMap*)newMap)->MakeBitmaps();
if (!err) newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
else
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
}
}
}
else
{
err = true;
sprintf(tempStr,"File %s is a current file and should be input as a universal mover.",shortFileName);
printNote(tempStr);
}
}
else if (IsTriCurFile (path))
{
TMap *newMap = 0;
TCurrentMover *newMover = CreateAndInitCurrentsMover (model->uMap,false,path,"TriCurFile",&newMap); // already have path
if (newMover)
{
TriCurMover *triCurMover = dynamic_cast<TriCurMover *>(newMover);
err = triCurMover -> SettingsDialog();
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
err = AddMoverToMap(newMap, timeFileChanged, newMover);
if(err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
/*err = ((PtCurMap*)newMap)->MakeBitmaps();
if(err)
{
newMap->Dispose(); delete newMap; newMap = 0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}*/
newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
}
}
}
else
{
err = true;
sprintf(tempStr,"File %s is a current file and should be input as a universal mover.",shortFileName);
printNote(tempStr);
}
}
else if (IsCATS3DFile (path))
{
TMap *newMap = 0;
char s[256],fileName[256];
TCurrentMover *newMover = 0;
Boolean timeFileChanged = false;
strcpy(s,path);
SplitPathFile (s, fileName);
strcat (nameStr, fileName);
newMover = CreateAndInitCurrentsMover (model->uMap,false,path,fileName,&newMap); // already have path
if (newMover)
{
TCATSMover3D *catsMover3D = dynamic_cast<TCATSMover3D *>(newMover);
err = CATSSettingsDialog (dynamic_cast<TCATSMover *>(newMover), model->uMap, &timeFileChanged);
if(err)
{
newMover->Dispose(); delete newMover; newMover = 0;
if (newMap) {newMap->Dispose(); delete newMap; newMap = 0;}
}
if(newMover && !err)
{
Boolean timeFileChanged = false;
if (!newMap)
{
err = AddMoverToMap (model->uMap, timeFileChanged, newMover);
}
else
{
err = model -> AddMap(newMap, 0);
if (err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
err = AddMoverToMap(newMap, timeFileChanged, newMover);
if(err)
{
newMap->Dispose(); delete newMap; newMap =0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}
/*err = ((PtCurMap*)newMap)->MakeBitmaps();
if(err)
{
newMap->Dispose(); delete newMap; newMap = 0;
newMover->Dispose(); delete newMover; newMover = 0;
return -1;
}*/
newMover->SetMoverMap(newMap);
if (model->ThereIsA3DMover(&arrowDepth)) InitAnalysisMenu(); // want to have it come and go?
}
}
}
}
else
{
err = true;
sprintf(tempStr,"File %s is not a recognizable map file.",shortFileName);
printError(tempStr);
}
if (!err)
{
model->NewDirtNotification();
}
return err;
}
OSErr NetCDFWindMoverCurv::TextRead(char *path, TMap **newMap, char *topFilePath) // don't want a map
{
// this code is for curvilinear grids
OSErr err = 0;
long i,j, numScanned, indexOfStart = 0;
int status, ncid, latIndexid, lonIndexid, latid, lonid, recid, timeid, numdims;
size_t latLength, lonLength, recs, t_len, t_len2;
float timeVal;
char recname[NC_MAX_NAME], *timeUnits=0, month[10];
char dimname[NC_MAX_NAME], s[256], topPath[256];
WORLDPOINTFH vertexPtsH=0;
float *lat_vals=0,*lon_vals=0,yearShift=0.;
static size_t timeIndex,ptIndex[2]={0,0};
static size_t pt_count[2];
Seconds startTime, startTime2;
double timeConversion = 1.;
char errmsg[256] = "",className[256]="";
char fileName[64],*modelTypeStr=0;
Point where;
OSType typeList[] = { 'NULL', 'NULL', 'NULL', 'NULL' };
MySFReply reply;
Boolean bTopFile = false, fIsNavy = false; // for now keep code around but probably don't need Navy curvilinear wind
//VelocityFH velocityH = 0;
char outPath[256];
if (!path || !path[0]) return 0;
strcpy(fPathName,path);
strcpy(s,path);
SplitPathFile (s, fileName);
strcpy(fFileName, fileName); // maybe use a name from the file
status = nc_open(path, NC_NOWRITE, &ncid);
//if (status != NC_NOERR) {err = -1; goto done;}
if (status != NC_NOERR)
{
#if TARGET_API_MAC_CARBON
err = ConvertTraditionalPathToUnixPath((const char *) path, outPath, kMaxNameLen) ;
status = nc_open(outPath, NC_NOWRITE, &ncid);
#endif
if (status != NC_NOERR) {err = -1; goto done;}
}
// check number of dimensions - 2D or 3D
status = nc_inq_ndims(ncid, &numdims);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_attlen(ncid,NC_GLOBAL,"generating_model",&t_len2);
if (status != NC_NOERR) {fIsNavy = false; /*goto done;*/}
else
{
fIsNavy = true;
// may only need to see keyword is there, since already checked grid type
modelTypeStr = new char[t_len2+1];
status = nc_get_att_text(ncid, NC_GLOBAL, "generating_model", modelTypeStr);
if (status != NC_NOERR) {fIsNavy = false; goto done;}
modelTypeStr[t_len2] = '\0';
strcpy(fFileName, modelTypeStr);
}
GetClassName(className);
if (!strcmp("NetCDF Wind",className))
SetClassName(fFileName); //first check that name is now the default and not set by command file ("NetCDF Wind")
//if (fIsNavy)
{
status = nc_inq_dimid(ncid, "time", &recid); //Navy
//if (status != NC_NOERR) {err = -1; goto done;}
if (status != NC_NOERR)
{ status = nc_inq_unlimdim(ncid, &recid); // issue of time not being unlimited dimension
if (status != NC_NOERR) {err = -1; goto done;}
}
}
/*else
{
status = nc_inq_unlimdim(ncid, &recid); // issue of time not being unlimited dimension
if (status != NC_NOERR) {err = -1; goto done;}
}*/
//if (fIsNavy)
status = nc_inq_varid(ncid, "time", &timeid);
if (status != NC_NOERR)
{
status = nc_inq_varid(ncid, "ProjectionHr", &timeid);
if (status != NC_NOERR) {err = -1; goto done;}
}
// if (status != NC_NOERR) {/*err = -1; goto done;*/timeid=recid;}
//if (!fIsNavy)
//status = nc_inq_attlen(ncid, recid, "units", &t_len); // recid is the dimension id not the variable id
//else // LAS has them in order, and time is unlimited, but variable/dimension names keep changing so leave this way for now
status = nc_inq_attlen(ncid, timeid, "units", &t_len);
if (status != NC_NOERR)
{
timeUnits = 0; // files should always have this info
timeConversion = 3600.; // default is hours
startTime2 = model->GetStartTime(); // default to model start time
//err = -1; goto done;
}
else
{
DateTimeRec time;
char unitStr[24], junk[10];
timeUnits = new char[t_len+1];
//if (!fIsNavy)
//status = nc_get_att_text(ncid, recid, "units", timeUnits); // recid is the dimension id not the variable id
//else
status = nc_get_att_text(ncid, timeid, "units", timeUnits);
if (status != NC_NOERR) {err = -1; goto done;}
timeUnits[t_len] = '\0'; // moved this statement before StringSubstitute, JLM 5/2/10
StringSubstitute(timeUnits, ':', ' ');
StringSubstitute(timeUnits, '-', ' ');
StringSubstitute(timeUnits, 'T', ' ');
StringSubstitute(timeUnits, 'Z', ' ');
numScanned=sscanf(timeUnits, "%s %s %hd %hd %hd %hd %hd %hd",
unitStr, junk, &time.year, &time.month, &time.day,
&time.hour, &time.minute, &time.second) ;
if (numScanned==5)
{time.hour = 0; time.minute = 0; time.second = 0; }
else if (numScanned==7) // has two extra time entries ??
time.second = 0;
else if (numScanned<8)
//else if (numScanned!=8)
{
//timeUnits = 0; // files should always have this info
//timeConversion = 3600.; // default is hours
//startTime2 = model->GetStartTime(); // default to model start time
err = -1; TechError("NetCDFWindMoverCurv::TextRead()", "sscanf() == 8", 0); goto done;
}
else
{
// code goes here, trouble with the DAYS since 1900 format, since converts to seconds since 1904
if (time.year ==1900) {time.year += 40; time.day += 1; /*for the 1900 non-leap yr issue*/ yearShift = 40.;}
DateToSeconds (&time, &startTime2); // code goes here, which start Time to use ??
if (!strcmpnocase(unitStr,"HOURS") || !strcmpnocase(unitStr,"HOUR"))
timeConversion = 3600.;
else if (!strcmpnocase(unitStr,"MINUTES") || !strcmpnocase(unitStr,"MINUTE"))
timeConversion = 60.;
else if (!strcmpnocase(unitStr,"SECONDS") || !strcmpnocase(unitStr,"SECOND"))
timeConversion = 1.;
else if (!strcmpnocase(unitStr,"DAYS") || !strcmpnocase(unitStr,"DAY"))
timeConversion = 24.*3600.;
}
}
if (fIsNavy)
{
status = nc_inq_dimid(ncid, "gridy", &latIndexid); //Navy
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_dimlen(ncid, latIndexid, &latLength);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_dimid(ncid, "gridx", &lonIndexid); //Navy
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_dimlen(ncid, lonIndexid, &lonLength);
if (status != NC_NOERR) {err = -1; goto done;}
// option to use index values?
status = nc_inq_varid(ncid, "grid_lat", &latid);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_varid(ncid, "grid_lon", &lonid);
if (status != NC_NOERR) {err = -1; goto done;}
}
else
{
for (i=0;i<numdims;i++)
{
if (i == recid) continue;
status = nc_inq_dimname(ncid,i,dimname);
if (status != NC_NOERR) {err = -1; goto done;}
if (!strncmpnocase(dimname,"X",1) || !strncmpnocase(dimname,"LON",3) || !strncmpnocase(dimname,"nx",2))
{
lonIndexid = i;
}
if (!strncmpnocase(dimname,"Y",1) || !strncmpnocase(dimname,"LAT",3) || !strncmpnocase(dimname,"ny",2))
{
latIndexid = i;
}
}
status = nc_inq_dimlen(ncid, latIndexid, &latLength);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_dimlen(ncid, lonIndexid, &lonLength);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_inq_varid(ncid, "LATITUDE", &latid);
if (status != NC_NOERR)
{
status = nc_inq_varid(ncid, "lat", &latid);
if (status != NC_NOERR)
{
status = nc_inq_varid(ncid, "latitude", &latid);
if (status != NC_NOERR) {err = -1; goto done;}
}
}
status = nc_inq_varid(ncid, "LONGITUDE", &lonid);
if (status != NC_NOERR)
{
status = nc_inq_varid(ncid, "lon", &lonid);
if (status != NC_NOERR)
{
status = nc_inq_varid(ncid, "longitude", &lonid);
if (status != NC_NOERR) {err = -1; goto done;}
}
}
}
pt_count[0] = latLength;
pt_count[1] = lonLength;
vertexPtsH = (WorldPointF**)_NewHandleClear(latLength*lonLength*sizeof(WorldPointF));
if (!vertexPtsH) {err = memFullErr; goto done;}
lat_vals = new float[latLength*lonLength];
lon_vals = new float[latLength*lonLength];
if (!lat_vals || !lon_vals) {err = memFullErr; goto done;}
status = nc_get_vara_float(ncid, latid, ptIndex, pt_count, lat_vals);
if (status != NC_NOERR) {err = -1; goto done;}
status = nc_get_vara_float(ncid, lonid, ptIndex, pt_count, lon_vals);
if (status != NC_NOERR) {err = -1; goto done;}
for (i=0;i<latLength;i++)
{
for (j=0;j<lonLength;j++)
{
//if (lat_vals[(latLength-i-1)*lonLength+j]==fill_value) // this would be an error
//lat_vals[(latLength-i-1)*lonLength+j]=0.;
//if (lon_vals[(latLength-i-1)*lonLength+j]==fill_value)
//lon_vals[(latLength-i-1)*lonLength+j]=0.;
INDEXH(vertexPtsH,i*lonLength+j).pLat = lat_vals[(latLength-i-1)*lonLength+j];
INDEXH(vertexPtsH,i*lonLength+j).pLong = lon_vals[(latLength-i-1)*lonLength+j];
}
}
fVertexPtsH = vertexPtsH;
status = nc_inq_dim(ncid, recid, recname, &recs);
if (status != NC_NOERR) {err = -1; goto done;}
if (recs<=0) {strcpy(errmsg,"No times in file. Error opening NetCDF wind file"); err = -1; goto done;}
fTimeHdl = (Seconds**)_NewHandleClear(recs*sizeof(Seconds));
if (!fTimeHdl) {err = memFullErr; goto done;}
for (i=0;i<recs;i++)
{
Seconds newTime;
// possible units are, HOURS, MINUTES, SECONDS,...
timeIndex = i;
//if (!fIsNavy)
//status = nc_get_var1_float(ncid, recid, &timeIndex, &timeVal); // recid is the dimension id not the variable id
//else
status = nc_get_var1_float(ncid, timeid, &timeIndex, &timeVal);
if (status != NC_NOERR) {err = -1; goto done;}
newTime = RoundDateSeconds(round(startTime2+timeVal*timeConversion));
//INDEXH(fTimeHdl,i) = startTime2+(long)(timeVal*timeConversion -yearShift*3600.*24.*365.25); // which start time where?
//if (i==0) startTime = startTime2+(long)(timeVal*timeConversion -yearShift*3600.*24.*365.25);
INDEXH(fTimeHdl,i) = newTime-yearShift*3600.*24.*365.25; // which start time where?
if (i==0) startTime = newTime-yearShift*3600.*24.*365.25;
}
if (model->GetStartTime() != startTime || model->GetModelTime()!=model->GetStartTime())
{
if (true) // maybe use NOAA.ver here?
{
short buttonSelected;
//buttonSelected = MULTICHOICEALERT(1688,"Do you want to reset the model start time to the first time in the file?",FALSE);
if(!gCommandFileRun) // also may want to skip for location files...
buttonSelected = MULTICHOICEALERT(1688,"Do you want to reset the model start time to the first time in the file?",FALSE);
else buttonSelected = 1; // TAP user doesn't want to see any dialogs, always reset (or maybe never reset? or send message to errorlog?)
switch(buttonSelected){
case 1: // reset model start time
//bTopFile = true;
model->SetModelTime(startTime);
model->SetStartTime(startTime);
model->NewDirtNotification(DIRTY_RUNBAR); // must reset the runbar
break;
case 3: // don't reset model start time
//bTopFile = false;
break;
case 4: // cancel
err=-1;// user cancel
goto done;
}
}
//model->SetModelTime(startTime);
//model->SetStartTime(startTime);
//model->NewDirtNotification(DIRTY_RUNBAR); // must reset the runbar
}
fNumRows = latLength;
fNumCols = lonLength;
status = nc_close(ncid);
if (status != NC_NOERR) {err = -1; goto done;}
//err = this -> SetInterval(errmsg);
//if(err) goto done;
// look for topology in the file
// for now ask for an ascii file, output from Topology save option
// need dialog to ask for file
//if (fIsNavy) // for now don't allow for wind files
{if (topFilePath[0]) {err = ReadTopology(topFilePath,newMap); goto done;}}
if (!gCommandFileRun)
{
short buttonSelected;
buttonSelected = MULTICHOICEALERT(1688,"Do you have an extended topology file to load?",FALSE);
switch(buttonSelected){
case 1: // there is an extended top file
bTopFile = true;
break;
case 3: // no extended top file
bTopFile = false;
break;
case 4: // cancel
err=-1;// stay at this dialog
goto done;
}
}
if(bTopFile)
{
#if TARGET_API_MAC_CARBON
mysfpgetfile(&where, "", -1, typeList,
(MyDlgHookUPP)0, &reply, M38c, MakeModalFilterUPP(STDFilter));
if (!reply.good)/* return USERCANCEL;*/
{
/*if (recs>0)
err = this -> ReadTimeData(indexOfStart,&velocityH,errmsg);
else {strcpy(errmsg,"No times in file. Error opening NetCDF file"); err = -1;}
if(err) goto done;*/
err = dynamic_cast<NetCDFWindMoverCurv *>(this)->ReorderPoints(newMap,errmsg);
//err = ReorderPoints(fStartData.dataHdl,newMap,errmsg); // if u, v input separately only do this once?
goto done;
}
else
strcpy(topPath, reply.fullPath);
#else
where = CenteredDialogUpLeft(M38c);
sfpgetfile(&where, "",
(FileFilterUPP)0,
-1, typeList,
(DlgHookUPP)0,
&reply, M38c,
(ModalFilterUPP)MakeUPP((ProcPtr)STDFilter, uppModalFilterProcInfo));
if (!reply.good)
{
/*if (recs>0)
err = this -> ReadTimeData(indexOfStart,&velocityH,errmsg);
else {strcpy(errmsg,"No times in file. Error opening NetCDF file"); err = -1;}
if(err) goto done;*/
err = dynamic_cast<NetCDFWindMoverCurv *>(this)->ReorderPoints(newMap,errmsg);
//err = ReorderPoints(fStartData.dataHdl,newMap,errmsg);
/*if (err)*/ goto done;
}
my_p2cstr(reply.fName);
#ifdef MAC
GetFullPath(reply.vRefNum, 0, (char *)reply.fName, topPath);
#else
strcpy(topPath, reply.fName);
#endif
#endif
strcpy (s, topPath);
err = ReadTopology(topPath,newMap);
goto done;
}
/*if (recs>0)
err = this -> ReadTimeData(indexOfStart,&velocityH,errmsg);
else {strcpy(errmsg,"No times in file. Error opening NetCDF wind file"); err = -1;}
if(err) goto done;*/
err = dynamic_cast<NetCDFWindMoverCurv *>(this)->ReorderPoints(newMap,errmsg);
//err = ReorderPoints(fStartData.dataHdl,newMap,errmsg);
done:
if (err)
{
printNote("Error opening NetCDF wind file");
if(fGrid)
{
fGrid ->Dispose();
delete fGrid;
fGrid = 0;
}
if(vertexPtsH) {DisposeHandle((Handle)vertexPtsH); vertexPtsH = 0; fVertexPtsH = 0;}
}
if (timeUnits) delete [] timeUnits;
if (lat_vals) delete [] lat_vals;
if (lon_vals) delete [] lon_vals;
if (modelTypeStr) delete [] modelTypeStr;
//if (velocityH) {DisposeHandle((Handle)velocityH); velocityH = 0;}
return err;
}
