Ignore:
Timestamp:
May 9, 2022, 12:35:40 PM (2 years ago)
Author:
dcugnet
Message:
  • Some variables are renamed or replaced by direct equivalents:
    • iso_indnum -> tracers(:)%iso_iName
    • niso_possibles -> niso
    • iqiso -> iqIsoPha ; index_trac -> itZonIso
    • ok_iso_verif -> isoCheck
    • ntraceurs_zone -> nzone ; ntraciso -> ntiso
    • qperemin -> min_qparent ; masseqmin -> min_qmass ; ratiomin -> min_ratio
  • Some renamed variables are only aliased with the older name (using USE <module>, ONLY: <oldName> => <newName>) in routines where they are repeated many times.
  • Few hard-coded indexes are now computed (examples: ilic, iso, ivap, irneb, iq_vap, iq_liq, iso_H2O, iso_HDO, iso_HTO, iso_O17, iso_O18).
  • The IF(isoCheck) test is now embedded in the check_isotopes_seq and check_isotopes_loc routines (lighter calling).
File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.F

    r4139 r4143  
    204204      TYPE(distrib),SAVE :: new_dist
    205205
    206       if (ok_iso_verif) then
    207          call check_isotopes(q0,ijb_u,ije_u,'leapfrog204: debut')
    208       endif !if (ok_iso_verif) then
     206      call check_isotopes(q0,ijb_u,ije_u,'leapfrog204: debut')
    209207     
    210208c$OMP MASTER
    … …  
    226224      itaufinp1 = itaufin +1
    227225
    228       if (ok_iso_verif) then
    229         call check_isotopes(q0,ijb_u,ije_u,'leapfrog 226')
    230       endif !if (ok_iso_verif) then
     226      call check_isotopes(q0,ijb_u,ije_u,'leapfrog 226')
    231227
    232228      itau = 0
    … …  
    243239      q=q0
    244240
    245       if (ok_iso_verif) then
    246         call check_isotopes(q,ijb_u,ije_u,'leapfrog 239')
    247       endif !if (ok_iso_verif) then
     241      call check_isotopes(q,ijb_u,ije_u,'leapfrog 239')
    248242     
    249243!      iday = day_ini+itau/day_step
    … …  
    324318      endif
    325319
    326         if (ok_iso_verif) then
    327            call check_isotopes(q,ijb_u,ije_u,'leapfrog 321')
    328         endif !if (ok_iso_verif) then
     320      call check_isotopes(q,ijb_u,ije_u,'leapfrog 321')
    329321
    330322#ifdef CPP_IOIPSL
    … …  
    406398
    407399
    408         if (ok_iso_verif) then
    409            call check_isotopes(q,ijb_u,ije_u,'leapfrog 400')
    410         endif !if (ok_iso_verif) then
     400         call check_isotopes(q,ijb_u,ije_u,'leapfrog 400')
    411401
    412402   2  CONTINUE ! Matsuno backward or leapfrog step begins here
    413403
    414404
    415         if (ok_iso_verif) then
    416            call check_isotopes(q,ijb_u,ije_u,'leapfrog 402')
    417         endif !if (ok_iso_verif) then
     405      call check_isotopes(q,ijb_u,ije_u,'leapfrog 402')
    418406
    419407c$OMP MASTER
    … …  
    497485
    498486
    499         if (ok_iso_verif) then
    500            call check_isotopes(q,ijb_u,ije_u,'leapfrog 471')
    501         endif !if (ok_iso_verif) then
     487      call check_isotopes(q,ijb_u,ije_u,'leapfrog 471')
    502488
    503489!ym  PAS D'AJUSTEMENT POUR LE MOMENT     
    … …  
    619605     
    620606     
    621         if (ok_iso_verif) then
    622            call check_isotopes(q,ijb_u,ije_u,'leapfrog 589')
    623         endif !if (ok_iso_verif) then
     607      call check_isotopes(q,ijb_u,ije_u,'leapfrog 589')
    624608     
    625609c-----------------------------------------------------------------------
    … …  
    684668      CALL geopot_loc  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
    685669       
    686         if (ok_iso_verif) then
    687            call check_isotopes(q,ijb_u,ije_u,'leapfrog 651')
    688         endif !if (ok_iso_verif) then
     670      call check_isotopes(q,ijb_u,ije_u,'leapfrog 651')
    689671     
    690672      call VTb(VTcaldyn)
    … …  
    725707c   -------------------------------------------------------------
    726708
    727         if (ok_iso_verif) then
    728            call check_isotopes(q,ijb_u,ije_u,
     709      call check_isotopes(q,ijb_u,ije_u,
    729710     &           'leapfrog 686: avant caladvtrac')
    730         endif !if (ok_iso_verif) then
    731711     
    732712      IF( forward. OR . leapf )  THEN
    … …  
    743723
    744724         !write(*,*) 'leapfrog 719'
    745          if (ok_iso_verif) then
    746            call check_isotopes(q,ijb_u,ije_u,
     725         call check_isotopes(q,ijb_u,ije_u,
    747726     &           'leapfrog 698: apres caladvtrac')
    748          endif !if (ok_iso_verif) then
    749727
    750728!      do j=1,nqtot
    … …  
    780758
    781759       !write(*,*) 'leapfrog 720'
    782         if (ok_iso_verif) then
    783            call check_isotopes(q,ijb_u,ije_u,'leapfrog 756')
    784         endif !if (ok_iso_verif) then
     760       call check_isotopes(q,ijb_u,ije_u,'leapfrog 756')
    785761
    786762       ! CRisi: pourquoi aller jusqu'à 2 et non pas jusqu'à nqtot??
    … …  
    790766
    791767       !write(*,*) 'leapfrog 724'       
    792         if (ok_iso_verif) then
    793            call check_isotopes(q,ijb_u,ije_u,'leapfrog 762')
    794         endif !if (ok_iso_verif) then
     768       call check_isotopes(q,ijb_u,ije_u,'leapfrog 762')
    795769 
    796770!       CALL FTRACE_REGION_END("integrd")
    … …  
    807781#endif   
    808782
    809         if (ok_iso_verif) then
    810            call check_isotopes(q,ijb_u,ije_u,'leapfrog 775')
    811         endif !if (ok_iso_verif) then
     783      call check_isotopes(q,ijb_u,ije_u,'leapfrog 775')
    812784
    813785c      do j=1,nqtot
    … …  
    11691141       ENDIF ! of IF( apphys )
    11701142
    1171         if (ok_iso_verif) then
    1172            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1132')
    1173         endif !if (ok_iso_verif) then
     1143       call check_isotopes(q,ijb_u,ije_u,'leapfrog 1132')
    11741144        !write(*,*) 'leapfrog 1134: iflag_phys=',iflag_phys
    11751145
    … …  
    12381208
    12391209cc$OMP END PARALLEL
    1240         if (ok_iso_verif) then
    1241            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1196')
    1242         endif !if (ok_iso_verif) then
     1210        call check_isotopes(q,ijb_u,ije_u,'leapfrog 1196')
    12431211
    12441212c-----------------------------------------------------------------------
    … …  
    14751443c              ENDIF
    14761444
    1477         if (ok_iso_verif) then
    1478            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1430')
    1479         endif !if (ok_iso_verif) then     
     1445       call check_isotopes(q,ijb_u,ije_u,'leapfrog 1430')
    14801446 
    14811447c   ********************************************************************
    … …  
    15671533      ENDIF
    15681534     
    1569         if (ok_iso_verif) then
    1570            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1509')
    1571         endif !if (ok_iso_verif) then
     1535      call check_isotopes(q,ijb_u,ije_u,'leapfrog 1509')
    15721536
    15731537      IF ( .NOT.purmats ) THEN
    … …  
    16561620            ENDIF
    16571621
    1658         if (ok_iso_verif) then
    1659            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1584')
    1660         endif !if (ok_iso_verif) then
     1622            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1584')
    16611623
    16621624c-----------------------------------------------------------------------
    … …  
    17011663            ENDIF ! of IF (itau.EQ.itaufin)
    17021664
    1703         if (ok_iso_verif) then
    1704            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1624')
    1705         endif !if (ok_iso_verif) then
     1665            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1624')
    17061666
    17071667c-----------------------------------------------------------------------
    … …  
    17411701
    17421702
    1743         if (ok_iso_verif) then
    1744            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1664')
    1745         endif !if (ok_iso_verif) then
     1703        call check_isotopes(q,ijb_u,ije_u,'leapfrog 1664')
    17461704
    17471705c       ........................................................
    … …  
    17881746
    17891747             
    1790         if (ok_iso_verif) then
    1791            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1698')
    1792         endif !if (ok_iso_verif) then 
     1748              call check_isotopes(q,ijb_u,ije_u,'leapfrog 1698')
    17931749
    17941750              IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
    … …  
    18571813
    18581814
    1859         if (ok_iso_verif) then
    1860            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1750')
    1861         endif !if (ok_iso_verif) then
     1815            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1750')
    18621816
    18631817      END IF ! of IF(.not.purmats)
Note: See TracChangeset for help on using the changeset viewer.