SURFEX v7.3: /home/dasprezs/EXPORT_v7_3/src/LIB/XRD38/FA/mt/faisc2

Go to the documentation of this file.
00001 C Jan-2011 P. Marguinaud Thread-safe FA
00002       SUBROUTINE FAISC2_MT (FA,  KREP, KRANGC )
00003       USE FA_MOD, ONLY : FA_COM
00004       USE PARKIND1, ONLY : JPRB
00005       USE YOMHOOK , ONLY : LHOOK, DR_HOOK
00006 C****
00007 C      Ce sous-programme initialise des tableaux "reference" de
00008 C      l'en-tete GRIB, section 2: les differents types de grille
00009 C      sont abordes (routine appelee une seule fois pour un cadre donne)
00010 C**
00011 C    Arguments : KREP   (Sortie) ==> Code-reponse du sous-programme;
00012 C                KRANGC (Entree) ==> Rang dans la table des cadres;
00013 C*
00014 C
00015 #include "precision.h"
00016 C
00017 C
00018       TYPE(FA_COM) :: FA
00019       INTEGER KREP, KRANGC
00020 C
00021       REAL ZPI, ZRAMDE, ZLATPRE, ZLATDER, ZLONPRE, ZLONDER
00022       REAL (KIND=JPDBLR), PARAMETER :: ONE = 1.0
00023 C
00024       INTEGER INLAT, INIVAU, INUMER, INIMES
00025 C
00026       LOGICAL LLMLAM
00027 C
00028 #include "facom2.h"
00029 #include "facom_mt.h"
00030 C**
00031 C     0.  -  INITIALISATIONS PREALABLES
00032 C-----------------------------------------------------------------------
00033 C
00034       REAL(KIND=JPRB) :: ZHOOK_HANDLE
00035       IF (LHOOK) CALL DR_HOOK('FAISC2_MT',0,ZHOOK_HANDLE)
00036       KREP=0
00037       IF (KRANGC.LE.0.OR.KRANGC.GT.FA%JPNXCA) THEN
00038         KREP=-66
00039         GOTO 1001
00040       ENDIF
00041       INLAT  = FA%NLATIT(KRANGC)
00042       INIVAU = FA%NNIVER(KRANGC)
00043       LLMLAM = FA%LIMLAM(KRANGC)
00044       ZPI = 2.*ASIN(1.)
00045 C Conversion des radians en 1/1000 de degre
00046       ZRAMDE = 180000./ZPI
00047 C
00048       IF (LLMLAM) GOTO 300
00049 C**
00050 C     1.  -  KSEC2 POUR LA REPRESENTATION SPECTRALE ARPEGE
00051 C-----------------------------------------------------------------------
00052 C
00053 C Type de representation de donnees
00054 C
00055 C FA%SSLAPO=sinus latitude du pole d'interet
00056 C (si=1, pole=pole N et pas de rotation)
00057 C FA%SCODIL=coeff de dilation (si =1, pas de dilatation)
00058       IF ((1.-FA%SSLAPO(KRANGC)).LE.1.E-10) THEN
00059         FA%NSEC2SP(1,KRANGC)=70
00060         IF (ABS(FA%SCODIL(KRANGC)-1.).LE.1.E-10) THEN
00061           FA%NSEC2SP(1,KRANGC)=50
00062         ENDIF
00063       ELSE
00064         FA%NSEC2SP(1,KRANGC)=80
00065         IF (ABS(FA%SCODIL(KRANGC)-1.).LE.1.E-10) THEN
00066           FA%NSEC2SP(1,KRANGC)=60
00067         ENDIF
00068       ENDIF
00069 C Troncature (3 fois la meme si triangulaire)
00070       FA%NSEC2SP(2,KRANGC) =FA%MTRONC(KRANGC)
00071       FA%NSEC2SP(3,KRANGC) =FA%MTRONC(KRANGC)
00072       FA%NSEC2SP(4,KRANGC) =FA%MTRONC(KRANGC)
00073 C Type de representation
00074       FA%NSEC2SP(5,KRANGC) =1
00075 C Mode de representation (2->complex packing)
00076       FA%NSEC2SP(6,KRANGC) =2
00077 C Reserves
00078       FA%NSEC2SP(7:11,KRANGC)=0
00079 C Nb de parametres pour la coord verticale
00080 C On prend ici le cas de la coordonnee hybride
00081 C mais le cas de la coord pression sera aisement
00082 C pris en compte + tard (KSEC2(12)=0).
00083       FA%NSEC2SP(12,KRANGC)=2*(INIVAU+1)
00084 C Latitude du pole sud de rotation
00085       FA%NSEC2SP(13,KRANGC)=0
00086 C Longitude du pole sud de rotation
00087       FA%NSEC2SP(14,KRANGC)=0
00088 C Lat et lon du pole d'etirement
00089       IF (FA%NTYPTR(KRANGC).GE.2) THEN
00090         FA%NSEC2SP(15,KRANGC)=ZRAMDE*ASIN(FA%SSLAPO(KRANGC))
00091         FA%NSEC2SP(16,KRANGC)=ZRAMDE*
00092      S    (SIGN(ONE,FA%SSLOPO(KRANGC))*ACOS(FA%SCLOPO(KRANGC)))
00093       ELSE
00094         FA%NSEC2SP(15,KRANGC)=0
00095         FA%NSEC2SP(16,KRANGC)=0
00096       ENDIF
00097 C Reserves
00098       FA%NSEC2SP(17:22,KRANGC)=0
00099 C**
00100 C     2.  -  KSEC2 POUR LA GRILLE DE GAUSS (ARPEGE)
00101 C-----------------------------------------------------------------------
00102 C
00103 C Type de representation de donnees
00104 C
00105 C FA%SSLAPO=sinus latitude du pole d'interet
00106 C (si=1, pole=pole N et pas de rotation)
00107 C FA%SCODIL=coeff de dilation (si =1, pas de dilatation)
00108       IF ((1.-FA%SSLAPO(KRANGC)).LE.1.E-10) THEN
00109         FA%NSEC2GG(1,KRANGC)=24
00110         IF (ABS(FA%SCODIL(KRANGC)-1.).LE.1.E-10) THEN
00111           FA%NSEC2GG(1,KRANGC)=4
00112         ENDIF
00113       ELSE
00114         FA%NSEC2GG(1,KRANGC)=34
00115         IF (ABS(FA%SCODIL(KRANGC)-1.).LE.1.E-10) THEN
00116           FA%NSEC2GG(1,KRANGC)=14
00117         ENDIF
00118       ENDIF
00119 C Nb de pts sur un parallele
00120       FA%NSEC2GG(2,KRANGC)=FA%NXLOPA(KRANGC)
00121 C Nb de pts sur une longitude
00122       FA%NSEC2GG(3,KRANGC)=INLAT
00123       ZLATPRE=ASIN(MAX(-ONE,MIN(ONE,FA%SINLAT(1,KRANGC))))
00124 C Latitude (1/1000 degre) du premier pt de grille
00125       FA%NSEC2GG(4,KRANGC)=ZRAMDE*ZLATPRE
00126 C Longitude (1/1000 degre) du premier pt de grille
00127       FA%NSEC2GG(5,KRANGC)=0
00128 C Flag pour la resolution (0->on ne donne pas l'increment)
00129       FA%NSEC2GG(6,KRANGC)=0
00130 C Latitude (1/1000 degre) du dernier pt de grille
00131       FA%NSEC2GG(7,KRANGC)=-FA%NSEC2GG(4,KRANGC)
00132 C Longitude (1/1000 degre) du dernier pt de grille.
00133 C (FA%NLOPAR(1,KRANGC)=nb de longitudes sur le 1er parallele)
00134       FA%NSEC2GG(8,KRANGC)=-360000/FA%NLOPAR(1,KRANGC)
00135 C Increment zonal (1/1000 degre)
00136 C Pas de sens ici.
00137       FA%NSEC2GG(9,KRANGC)=0
00138 C Nb de paralleles entre le pole et l'equateur
00139       FA%NSEC2GG(10,KRANGC)=(INLAT+1)/2
00140 C Flag pour le mode de balayage
00141       FA%NSEC2GG(11,KRANGC)=0
00142 C Nombre de parametres pour la coord. verticale
00143       FA%NSEC2GG(12,KRANGC)=0
00144 C Latitude du pole sud de rotation
00145       FA%NSEC2GG(13,KRANGC)=0
00146 C Longitude du pole sud de rotation
00147       FA%NSEC2GG(14,KRANGC)=0
00148 C Latitude et longitude du pole d'etirement
00149       IF (FA%NTYPTR(KRANGC).GE.2) THEN
00150         FA%NSEC2GG(15,KRANGC)=ZRAMDE*ASIN(FA%SSLAPO(KRANGC))
00151         FA%NSEC2GG(16,KRANGC)=ZRAMDE*
00152      S     (SIGN(ONE,FA%SSLOPO(KRANGC))*ACOS(FA%SCLOPO(KRANGC)))
00153       ELSE
00154         FA%NSEC2GG(15,KRANGC)=0
00155         FA%NSEC2GG(16,KRANGC)=0
00156       ENDIF
00157 C Flag:      0 -> grille reguliere,  1 -> grille reduite
00158       IF (FA%NLOPAR(1,KRANGC)==FA%NLOPAR((1+INLAT)/2,KRANGC)) THEN
00159         FA%NSEC2GG(17,KRANGC)=0
00160       ELSE
00161         FA%NSEC2GG(17,KRANGC)=1
00162       ENDIF
00163 C Flag:      0 -> Terre ronde     , 64 -> Terre ellipsoide
00164       FA%NSEC2GG(18,KRANGC)=0
00165 C Flag sur les composantes des vecteurs (0->geographique, 8->grille)
00166       FA%NSEC2GG(19,KRANGC)=0
00167 C Reserves
00168       FA%NSEC2GG(20:22,KRANGC)=0
00169 C Pour les grilles reduites, nb de points sur chaque parallele
00170       FA%NSEC2GG(23:22+(1+INLAT)/2,KRANGC)=
00171      S                 FA%NLOPAR(1:(1+INLAT)/2,KRANGC)
00172       FA%NSEC2GG(23-MOD(INLAT,2)+(1+INLAT)/2:22+INLAT,KRANGC)=
00173      S                               FA%NLOPAR((1+INLAT)/2:1:-1,KRANGC)
00174       GOTO 600
00175 C**
00176 C     3.  -  KSEC2 POUR LA GRILLE LAT-LON (CAS FULL-POS, ARPEGE OU ALADIN)
00177 C-------------------------------------------------------------------------
00178 C
00179 300   CONTINUE
00180 C TEST POUR NEW EGGX
00181       IF (FA%SINLAT(1,KRANGC) .GE. 0) THEN
00182 C OLD EGGX
00183       ZLATPRE=FA%SINLAT(7,KRANGC)
00184       ZLONPRE=FA%SINLAT(4,KRANGC)
00185       ZLATDER=FA%SINLAT(5,KRANGC)
00186       ZLONDER=FA%SINLAT(6,KRANGC)
00187 C Type de representation de donnees
00188       ELSE
00189 C NEW EGGX
00190       ZLATPRE=FA%SINLAT(16,KRANGC)
00191       ZLONPRE=FA%SINLAT(13,KRANGC)
00192       ZLATDER=FA%SINLAT(14,KRANGC)
00193       ZLONDER=FA%SINLAT(15,KRANGC)
00194       ENDIF
00195 C
00196 C Type de representation de donnees
00197       FA%NSEC2LL(1,KRANGC)=0
00198 C Nb de pts sur un parallele
00199       FA%NSEC2LL(2,KRANGC)=FA%NXLOPA(KRANGC)
00200 C Nb de pts sur une longitude
00201       FA%NSEC2LL(3,KRANGC)=INLAT
00202 C Latitude (1/1000 degre) du premier pt de grille
00203       FA%NSEC2LL(4,KRANGC)=NINT(ZRAMDE*ZLATPRE)
00204 C Longitude (1/1000 degre) du premier pt de grille
00205       FA%NSEC2LL(5,KRANGC)=NINT(ZRAMDE*ZLONPRE)
00206       IF (FA%NSEC2LL(5,KRANGC).GT.180000)
00207      S                      FA%NSEC2LL(5,KRANGC)=FA%NSEC2LL(5,KRANGC)-360000
00208 C Flag pour la resolution (128->on donne l'increment: grille reguliere)
00209       FA%NSEC2LL(6,KRANGC)=128
00210 C Latitude (1/1000 degre) du dernier pt de grille
00211       FA%NSEC2LL(7,KRANGC)=NINT(ZRAMDE*ZLATDER)
00212 C Longitude (1/1000 degre) du dernier pt de grille
00213       FA%NSEC2LL(8,KRANGC)=NINT(ZRAMDE*ZLONDER)
00214       IF (FA%NSEC2LL(8,KRANGC).GT.180000)
00215      S                      FA%NSEC2LL(8,KRANGC)=FA%NSEC2LL(8,KRANGC)-360000
00216 C Increment zonal (1/1000 degre)
00217       IF (ZLONPRE.GT.ZLONDER) THEN
00218         FA%NSEC2LL(9,KRANGC)=
00219      S    NINT((ZLONDER+2.*ZPI-ZLONPRE)*ZRAMDE/(FA%NXLOPA(KRANGC)-1))
00220       ELSE
00221         FA%NSEC2LL(9,KRANGC)=
00222      S        NINT((ZLONDER-ZLONPRE)*ZRAMDE/(FA%NXLOPA(KRANGC)-1))
00223       ENDIF
00224 C Increment meridien (1/1000 degre)
00225       FA%NSEC2LL(10,KRANGC)=
00226      S          NINT((ZLATPRE-ZLATDER)*ZRAMDE/(INLAT-1))
00227 C Flag pour le mode de balayage: W->E et S->N = 64; W->E et N->S = 0
00228 C Full-Pos produit des champs lat-lon ranges S->N pour ARP et ALD.
00229 C Or la BDAP attend un rangt N->S pour les grilles lat-lon.
00230 C FA renverse donc les champs issus de Full-Pos avant codage GRIBEX.
00231 C 
00232       FA%NSEC2LL(11,KRANGC)=0
00233 C Nombre de parametres pour la coord. verticale
00234       FA%NSEC2LL(12,KRANGC)=0
00235 C Latitude du pole sud de rotation
00236       FA%NSEC2LL(13,KRANGC)=0
00237 C Longitude du pole sud de rotation
00238       FA%NSEC2LL(14,KRANGC)=0
00239 C Latitude et longitude du pole d'etirement
00240       FA%NSEC2LL(15,KRANGC)=0
00241       FA%NSEC2LL(16,KRANGC)=0
00242 C Flag:      0 -> grille reguliere,  1 -> grille reduite
00243       FA%NSEC2LL(17,KRANGC)=0
00244 C Flag:      0 -> Terre ronde     , 64 -> Terre ellipsoide
00245       FA%NSEC2LL(18,KRANGC)=0
00246 C Flag sur les composantes des vecteurs (0->geographique, 8->grille)
00247       FA%NSEC2LL(19,KRANGC)=0
00248 C Reserves
00249       FA%NSEC2LL(20:22,KRANGC)=0
00250 C**
00251 C     4.  -  KSEC2 POUR LA GRILLE LAT-LON QUASI-REGULIERE ALADIN
00252 C            (en fait, tenue de camouflage pour les coeff spectraux
00253 C             que l'on va ranger en balayant le 1/4 de l'ellipse
00254 C             verticalement: axe X=axe n (nb d'onde meridien) et
00255 C             axe Y=axe m (nb d'onde zonal) afin de suivre le rangt
00256 C             dans le modele. Seuls les coeff spectraux qui seront
00257 C             compactes sont stockes sur la grille lat-lon, soit
00258 C             tous sauf ceux des axes et ceux inclus dans le carre
00259 C             delimite par la ss-tronc de non-compactage).
00260 C-----------------------------------------------------------------------
00261 C
00262 C Type de representation de donnees
00263       FA%NSEC2AL(1,KRANGC)=0
00264 C Nb de pts sur un parallele: valeur manquante
00265       FA%NSEC2AL(2,KRANGC)=2**16 -1
00266 C Nb de pts sur une longitude: nombre d'onde zonal max -1
00267 C associe au nombre d'onde meridien 1 (les CSP sur les axes sont
00268 C extraits des champs de CSP puisque non compactes)
00269       FA%NSEC2AL(3,KRANGC)=(FA%NOZPAR(6,KRANGC)-
00270      S                      FA%NOZPAR(5,KRANGC)+1)/4 -1
00271 C Latitude (1/1000 degre) du premier pt de grille: valeur bidon
00272       FA%NSEC2AL(4,KRANGC)=0
00273 C Longitude (1/1000 degre) du premier pt de grille: valeur bidon
00274       FA%NSEC2AL(5,KRANGC)=0
00275 C Flag pour la resolution (128->on donne l'increment: grille reguliere)
00276       FA%NSEC2AL(6,KRANGC)=0
00277 C Latitude (1/1000 degre) du dernier pt de grille: valeur bidon
00278       FA%NSEC2AL(7,KRANGC)=40000
00279 C Longitude (1/1000 degre) du dernier pt de grille: valeur bidon
00280       FA%NSEC2AL(8,KRANGC)=40000
00281 C Increment zonal (1/1000 degre)
00282       FA%NSEC2AL(9,KRANGC)=2**16 -1
00283 C Increment meridien (1/1000 degre): deduit des valeurs bidon
00284       FA%NSEC2AL(10,KRANGC)=(FA%NSEC2AL(7,KRANGC)-FA%NSEC2AL(4,KRANGC))/
00285      S                   (FA%NSEC2AL(3,KRANGC)-1)
00286 C Flag pour le mode de balayage
00287       FA%NSEC2AL(11,KRANGC)=0
00288 C Nombre de parametres pour la coord. verticale
00289       FA%NSEC2AL(12,KRANGC)=0
00290 C Latitude du pole sud de rotation
00291       FA%NSEC2AL(13,KRANGC)=0
00292 C Longitude du pole sud de rotation
00293       FA%NSEC2AL(14,KRANGC)=0
00294 C Latitude et longitude du pole d'etirement
00295       FA%NSEC2AL(15,KRANGC)=0
00296       FA%NSEC2AL(16,KRANGC)=0
00297 C Flag:      0 -> grille reguliere,  1 -> grille reduite
00298       FA%NSEC2AL(17,KRANGC)=1
00299 C Flag:      0 -> Terre ronde     , 64 -> Terre ellipsoide
00300       FA%NSEC2AL(18,KRANGC)=0
00301 C Flag sur les composantes des vecteurs (0->geographique, 8->grille)
00302       FA%NSEC2AL(19,KRANGC)=0
00303 C Reserves
00304       FA%NSEC2AL(20:22,KRANGC)=0
00305 C Les valeurs (22+1:22+FA%MTRONC(KRANGC)-1) representant les nb de pts
00306 C le long de chaque parallele (ici, le nb de coeff spectraux
00307 C pour un meme n (et -n), excepte le triangle et les axes non
00308 C compactes) dependent de la ss-troncature qui depend du fichier
00309 C et ne seront donc pas stockes dans le tableau FA%NSEC2AL qui
00310 C depend du cadre. Le tableau FA%NSC2ALF(FA%JPXTRO-1,FA%JPNXFA) les
00311 C contiendra.
00312 C
00313 C**
00314 C     5.  -  KSEC2 POUR LA GRILLE LAMBERT CONFORME (CAS GENERAL ALADIN)
00315 C-------------------------------------------------------------------------
00316 C
00317 C Type de representation de donnees
00318       FA%NSEC2LA(1,KRANGC)=3
00319 C Nb de pts sur un parallele
00320       FA%NSEC2LA(2,KRANGC)=FA%NXLOPA(KRANGC)
00321 C Nb de pts sur une longitude
00322       FA%NSEC2LA(3,KRANGC)=INLAT
00323 C
00324 C Les parametres communs sont regroupes
00325 C Flag pour la resolution (128->on donne l'increment: grille reguliere)
00326       FA%NSEC2LA(6,KRANGC)=128
00327 C Reserve
00328       FA%NSEC2LA(8,KRANGC)=0
00329 C Flag pour le mode de balayage: W->E et S->N = 64; W->E et N->S = 0
00330       FA%NSEC2LA(11,KRANGC)=64
00331 C Nombre de parametres pour la coord. verticale
00332       FA%NSEC2LA(12,KRANGC)=0
00333 C Latitude (1/1000 degre) du premier pt de grille
00334       FA%NSEC2LA(4,KRANGC)=NINT(ZRAMDE*ZLATPRE)
00335 C Longitude (1/1000 degre) du premier pt de grille
00336       FA%NSEC2LA(5,KRANGC)=NINT(ZRAMDE*ZLONPRE)
00337       IF (FA%NSEC2LA(5,KRANGC).GT.180000)
00338      S                      FA%NSEC2LA(5,KRANGC)=FA%NSEC2LA(5,KRANGC)-360000
00339 C
00340 C TEST POUR OLD/NEW EGGX
00341       IF (FA%SINLAT(1,KRANGC) .GE. 0) THEN
00342 C Old EGGX
00343 C Orientation de la grille
00344       FA%NSEC2LA(7,KRANGC)=NINT(ZRAMDE*FA%SINLAT(8,KRANGC))
00345       IF (FA%NSEC2LA(7,KRANGC).GT.180000)
00346      S                      FA%NSEC2LA(7,KRANGC)=FA%NSEC2LA(7,KRANGC)-360000
00347 C Dimension de la maille dans la direction X
00348       FA%NSEC2LA(9,KRANGC)=NINT(FA%SINLAT(15,KRANGC))
00349 C Dimension de la maille dans la direction Y
00350       FA%NSEC2LA(10,KRANGC)=NINT(FA%SINLAT(16,KRANGC))
00351 C Flag pour le centre de projection
00352 C (0: le pole Nord est sur le plan de projection
00353 C  et 1 seul centre de projection est utilise;
00354 C  128: idem sauf que c'est le pole Sud)
00355       IF (FA%SINLAT(9,KRANGC).GE.0) THEN
00356         FA%NSEC2LA(13,KRANGC)=0
00357       ELSE
00358         FA%NSEC2LA(13,KRANGC)=128
00359       ENDIF
00360 C Premiere latitude depuis le pole ou le cone coupe la sphere
00361       FA%NSEC2LA(14,KRANGC)=NINT(ZRAMDE*FA%SINLAT(9,KRANGC))
00362 C Deuxieme latitude depuis le pole ou le cone coupe la sphere
00363 C Dans Aladin, le plan de projection est rarement secant (cela
00364 C releve plus d'un domaine mal defini que d'un choix) et cette
00365 C possibilite va disparaitre bientot. Comme le calcul de cette
00366 C seconde latitude n'est pas aise (pb de convergence), on va
00367 C declarer la grille tangente! mais avec un WARNING...
00368       FA%NSEC2LA(15,KRANGC)=FA%NSEC2LA(14,KRANGC)
00369       IF (ABS(FA%SINLAT(10,KRANGC)-SIN(FA%SINLAT(9,KRANGC))).GT.1.E-10 
00370      S    .AND. FA%LFAMOP) THEN
00371         WRITE (UNIT=FA%NULOUT,FMT=*)
00372      S          ' FAISC2: WARNING !! La grille Lambert coupe en fait',
00373      S          ' la sphere, mais sera consideree comme tangente'
00374       ENDIF
00375       ELSE
00376 C NEW EGGX
00377 C Orientation de la grille
00378       FA%NSEC2LA(7,KRANGC)=NINT(ZRAMDE*FA%SINLAT(3,KRANGC))
00379       IF (FA%NSEC2LA(7,KRANGC).GT.180000)
00380      S                      FA%NSEC2LA(7,KRANGC)=FA%NSEC2LA(7,KRANGC)-360000
00381 C Dimension de la maille dans la direction X
00382       FA%NSEC2LA(9,KRANGC)=NINT(FA%SINLAT(7,KRANGC))
00383 C Dimension de la maille dans la direction Y
00384       FA%NSEC2LA(10,KRANGC)=NINT(FA%SINLAT(8,KRANGC))
00385 C Flag pour le centre de projection
00386 C (0: le pole Nord est sur le plan de projection
00387 C  et 1 seul centre de projection est utilise;
00388 C  128: idem sauf que c'est le pole Sud)
00389       IF (FA%SINLAT(4,KRANGC).GE.0) THEN
00390         FA%NSEC2LA(13,KRANGC)=0
00391       ELSE
00392         FA%NSEC2LA(13,KRANGC)=128
00393       ENDIF
00394 C Premiere latitude depuis le pole ou le cone coupe la sphere
00395       FA%NSEC2LA(14,KRANGC)=NINT(ZRAMDE*FA%SINLAT(4,KRANGC))
00396 C NEW EGGX toujours tangent
00397       FA%NSEC2LA(15,KRANGC)=FA%NSEC2LA(14,KRANGC)
00398       ENDIF
00399 
00400 C Reserve
00401       FA%NSEC2LA(16,KRANGC)=0
00402 C Flag:      0 -> grille reguliere
00403       FA%NSEC2LA(17,KRANGC)=0
00404 C Flag:      0 -> Terre ronde     , 64 -> Terre ellipsoide
00405       FA%NSEC2LA(18,KRANGC)=0
00406 C Flag sur les composantes des vecteurs (0->geographique, 8->grille)
00407       FA%NSEC2LA(19,KRANGC)=8
00408 C Latitude du pole sud
00409       FA%NSEC2LA(20,KRANGC)=0
00410 C Longitude du pole sud
00411       FA%NSEC2LA(21,KRANGC)=0
00412 C Reserve
00413       FA%NSEC2LA(22,KRANGC)=0
00414 C**
00415 C     6.  -  PARTIE REELLE DE LA SECTION 2 DE GRIBEX
00416 C-----------------------------------------------------------------------
00417 C
00418 600   CONTINUE
00419 C Angle de rotation
00420       FA%XSEC2(1,KRANGC)=0.
00421 C Coefficient d'etirement
00422       FA%XSEC2(2,KRANGC)=FA%SCODIL(KRANGC)
00423 C Reserve
00424       FA%XSEC2(3:10,KRANGC)=0.
00425 C Parametres pour la coordonnee verticale
00426       FA%XSEC2(11:11+INIVAU,KRANGC)=FA%SFOHYB(1,0:INIVAU,KRANGC)*
00427      S                           FA%SPREFE(KRANGC)
00428       FA%XSEC2(12+INIVAU:12+2*INIVAU,KRANGC)=
00429      S     FA%SFOHYB(2,0:INIVAU,KRANGC)
00430 C**
00431 C    10.  -  PHASE TERMINALE : MESSAGERIE, AVEC "ABORT" EVENTUEL,
00432 C            VIA LE SOUS-PROGRAMME "FAIPAR" .
00433 C-----------------------------------------------------------------------
00434 C
00435  1001 CONTINUE
00436 C
00437       IF (FA%LFAMOP.OR.KREP.NE.0) THEN
00438         INIMES=2
00439         CLNSPR='FAISC2'
00440         INUMER=FA%JPNIIL
00441 C
00442         WRITE (UNIT=CLMESS,FMT='(''KREP='',I4,'', KRANGC='',I4)')
00443      S     KREP, KRANGC
00444         CALL FAIPAR_MT (FA, INUMER,INIMES,KREP,.FALSE.,CLMESS,
00445      S                  CLNSPR,CLNSPR,.FALSE.)
00446       ENDIF
00447 C
00448       IF (LHOOK) CALL DR_HOOK('FAISC2_MT',1,ZHOOK_HANDLE)
00449       END
00450