gridspec_api/mosaic/tst_two_tiles_print_polytopes.c

/*
 * Test the definition of a structured grid.
 *
 * "$Id: tst_two_tiles_print_polytopes.c 507 2011-02-14 21:27:52Z dkindig $"
 * */

#include <stdio.h>
#include <stdlib.h>

#include <libcf_src.h>
#include <netcdf.h>
#include <nccf_grid.h>
#include <nccf_coord.h>
#include <nccf_contacts.h>
#include <nccf_mosaic.h>
#include <nccf_varObj.h>
#include <nccf_handle_error.h>

int main(int argc,char **argv){

  /* Initialize the the cflists */
  CFLIST_STRUCTURED_GRID = NULL;
  CFLIST_CONTACTS = NULL;

  /* Variables for the grid */
  // Each Id goes with a ncid
  int contactId, mosaicid, i, j;
  int status, ncid;
  int iGrid; // Some counters
  int ndims = 2;         // Number of x, y dims
  int ngrids = 2;
  int nlon = 10;         // number of x faces, gives nlon + 1 edges
  int dims[] = {nlon, nlon};
  double period[] = {0., 0.};

  int save = 1;

  nc_type nc_mode = NC_CLOBBER;

  char *name = "tst_two_tile";
  char grid_filename[STRING_SIZE] = "tst_two_tiles_contacts.nc";
  char mosaic_filename[STRING_SIZE] = "tst_two_tiles_mosaic.nc";
  char gridname[STRING_SIZE];
  char filename[STRING_SIZE];
  const char *dimnames[] = {"ni", "nj"};

  double *lon1, *lat1, *lon2, *lat2, dlon, dlat;
  const double lonMin =   0.0;
  const double lonMax = 360.0;
  const double latMin = -90.0;
  const double latMax =  90.0;

  lon1 = ( double * )malloc( sizeof( double ) * nlon * nlon );
  lat1 = ( double * )malloc( sizeof( double ) * nlon * nlon );
  lon2 = ( double * )malloc( sizeof( double ) * nlon * nlon );
  lat2 = ( double * )malloc( sizeof( double ) * nlon * nlon );

  dlat = ( latMax - latMin )/( nlon - 1 );
  dlon = ( lonMax - lonMin )/( nlon - 1 )/2;

  for( j = 0; j < nlon; j++){
    for( i = 0; i < nlon; i++ ){
      lon1[j + nlon * i] = lonMin + j * dlon;
      lon2[j + nlon * i] = lonMin + j * dlon + 180;
      lat1[j + nlon * i] = latMin + i * dlat;
      lat2[j + nlon * i] = latMin + i * dlat;
    }
  }
  /* Ccoordinate ids - written over for each tile. */
//  coordIds = (int *) malloc( ndims * sizeof( int ));

  int coordIds[ndims];
//  gridid = ( int* )malloc( ngrids * sizeof( int ));
  int gridid[ndims];

/* Define a grid */
  ngrids = 2;
/* Calculate the coordinates for the current tile */

/* Define tile 0 */
  if(( status = nccf_def_lon_coord( ndims, dims, dimnames, lon1, save,
                                         &coordIds[0] ))) ERR;
  if(( status = nccf_def_lat_coord( ndims, dims, dimnames, lat1, save,
                                         &coordIds[1] ))) ERR;

  iGrid = 0;
  sprintf( gridname, "%s%d", name, iGrid + 1 );
  sprintf( filename, "%s.nc", gridname );
  if(( status = nccf_def_struct_grid( coordIds,
               gridname, &gridid[iGrid] ))) ERR;

/* Write tile 0 */
  nc_create( filename, nc_mode, &ncid );             
  nccf_put_struct_grid( ncid, filename, gridid[iGrid] );
  nc_close( ncid );

/* Define tile 1 */
  if(( status = nccf_def_lon_coord( ndims, dims, dimnames, lon2, save,
                                         &coordIds[0] ))) ERR;
  if(( status = nccf_def_lat_coord( ndims, dims, dimnames, lat2, save,
                                         &coordIds[1] ))) ERR;

  iGrid = 1;
  sprintf( gridname, "%s%d", name, iGrid + 1 );
  sprintf( filename, "%s.nc", gridname );
  if(( status = nccf_def_struct_grid(coordIds,
               gridname, &gridid[iGrid] ))) ERR;
  
/* Write tile 1 */
  nc_create( filename, nc_mode, &ncid );             
  nccf_put_struct_grid( ncid, filename, gridid[iGrid] );
  nc_close( ncid );

/* Define the contacts */
  if(( status = nccf_def_contacts( ngrids, gridid, gridname,
                                period, &contactId ))) ERR;

/* Write the file */
  if(( status = nc_create( grid_filename, nc_mode, &ncid ))) ERR;
  if(( status = nccf_put_contacts( ncid, filename, contactId ))) ERR;
  if(( status = nc_close( ncid ))) ERR;

/* Define the mosaic */  
  if(( status = nccf_def_mosaic( contactId, name, &mosaicid )));

/* Write the mosaic */
  if(( status = nc_create( mosaic_filename, nc_mode, &ncid ))) ERR;
  if(( status = nccf_put_mosaic( ncid, mosaicid ))) ERR;
  if(( status = nc_close( ncid ))) ERR;

  if(( status = nccf_print_contacts_as_polytopes(contactId,"polytopes" ))) ERR;

/* Free the mosaic */
  if(( status = nccf_free_mosaic( mosaicid )));

/* Free the contacts */
  if(( status = nccf_free_contacts( contactId ))) ERR;

/* Free the grid and the coordinates */
  for( iGrid = 0; iGrid < ngrids; iGrid++ ){
    if(( status = nccf_get_struct_grid_coordIds( gridid[iGrid],
        coordIds))) ERR;
    if(( status = nccf_free_struct_grid( gridid[iGrid] ))) ERR;
    for (i = 0; i < ndims; ++i) {
      if((status = nccf_free_coord(coordIds[i]))) ERR;
    }
  }


  free( lon1 );
  free( lon2 );
  free( lat1 );
  free( lat2 );
//  free(coordIds);
//  free(gridid);

/* Close the file*/
  return 0;
}