CALL SEQUENCE...
       calpres   GPST GPSLNSP PHIS DATAIN LNPM LNPT LNPHT \ 
                 INPUT OUTPUT [-a VAL1 -b VAL2]
 
 
      DESCRIPTION...
       calpres - CALCUL DU LOG DE LA PRESSION HYDROSTATIQUE SUR
                 LES SURFACES DES NIVEAUX GAL-CHEN
 
 										
      AUTHORS  - G. Bergeron et C. Thurre, JUILLET 01/1993
      LAST REVISION: $Header: First-release @ 2015-01-07 16:36:15 -0500 (Bernard Dugas) $
 
 
      BUT     -  CE PROGRAMME EVALUE LA PRESSION HYDROSTATIQUE (PI)
                 CORRESPONDANTE POUR LES POINTS DE GRILLE SITUES SUR
                 TOUS LES NIVEAUX GAL-CHEN.
 										
                 NOUS GENERONS DEUX ENSEMBLES DE RESULTATS:
                 LE PREMIER EST ASSOCIE AUX NIVEAUX MOMENTUM (LNPM).
                 LE SECOND EST ASSOCIE AUX NIVEAUX THERMODYNAMIQUE (LNPT).
                 NOUS CALCULONS EGALEMENT LE LOG DE LA PRESSION AU NIVEAU
                 DU TOIT (LNPHT).
 
                 POUR CHAQUE POINT DE GRILLE, NOUS INTEGRONS L'EQUATION
                 HYDROSTATIQUE DLN(PI)=-(G/RT)DZ POUR OBTENIR LA VARIABLE
                 LN(PI). LA PRESSSION HYDROSTATIQUE EST REPRESENTEE PAR PI
                 ET DZ REPRESENTE LA HAUTEUR PHYSIQUE ENTRE DEUX NIVEAUX
                 GAL-CHEN DE MEME TYPE. UNE FOIS QUE NOUS CONNAISSONS LN(PI)
                 POUR LES NIVEAUX DE TYPE MOMENTUM, NOUS EFFECTUONS UNE
                 MOYENNE VERTICALE POUR OBTENIR LES LN(PI) ASSOCIES A
                 L'ENSEMBLE DES POINTS DE GRILLE SITUEES SUR LES NIVEAUX
                 DE TYPE THERMODYNAMIQUE.
 
 
      FICHIERS D'ENTRES...
       GPST    = TEMPERATURE SUR DES NIVEAUX GAL-CHEN (NOM: 'TEMP' OU 'TT')
       GPSLNSP = LOG DE LA PRESSION DE SURFACE        (NOM: 'LNSP' OU 'LP')
       PHIS    = GEOPOTENTIEL AU SOL                  (NOM: 'PHIS' OU 'MT')
       DATAIN  = LISTE DES PARAMETRES PROPRES A L'INTEGRATION (optionel)
                                                      (NOM: 'DATA' OU 'XX')
 
      FICHIERS DE SORTIES...
       LNPM    = LN(PI) CORRESPONDANT AUX POINTS DE GRILLE SITUES SUR DES
                 NIVEAUX DE TYPE MOMENTUM             (NOM: 'LNPM' OU 'PM').
       LNPT    = LN(PI) CORRESPONDANT AUX POINTS DE GRILLE SITUES SUR DES
                 NIVEAUX DE TYPE THERMODYNAMIQUE      (NOM: 'LNPT' OU 'PT')
       LNPHT   = LN(PI) CORRESPONDANT AUX POINTS DE GRILLE SITUES AU TOIT
                 DU MODELE                            (NOM: 'LNPH' OU 'PH')
 
      NOTONS QUE SI GPST CONTIENT DES 'TEMP', GPSLNSP, PHIS ET DATAIN
      DEVRONT CONTENIR DES 'LNSP', 'PHIS' ET 'DATA', RESPECTIVEMENT. 
      ET LES FICHIERS DE SORTIES RESPECTENT LA MEME CONVENTION. LES
      NOMS A DEUX CARACTERES SONT REQUIS POUR LES FICHIERS DE TYPE
      CMC/RPN 1989. ILS NE LE SONT PLUS DEPUIS LA VERSION 2000.
 
 
      PARAMETRES D'ENTREE (FACULTATIFS, SI LE FICHIER DATAIN EXISTE)...
 
        HTOIT  (a) = HAUTEUR DU TOIT DU MODELE (OU BIEN DATAIN(8))
        TMOYEN (b) = TEMPERATURE MOYENNE UTILISEE POUR LE CALCUL
                     DES HAUTEURS PSEUDO-SIGMA (OU BIEN DATAIN(11))
 
        READ(5,5011,END=901) HTOIT,TMOYEN
 5011   FORMAT(10X,2F10.1)
 
 
      EXEMPLE DE CARTE D'ENTRE...
 
 CALPRES      38770.      250.
 
         0123456789012345678901234567890123456789012345678901234567890
         1         2         3         4         5         6         7
 
 
      EXIT CONDITIONS...
        0 NO PROBLEMS
        1 PROBLEM READING THE DIRECTIVE CARD
        2 TOO MANY OR NOT ENOUGH VARTICAL LEVELS IN GPST
        3 UNABLE TO READ PHIS
        4 PHIS AND GPST NOT THE SAME SIZE
        5 UNABLE TO READ FIRST SET IN GPST
        6 FILE GPST CONTAINS ANOTHER VARIABLE THAN 'TEMP'
        7 VARIABLE NUMBER OF VERTICAL LEVELS IN GPST
        8 MORE THAN ONE SIZE IN GPST
        9 PROBLEMS WHILE READING GPSLNSP
       10 GPSLNSP AND GSPT NOT THE SAME SIZE
       11 UNABLE TO DETERMINE SIZE AND VERTICAL LEVELS OF GSPT
       12 GPST DOES NOT CONTAIN GRIDDED DATA
       13 GPST DOES NOT CONTAIN GAL-CHEN LEVELS