ubuntu-fr

grdraster 1 -R-66/-60/30/35 -Gbermuda.nc -V

wget "http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO5/TOPO/ETOPO5/ETOPO5.DOS"

sudo mv ETOPO5.DOS $GMTHOME/share/dbase/etopo5.i2

wget "http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO2/ETOPO2v2-2006/ETOPO2v2g/raw_binary/ETOPO2v2g_f4_LSB.zip"
wget "http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/bedrock/grid_registered/binary_float/etopo1_bed_g.zip"
wget "http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/ice_surface/grid_registered/binary_float/etopo1_ice_g.zip"

unzip ETOPO2v2g_f4_LSB.zip
unzip etopo1_bed_g.zip
unzip etopo1_ice_g.zip

sudo ./gmtf4toi4 ETOPO2v2g_f4_LSB.flt $GMTHOME/share/dbase/etopo2.i4
sudo ./gmtf4toi4 etopo1_bed_g.flt $GMTHOME/share/dbase/etopo1bed.i4
sudo ./gmtf4toi4 etopo1_ice_g.flt $GMTHOME/share/dbase/etopo1ice.i4

cp $GMTHOME/share/dbase/grdraster.info $GMTHOME/share/dbase/grdraster.info.bak
sudo gedit $GMTHOME/share/dbase/grdraster.info

12 "ETOPO2 global topography"    "m"    -R-180/180/-90/90        -I2m        GG l 1        0    none    etopo2.i4
13 "ETOPO1 global topo Bedrock"    "m"    -R-180:0.5/180/-90:0.5/90    -I1m        GG l 1        0    none    etopo1bed.i4
14 "ETOPO1 global topo Ice"    "m"    -R-180:0.5/180/-90:0.5/90    -I1m        GG l 1        0    none    etopo1ice.i4

#include <fcntl.h>
#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <string.h>
#include <sys/stat.h>
#include <sys/mman.h>

#define BLOCK 20000

#define FORMAT long int


long tailleFichier(char *fic) {
    struct stat etat_fichier;
    if (stat(fic, & etat_fichier) !=0) {
        perror("stat");
        exit(EXIT_FAILURE);
        return 0;
    }
    return etat_fichier.st_size;
}


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

    int fichierHdl;    
    FORMAT *converti;
    FILE* sortie;
    float * projection;
    long taille_fichier;

    //***************************
    // Verification des arguments
    //***************************
    printf("Tool to convert data stored as float (4 bytes) to data stored as long signed integer (4 bytes)\n08/2008 - C. MEICHEL\n");

    if (argc<3) {
        printf("Syntaxe : %s fichier_a_lire fichier_de_sortie\n", argv[0]);
        exit(EXIT_FAILURE);
    }


    //*******************************************
    // Mise en memoire du fichier (lecture seule)
    //*******************************************

    if ((fichierHdl=open(argv[1], 0)) <0) {
        perror("open");
        exit(EXIT_FAILURE);
    }

    taille_fichier = tailleFichier(argv[1]);
    projection = (float *) mmap(NULL, taille_fichier, PROT_READ, MAP_SHARED, fichierHdl, 0);
    if (projection == (float *) MAP_FAILED) {
        perror("mmap");
        exit(EXIT_FAILURE);
    }

    close(fichierHdl);

    if ((sortie = fopen(argv[2], "w")) == NULL) {
        return 0;
    }
    converti = (FORMAT*) malloc(BLOCK*sizeof(FORMAT));


    //**********************
    // Traitement du fichier
    //**********************
    
    size_t nbblock = taille_fichier / (BLOCK*sizeof(float));
    size_t k;
    size_t curseur = 0;
    size_t n;
    printf("%d blocks (%d values) to write\n", nbblock, BLOCK);
    for(k=0; k<nbblock; k++) {
        for(n=0; n<BLOCK; n++) {
            converti[n] = (FORMAT) projection[curseur++];
        }
        printf("%d out of 100 converted\r", k*100/(nbblock));
        fwrite(converti, (BLOCK)*sizeof(FORMAT), 1, sortie);
    }
    for(n=0; n<(taille_fichier % (BLOCK*sizeof(float)))/sizeof(FORMAT); n++) 
        converti[n] = (FORMAT) projection[curseur++];
    fwrite(converti, taille_fichier % (BLOCK*sizeof(FORMAT)), 1, sortie);


    //******************************************
    // Fin du code et restauration de la memoire
    //******************************************
    fclose(sortie);
    free(converti);
    munmap((void *) projection, taille_fichier);
    return 0;
}

gcc -o gmtf4toi4 gmtf4toi4.c