Skip to content
Snippets Groups Projects
caresolv.f90 173 KiB
Newer Older
!     ######spl
      MODULE MODI_CARESOLV
!     ####################
!
INTERFACE
!
SUBROUTINE CARESOLV(HCARIN)
CHARACTER(LEN=*) :: HCARIN
END SUBROUTINE CARESOLV
!
END INTERFACE
!
END MODULE MODI_CARESOLV
!     ######spl
      SUBROUTINE CARESOLV(HCARIN)
!     ###########################
!
!!****  *CARESOLV* - 
!!
!!    PURPOSE
!!    -------
!      
!
!!**  METHOD
!!    ------
!!     
!!     N.A.
!!
!!    EXTERNAL
!!    --------
!!      None
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      Module MODD_CST  : declares model physical constants
!!
!!      Module MODN_NCAR : defines NAM_DIRTRA_POS namelist 
!!                         (former NCAR common)
!!
!!       NIOFFD     : Label normalisation (=0 none, =/=0 active)
!!       NULBLL     : Nb of contours between 2 labelled contours
!!       NIOFFM     : =0    --> message at picture bottom
!!                    =/= 0 --> no message
!!       NIOFFP     : Special point value detection
!!                    (=0 none, =/=0 active)
!!       NHI        : Extrema detection
!!                    (=0 --> H+L, <0 nothing)
!!       NINITA     : For streamlimes
!!       NINITB     : Not yet implemented
!!       NIGRNC     : Not yet implemented
!!       NDOT       : Line style
!!                    (=0|1|1023|65535 --> solid lines;
!!                    <0 --> solid lines for positive values and
!!                    dotted lines(ABS(NDOT))for negative values;
!!                    >0 --> dotted lines(ABS(NDOT)) )
!!       NIFDC      : Coastline data style (0 none, 1 NCAR, 2 IGN)
!!       NLPCAR     : Number of land-mark points to be plotted
!!       NIMNMX     : Contour selection option
!!                    (=-1 Min, max and inc. automatically set;
!!                    =0 Min, max automatically set; inc. given;
!!                    >0 Min, max, inc. given by user)
!!       NISKIP     : Rate for drawing velocity vectors
!!       CTYPHOR    : Horizontal cross-section type
!!                    (='K' --> model level section;
!!                     ='Z' --> constant-altitude section;
!!                     ='P' --> isobar section (planned)
!!                     ='T' --> isentrope section (planned)
!!       XSPVAL     : Special value
!!       XSIZEL     : Label size
!!       XLATCAR, XLONCAR :  Lat. and Long. of land-mark points
!!       LXY        : If =.TRUE., plots  a grid-mesh stencil background
!!       LXZ        : If =.TRUE., plots  a model-level stencil background 
!!
!!      Module MODN_PARA  : Defines NAM_DOMAIN_POS namelist 
!!                          (former PARA common)
!!
!!       XIDEBCOU, XJDEBCOU : Origin of a vertical cross-section
!!                            in cartesian (or conformal) real values
!!       XHMIN      : Altitude of the vert. cross-section
!!                    bottom (in meters above sea-level)
!!       XHMAX      : Altitude of the vert. cross-section
!!                    top (in meters above sea-level)
!!
!!      Module MODD_ALLVAR 
!!
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      J. Duron    * Laboratoire d'Aerologie *
!!
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       06/06/94
!!      Updated   PM   02/12/94
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
#ifdef NAGf95
USE F90_UNIX  ! for FLUSH and GETENV
#endif
USE MODD_RESOLVCAR
USE MODE_GRIDPROJ
USE MODD_SUPER
USE MODD_EXPR 
USE MODD_RADAR 
USE MODD_MASK3D
USE MODN_NCAR
USE MODN_PARA
USE MODI_RESOLVT
USE MODI_RESOLVN
USE MODI_RESOLVON
USE MODI_RESOLVP
USE MODI_RESOLVK
USE MODI_RESOLVZ
USE MODI_RESOLVX
USE MODI_CAREAL
USE MODI_CARINT
USE MODI_RESOLVI
USE MODI_RESOLVIARRAY
USE MODI_RESOLVL
USE MODI_CARMEMORY
USE MODI_RESOLVXISOLEV
USE MODD_FILES_DIACHRO
USE MODD_DIM1
USE MODD_GRID1
USE MODD_GRID, ONLY: XLONORI,XLATORI
USE MODD_HACH
USE MODD_PVT
USE MODD_COORD
USE MODD_DEFCV
USE MODD_TIT
USE MODD_RSISOCOL
USE MODD_MEMCV
USE MODD_CTL_AXES_AND_STYL
USE MODD_PARAMETERS
USE MODD_ALLOC_FORDIACHRO
USE MODI_LOADMNMX_FT_PVKT
USE MODI_LOADMNMXINT_ISO
USE MODI_LOADXISOLEVP
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
USE MODD_TRAJ3D
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
USE MODI_WRITEDIR

IMPLICIT NONE
!
!*       0.1   Dummy arguments
!              ---------------

CHARACTER(LEN=*) :: HCARIN
!CHARACTER(LEN=32) :: HGRP
!
!*       0.1   Local variables
!              ---------------
!
CHARACTER(LEN=80)  :: YCAR80
!CHARACTER(LEN=LEN(HCARIN)+240) :: YCARIN, YCAR
!CHARACTER(LEN=800) :: YCAROUT
CHARACTER(LEN=LEN(HCARIN)) :: YCARIN, YCAROUT, YCAR
!
INTEGER   ::   INDPARTIEL
INTEGER   ::   ILENC, ILENGRP, ILENC2, IETOILE, ILOG, ILEN, INDEXPR
INTEGER   ::   INDP, INDT, INDK, INDZ, INDCV, INDPV, INDPVT, INDPH, INDON, &
               INDTOT, INDMIN, INDFI, INDN, INDID, INDPVKT, INDPXT, INDPYT
INTEGER   ::   INDNDOMAINL, INDNDOMAINR, INDNDOMAINB, INDNDOMAINT
INTEGER   ::   INDNSZLBX, INDNSZLBY, INDTMP
INTEGER   ::   INDPVKT1, INDZT, INDXT, INDYT, INDXY, INDZTPVKT1, INDXYZ
INTEGER   ::   INDFT, INDFT1, INDMASK, INDMASKCUM, INDMASKSUM
INTEGER   ::   INDTK, INDPR, INDRS, INDRS1, INDEV, INDUMVMPV
INTEGER   ::   INDIINF, INDJINF, INDISUP, INDJSUP, INDIM, INDNZSTR, INDNARSTR
INTEGER   ::   INDIDEBCOU, INDJDEBCOU, INDXIDEBCOU, INDXJDEBCOU, INDNLANGLE, &
               INDNLMAX, INDXHMIN, INDXHMAX, INDXISOMIN, INDXISOMAX
INTEGER   ::   INDXPMIN, INDXPMAX, INDXPINT, INDXSSP, INDXLWSTR, INDXARLSTR
INTEGER   ::   INDXISOMIN_, INDXISOMAX_, INDXDIAINT_, INDXANGULVT
INTEGER   ::   INDXDIAINT, INDLXY, INDLXZ, INDLISO, INDLMINMAX, INDATFILE,  &
               INDLCOLAREA, INDLCOLINE, INDLISOWHI, INDLCOLBR, INDLCOLAREASEL, &
               INDLSPVALT, INDLSEGM, INDLPRESY, INDLSPSECT, INDLEGVECT , &
               INDLINTERPTOP, INDLCOLISONE, INDLCOLRSONE, INDLCOLRS1ONE,   &
               INDLCOLINESEL, INDLTABCOLDEF, INDVISU, INDNOVISU, INDXSIZEL, &
               INDLMNMXUSER, INDLCOLUSER, INDLVECTMNMX, INDLANIMK, INDLANIMT, &
               INDLMNMXLOC, INDLULMVTMOLD, INDLTITFTUSER, INDLFMTAXEX,  &
               INDLTIMEUSER, INDXTIMEMIN, INDXTIMEMAX, INDLSTREAM, &
               INDLNOUVRS, INDLMYHEURX, INDNHEURXLBL, INDNHEURXGRAD
INTEGER   ::   INDLINZEROPV, INDLBLFT1SUP, INDL24H
INTEGER   ::   INDNVERBIA, INDLISOWHI2, INDLISOWHI3, INDLFMTAXEY
INTEGER   ::   INDLINVWB , INDLGEOG, INDLMASK3D, INDMSKTOP, INDSV3
INTEGER   ::   INDLINVPTIR, INDLDOMAIN, INDLNOLABELX, INDLNOLABELY
INTEGER   ::   INDLNOLBLBAR
INTEGER   ::   INDLMASK3D_XY,INDLMASK3D_XZ,INDLMASK3D_YZ,INDLXYZ00
INTEGER   ::   INDLINDSP, INDLOGNEP, INDLTABCOLDEF2, INDLCONT, INDLRELIEF, &
               INDLCONV2XY, INDLINDAX, INDLCHREEL,INDLCOLUSERUV,INDL2CONT, &
               INDLCONVG2MASS
INTEGER   ::   INDLSPLO, INDSPO, INDOSPLO, INDPHALO, INDPHAO, INDLFTBAUTO, &
               INDLFT1BAUTO
INTEGER   ::   INDLPRINT, INDLPOINTG, INDL2DBX, INDL2DBY, INDLXYO, INDLPRINTXY
INTEGER   ::   INDLPRDAT, INDLINTERPOLSTR, INDL3D
INTEGER   ::   INDLVPTUSER, INDLVPTVUSER, INDLVPTPVUSER, INDLXABSC, INDLXMINTOP
INTEGER   ::   INDLVPTXYUSER, INDLFACTIMP, INDLFACTAXEX, INDLFACTAXEY
INTEGER   ::   INDLAXEXUSER, INDLAXEYUSER
INTEGER   ::   INDLHACH1, INDLHACH2, INDLHACH3, INDLHACH4, INDLHACHSEL
INTEGER   ::   INDLGREY
INTEGER   ::   INDLHEURX
INTEGER   ::   INDLABEL1, INDLBLUSER1, INDLBLUSER2, INDLBLUSER3, INDLBLUSER4
INTEGER   ::   INDXFACTAXEX, INDXFACTAXEY, INDXAXEXUSERD, INDXAXEYUSERD
INTEGER   ::   INDXAXEXUSERF, INDXAXEYUSERF
INTEGER   ::   INDXVPTL, INDXVPTR, INDXVPTB, INDXVPTT
INTEGER   ::   INDXVPTVL, INDXVPTVR, INDXVPTVB, INDXVPTVT
INTEGER   ::   INDXVPTPVL, INDXVPTPVR, INDXVPTPVB, INDXVPTPVT
INTEGER   ::   INDXVPTXYL, INDXVPTXYR, INDXVPTXYB, INDXVPTXYT
INTEGER   ::   INDXISOLEV, INDXFTMIN, INDXFTMAX, INDXPVKTMIN, INDXPVKTMAX
INTEGER   ::   INDXFT_ADTIM1, INDXFT_ADTIM2, INDXFT_ADTIM3, INDXFT_ADTIM4
INTEGER   ::   INDXFT_ADTIM5, INDXFT_ADTIM6, INDXFT_ADTIM7, INDXFT_ADTIM8
INTEGER   ::   INDXFT1_ADTIM1, INDXFT1_ADTIM2, INDXFT1_ADTIM3, INDXFT1_ADTIM4
INTEGER   ::   INDXFT1_ADTIM5, INDXFT1_ADTIM6, INDXFT1_ADTIM7, INDXFT1_ADTIM8
INTEGER   ::   INDXISOLEV_ , INDXPARCOLUV
INTEGER   ::   INDXFT1MIN, INDXFT1MAX, INDXFT1MIN_, INDXFT1MAX_
INTEGER   ::   INDXVARMIN, INDXVARMAX, INDXZTMIN, INDXZTMAX
INTEGER   ::   INDXLATCAR, INDXLONCAR, INDXXL, INDXXH, INDXYL, INDXYH
INTEGER   ::   INDXICAR, INDXJCAR
INTEGER   ::   INDCNOMCAR, INDCSYMCAR, INDXPOSNOM, INDXSZNOM, INDXSZSYM
INTEGER   ::   INDXPOSXVARNPV1TOP,INDXPOSYVARNPV1TOP
INTEGER   ::   INDXPOSXVARNPV5BOT,INDXPOSYVARNPV5BOT
INTEGER   ::   INDXSZVARNPVTOP,INDXSZVARNPVBOT
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! 19/12/2008 : modification pour controler la taille et le format des labels !!
!! pour les retrotrajectoires                                                 !!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
INTEGER   ::   INDNSZRTRAJ,INDLFMTRTRAJ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
INTEGER   ::   INDXLWPH1, INDXLWPH2, INDXLWPH3, INDXLWPH4
INTEGER   ::   INDXLWPH5, INDXLWPH6, INDXLWPH7, INDXLWPH8
INTEGER   ::   INDXSZTITYT, INDXSZTITYM, INDXSZTITYB
INTEGER   ::   INDXSZTITT1, INDXSZTITT2, INDXSZTITT3
INTEGER   ::   INDXSZTITB1, INDXSZTITB2, INDXSZTITB3
INTEGER   ::   INDXSZTITVAR1, INDXSZTITVAR2, INDXSZTITVAR3, INDXSZTITVAR4
INTEGER   ::   INDXSZTITVAR5, INDXSZTITVAR6, INDXSZTITVAR7, INDXSZTITVAR8
INTEGER   ::   INDXPOSTITYT, INDXPOSTITYM, INDXPOSTITYB
INTEGER   ::   INDXPOSTITT1, INDXPOSTITT2, INDXPOSTITT3
INTEGER   ::   INDXPOSTITB1, INDXPOSTITB2, INDXPOSTITB3
INTEGER   ::   INDXPOSTITVAR1, INDXPOSTITVAR2, INDXPOSTITVAR3, INDXPOSTITVAR4
INTEGER   ::   INDXPOSTITVAR5, INDXPOSTITVAR6, INDXPOSTITVAR7, INDXPOSTITVAR8
INTEGER   ::   INDXYPOSTITYT, INDXYPOSTITYM, INDXYPOSTITYB
INTEGER   ::   INDXYPOSTITT1, INDXYPOSTITT2, INDXYPOSTITT3
INTEGER   ::   INDXYPOSTITB1, INDXYPOSTITB2, INDXYPOSTITB3
INTEGER   ::   INDXYPOSTITVAR1, INDXYPOSTITVAR2, INDXYPOSTITVAR3,INDXYPOSTITVAR4
INTEGER   ::   INDXYPOSTITVAR5, INDXYPOSTITVAR6, INDXYPOSTITVAR7,INDXYPOSTITVAR8
INTEGER   ::   INDICOLNOM, INDICOLSYM
INTEGER   ::   INDXZL, INDXZH
INTEGER   ::   INDNIFDC, INDNIGRNC, INDNDOT, INDNISKIP, INDNISKIPVX, INDNISKIPVY
INTEGER   ::   INDNIOFFD, INDNULBLL, INDNHI, INDNIMNMX, INDNPROFILE, INDNSD
INTEGER   ::   INDMINUS, INDPLUS, INDNFT1ITVXMJ, INDNFT1ITVXMN
INTEGER   ::   INDNFT1ITVYMJ, INDNFT1ITVYMN, INDNINDCOLUV
INTEGER   ::   INDNFTITVXMJ, INDNFTITVXMN, INDNFTITVYMJ, INDNFTITVYMN
INTEGER   ::   INDNCHITVXMJ, INDNCHITVXMN, INDNCHITVYMJ, INDNCHITVYMN
INTEGER   ::   INDNCHPCITVXMJ, INDNCHPCITVXMN, INDNCHPCITVYMJ, INDNCHPCITVYMN
INTEGER   ::   INDNCVITVXMJ, INDNCVITVXMN, INDNCVITVYMJ, INDNCVITVYMN
INTEGER   ::   INDNPVITVXMJ, INDNPVITVXMN, INDNPVITVYMJ, INDNPVITVYMN
INTEGER   ::   INDNMASKITVXMJ, INDNMASKITVXMN, INDNMASKITVYMJ, INDNMASKITVYMN
INTEGER   ::   INDNXYITVXMJ, INDNXYITVXMN, INDNXYITVYMJ, INDNXYITVYMN
INTEGER   ::   INDXPVMIN, INDXPVMAX, INDXPVMINT, INDXPVMAXT
INTEGER   ::   INDXLWV, INDXLWDEF, INDXLWVDEF, INDXLWCONT, INDLMARKER
INTEGER   ::   INDXLWFTALL, INDXLWSEGM
INTEGER   ::   INDXLW, INDXLW1, INDXLW2, INDXLW3, INDXLW4, INDXLWDOMAIN
INTEGER   ::   INDXLWPV1, INDXLWPV2, INDXLWPV3, INDXLWPV4, INDXLWPV5
INTEGER   ::   INDXLWPV6, INDXLWPV7, INDXLWPV8, INDXLWTRACECV
INTEGER   ::   INDXLWPV9, INDXLWPV10
INTEGER   ::   INDXLWPV11, INDXLWPV12, INDXLWPV13, INDXLWPV14, INDXLWPV15
INTEGER   ::   INDXSTYLPV1, INDXSTYLPV2, INDXSTYLPV3, INDXSTYLPV4, INDXSTYLPV5
INTEGER   ::   INDXSTYLPV6, INDXSTYLPV7, INDXSTYLPV8
INTEGER   ::   INDXSTYLPV9, INDXSTYLPV10
INTEGER   ::   INDXSTYLPV11, INDXSTYLPV12, INDXSTYLPV13, INDXSTYLPV14, INDXSTYLPV15
INTEGER   ::   INDXAMX, INDXVHC, INDXVRL, INDXVLC, INDXVRLPH, INDXVHCPH
INTEGER   ::   INDNIRS, INDNJRS, INDXIRS, INDXJRS, INDXSPVAL, INDXSPVALT
INTEGER   ::   INDLCOLZERO, INDNCOLZERO, INDLFT1STYLUSER, INDLFTSTYLUSER
INTEGER   ::   INDNCOLSEGM, INDLFT3C, INDLFT4C, INDLFTCLIP
INTEGER   ::   INDNCOLUV1, INDNCOLUV2, INDNCOLUV3, INDNCOLUV4, INDNCOLUV5
INTEGER   ::   INDNCOLISONE1, INDNCOLISONE2, INDNCOLISONE3,INDNCOLISONE4
INTEGER   ::   INDNCOLISONE5, INDNCOLRSONE, INDNCOLRS1ONE1,INDNCOLRS1ONE2
INTEGER   ::   INDNCOLRS1ONE3, INDNCOLRS1ONE4,INDNCOLRS1ONE5
INTEGER   ::   INDXSZTITXL, INDXSZTITXM, INDXSZTITXR
INTEGER   ::   INDLDEFCV2, INDLDEFCV2LL, INDLDEFCV2IND
INTEGER   ::   INDXIDEBCV, INDXIFINCV, INDXJDEBCV, INDXJFINCV
INTEGER   ::   INDXIDEBCVLL, INDXIFINCVLL, INDXJDEBCVLL, INDXJFINCVLL
INTEGER   ::   INDNIDEBCV, INDNIFINCV, INDNJDEBCV, INDNJFINCV
INTEGER   ::   INDLSYMB, INDLTEXTG, INDLTEXTIT, INDLSYMBTEXTG, INDLSTI
INTEGER   ::   INDLTRACECV, INDLM5S3, INDLCVZOOM, INDLDILW, INDLVST
INTEGER   ::   INDLVSUPSCA, INDLXYWINCUR, INDLXYNVARTOP, INDLXYSTYLTOP
INTEGER   ::   INDLRADAR, INDXLATRAD1,INDXLATRAD2,INDXLATRAD3,INDXLATRAD4
INTEGER   ::   INDLRADIST, INDLRADRAY
INTEGER   ::   INDXLONRAD1,INDXLONRAD2,INDXLONRAD3,INDXLONRAD4
INTEGER   ::   INDXPORTRAD1,INDXPORTRAD2,INDXPORTRAD3,INDXPORTRAD4
INTEGER   ::   INDXLWRAD1,INDXLWRAD2,INDXLWRAD3,INDXLWRAD4
INTEGER   ::   INDCSYMRAD1,INDCSYMRAD2,INDCSYMRAD3,INDCSYMRAD4
INTEGER   ::   INDXISOREF,INDXISOREF_,INDLSPOT
INTEGER   ::   INDLFT1LUSER,INDNFT1STY1,INDNFT1STY2,INDNFT1STY3,INDNFT1STY4
INTEGER   ::   INDNFT1STY5,INDNFT1STY6,INDNFT1STY7,INDNFT1STY8,INDNFT1STY9
INTEGER   ::   INDNFT1STY10,INDNFT1STY11,INDNFT1STY12,INDNFT1STY13,INDNFT1STY14
INTEGER   ::   INDNFT1STY15
INTEGER   ::   INDNFT1COL1,INDNFT1COL2,INDNFT1COL3,INDNFT1COL4,INDNFT1COL5
INTEGER   ::   INDNFT1COL6,INDNFT1COL7,INDNFT1COL8,INDNFT1COL9,INDNFT1COL10
INTEGER   ::   INDNFT1COL11,INDNFT1COL12,INDNFT1COL13,INDNFT1COL14,INDNFT1COL15
INTEGER   ::   INDXFT1LW1,INDXFT1LW2,INDXFT1LW3,INDXFT1LW4,INDXFT1LW5
INTEGER   ::   INDXFT1LW6,INDXFT1LW7,INDXFT1LW8,INDXFT1LW9,INDXFT1LW10
INTEGER   ::   INDXFT1LW11,INDXFT1LW12,INDXFT1LW13,INDXFT1LW14,INDXFT1LW15
INTEGER   ::   INDCFT1TIT1,INDCFT1TIT2,INDCFT1TIT3,INDCFT1TIT4,INDCFT1TIT5
INTEGER   ::   INDCFT1TIT6,INDCFT1TIT7,INDCFT1TIT8,INDCFT1TIT9,INDCFT1TIT10
INTEGER   ::   INDCFT1TIT11,INDCFT1TIT12,INDCFT1TIT13,INDCFT1TIT14,INDCFT1TIT15
INTEGER   ::   INDXVPTFT1L,INDXVPTFT1R,INDXVPTFT1B,INDXVPTFT1T
INTEGER   ::   INDLVPTFT1USER, INDLVARNPVUSER, INDNSTYLINZEROPV
INTEGER   ::   INDCVARNPV1,INDCVARNPV2,INDCVARNPV3,INDCVARNPV4,INDCVARNPV5
INTEGER   ::   INDCVARNPV6,INDCVARNPV7,INDCVARNPV8,INDCVARNPV9,INDCVARNPV10
INTEGER   ::   INDCVARNPV11,INDCVARNPV12,INDCVARNPV13,INDCVARNPV14,INDCVARNPV15
INTEGER   ::   INDCVARNPH1,INDCVARNPH2,INDCVARNPH3,INDCVARNPH4,INDCVARNPH5
INTEGER   ::   INDCVARNPH6,INDCVARNPH7,INDCVARNPH8,INDLVARNPHUSER
INTEGER   ::   INDNPHCOL1,INDNPHCOL2,INDNPHCOL3,INDNPHCOL4,INDNPHCOL5
INTEGER   ::   INDNPHCOL6,INDNPHCOL7,INDNPHCOL8
INTEGER   ::   INDNPHSTY1,INDNPHSTY2,INDNPHSTY3,INDNPHSTY4,INDNPHSTY5
INTEGER   ::   INDNPHSTY6,INDNPHSTY7,INDNPHSTY8
INTEGER   ::   INDLPHCOLUSER,INDLPHSTYUSER
Gaelle Tanguy's avatar
Gaelle Tanguy committed
INTEGER   ::   INDLPATCH
#ifdef RHODES
INTEGER          :: ISTAF
#endif
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
INTEGER   ::   INDXXPART,INDXYPART,INDXZPART,INDLTRAJ3D,INDLFLUX3D
INTEGER   ::   INDLTRAJ_GROUP
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
INTEGER   ::   INBMIN, J, JJ, JM
INTEGER   ::   ISTA
INTEGER   ::   INBV, IND9999
INTEGER   ::   INDQ1,INDQ2
INTEGER,DIMENSION(30,100) :: IIMIN
!!!!!!!!!!!!JOEL!!!!!!!!!!
Gaelle Tanguy's avatar
Gaelle Tanguy committed
INTEGER,DIMENSION(603)     :: IT
346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000
REAL,DIMENSION(100) :: ZISOLEV
REAL      ::   ZISO, ZX, ZY
LOGICAl   ::   GXI=.FALSE., GXJ=.FALSE.
!!! NOVEMBRE 2009 G. TANGUY
INTEGER :: INDL90TITYT,INDL90TITYM,INDL90TITYB
!------------------------------------------------------------------------------
!
YCARIN=' '
YCARIN = ADJUSTL(HCARIN)
YCARIN = ADJUSTL(YCARIN)
ILENC = LEN_TRIM(YCARIN)
if(nverbia >0)then
print *,' **entree caresolv YCARIN(1:LEN_TRIM(YCARIN)) ',YCARIN(1:ILENC)
endif
!
INDID = INDEX(YCARIN,'_IDEM')
IF(INDID /= 0)THEN
  CALL CARMEMORY(YCAR,2)
  YCAR=ADJUSTL(YCAR)
! print *,' AP CARMEMORY YCAR ',YCAR(1:LEN_TRIM(YCAR))
  ILENGRP=0
! ILENGRP=INDEX(YCAR,'_')-1
  ILENC2 = LEN_TRIM(YCAR)
  YCARIN(INDID:INDID+ILENC2-ILENGRP-1)=YCAR(ILENGRP+1:ILENC2)
  ILENC = LEN_TRIM(YCARIN)

  DO J=1,20
    IF((INDID-1+ILENC2-ILENGRP) >= ILENC)THEN
      EXIT
    ELSE
      YCARIN(ILENC:ILENC)=' '
      ILENC=ILENC-1
    ENDIF
  ENDDO

ENDIF

INDNSZLBX=INDEX(YCARIN,'NSZLBX')
INDNSZLBY=INDEX(YCARIN,'NSZLBY')
INDIINF=INDEX(YCARIN,'NIINF') 
INDISUP=INDEX(YCARIN,'NISUP')
INDJINF=INDEX(YCARIN,'NJINF')
INDJSUP=INDEX(YCARIN,'NJSUP')
INDNHEURXLBL=INDEX(YCARIN,'NHEURXLBL')
INDNHEURXGRAD=INDEX(YCARIN,'NHEURXGRAD')
INDIDEBCOU=INDEX(YCARIN,'NIDEBCOU') 
INDJDEBCOU=INDEX(YCARIN,'NJDEBCOU') 
INDXTIMEMIN=INDEX(YCARIN,'XTIMEMIN')
INDXTIMEMAX=INDEX(YCARIN,'XTIMEMAX')
INDXIDEBCOU=INDEX(YCARIN,'XIDEBCOU') 
INDXJDEBCOU=INDEX(YCARIN,'XJDEBCOU') 
INDXSSP=INDEX(YCARIN,'XSSP')
INDXARLSTR=INDEX(YCARIN,'XARLSTR')
INDXLWSTR=INDEX(YCARIN,'XLWSTR')
INDXANGULVT=INDEX(YCARIN,'XANGULVT')
INDXFACTAXEX=INDEX(YCARIN,'XFACTAXEX')
INDXFACTAXEY=INDEX(YCARIN,'XFACTAXEY')
INDXAXEXUSERD=INDEX(YCARIN,'XAXEXUSERD')
INDXAXEYUSERD=INDEX(YCARIN,'XAXEYUSERD')
INDXAXEXUSERF=INDEX(YCARIN,'XAXEXUSERF')
INDXAXEYUSERF=INDEX(YCARIN,'XAXEYUSERF')
INDXSIZEL=INDEX(YCARIN,'XSIZEL')
INDXSZTITXL=INDEX(YCARIN,'XSZTITXL')
INDXSZTITXM=INDEX(YCARIN,'XSZTITXM')
INDXSZTITXR=INDEX(YCARIN,'XSZTITXR')
INDXSZTITYT=INDEX(YCARIN,'XSZTITYT')
INDXSZTITYM=INDEX(YCARIN,'XSZTITYM')
INDXSZTITYB=INDEX(YCARIN,'XSZTITYB')

INDXSZTITT1=INDEX(YCARIN,'XSZTITT1')
INDXSZTITT2=INDEX(YCARIN,'XSZTITT2')
INDXSZTITT3=INDEX(YCARIN,'XSZTITT3')
INDXPOSTITYT=INDEX(YCARIN,'XPOSTITYT')
INDXPOSTITYM=INDEX(YCARIN,'XPOSTITYM')
INDXPOSTITYB=INDEX(YCARIN,'XPOSTITYB')
INDXPOSTITT1=INDEX(YCARIN,'XPOSTITT1')
INDXPOSTITT2=INDEX(YCARIN,'XPOSTITT2')
INDXPOSTITT3=INDEX(YCARIN,'XPOSTITT3')
INDXYPOSTITT1=INDEX(YCARIN,'XYPOSTITT1')
INDXYPOSTITT2=INDEX(YCARIN,'XYPOSTITT2')
INDXYPOSTITT3=INDEX(YCARIN,'XYPOSTITT3')

INDXSZTITB1=INDEX(YCARIN,'XSZTITB1')
INDXSZTITB2=INDEX(YCARIN,'XSZTITB2')
INDXSZTITB3=INDEX(YCARIN,'XSZTITB3')
INDXPOSTITB1=INDEX(YCARIN,'XPOSTITB1')
INDXPOSTITB2=INDEX(YCARIN,'XPOSTITB2')
INDXPOSTITB3=INDEX(YCARIN,'XPOSTITB3')
INDXYPOSTITYT=INDEX(YCARIN,'XYPOSTITYT')
INDXYPOSTITYM=INDEX(YCARIN,'XYPOSTITYM')
INDXYPOSTITYB=INDEX(YCARIN,'XYPOSTITYB')
INDXYPOSTITB1=INDEX(YCARIN,'XYPOSTITB1')
INDXYPOSTITB2=INDEX(YCARIN,'XYPOSTITB2')
INDXYPOSTITB3=INDEX(YCARIN,'XYPOSTITB3')
INDXSZTITVAR1=INDEX(YCARIN,'XSZTITVAR1')
INDXSZTITVAR2=INDEX(YCARIN,'XSZTITVAR2')
INDXSZTITVAR3=INDEX(YCARIN,'XSZTITVAR3')
INDXSZTITVAR4=INDEX(YCARIN,'XSZTITVAR4')
INDXSZTITVAR5=INDEX(YCARIN,'XSZTITVAR5')
INDXSZTITVAR6=INDEX(YCARIN,'XSZTITVAR6')
INDXSZTITVAR7=INDEX(YCARIN,'XSZTITVAR7')
INDXSZTITVAR8=INDEX(YCARIN,'XSZTITVAR8')
INDXPOSTITVAR1=INDEX(YCARIN,'XPOSTITVAR1')
INDXPOSTITVAR2=INDEX(YCARIN,'XPOSTITVAR2')
INDXPOSTITVAR3=INDEX(YCARIN,'XPOSTITVAR3')
INDXYPOSTITVAR1=INDEX(YCARIN,'XYPOSTITVAR1')
INDXYPOSTITVAR2=INDEX(YCARIN,'XYPOSTITVAR2')
INDXYPOSTITVAR3=INDEX(YCARIN,'XYPOSTITVAR3')
INDXPOSTITVAR4=INDEX(YCARIN,'XPOSTITVAR4')
INDXPOSTITVAR5=INDEX(YCARIN,'XPOSTITVAR5')
INDXPOSTITVAR6=INDEX(YCARIN,'XPOSTITVAR6')
INDXYPOSTITVAR4=INDEX(YCARIN,'XYPOSTITVAR4')
INDXYPOSTITVAR5=INDEX(YCARIN,'XYPOSTITVAR5')
INDXYPOSTITVAR6=INDEX(YCARIN,'XYPOSTITVAR6')
INDXPOSTITVAR7=INDEX(YCARIN,'XPOSTITVAR7')
INDXPOSTITVAR8=INDEX(YCARIN,'XPOSTITVAR8')
INDXYPOSTITVAR7=INDEX(YCARIN,'XYPOSTITVAR7')
INDXYPOSTITVAR8=INDEX(YCARIN,'XYPOSTITVAR8')
!*JD*Mars 2009
INDXPOSXVARNPV1TOP=INDEX(YCARIN,'XPOSXVARNPV1TOP')
INDXPOSYVARNPV1TOP=INDEX(YCARIN,'XPOSYVARNPV1TOP')
INDXPOSXVARNPV5BOT=INDEX(YCARIN,'XPOSXVARNPV5BOT')
INDXPOSYVARNPV5BOT=INDEX(YCARIN,'XPOSYVARNPV5BOT')
INDXSZVARNPVTOP=INDEX(YCARIN,'XSZVARNPVTOP')
INDXSZVARNPVBOT=INDEX(YCARIN,'XSZVARNPVBOT')
!*JD*Mars 2009

INDXIDEBCVLL=INDEX(YCARIN,'XIDEBCVLL')
INDXIDEBCV=0
IF(INDXIDEBCVLL == 0)THEN
  INDXIDEBCV=INDEX(YCARIN,'XIDEBCV')
ENDIF
INDXJDEBCVLL=INDEX(YCARIN,'XJDEBCVLL')
INDXJDEBCV=0
IF(INDXJDEBCVLL == 0)THEN
  INDXJDEBCV=INDEX(YCARIN,'XJDEBCV')
ENDIF
INDXIFINCVLL=INDEX(YCARIN,'XIFINCVLL')
INDXIFINCV=0
IF(INDXIFINCVLL == 0)THEN
  INDXIFINCV=INDEX(YCARIN,'XIFINCV')
ENDIF
INDXJFINCVLL=INDEX(YCARIN,'XJFINCVLL')
INDXJFINCV=0
IF(INDXJFINCVLL == 0)THEN
  INDXJFINCV=INDEX(YCARIN,'XJFINCV')
ENDIF
INDNIDEBCV=INDEX(YCARIN,'NIDEBCV')
INDNJDEBCV=INDEX(YCARIN,'NJDEBCV')
INDNIFINCV=INDEX(YCARIN,'NIFINCV')
INDNJFINCV=INDEX(YCARIN,'NJFINCV')
INDXAMX=INDEX(YCARIN,'XAMX')
INDXVHC=INDEX(YCARIN,'XVHC=')
IF(INDXVHC == 0)THEN
  INDXVHC=INDEX(YCARIN,'XVHC =')
ENDIF
INDXVHCPH=INDEX(YCARIN,'XVHCPH')
INDXVLC=INDEX(YCARIN,'XVLC')
INDXVRL=INDEX(YCARIN,'XVRL=')
IF(INDXVRL == 0)THEN
  INDXVRL=INDEX(YCARIN,'XVRL =')
ENDIF
INDXVRLPH=INDEX(YCARIN,'XVRLPH')
INDVISU=INDEX(YCARIN,'VISU')
INDNOVISU=INDEX(YCARIN,'NOVISU')
INDNVERBIA=INDEX(YCARIN,'NVERBIA')
INDXPMIN=INDEX(YCARIN,'XPMIN')
INDXPMAX=INDEX(YCARIN,'XPMAX')
INDXPINT=INDEX(YCARIN,'XPINT')
INDXHMIN=INDEX(YCARIN,'XHMIN')
INDXHMAX=INDEX(YCARIN,'XHMAX')
INDXLATRAD1=INDEX(YCARIN,'XLATRAD1')
INDXLATRAD2=INDEX(YCARIN,'XLATRAD2')
INDXLATRAD3=INDEX(YCARIN,'XLATRAD3')
INDXLATRAD4=INDEX(YCARIN,'XLATRAD4')
INDXLONRAD1=INDEX(YCARIN,'XLONRAD1')
INDXLONRAD2=INDEX(YCARIN,'XLONRAD2')
INDXLONRAD3=INDEX(YCARIN,'XLONRAD3')
INDXLONRAD4=INDEX(YCARIN,'XLONRAD4')
INDXSPVALT=INDEX(YCARIN,'XSPVALT')
INDXSPVAL=INDEX(YCARIN,'XSPVAL=')
IF(INDXSPVAL == 0)THEN
  INDXSPVAL=INDEX(YCARIN,'XSPVAL =')
  IF(INDXSPVAL == 0)THEN
    INDXSPVAL=INDEX(YCARIN,'XSPVAL  =')
  ENDIF
ENDIF
INDXPORTRAD1=INDEX(YCARIN,'XPORTRAD1')
INDXPORTRAD2=INDEX(YCARIN,'XPORTRAD2')
INDXPORTRAD3=INDEX(YCARIN,'XPORTRAD3')
INDXPORTRAD4=INDEX(YCARIN,'XPORTRAD4')
INDXLWRAD1=INDEX(YCARIN,'XLWRAD1')
INDXLWRAD2=INDEX(YCARIN,'XLWRAD2')
INDXLWRAD3=INDEX(YCARIN,'XLWRAD3')
INDXLWRAD4=INDEX(YCARIN,'XLWRAD4')
INDXISOMIN_=INDEX(YCARIN,'XISOMIN_')
INDXISOMAX_=INDEX(YCARIN,'XISOMAX_')
INDXISOMIN=INDEX(YCARIN,'XISOMIN')
INDXISOMAX=INDEX(YCARIN,'XISOMAX')
IF(INDXISOMIN_ /= 0)THEN
  INDXISOMIN=0
ENDIF
IF(INDXISOMAX_ /= 0)THEN
  INDXISOMAX=0
ENDIF
INDNLANGLE=INDEX(YCARIN,'NLANGLE')
INDNLMAX=INDEX(YCARIN,'NLMAX')
INDNZSTR=INDEX(YCARIN,'NZSTR')
INDNARSTR=INDEX(YCARIN,'NARSTR')
INDNIOFFD=INDEX(YCARIN,'NIOFFD')
INDNSD=INDEX(YCARIN,'NSD')
INDNULBLL=INDEX(YCARIN,'NULBLL')
INDNDOT=INDEX(YCARIN,'NDOT')
INDNISKIP=INDEX(YCARIN,'NISKIP=')
IF(INDNISKIP == 0)THEN
  INDNISKIP=INDEX(YCARIN,'NISKIP =')
ENDIF
INDNISKIPVX=INDEX(YCARIN,'NISKIPVX')
INDNISKIPVY=INDEX(YCARIN,'NISKIPVY')
INDNHI=INDEX(YCARIN,'NHI')
INDNIMNMX=INDEX(YCARIN,'NIMNMX')
INDNFT1ITVXMJ=INDEX(YCARIN,'NFT1ITVXMJ')
INDNFT1ITVXMN=INDEX(YCARIN,'NFT1ITVXMN')
INDNFT1ITVYMJ=INDEX(YCARIN,'NFT1ITVYMJ')
INDNFT1ITVYMN=INDEX(YCARIN,'NFT1ITVYMN')
INDNFTITVXMJ=INDEX(YCARIN,'NFTITVXMJ')
INDNFTITVXMN=INDEX(YCARIN,'NFTITVXMN')
INDNFTITVYMJ=INDEX(YCARIN,'NFTITVYMJ')
INDNFTITVYMN=INDEX(YCARIN,'NFTITVYMN')
INDNCHITVXMJ=INDEX(YCARIN,'NCHITVXMJ')
INDNCHITVXMN=INDEX(YCARIN,'NCHITVXMN')
INDNCHITVYMJ=INDEX(YCARIN,'NCHITVYMJ')
INDNCHITVYMN=INDEX(YCARIN,'NCHITVYMN')
INDNCHPCITVXMJ=INDEX(YCARIN,'NCHPCITVXMJ')
INDNCHPCITVXMN=INDEX(YCARIN,'NCHPCITVXMN')
INDNCHPCITVYMJ=INDEX(YCARIN,'NCHPCITVYMJ')
INDNCHPCITVYMN=INDEX(YCARIN,'NCHPCITVYMN')
INDNCVITVXMJ=INDEX(YCARIN,'NCVITVXMJ')
INDNCVITVXMN=INDEX(YCARIN,'NCVITVXMN')
INDNCVITVYMJ=INDEX(YCARIN,'NCVITVYMJ')
INDNCVITVYMN=INDEX(YCARIN,'NCVITVYMN')
INDNPVITVXMJ=INDEX(YCARIN,'NPVITVXMJ')
INDNPVITVXMN=INDEX(YCARIN,'NPVITVXMN')
INDNPVITVYMJ=INDEX(YCARIN,'NPVITVYMJ')
INDNPVITVYMN=INDEX(YCARIN,'NPVITVYMN')
INDNXYITVXMJ=INDEX(YCARIN,'NXYITVXMJ')
INDNXYITVXMN=INDEX(YCARIN,'NXYITVXMN')
INDNXYITVYMJ=INDEX(YCARIN,'NXYITVYMJ')
INDNXYITVYMN=INDEX(YCARIN,'NXYITVYMN')
INDNMASKITVXMJ=INDEX(YCARIN,'NMASKITVXMJ')
INDNMASKITVXMN=INDEX(YCARIN,'NMASKITVXMN')
INDNMASKITVYMJ=INDEX(YCARIN,'NMASKITVYMJ')
INDNMASKITVYMN=INDEX(YCARIN,'NMASKITVYMN')
INDXDIAINT_=INDEX(YCARIN,'XDIAINT_')
INDXDIAINT=INDEX(YCARIN,'XDIAINT')
IF(INDXDIAINT_ /= 0)THEN
  INDXDIAINT=0
ENDIF
!!!!!!!!!!
INDLVARNPVUSER=INDEX(YCARIN,'LVARNPVUSER')
INDLVPTFT1USER=INDEX(YCARIN,'LVPTFT1USER')
INDLVARNPHUSER=INDEX(YCARIN,'LVARNPHUSER')
!!!!!!!!!!
INDLM5S3=INDEX(YCARIN,'LM5S3')
INDLSYMBTEXTG=INDEX(YCARIN,'LSYMBTEXTG')
INDLCVZOOM=INDEX(YCARIN,'LCVZOOM')
INDLVST=INDEX(YCARIN,'LVST')
INDLDILW=INDEX(YCARIN,'LDILW')
INDLXYNVARTOP=INDEX(YCARIN,'LXYNVARTOP')
INDLXYSTYLTOP=INDEX(YCARIN,'LXYSTYLTOP')
INDLXYWINCUR=INDEX(YCARIN,'LXYWINCUR')
INDLVSUPSCA=INDEX(YCARIN,'LVSUPSCA')
INDLSYMB=INDEX(YCARIN,'LSYMB=')
IF(INDLSYMB == 0)THEN
  INDLSYMB=INDEX(YCARIN,'LSYMB =')
ENDIF
INDLTEXTG=INDEX(YCARIN,'LTEXTG')
INDLTEXTIT=INDEX(YCARIN,'LTEXTIT')
INDLTRACECV=INDEX(YCARIN,'LTRACECV')
INDLSTI=INDEX(YCARIN,'LSTI')
INDLSEGM=INDEX(YCARIN,'OLSEGM')
IF(INDLSEGM == 0)THEN
  INDLSEGM=INDEX(YCARIN,'LSEGM')
ELSE
  INDLSEGM=0
ENDIF
INDLXY=INDEX(YCARIN,'LXY=')
IF(INDLXY == 0)THEN
  INDLXY=INDEX(YCARIN,'LXY =')
ENDIF
INDLXZ=INDEX(YCARIN,'LXZ')
INDLVECTMNMX=INDEX(YCARIN,'LVECTMNMX')
INDLISO=INDEX(YCARIN,'LISO=')
IF(INDLISO == 0)THEN
  INDLISO=INDEX(YCARIN,'LISO =')
ENDIF
INDLANIMK=INDEX(YCARIN,'LANIMK')
INDLANIMT=INDEX(YCARIN,'LANIMT')
INDLMINMAX=INDEX(YCARIN,'LMINMAX')
INDATFILE=INDEX(YCARIN,'LDATFILE')
INDLINTERPTOP=INDEX(YCARIN,'LINTERPTOP=')
INDLSPSECT=INDEX(YCARIN,'LSPSECT')
INDLSPVALT=INDEX(YCARIN,'LSPVALT')
INDLPRESY=INDEX(YCARIN,'LPRESY')
INDLEGVECT=INDEX(YCARIN,'LEGVECT')
INDLSTREAM=INDEX(YCARIN,'LSTREAM')
INDLINTERPOLSTR=INDEX(YCARIN,'LINTERPOLSTR')
INDLNOLBLBAR=INDEX(YCARIN,'LNOLBLBAR')
INDLNOLABELX=INDEX(YCARIN,'LNOLABELX')
INDLNOLABELY=INDEX(YCARIN,'LNOLABELY')
INDLNOUVRS=INDEX(YCARIN,'LNOUVRS')
INDLMYHEURX=INDEX(YCARIN,'LMYHEURX')
INDLRADAR=INDEX(YCARIN,'LRADAR')
INDLRADIST=INDEX(YCARIN,'LRADIST')
INDLRADRAY=INDEX(YCARIN,'LRADRAY')
INDLCOLAREA=INDEX(YCARIN,'LCOLAREA=')
IF(INDLCOLAREA == 0)THEN
  INDLCOLAREA=INDEX(YCARIN,'LCOLAREA =')
ENDIF
INDLFT1LUSER=INDEX(YCARIN,'LFT1LUSER')
INDLFTBAUTO=INDEX(YCARIN,'LFTBAUTO')
INDLFT1BAUTO=INDEX(YCARIN,'LFT1BAUTO')
INDLTITFTUSER=INDEX(YCARIN,'LTITFTUSER')
INDLPHCOLUSER=INDEX(YCARIN,'LPHCOLUSER')
INDLPHSTYUSER=INDEX(YCARIN,'LPHSTYUSER')
INDLXABSC=INDEX(YCARIN,'LXABSC')
INDLXMINTOP=INDEX(YCARIN,'LXMINTOP')
INDLABEL1=INDEX(YCARIN,'LABEL1')
INDLDEFCV2LL=INDEX(YCARIN,'LDEFCV2LL')
INDLDEFCV2IND=INDEX(YCARIN,'LDEFCV2IND')
INDLDEFCV2=0
INDLDEFCV2=INDEX(YCARIN,'LDEFCV2=')
IF(INDLDEFCV2 == 0)THEN
  INDLDEFCV2=INDEX(YCARIN,'LDEFCV2 =')
ENDIF
IF(INDLDEFCV2 == 0)THEN
  INDLDEFCV2=INDEX(YCARIN,'LDEFCV2  =')
ENDIF
IF(NVERBIA > 0)THEN
print *,' INDLDEFCV2LL,INDLDEFCV2IND,INDLDEFCV2 ',INDLDEFCV2LL,INDLDEFCV2IND,INDLDEFCV2
ENDIF
INDLBLUSER1=INDEX(YCARIN,'LBLUSER1')
INDLBLUSER2=INDEX(YCARIN,'LBLUSER2')
INDLBLUSER3=INDEX(YCARIN,'LBLUSER3')
INDLBLUSER4=INDEX(YCARIN,'LBLUSER4')
INDLINDSP=INDEX(YCARIN,'LINDSP')
INDLINDAX=INDEX(YCARIN,'LINDAX')
INDLCOLUSERUV=INDEX(YCARIN,'LCOLUSERUV')
INDLTIMEUSER=INDEX(YCARIN,'LTIMEUSER')
INDLCHREEL=INDEX(YCARIN,'LCHREEL')
INDLOGNEP=INDEX(YCARIN,'LOGNEP')
INDLCOLISONE=INDEX(YCARIN,'LCOLISONE')
INDLCOLRSONE=INDEX(YCARIN,'LCOLRSONE')
INDLCOLRS1ONE=INDEX(YCARIN,'LCOLRS1ONE')
INDLCOLINE=INDEX(YCARIN,'LCOLINE')
INDL24H=INDEX(YCARIN,'L24H')
INDLCONT=INDEX(YCARIN,'LCONT')
INDL2CONT=INDEX(YCARIN,'L2CONT')
INDLRELIEF=INDEX(YCARIN,'LRELIEF')
INDLCONV2XY=INDEX(YCARIN,'LCONV2XY')
INDLCONVG2MASS=INDEX(YCARIN,'LCONVG2MASS')
INDLCOLZERO=INDEX(YCARIN,'LCOLZERO')
INDL3D=INDEX(YCARIN,'L3D')
INDLMARKER=INDEX(YCARIN,'LMARKER')
INDLPRDAT=INDEX(YCARIN,'LPRDAT')
INDLPRINTXY=INDEX(YCARIN,'LPRINTXY')
INDLPRINT=INDEX(YCARIN,'LPRINT')
IF((INDLPRINTXY == 0 .AND. INDLPRINT == 0) .OR. (INDLPRINTXY /= INDLPRINT))THEN
ELSE
 INDLPRINT=INDEX(YCARIN(INDLPRINTXY+1:LEN_TRIM(YCARIN)),'LPRINT')
 IF(INDLPRINT /= 0)THEN
   INDLPRINT=INDLPRINT+INDLPRINTXY
 ENDIF
ENDIF
INDLPOINTG=INDEX(YCARIN,'LPOINTG')
INDLXYO=INDEX(YCARIN,'LXYO')
INDL2DBX=INDEX(YCARIN,'L2DBX')
INDL2DBY=INDEX(YCARIN,'L2DBY')
INDLTIMEUSER=INDEX(YCARIN,'LTIMEUSER=')
INDLCOLUSER=INDEX(YCARIN,'LCOLUSER=')
IF(INDLCOLUSER == 0)THEN
  INDLCOLUSER=INDEX(YCARIN,'LCOLUSER =')
ENDIF
IF(INDLCOLUSER == 0)THEN
  INDLCOLUSER=INDEX(YCARIN,'LCOLUSER  =')
ENDIF
INDLISOWHI2=INDEX(YCARIN,'LISOWHI2')
INDLISOWHI3=INDEX(YCARIN,'LISOWHI3')
INDLISOWHI=INDEX(YCARIN,'LISOWHI=')
IF(INDLISOWHI == 0)THEN
  INDLISOWHI=INDEX(YCARIN,'LISOWHI =')
ENDIF
IF(INDLISOWHI == 0)THEN
  INDLISOWHI=INDEX(YCARIN,'LISOWHI  =')
ENDIF
IF(INDLISOWHI == 0)THEN
  INDLISOWHI=INDEX(YCARIN,'LISOWHI   =')
ENDIF
INDLCOLBR=INDEX(YCARIN,'LCOLBR')
INDLCOLAREASEL=INDEX(YCARIN,'LCOLAREASEL')
INDLCOLINESEL=INDEX(YCARIN,'LCOLINESEL')
INDLINVWB=INDEX(YCARIN,'LINVWB')
INDLINVPTIR=INDEX(YCARIN,'LINVPTIR')
INDLDOMAIN=INDEX(YCARIN,'LDOMAIN')
INDLGEOG=INDEX(YCARIN,'LGEOG')
INDLBLFT1SUP=INDEX(YCARIN,'LBLFT1SUP')
INDLMASK3D_XY=INDEX(YCARIN,'LMASK3D_XY')
INDLMASK3D_XZ=INDEX(YCARIN,'LMASK3D_XZ')
INDLMASK3D_YZ=INDEX(YCARIN,'LMASK3D_YZ')
INDLMASK3D=INDEX(YCARIN,'LMASK3D=')
IF(INDLMASK3D == 0)THEN
  INDLMASK3D=INDEX(YCARIN,'LMASK3D =')
ENDIF
IF(INDLMASK3D == 0)THEN
  INDLMASK3D=INDEX(YCARIN,'LMASK3D  =')
ENDIF
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
INDLTRAJ3D=INDEX(YCARIN,'LTRAJ3D=')
IF(INDLTRAJ3D == 0)THEN
  INDLTRAJ3D=INDEX(YCARIN,'LTRAJ3D =')
ENDIF
IF(INDLTRAJ3D == 0)THEN
  INDLTRAJ3D=INDEX(YCARIN,'LTRAJ3D  =')
ENDIF
!
INDLFLUX3D=INDEX(YCARIN,'LFLUX3D=')
IF(INDLFLUX3D == 0)THEN
  INDLFLUX3D=INDEX(YCARIN,'LFLUX3D =')
ENDIF
IF(INDLFLUX3D == 0)THEN
  INDLFLUX3D=INDEX(YCARIN,'LFLUX3D  =')
ENDIF
!
INDLTRAJ_GROUP=INDEX(YCARIN,'LTRAJ_GROUP=')
IF(INDLTRAJ_GROUP == 0)THEN
  INDLTRAJ_GROUP=INDEX(YCARIN,'LTRAJ_GROUP =')
ENDIF
IF(INDLTRAJ_GROUP == 0)THEN
  INDLTRAJ_GROUP=INDEX(YCARIN,'LTRAJ_GROUP  =')
ENDIF
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
!IF(INDLMASK3D /= 0)THEN
! LMASK3D_XY=.TRUE.
! LMASK3D_XZ=.TRUE.
! LMASK3D_YZ=.TRUE.
!ENDIF
INDLXYZ00=INDEX(YCARIN,'LXYZ00')
INDLFT3C=INDEX(YCARIN,'LFT3C')
INDLFT4C=INDEX(YCARIN,'LFT4C')
INDLFTCLIP=INDEX(YCARIN,'LFTCLIP')
INDLFT1STYLUSER=INDEX(YCARIN,'LFT1STYLUSER')
INDLFTSTYLUSER=INDEX(YCARIN,'LFTSTYLUSER')
INDLHEURX=INDEX(YCARIN,'LHEURX')
INDLHACH1=INDEX(YCARIN,'LHACH1')
INDLHACH2=INDEX(YCARIN,'LHACH2')
INDLHACH3=INDEX(YCARIN,'LHACH3')
INDLHACH4=INDEX(YCARIN,'LHACH4')
INDLHACHSEL=INDEX(YCARIN,'LHACHSEL')
INDLGREY=INDEX(YCARIN,'LGREY')
INDLTABCOLDEF2=INDEX(YCARIN,'LTABCOLDEF2')
IF(INDLTABCOLDEF2 == 0)THEN
INDLTABCOLDEF=INDEX(YCARIN,'LTABCOLDEF')
ELSE
LTABCOLDEF=.TRUE.
ENDIF
INDLTABCOLDEF=INDEX(YCARIN,'LTABCOLDEF=')
IF(INDLTABCOLDEF == 0)THEN
  INDLTABCOLDEF=INDEX(YCARIN,'LTABCOLDEF =')
ENDIF
INDLMNMXUSER=INDEX(YCARIN,'LMNMXUSER')
!Mars 2009
INDLINZEROPV=INDEX(YCARIN,'LINZEROPV')
INDNSTYLINZEROPV=INDEX(YCARIN,'NSTYLINZEROPV')
IF(INDNSTYLINZEROPV /= 0 .AND. (ABS(INDNSTYLINZEROPV-INDLINZEROPV) == 4))THEN
  INDLINZEROPV=0
ENDIF
!Mars 2009
INDLMNMXLOC=INDEX(YCARIN,'LMNMXLOC')
INDLULMVTMOLD=INDEX(YCARIN,'LULMVTMOLD')
INDLVPTUSER=INDEX(YCARIN,'LVPTUSER')
INDLVPTVUSER=INDEX(YCARIN,'LVPTVUSER')
INDLVPTPVUSER=INDEX(YCARIN,'LVPTPVUSER')
INDLVPTXYUSER=INDEX(YCARIN,'LVPTXYUSER')
INDLFACTIMP=INDEX(YCARIN,'LFACTIMP')
INDLFACTAXEX=INDEX(YCARIN,'LFACTAXEX')
INDLFACTAXEY=INDEX(YCARIN,'LFACTAXEY')
INDLAXEYUSER=INDEX(YCARIN,'LAXEYUSER')
INDLAXEXUSER=INDEX(YCARIN,'LAXEXUSER')
INDLFMTAXEX=INDEX(YCARIN,'LFMTAXEX')
INDLFMTAXEY=INDEX(YCARIN,'LFMTAXEY')
INDLSPOT=INDEX(YCARIN,'LSPOT')
INDNIFDC=INDEX(YCARIN,'NIFDC')
INDNDOMAINL=INDEX(YCARIN,'NDOMAINL')
INDNDOMAINR=INDEX(YCARIN,'NDOMAINR')
INDNDOMAINB=INDEX(YCARIN,'NDOMAINB')
INDNDOMAINT=INDEX(YCARIN,'NDOMAINT')
INDNIGRNC=INDEX(YCARIN,'NIGRNC')
INDNPROFILE=INDEX(YCARIN,'PROFILE')
INDXPVMIN=INDEX(YCARIN,'XPVMIN_')
INDXPVMAX=INDEX(YCARIN,'XPVMAX_')
INDXPVMINT=INDEX(YCARIN,'XPVMINT')
INDXPVMAXT=INDEX(YCARIN,'XPVMAXT')
INDXLWFTALL=INDEX(YCARIN,'XLWFTALL')
INDXLWSEGM=INDEX(YCARIN,'XLWSEGM')
INDXLWDOMAIN=INDEX(YCARIN,'XLWDOMAIN')
INDXLWTRACECV=INDEX(YCARIN,'XLWTRACECV')
INDXLWV=INDEX(YCARIN,'XLWV=')
IF(INDXLWV == 0)THEN
  INDXLWV=INDEX(YCARIN,'XLWV =')
ENDIF
IF(INDXLWV == 0)THEN
  INDXLWV=INDEX(YCARIN,'XLWV  =')
ENDIF
IF(INDXLWV == 0)THEN
  INDXLWV=INDEX(YCARIN,'XLWV   =')
ENDIF
INDXLW=INDEX(YCARIN,'XLW=')
IF(INDXLW == 0)THEN
  INDXLW=INDEX(YCARIN,'XLW =')
ENDIF
IF(INDXLW == 0)THEN
  INDXLW=INDEX(YCARIN,'XLW  =')
ENDIF
IF(INDXLW == 0)THEN
  INDXLW=INDEX(YCARIN,'XLW   =')
ENDIF
INDXLW1=INDEX(YCARIN,'XLW1')
INDXLW2=INDEX(YCARIN,'XLW2')
INDXLW3=INDEX(YCARIN,'XLW3')
INDXLW4=INDEX(YCARIN,'XLW4')
INDXLWDEF=INDEX(YCARIN,'XLWDEF')
INDXLWVDEF=INDEX(YCARIN,'XLWVDEF')
INDXLWCONT=INDEX(YCARIN,'XLWCONT')
INDXLWPV1=INDEX(YCARIN,'XLWPV1')
INDXLWPV2=INDEX(YCARIN,'XLWPV2')
INDXLWPV3=INDEX(YCARIN,'XLWPV3')
INDXLWPV4=INDEX(YCARIN,'XLWPV4')
INDXLWPV5=INDEX(YCARIN,'XLWPV5')
INDXLWPV6=INDEX(YCARIN,'XLWPV6')
INDXLWPV7=INDEX(YCARIN,'XLWPV7')
INDXLWPV8=INDEX(YCARIN,'XLWPV8')
INDXLWPV9=INDEX(YCARIN,'XLWPV9')
INDXLWPV10=INDEX(YCARIN,'XLWPV10')
INDXLWPV11=INDEX(YCARIN,'XLWPV11')
INDXLWPV12=INDEX(YCARIN,'XLWPV12')
INDXLWPV13=INDEX(YCARIN,'XLWPV13')
INDXLWPV14=INDEX(YCARIN,'XLWPV14')
INDXLWPV15=INDEX(YCARIN,'XLWPV15')
INDXSTYLPV1=INDEX(YCARIN,'XSTYLPV1')
INDXSTYLPV2=INDEX(YCARIN,'XSTYLPV2')
INDXSTYLPV3=INDEX(YCARIN,'XSTYLPV3')
INDXSTYLPV4=INDEX(YCARIN,'XSTYLPV4')
INDXSTYLPV5=INDEX(YCARIN,'XSTYLPV5')
INDXSTYLPV6=INDEX(YCARIN,'XSTYLPV6')
INDXSTYLPV7=INDEX(YCARIN,'XSTYLPV7')
INDXSTYLPV8=INDEX(YCARIN,'XSTYLPV8')
INDXSTYLPV9=INDEX(YCARIN,'XSTYLPV9')
INDXSTYLPV10=INDEX(YCARIN,'XSTYLPV10')
INDXSTYLPV11=INDEX(YCARIN,'XSTYLPV11')
INDXSTYLPV12=INDEX(YCARIN,'XSTYLPV12')
INDXSTYLPV13=INDEX(YCARIN,'XSTYLPV13')
INDXSTYLPV14=INDEX(YCARIN,'XSTYLPV14')
INDXSTYLPV15=INDEX(YCARIN,'XSTYLPV15')
INDXPARCOLUV=INDEX(YCARIN,'XPARCOLUV')
INDXISOLEV_=INDEX(YCARIN,'XISOLEV_')
INDXISOLEV=INDEX(YCARIN,'XISOLEV')
IF(INDXISOLEV_ /= 0)THEN
  INDXISOLEV=0
ENDIF
INDXICAR=INDEX(YCARIN,'XICAR')
INDXJCAR=INDEX(YCARIN,'XJCAR')
INDXLWPH1=INDEX(YCARIN,'XLWPH1')
INDXLWPH2=INDEX(YCARIN,'XLWPH2')
INDXLWPH3=INDEX(YCARIN,'XLWPH3')
INDXLWPH4=INDEX(YCARIN,'XLWPH4')
INDXLWPH5=INDEX(YCARIN,'XLWPH5')
INDXLWPH6=INDEX(YCARIN,'XLWPH6')
INDXLWPH7=INDEX(YCARIN,'XLWPH7')
INDXLWPH8=INDEX(YCARIN,'XLWPH8')
INDXLATCAR=INDEX(YCARIN,'XLATCAR')
INDXLONCAR=INDEX(YCARIN,'XLONCAR')
INDCNOMCAR=INDEX(YCARIN,'CNOMCAR')
INDCSYMCAR=INDEX(YCARIN,'CSYMCAR')
INDCSYMRAD1=INDEX(YCARIN,'CSYMRAD1')
INDCSYMRAD2=INDEX(YCARIN,'CSYMRAD2')
INDCSYMRAD3=INDEX(YCARIN,'CSYMRAD3')
INDCSYMRAD4=INDEX(YCARIN,'CSYMRAD4')
INDXPOSNOM=INDEX(YCARIN,'XPOSNOM')
INDXSZNOM=INDEX(YCARIN,'XSZNOM')
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! 19/12/2008 : modification pour controler la taille et le format des labels !!
!! pour les retrotrajectoires                                                 !!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
INDNSZRTRAJ=INDEX(YCARIN,'NSZRTRAJ')
INDLFMTRTRAJ=INDEX(YCARIN,'LFMTRTRAJ')
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
INDXSZSYM=INDEX(YCARIN,'XSZSYM')
INDICOLNOM=INDEX(YCARIN,'ICOLNOM')
INDICOLSYM=INDEX(YCARIN,'ICOLSYM')
INDNINDCOLUV=INDEX(YCARIN,'NINDCOLUV')
INDXXL=INDEX(YCARIN,'XXL')
INDXXH=INDEX(YCARIN,'XXH')
INDXVPTXYL=INDEX(YCARIN,'XVPTXYL')
INDXYL=0
IF(INDXVPTXYL == 0)THEN
  INDXYL=INDEX(YCARIN,'XYL')
! Aout 2001 Pour tenir compte de EMIS_XYLE (xylene chimie)
  IF(INDXYL /= 0)THEN
   INDTMP=0
   INDTMP=INDEX(YCARIN,'_XYLE')
   IF(INDTMP /= 0)THEN
    INDXYL=0
   ENDIF
  ENDIF
! Aout 2001 Pour tenir compte de EMIS_XYLE (xylene chimie)
ENDIF
INDXYH=INDEX(YCARIN,'XYH')
INDXZL=INDEX(YCARIN,'XZL')
INDXZH=INDEX(YCARIN,'XZH')
INDXFT_ADTIM1=INDEX(YCARIN,'XFT_ADTIM1')
INDXFT_ADTIM2=INDEX(YCARIN,'XFT_ADTIM2')
INDXFT_ADTIM3=INDEX(YCARIN,'XFT_ADTIM3')
INDXFT_ADTIM4=INDEX(YCARIN,'XFT_ADTIM4')
INDXFT_ADTIM5=INDEX(YCARIN,'XFT_ADTIM5')
INDXFT_ADTIM6=INDEX(YCARIN,'XFT_ADTIM6')
INDXFT_ADTIM7=INDEX(YCARIN,'XFT_ADTIM7')
INDXFT_ADTIM8=INDEX(YCARIN,'XFT_ADTIM8')
INDXFT1_ADTIM1=INDEX(YCARIN,'XFT1_ADTIM1')
INDXFT1_ADTIM2=INDEX(YCARIN,'XFT1_ADTIM2')
INDXFT1_ADTIM3=INDEX(YCARIN,'XFT1_ADTIM3')
INDXFT1_ADTIM4=INDEX(YCARIN,'XFT1_ADTIM4')
INDXFT1_ADTIM5=INDEX(YCARIN,'XFT1_ADTIM5')
INDXFT1_ADTIM6=INDEX(YCARIN,'XFT1_ADTIM6')
INDXFT1_ADTIM7=INDEX(YCARIN,'XFT1_ADTIM7')
INDXFT1_ADTIM8=INDEX(YCARIN,'XFT1_ADTIM8')
INDXFTMIN=INDEX(YCARIN,'XFTMIN')
INDXFTMAX=INDEX(YCARIN,'XFTMAX')
INDXFT1MIN=INDEX(YCARIN,'XFT1MIN')
INDXFT1MAX=INDEX(YCARIN,'XFT1MAX')
print *,' ***CARESOLV INDXFT1MIN,INDXFT1MAX ',INDXFT1MIN,INDXFT1MAX
INDXFT1MIN_=INDEX(YCARIN,'XFT1MIN_')
INDXFT1MAX_=INDEX(YCARIN,'XFT1MAX_')
IF(INDXFT1MIN_ /= 0)THEN
  INDXFT1MIN=0
ENDIF
IF(INDXFT1MAX_ /= 0)THEN
  INDXFT1MAX=0
ENDIF
!!!!!!!!!!!!JOEL!!!!!!!!!!
!!!!!!!!!!!!JOEL!!!!!!!!!!
INDXXPART=INDEX(YCARIN,'XXPART')
INDXYPART=INDEX(YCARIN,'XYPART')
INDXZPART=INDEX(YCARIN,'XZPART')