Changeset 3945 for LMDZ6


Ignore:
Timestamp:
Jun 18, 2021, 5:24:43 PM (3 years ago)
Author:
Laurent Fairhead
Message:

Resolving issue in infotrac concerning tracers from INCA and REPROOBUS that
need to be treated before isotopes.
AC

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3d_common/infotrac.F90

    r3924 r3945  
    287287   
    288288!jyg<
    289 ! Transfert number of tracers to Reprobus
    290 !!    IF (type_trac == 'repr') THEN
    291 !!#ifdef REPROBUS
    292 !!       CALL Init_chem_rep_trac(nbtr)
    293 !!#endif
    294 !!    ENDIF
    295 !>jyg
    296289       
    297290!
     
    300293    ALLOCATE(tnom_0(nqtrue), hadv(nqtrue), vadv(nqtrue),tnom_transp(nqtrue))
    301294
    302 !
    303 !jyg<
    304 !!    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), solsym(nbtr))
    305 !!    conv_flg(:) = 1 ! convection activated for all tracers
    306 !!    pbl_flg(:)  = 1 ! boundary layer activated for all tracers
    307 !>jyg
    308295
    309296!-----------------------------------------------------------------------
     
    334321!---------------------------------------------------------------------
    335322    IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag' .OR. type_trac == 'co2i') THEN
     323
    336324       IF(ierr.EQ.0) THEN
     325
    337326          ! Continue to read tracer.def
    338327          DO iq=1,nqtrue
     
    378367
    379368          ENDDO!DO iq=1,nqtrue
     369
    380370          CLOSE(90) 
    381371
    382372       ELSE ! Without tracer.def, set default values
     373
    383374         if (planet_type=="earth") then
    384375          ! for Earth, default is to have 4 tracers
     
    405396          tnom_transp(1) = 'dummy'
    406397         endif ! of if (planet_type=="earth")
     398
    407399       ENDIF
    408400       
     
    470462!
    471463
     464! Transfert number of tracers to Reprobus
     465    IF (type_trac == 'repr') THEN
     466#ifdef REPROBUS
     467       CALL Init_chem_rep_trac(nbtr,nqo,tnom_0)
     468#endif
     469    ENDIF
     470!
     471! Allocate variables depending on nbtr
     472!
     473    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), solsym(nbtr))
     474    conv_flg(:) = 1 ! convection activated for all tracers
     475    pbl_flg(:)  = 1 ! boundary layer activated for all tracers
     476
     477    IF (type_trac == 'inca' .OR. type_trac == 'inco') THEN   ! config_inca='aero' ou 'chem'
     478!>jyg
     479! le module de chimie fournit les noms des traceurs
     480! et les schemas d'advection associes. excepte pour ceux lus
     481! dans traceur.def
     482       IF (ierr .eq. 0) then
     483
     484          DO iq=1,nqo+nqCO2
     485
     486             write(*,*) 'infotrac 237: iq=',iq
     487             ! CRisi: ajout du nom du fluide transporteur
     488             ! mais rester retro compatible
     489             READ(90,'(I2,X,I2,X,A)',IOSTAT=IOstatus) hadv(iq),vadv(iq),tchaine
     490             write(lunout,*) 'iq,hadv(iq),vadv(iq)=',iq,hadv(iq),vadv(iq)
     491             write(lunout,*) 'tchaine=',trim(tchaine)
     492             write(*,*) 'infotrac 238: IOstatus=',IOstatus
     493             if (IOstatus.ne.0) then
     494                CALL abort_gcm('infotrac_init','Pb dans la lecture de traceur.def',1)
     495             endif
     496             ! Y-a-t-il 1 ou 2 noms de traceurs? -> On regarde s'il y a un
     497             ! espace ou pas au milieu de la chaine.
     498             continu=.true.
     499             nouveau_traceurdef=.false.
     500             iiq=1
     501
     502             do while (continu)
     503                if (tchaine(iiq:iiq).eq.' ') then
     504                  nouveau_traceurdef=.true.
     505                  continu=.false.
     506                else if (iiq.lt.LEN_TRIM(tchaine)) then
     507                  iiq=iiq+1
     508                else
     509                  continu=.false.
     510                endif
     511             enddo
     512
     513             write(*,*) 'iiq,nouveau_traceurdef=',iiq,nouveau_traceurdef
     514
     515             if (nouveau_traceurdef) then
     516                write(lunout,*) 'C''est la nouvelle version de traceur.def'
     517                tnom_0(iq)=tchaine(1:iiq-1)
     518                tnom_transp(iq)=tchaine(iiq+1:30)
     519             else
     520                write(lunout,*) 'C''est l''ancienne version de traceur.def'
     521                write(lunout,*) 'On suppose que les traceurs sont tous d''air'
     522                tnom_0(iq)=tchaine
     523                tnom_transp(iq) = 'air'
     524             endif
     525
     526             write(lunout,*) 'tnom_0(iq)=<',trim(tnom_0(iq)),'>'
     527             write(lunout,*) 'tnom_transp(iq)=<',trim(tnom_transp(iq)),'>'
     528
     529          ENDDO  !DO iq=1,nqo
     530          CLOSE(90) 
     531
     532       ELSE  !! if traceur.def doesn't exist
     533
     534          tnom_0(1)='H2Ov'
     535          tnom_transp(1) = 'air'
     536          tnom_0(2)='H2Ol'
     537          tnom_transp(2) = 'air'
     538          hadv(1) = 10
     539          hadv(2) = 10
     540          vadv(1) = 10
     541          vadv(2) = 10
     542
     543       ENDIF
     544 
     545#ifdef INCA
     546       CALL init_transport( &
     547            hadv_inca, &
     548            vadv_inca, &
     549            conv_flg_inca, &
     550            pbl_flg_inca,  &
     551            solsym_inca)
     552
     553       conv_flg(1+nqCO2:nbtr) = conv_flg_inca
     554       pbl_flg(1+nqCO2:nbtr) = pbl_flg_inca
     555       solsym(1+nqCO2:nbtr) = solsym_inca
     556
     557       IF (type_trac == 'inco') THEN
     558          conv_flg(1:nqCO2) = 1
     559          pbl_flg(1:nqCO2) = 1
     560          solsym(1:nqCO2) = 'CO2'
     561       ENDIF
     562#endif
     563
     564!jyg<
     565       DO iq = nqo+nqCO2+1, nqtrue
     566          hadv(iq) = hadv_inca(iq-nqo-nqCO2)
     567          vadv(iq) = vadv_inca(iq-nqo-nqCO2)
     568          tnom_0(iq)=solsym_inca(iq-nqo-nqCO2)
     569          tnom_transp(iq) = 'air'
     570       END DO
     571
     572    ENDIF ! (type_trac == 'inca' or 'inco')
     573
    472574!-----------------------------------------------------------------------
    473575!
     
    575677    WRITE(lunout,*) trim(modname),': Information stored in infotrac :'
    576678    WRITE(lunout,*) trim(modname),': iadv  niadv tname  ttext :'
     679
    577680    DO iq=1,nqtot
    578        WRITE(lunout,*) iadv(iq),niadv(iq),&
    579        ' ',trim(tname(iq)),' ',trim(ttext(iq))
     681       WRITE(lunout,*) iadv(iq),niadv(iq), ' ',trim(tname(iq)),' ',trim(ttext(iq))
    580682    END DO
    581683
     
    720822!    endif !if (ntraciso.gt.0) then
    721823
    722 
    723 ! Transfert number of tracers to Reprobus
    724     IF (type_trac == 'repr') THEN
    725 #ifdef REPROBUS
    726        CALL Init_chem_rep_trac(nbtr,nqo,tnom_0)
    727 #endif
    728     ENDIF
    729 !
    730 ! Allocate variables depending on nbtr
    731 !
    732     ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), solsym(nbtr))
    733     conv_flg(:) = 1 ! convection activated for all tracers
    734     pbl_flg(:)  = 1 ! boundary layer activated for all tracers
    735 !
    736 !!    ELSE  ! type_trac=inca : config_inca='aero' ou 'chem'
    737 !
    738     IF (type_trac == 'inca' .OR. type_trac == 'inco') THEN   ! config_inca='aero' ou 'chem'
    739 !>jyg
    740 ! le module de chimie fournit les noms des traceurs
    741 ! et les schemas d'advection associes. excepte pour ceux lus
    742 ! dans traceur.def
    743        IF (ierr .eq. 0) then
    744           DO iq=1,nqo+nqCO2
    745 
    746              write(*,*) 'infotrac 237: iq=',iq
    747              ! CRisi: ajout du nom du fluide transporteur
    748              ! mais rester retro compatible
    749              READ(90,'(I2,X,I2,X,A)',IOSTAT=IOstatus) hadv(iq),vadv(iq),tchaine
    750              write(lunout,*) 'iq,hadv(iq),vadv(iq)=',iq,hadv(iq),vadv(iq)
    751              write(lunout,*) 'tchaine=',trim(tchaine)
    752              write(*,*) 'infotrac 238: IOstatus=',IOstatus
    753              if (IOstatus.ne.0) then
    754                 CALL abort_gcm('infotrac_init','Pb dans la lecture de traceur.def',1)
    755              endif
    756              ! Y-a-t-il 1 ou 2 noms de traceurs? -> On regarde s'il y a un
    757              ! espace ou pas au milieu de la chaine.
    758              continu=.true.
    759              nouveau_traceurdef=.false.
    760              iiq=1
    761              do while (continu)
    762                 if (tchaine(iiq:iiq).eq.' ') then
    763                   nouveau_traceurdef=.true.
    764                   continu=.false.
    765                 else if (iiq.lt.LEN_TRIM(tchaine)) then
    766                   iiq=iiq+1
    767                 else
    768                   continu=.false.
    769                 endif
    770              enddo
    771              write(*,*) 'iiq,nouveau_traceurdef=',iiq,nouveau_traceurdef
    772              if (nouveau_traceurdef) then
    773                 write(lunout,*) 'C''est la nouvelle version de traceur.def'
    774                 tnom_0(iq)=tchaine(1:iiq-1)
    775                 tnom_transp(iq)=tchaine(iiq+1:30)
    776              else
    777                 write(lunout,*) 'C''est l''ancienne version de traceur.def'
    778                 write(lunout,*) 'On suppose que les traceurs sont tous d''air'
    779                 tnom_0(iq)=tchaine
    780                 tnom_transp(iq) = 'air'
    781              endif
    782              write(lunout,*) 'tnom_0(iq)=<',trim(tnom_0(iq)),'>'
    783              write(lunout,*) 'tnom_transp(iq)=<',trim(tnom_transp(iq)),'>'
    784 
    785           ENDDO  !DO iq=1,nqo
    786           CLOSE(90) 
    787        ELSE  !! if traceur.def doesn't exist
    788           tnom_0(1)='H2Ov'
    789           tnom_transp(1) = 'air'
    790           tnom_0(2)='H2Ol'
    791           tnom_transp(2) = 'air'
    792           hadv(1) = 10
    793           hadv(2) = 10
    794           vadv(1) = 10
    795           vadv(2) = 10
    796        ENDIF
    797  
    798 #ifdef INCA
    799        CALL init_transport( &
    800             hadv_inca, &
    801             vadv_inca, &
    802             conv_flg_inca, &
    803             pbl_flg_inca,  &
    804             solsym_inca)
    805 
    806        conv_flg(1+nqCO2:nbtr) = conv_flg_inca
    807        pbl_flg(1+nqCO2:nbtr) = pbl_flg_inca
    808        solsym(1+nqCO2:nbtr) = solsym_inca
    809 
    810        IF (type_trac == 'inco') THEN
    811           conv_flg(1:nqCO2) = 1
    812           pbl_flg(1:nqCO2) = 1
    813           solsym(1:nqCO2) = 'CO2'
    814        ENDIF
    815 #endif
    816 
    817 !jyg<
    818        DO iq = nqo+nqCO2+1, nqtrue
    819           hadv(iq) = hadv_inca(iq-nqo-nqCO2)
    820           vadv(iq) = vadv_inca(iq-nqo-nqCO2)
    821           tnom_0(iq)=solsym_inca(iq-nqo-nqCO2)
    822           tnom_transp(iq) = 'air'
    823        END DO
    824 
    825     ENDIF ! (type_trac == 'inca' or 'inco')
    826 
    827824!-----------------------------------------------------------------------
    828825! Finalize :
Note: See TracChangeset for help on using the changeset viewer.