Ignore:
Timestamp:
Jul 23, 2024, 3:29:36 PM (8 weeks ago)
Author:
abarral
Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/isotrac_routines_mod.F90

    r5101 r5103  
    1414
    1515
    16          subroutine uncompress_commun_zone(ncas, cas, &
     16         SUBROUTINE uncompress_commun_zone(ncas, cas, &
    1717         xtp_cas,xtp,xtwater_cas,xtwater,xtevap_cas,xtevap, &
    1818                 ncum,izone)
    … …  
    5151         enddo !do il=1,ncas
    5252
    53          end subroutine uncompress_commun_zone
    54 
    55 
    56          subroutine uncompress_commun_zone_revap(ncas, cas, &
     53         END SUBROUTINE uncompress_commun_zone
     54
     55
     56         SUBROUTINE uncompress_commun_zone_revap(ncas, cas, &
    5757                xtp_cas,xtp,xtwater_cas,xtwater, &
    5858                xtevap_cas,xtevap, &
    … …  
    159159            if (iso_HDO.gt.0) then
    160160                if (xtevap_cas(iso_eau,il).gt.ridicule_evap) then
    161                     call iso_verif_aberrant_choix( &
     161                    CALL iso_verif_aberrant_choix( &
    162162                       xtevap_cas(iso_HDO,il),xtevap_cas(iso_eau,il), &
    163163                       ridicule_trac,deltalimtrac*2,'compress 344a')
    … …  
    165165                ieau=index_trac(izone_revap,iso_eau)
    166166                iHDO=index_trac(izone_revap,iso_HDO)
    167                 call iso_verif_aberrant_choix(xtevap(iHDO,cas(il)), &
     167                CALL iso_verif_aberrant_choix(xtevap(iHDO,cas(il)), &
    168168                    xtevap(ieau,cas(il)),ridicule_trac,deltalimtrac*2, &
    169169                    'compress 344b')
    … …  
    212212!     :           'sortie de uncompress_commun_zone_revap'
    213213         
    214          end subroutine uncompress_commun_zone_revap
    215 
    216 
    217 
    218          subroutine compress_cond_facftmr_zone( &
     214         END SUBROUTINE uncompress_commun_zone_revap
     215
     216
     217
     218         SUBROUTINE compress_cond_facftmr_zone( &
    219219          ncas,  cas, &
    220220          Eqi_prime_cas,Eqi_prime, &
    … …  
    265265            else !if (qp_avantevap_cas(cas(il)).gt.0.0) then
    266266#ifdef ISOVERIF
    267                 call iso_verif_egalite_choix( &
     267                CALL iso_verif_egalite_choix( &
    268268                         (Eqi_prime(cas(il))),0.0, &
    269269                         'compress_stewart 495',errmax,errmaxrel)
    … …  
    278278            else !if (Pqisup(cas(il)).gt.0.0) then
    279279#ifdef ISOVERIF
    280                 call iso_verif_egalite_choix(water(cas(il)),0.0, &
     280                CALL iso_verif_egalite_choix(water(cas(il)),0.0, &
    281281                         'compress_stewart 507',errmax,errmaxrel)
    282282#endif
    … …  
    288288
    289289!#ifdef ISOVERIF
    290 !            call iso_verif_noNaN(water_cas(il),'compress_stewart 518')
     290!            CALL iso_verif_noNaN(water_cas(il),'compress_stewart 518')
    291291!            evap_cas(il)=evap(cas(il))
    292292!     :        *(xtp_avantevap(ieau,cas(il))/qp_avantevap(cas(il)))
    … …  
    301301          enddo
    302302
    303          end subroutine compress_cond_facftmr_zone
    304 
    305          subroutine compress_cond_nofftmr_zone( &
     303         END SUBROUTINE compress_cond_facftmr_zone
     304
     305         SUBROUTINE compress_cond_nofftmr_zone( &
    306306          ncas,  cas, &
    307307          Eqi_prime_cas,Eqi_prime, &
    … …  
    363363            else !if (Pqisup(cas(il)).gt.0.0) then
    364364#ifdef ISOVERIF
    365                 call iso_verif_egalite_choix(water(cas(il)),0.0, &
     365                CALL iso_verif_egalite_choix(water(cas(il)),0.0, &
    366366                         'compress_stewart 654',errmax,errmaxrel)
    367367#endif
    … …  
    391391          enddo
    392392
    393          end subroutine compress_cond_nofftmr_zone
    394 
    395          subroutine compress_noevap_zone( &
     393         END SUBROUTINE compress_cond_nofftmr_zone
     394
     395         SUBROUTINE compress_noevap_zone( &
    396396          ncas,  cas, &
    397397          Pqisup_cas,Pqisup,  &
    … …  
    436436                water_cas(il)=0.0
    437437#ifdef ISOVERIF                 
    438                 call iso_verif_egalite_choix(water(cas(il)), &
     438                CALL iso_verif_egalite_choix(water(cas(il)), &
    439439                     0.0,'compress_stewart 709',errmax,errmaxrel)
    440440#endif               
    … …  
    453453          enddo
    454454
    455          end subroutine compress_noevap_zone
    456 
    457          subroutine compress_evap_liq_zone(iflag_con,ncas,   &
     455         END SUBROUTINE compress_noevap_zone
     456
     457         SUBROUTINE compress_evap_liq_zone(iflag_con,ncas,   &
    458458          cas, &
    459459          Pqisup_cas,Pqisup,  &
    … …  
    521521            else
    522522#ifdef ISOVERIF
    523               call iso_verif_egalite(( &
     523              CALL iso_verif_egalite(( &
    524524                         Eqi_prime(cas(il))),0.0, &
    525525                         'compress_stewart 979')
    … …  
    534534            else !if (Pqisup(cas(il)).gt.0.0) then
    535535#ifdef ISOVERIF
    536                 call iso_verif_egalite_choix(water(cas(il)),0.0, &
     536                CALL iso_verif_egalite_choix(water(cas(il)),0.0, &
    537537                         'compress_stewart 507',errmax,errmaxrel)
    538538#endif
    … …  
    593593          else ! if (Pqisup(cas(il).gt.0.0) then
    594594#ifdef ISOVERIF
    595                call iso_verif_egalite((Pqiinf(cas(il))), &
     595               CALL iso_verif_egalite((Pqiinf(cas(il))), &
    596596                 0.0,'compress_stewart 1047a')
    597                call iso_verif_egalite(( &
     597               CALL iso_verif_egalite(( &
    598598                 Eqi_prime(cas(il))),0.0,'compress_stewart 1047b')
    599                call iso_verif_egalite(( &
     599               CALL iso_verif_egalite(( &
    600600                 Pqiinf_par(cas(il))),0.0,'compress_stewart 1047c')
    601                call iso_verif_egalite((Eqi_par(cas(il))), &
     601               CALL iso_verif_egalite((Eqi_par(cas(il))), &
    602602                 0.0,'compress_stewart 1047d')
    603603#endif             
    … …  
    622622#endif
    623623
    624          end subroutine compress_evap_liq_zone
    625 
    626          subroutine compress_evap_glace_zone(iflag_con, &
     624         END SUBROUTINE compress_evap_liq_zone
     625
     626         SUBROUTINE compress_evap_glace_zone(iflag_con, &
    627627          ncas, cas, &
    628628          water_cas,water,     &
    … …  
    685685            else
    686686#ifdef ISOVERIF
    687                 call iso_verif_egalite(( &
     687                CALL iso_verif_egalite(( &
    688688                         Eqi_prime(cas(il))),0.0, &
    689689                         'compress_stewart 979b')
    … …  
    699699            else !if (Pqisup(cas(il)).gt.0.0) then
    700700#ifdef ISOVERIF
    701                 call iso_verif_egalite_choix(water(cas(il)),0.0, &
     701                CALL iso_verif_egalite_choix(water(cas(il)),0.0, &
    702702                         'compress_stewart 507',errmax,errmaxrel)
    703703#endif
    … …  
    761761
    762762#ifdef ISOVERIF
    763               call iso_verif_egalite((Pqiinf(cas(il))), &
     763              CALL iso_verif_egalite((Pqiinf(cas(il))), &
    764764                 0.0,'compress_stewart 1347a')
    765                call iso_verif_egalite(( &
     765               CALL iso_verif_egalite(( &
    766766                 Eqi_prime(cas(il))),0.0,'compress_stewart 1347b')
    767                call iso_verif_egalite(( &
     767               CALL iso_verif_egalite(( &
    768768                 Pqiinf_par(cas(il))),0.0,'compress_stewart 1347c')
    769                call iso_verif_egalite((Eqi_par(cas(il))), &
     769               CALL iso_verif_egalite((Eqi_par(cas(il))), &
    770770                 0.0,'compress_stewart 1347d')
    771771#endif             
    … …  
    790790#endif             
    791791
    792          end subroutine compress_evap_glace_zone
    793 
    794          subroutine uncompress_ilp_zone( &
     792         END SUBROUTINE compress_evap_glace_zone
     793
     794         SUBROUTINE uncompress_ilp_zone( &
    795795             ncas,cas, &
    796796             zxtrfln_cas,zxt_cas,zxtrfl,zxtrfln,zxt,klon, &
    … …  
    867867              xtrevap_tag(ixt,il)=0.0
    868868#ifdef ISOVERIF
    869               call iso_verif_positif_choix(zxt(ixt,cas(il)), &
     869              CALL iso_verif_positif_choix(zxt(ixt,cas(il)), &
    870870                    0.0,'compress 1508')
    871871#endif             
    … …  
    876876!     :            zxt(iso_eau:ntraciso:3,cas(1))
    877877
    878          end subroutine uncompress_ilp_zone
    879 
    880          subroutine compress_ilp_evap_liq_zone( &
     878         END SUBROUTINE uncompress_ilp_zone
     879
     880         SUBROUTINE compress_ilp_evap_liq_zone( &
    881881             ncas,cas, &
    882882             zxt_cas,zxt, &
    … …  
    947947!     :              il,cas(il),zrfln(cas(il)),zrfl_ancien(cas(il)),
    948948!     :              zqev_diag(cas(il))   
    949               call iso_verif_egalite(zqev_diag(cas(il)),0.0, &
     949              CALL iso_verif_egalite(zqev_diag(cas(il)),0.0, &
    950950                 'compress_stewart 1591a')
    951               call iso_verif_egalite(zrfln(cas(il)),0.0, &
     951              CALL iso_verif_egalite(zrfln(cas(il)),0.0, &
    952952                 'compress_stewart 1591b')
    953953#endif               
    … …  
    961961         enddo !do il=1,ncas
    962962
    963          end subroutine compress_ilp_evap_liq_zone   
    964 
    965 
    966          subroutine compress_ilp_evap_glace_zone( &
     963         END SUBROUTINE  compress_ilp_evap_liq_zone
     964
     965
     966         SUBROUTINE compress_ilp_evap_glace_zone( &
    967967             ncas,cas, &
    968968             zxt_cas,zxt, &
    … …  
    10221022          else !if (zrfl_ancien(cas(il)).gt.0.0) then
    10231023#ifdef ISOVERIF 
    1024               call iso_verif_egalite(zqev_diag(cas(il)),0.0, &
     1024              CALL iso_verif_egalite(zqev_diag(cas(il)),0.0, &
    10251025                 'compress_stewart 1791a')
    1026               call iso_verif_egalite(zrfln(cas(il)),0.0, &
     1026              CALL iso_verif_egalite(zrfln(cas(il)),0.0, &
    10271027                 'compress_stewart 1791b')
    10281028#endif               
    … …  
    10321032         enddo
    10331033
    1034          end subroutine compress_ilp_evap_glace_zone   
    1035 
    1036 
    1037       subroutine ajoute_revap(ncas,cas, &
     1034         END SUBROUTINE  compress_ilp_evap_glace_zone
     1035
     1036
     1037      SUBROUTINE ajoute_revap(ncas,cas, &
    10381038          klon,izone,zxt,xtrevap_tag)
    10391039
    … …  
    10661066      do i=1,ncas
    10671067        do ixt=1,ntraciso
    1068             call iso_verif_positif_choix(zxt(ixt,cas(i)),0.0, &
     1068            CALL iso_verif_positif_choix(zxt(ixt,cas(i)),0.0, &
    10691069                         'ajoute_revap 29')
    10701070        enddo
    … …  
    10951095      do i=1,ncas
    10961096        do ixt=1,ntraciso
    1097             call iso_verif_positif_choix(zxt(ixt,cas(i)),0.0, &
     1097            CALL iso_verif_positif_choix(zxt(ixt,cas(i)),0.0, &
    10981098                         'ajoute_revap 40')
    10991099        enddo
    … …  
    11021102
    11031103      return
    1104       end subroutine ajoute_revap
     1104      END SUBROUTINE ajoute_revap
    11051105
    11061106
    … …  
    11281128      !logical is_in_triangle
    11291129     
    1130       is_in_bassin=.false.
     1130      is_in_bassin=.FALSE.
    11311131#ifdef ISOVERIF     
    11321132      write(*,*) 'is_in_basin 84: entree,bassin=',bassin
    … …  
    11381138          if (is_in_rectangle(lon,lat,-67.0,28.0,20.0,-45.0)) then
    11391139            ! boite sud
    1140             is_in_bassin=.true.
     1140            is_in_bassin=.TRUE.
    11411141            return
    11421142          endif
    11431143          if (is_in_rectangle(lon,lat,-100.0,40.0,-5.3,28.0)) then
    11441144            ! ouest gibraltar   
    1145             is_in_bassin=.true.
     1145            is_in_bassin=.TRUE.
    11461146            return
    11471147          endif
    11481148          if (is_in_rectangle(lon,lat,-100.0,48.0,0.0,40.0)) then
    11491149            ! Ouest France
    1150             is_in_bassin=.true.
     1150            is_in_bassin=.TRUE.
    11511151            return
    11521152          endif
    11531153          if (is_in_rectangle(lon,lat,-90.0,80.0,10.0,46.0)) then
    11541154            ! Atlantic Nord
    1155             is_in_bassin=.true.
     1155            is_in_bassin=.TRUE.
    11561156            return
    11571157          endif
    … …  
    11591159                 -62.0,0.0,-62.0,30.0,-112.0,30.0)) then
    11601160            ! golfe du Mexique
    1161             is_in_bassin=.true.
     1161            is_in_bassin=.TRUE.
    11621162            return
    11631163          endif
    … …  
    11681168#endif         
    11691169          if (is_in_rectangle(lon,lat,0.0,48.0,45.0,29.0)) then
    1170             is_in_bassin=.true.
     1170            is_in_bassin=.TRUE.
    11711171            return
    11721172          endif
    11731173          if (is_in_rectangle(lon,lat,-5.3,42.0,45.0,29.0)) then
    1174             is_in_bassin=.true.
     1174            is_in_bassin=.TRUE.
    11751175            return
    11761176          endif
    … …  
    11811181#endif           
    11821182          if (is_in_rectangle(lon,lat,20.0,30.0,110.0,-45.0)) then
    1183             is_in_bassin=.true.
     1183            is_in_bassin=.TRUE.
    11841184            return
    11851185          endif
    … …  
    11871187                 90.0,30.0,90.0,-45.0,150.0,-45.0)) then
    11881188            ! Ouest Australie
    1189             is_in_bassin=.true.
     1189            is_in_bassin=.TRUE.
    11901190            return
    11911191          endif   
    … …  
    11961196#endif           
    11971197          if (is_in_rectangle(lon,lat,20.0,0.0,120.0,-45.0)) then
    1198             is_in_bassin=.true.
     1198            is_in_bassin=.TRUE.
    11991199            return
    12001200          endif
    … …  
    12051205#endif           
    12061206          if (is_in_rectangle(lon,lat,20.0,30.0,76.0,0.0)) then
    1207             is_in_bassin=.true.
     1207            is_in_bassin=.TRUE.
    12081208            return
    12091209          endif
    … …  
    12141214#endif           
    12151215          if (is_in_rectangle(lon,lat,76.0,30.0,110.0,0.0)) then
    1216             is_in_bassin=.true.
     1216            is_in_bassin=.TRUE.
    12171217            return
    12181218          endif         
    … …  
    12241224         if (is_in_rectangle(lon,lat,-180.0,80.0,-100.0,-45.0)) then
    12251225             ! pacifique Est
    1226             is_in_bassin=.true.
     1226            is_in_bassin=.TRUE.
    12271227            return
    12281228          endif
    12291229          if (is_in_rectangle(lon,lat,110.0,80.0,180.0,28.0)) then
    12301230            ! Pacifique Nord Ouest
    1231             is_in_bassin=.true.
     1231            is_in_bassin=.TRUE.
    12321232            return
    12331233          endif
    12341234          if (is_in_rectangle(lon,lat,120.0,80.0,180.0,-45.0)) then
    12351235            ! Pacifique central Sud
    1236             is_in_bassin=.true.
     1236            is_in_bassin=.TRUE.
    12371237            return
    12381238          endif
    … …  
    12401240                  90.0,28.0,150.0,-45.0,150.0,28.0)) then
    12411241            ! Pacifque Sud Ouest
    1242             is_in_bassin=.true.
     1242            is_in_bassin=.TRUE.
    12431243            return
    12441244          endif
    … …  
    12461246                  -62.0,0.0,-112.0,30.0,-112.0,0.0)) then
    12471247            ! Ouest Amérique centrale
    1248             is_in_bassin=.true.
     1248            is_in_bassin=.TRUE.
    12491249            return
    12501250          endif
    12511251          if (is_in_rectangle(lon,lat,-180.0,0.0,-67.0,-45.0)) then
    12521252            ! Ouest Chili
    1253             is_in_bassin=.true.
     1253            is_in_bassin=.TRUE.
    12541254            return
    12551255          endif
    … …  
    12601260#endif           
    12611261          if (lat.lt.-45.0+0.2) then
    1262                 is_in_bassin=.true.
     1262                is_in_bassin=.TRUE.
    12631263                return   
    12641264          endif 
    … …  
    12691269#endif           
    12701270          if (abs(lat).gt.35.0) then
    1271                 is_in_bassin=.true.
     1271                is_in_bassin=.TRUE.
    12721272                return   
    12731273          endif
    … …  
    12781278#endif           
    12791279          if (abs(lat).lt.15.0) then
    1280                 is_in_bassin=.true.
     1280                is_in_bassin=.TRUE.
    12811281                return   
    12821282          endif
    … …  
    12881288          if ((abs(lat).ge.15.0).and. &
    12891289                 (abs(lat).le.35.0)) then
    1290                 is_in_bassin=.true.
     1290                is_in_bassin=.TRUE.
    12911291                return   
    12921292          endif   
    … …  
    13141314      end function is_in_bassin
    13151315
    1316       subroutine find_bassin(lat,lon,bassin)
     1316      SUBROUTINE find_bassin(lat,lon,bassin)
    13171317      use isotrac_mod, only: izone_poubelle,ntraceurs_zone=>ntiso,option_traceurs, &
    13181318&        bassin_map
    … …  
    13301330      !logical is_in_bassin
    13311331
    1332       continu=.true.
     1332      continu=.TRUE.
    13331333      bassin=1
    13341334     
    … …  
    13421342!#endif       
    13431343         if (is_in_bassin(lat,lon,bassin)) then
    1344                 continu=.false.
     1344                continu=.FALSE.
    13451345#ifdef ISOVERIF                 
    13461346                write(*,*) 'find_bassin 173: trouve: bassin=',bassin
    … …  
    13531353         endif
    13541354         if (bassin.eq.izone_poubelle) then
    1355                 continu=.false.
     1355                continu=.FALSE.
    13561356                bassin=izone_poubelle
    13571357!#ifdef ISOVERIF                 
    … …  
    13641364       ! donc bassin<=ntraceurs_zone-1
    13651365#ifdef ISOVERIF
    1366        call iso_verif_positif(float(ntraceurs_zone-1-bassin), &
     1366       CALL iso_verif_positif(float(ntraceurs_zone-1-bassin), &
    13671367                 'find_bassin 195')
    13681368#endif
    13691369
    13701370        return
    1371         end subroutine find_bassin
    1372 
    1373         subroutine initialise_bassins_boites(presnivs)
     1371        END SUBROUTINE find_bassin
     1372
     1373        SUBROUTINE initialise_bassins_boites(presnivs)
    13741374        USE dimphy, only: klev
    13751375        USE geometry_mod, ONLY: longitude_deg, latitude_deg
    … …  
    13831383       if (option_traceurs.eq.3) then
    13841384           ! initialisation de bassin_map
    1385            call bassin_map_init(latitude_deg,longitude_deg,bassin_map)
     1385           CALL bassin_map_init(latitude_deg,longitude_deg,bassin_map)
    13861386       else if (option_traceurs.eq.20) then
    13871387           ! initialisation de bassin_map selon < ou > 35° lat
    13881388           write(*,*) 'physiq 1681: init de la map pour tag 20'
    1389            call bassin_map_init_opt20(latitude_deg,bassin_map)
     1389           CALL bassin_map_init_opt20(latitude_deg,bassin_map)
    13901390       else if (option_traceurs.eq.5) then
    1391            call boite_AMMA_init(latitude_deg,longitude_deg,presnivs,boite_map)
     1391           CALL boite_AMMA_init(latitude_deg,longitude_deg,presnivs,boite_map)
    13921392       else if (option_traceurs.eq.21) then
    1393            call boite_UT_extra_init(latitude_deg,longitude_deg,presnivs,boite_map)
     1393           CALL boite_UT_extra_init(latitude_deg,longitude_deg,presnivs,boite_map)
    13941394       endif
    13951395
    13961396        return
    1397         end subroutine initialise_bassins_boites
    1398 
    1399         subroutine bassin_map_init(lat,lon,bassin_map)
     1397        END SUBROUTINE initialise_bassins_boites
     1398
     1399        SUBROUTINE bassin_map_init(lat,lon,bassin_map)
    14001400        USE dimphy, only: klon
    14011401#ifdef ISOVERIF
    … …  
    14131413
    14141414        do i=1,klon
    1415              call find_bassin(lat(i),lon(i),bassin_map(i))
     1415             CALL find_bassin(lat(i),lon(i),bassin_map(i))
    14161416#ifdef ISOVERIF             
    14171417             write(*,*) 'init 233: i,lat,lon,bassin=',i,lat(i),lon(i), &
    … …  
    14211421
    14221422        return
    1423         end subroutine bassin_map_init
     1423        END SUBROUTINE bassin_map_init
    14241424
    14251425        function is_in_rectangle(x,y,x1,y1,x2,y2)
    … …  
    14401440        if ((x-x2.lt.0.1).and.(x-x1.gt.-0.1).and. &
    14411441                  (y-y1.lt.0.1).and.(y-y2.gt.-0.1)) then
    1442           is_in_rectangle=.true.
     1442          is_in_rectangle=.TRUE.
    14431443        else
    1444           is_in_rectangle=.false.
     1444          is_in_rectangle=.FALSE.
    14451445        endif
    14461446!#ifdef ISOVERIF       
    … …  
    14741474        det2=(x2-x)*(y3-y)-(y2-y)*(x3-x)
    14751475        det3=(x3-x)*(y1-y)-(y3-y)*(x1-x)
    1476         is_in_triangle=.false.
     1476        is_in_triangle=.FALSE.
    14771477        if ((det1*det2.gt.0.0).and.(det2*det3.gt.0.0)) then
    1478           is_in_triangle=.true.
     1478          is_in_triangle=.TRUE.
    14791479        else
    1480           is_in_triangle=.false.     
     1480          is_in_triangle=.FALSE.
    14811481        endif
    14821482!#ifdef ISOVERIF       
    … …  
    14881488
    14891489
    1490         subroutine isotrac_recolorise_tmin(xt,t)
     1490        SUBROUTINE isotrac_recolorise_tmin(xt,t)
    14911491        USE dimphy, only: klon, klev
    14921492        USE isotrac_mod, only: zone_temp,nzone_temp
    … …  
    15641564        do k=1,klev
    15651565          do i=1,klon
    1566             call iso_verif_traceur(xt(1,i,k),'recolorise 403')
     1566            CALL iso_verif_traceur(xt(1,i,k),'recolorise 403')
    15671567          enddo !do i=1,klon
    15681568        enddo !do k=1,klev
    … …  
    15701570
    15711571        return
    1572         end subroutine isotrac_recolorise_tmin
    1573 
    1574         subroutine isotrac_recolorise_tmin_sfrev(xt,t)
     1572        END SUBROUTINE isotrac_recolorise_tmin
     1573
     1574        SUBROUTINE isotrac_recolorise_tmin_sfrev(xt,t)
    15751575        USE dimphy, only: klon,klev
    15761576        USE isotrac_mod, only: nzone_temp,zone_temp
    … …  
    16191619        do k=1,klev
    16201620          do i=1,klon
    1621             call iso_verif_traceur(xt(1,i,k),'recolorise 594')
     1621            CALL iso_verif_traceur(xt(1,i,k),'recolorise 594')
    16221622          enddo !do i=1,klon
    16231623        enddo !do k=1,klev
    … …  
    16251625
    16261626        return
    1627         end subroutine isotrac_recolorise_tmin_sfrev
    1628 
    1629 
    1630         subroutine isotrac_recolorise_saturation(xt,rh,lat,pres)
     1627        END SUBROUTINE isotrac_recolorise_tmin_sfrev
     1628
     1629
     1630        SUBROUTINE isotrac_recolorise_saturation(xt,rh,lat,pres)
    16311631        USE dimphy, only: klon,klev
    16321632#ifdef ISOVERIF
    … …  
    16561656        do k=1,klev
    16571657          do i=1,klon
    1658             call iso_verif_traceur(xt(1,i,k),'recolorise 612')
     1658            CALL iso_verif_traceur(xt(1,i,k),'recolorise 612')
    16591659          enddo !do i=1,klon
    16601660        enddo !do k=1,klev       
    … …  
    16841684        do k=1,klev
    16851685          do i=1,klon
    1686             call iso_verif_traceur(xt(1,i,k),'recolorise 637')
     1686            CALL iso_verif_traceur(xt(1,i,k),'recolorise 637')
    16871687          enddo !do i=1,klon
    16881688        enddo !do k=1,klev
    … …  
    16901690
    16911691        return
    1692         end subroutine isotrac_recolorise_saturation
    1693 
    1694         subroutine isotrac_recolorise_boite(xt,boite_map)
     1692        END SUBROUTINE isotrac_recolorise_saturation
     1693
     1694        SUBROUTINE isotrac_recolorise_boite(xt,boite_map)
    16951695        USE dimphy, only: klon,klev
    16961696#ifdef ISOVERIF
    … …  
    16991699        implicit none
    17001700
    1701         ! subroutine écrite à la base pour tagguer 3 boites AMMA.
     1701        ! SUBROUTINE écrite à la base pour tagguer 3 boites AMMA.
    17021702        ! Mais ça peut être générique, selon comment est initialisée boite_map
    17031703
    … …  
    17361736        do k=1,klev
    17371737          do i=1,klon
    1738             call iso_verif_traceur(xt(1,i,k),'recolorise 514')
     1738            CALL iso_verif_traceur(xt(1,i,k),'recolorise 514')
    17391739          enddo !do i=1,klon
    17401740        enddo !do k=1,klev
    … …  
    17421742
    17431743        return
    1744         end subroutine isotrac_recolorise_boite
    1745 
    1746         subroutine isotrac_recolorise_extra(xt,rlat)
     1744        END SUBROUTINE isotrac_recolorise_boite
     1745
     1746        SUBROUTINE isotrac_recolorise_extra(xt,rlat)
    17471747        USE dimphy, only: klon,klev
    17481748        usE isotrac_mod, only: lim_tag20,izone_trop,izone_extra
    … …  
    17521752        implicit none
    17531753
    1754         ! subroutine écrite pour l'option de taggage 20
     1754        ! SUBROUTINE écrite pour l'option de taggage 20
    17551755        ! permet de retagguer la vapeur tropicale en vapeur
    17561756        ! extratropicale dès qu'elle atteint 35° de latitude
    … …  
    17831783        do k=1,klev
    17841784          do i=1,klon
    1785             call iso_verif_traceur(xt(1,i,k),'recolorise 741')
     1785            CALL iso_verif_traceur(xt(1,i,k),'recolorise 741')
    17861786          enddo !do i=1,klon
    17871787        enddo !do k=1,klev
    … …  
    17891789
    17901790        return
    1791         end subroutine isotrac_recolorise_extra
    1792 
    1793         subroutine isotrac_recolorise_conv(xt,rlat,presnivs,rain_con)
     1791        END SUBROUTINE isotrac_recolorise_extra
     1792
     1793        SUBROUTINE isotrac_recolorise_conv(xt,rlat,presnivs,rain_con)
    17941794        USE dimphy, only: klon,klev
    17951795        use isotrac_mod, only: lim_precip_tag22, &
    … …  
    18001800        implicit none
    18011801
    1802         ! subroutine écrite pour l'option de taggage 20
     1802        ! SUBROUTINE écrite pour l'option de taggage 20
    18031803        ! permet de retagguer la vapeur tropicale en vapeur
    18041804        ! extratropicale dès qu'elle atteint 35° de latitude
    … …  
    18521852        do k=1,klev
    18531853          do i=1,klon
    1854             call iso_verif_traceur(xt(1,i,k),'recolorise 741')
     1854            CALL iso_verif_traceur(xt(1,i,k),'recolorise 741')
    18551855          enddo !do i=1,klon
    18561856        enddo !do k=1,klev
    … …  
    18581858
    18591859        return
    1860         end subroutine isotrac_recolorise_conv
    1861 
    1862 
    1863         subroutine boite_AMMA_init(lat,lon,presnivs,boite_map)
     1860        END SUBROUTINE isotrac_recolorise_conv
     1861
     1862
     1863        SUBROUTINE boite_AMMA_init(lat,lon,presnivs,boite_map)
    18641864        USE dimphy, only: klon,klev
    18651865#ifdef ISOVERIF
    … …  
    19111911
    19121912        return
    1913         end subroutine boite_AMMA_init
     1913        END SUBROUTINE boite_AMMA_init
    19141914
    19151915       
    1916         subroutine boite_UT_extra_init(lat,lon,presnivs,boite_map)
     1916        SUBROUTINE boite_UT_extra_init(lat,lon,presnivs,boite_map)
    19171917        USE dimphy, only: klon,klev
    19181918        use isotrac_mod, only: izone_extra,izone_trop
    … …  
    19491949
    19501950        return
    1951         end subroutine boite_UT_extra_init
     1951        END SUBROUTINE boite_UT_extra_init
    19521952
    19531953
    … …  
    20082008                find_index=nzone_pres
    20092009          else !if (t(i,k).ge.zone_temp1) then
    2010                 continu=.true.
     2010                continu=.TRUE.
    20112011                find_index=2
    20122012                do while (continu)
    20132013                  if (pres.ge.zone_pres(find_index)) then
    2014                      continu=.false.
     2014                     continu=.FALSE.
    20152015                     ! c'est izone_temp, zone trouvée
    20162016                  else   
    … …  
    20482048        end function index_zone_latpres
    20492049
    2050         subroutine iso_recolorise_condensation(qt,cond, &
     2050        SUBROUTINE iso_recolorise_condensation(qt,cond, &
    20512051                 xt,zxtcond,tcond,ep,xtres, &
    20522052                 seuil_in)
    … …  
    20912091#ifdef ISOVERIF
    20922092           do ixt=1,ntraciso
    2093               call iso_verif_positif(xt(ixt)-zxtcond(ixt), &
     2093              CALL iso_verif_positif(xt(ixt)-zxtcond(ixt), &
    20942094                 'iso_trac 898')
    20952095           enddo  !do ixt=1,ntraciso 
    2096            call iso_verif_traceur_justmass(xt, &
     2096           CALL iso_verif_traceur_justmass(xt, &
    20972097                'iso_trac_routines 906')
    2098            call iso_verif_traceur_justmass(zxtcond, &
     2098           CALL iso_verif_traceur_justmass(zxtcond, &
    20992099                'iso_trac_routines 908')
    21002100#endif           
    21012101          ! bidouille
    21022102          if (bidouille_anti_divergence) then
    2103               call iso_verif_traceur_jbidouille(xt)
    2104               call iso_verif_traceur_jbidouille(zxtcond) 
     2103              CALL iso_verif_traceur_jbidouille(xt)
     2104              CALL iso_verif_traceur_jbidouille(zxtcond)
    21052105          endif
    21062106
    … …  
    21152115#ifdef ISOVERIF
    21162116        do ixt=1,ntraciso
    2117           call iso_verif_positif(xtres(ixt), &
     2117          CALL iso_verif_positif(xtres(ixt), &
    21182118                    'iso_trac_routines 921')
    21192119        enddo
    … …  
    21862186#ifdef ISOVERIF
    21872187        do ixt=1,ntraciso
    2188           call iso_verif_positif(xtres(ixt), &
     2188          CALL iso_verif_positif(xtres(ixt), &
    21892189                    'iso_trac_routines 940')
    21902190        enddo
    … …  
    22372237#ifdef ISOVERIF             
    22382238        do ixt=1,ntraciso
    2239           call iso_verif_positif(xtres(ixt),'iso_trac_routines 953')
     2239          CALL iso_verif_positif(xtres(ixt),'iso_trac_routines 953')
    22402240        enddo
    22412241        if (nzone_temp.ge.5) then
    … …  
    22552255
    22562256            return
    2257             end subroutine iso_recolorise_condensation
    2258 
    2259         subroutine bassin_map_init_opt20(lat,bassin_map)
     2257            END SUBROUTINE iso_recolorise_condensation
     2258
     2259        SUBROUTINE bassin_map_init_opt20(lat,bassin_map)
    22602260        USE dimphy, only: klon
    22612261        use isotrac_mod, only: izone_cont,izone_trop,lim_tag20
    … …  
    22822282       
    22832283        return
    2284         end subroutine bassin_map_init_opt20
    2285 
    2286         subroutine isotrac_recolorise_general(xt_seri,t_seri,zx_rh,presnivs)
     2284        END SUBROUTINE bassin_map_init_opt20
     2285
     2286        SUBROUTINE isotrac_recolorise_general(xt_seri,t_seri,zx_rh,presnivs)
    22872287        USE geometry_mod, ONLY: latitude_deg
    22882288        USE dimphy, only: klon,klev
    … …  
    22972297
    22982298              if (option_traceurs.eq.4) then
    2299          call isotrac_recolorise_tmin(xt_seri,t_seri)
     2299         CALL isotrac_recolorise_tmin(xt_seri,t_seri)
    23002300      elseif ((option_traceurs.eq.5).or. &
    23012301                  (option_traceurs.eq.21)) then
    2302         call isotrac_recolorise_boite(xt_seri,boite_map)
     2302        CALL isotrac_recolorise_boite(xt_seri,boite_map)
    23032303      elseif (option_traceurs.eq.13) then
    2304         call isotrac_recolorise_tmin_sfrev(xt_seri,t_seri)
     2304        CALL isotrac_recolorise_tmin_sfrev(xt_seri,t_seri)
    23052305      elseif (option_traceurs.eq.14) then
    2306         call isotrac_recolorise_saturation(xt_seri,zx_rh,latitude_deg,presnivs)
     2306        CALL isotrac_recolorise_saturation(xt_seri,zx_rh,latitude_deg,presnivs)
    23072307      elseif (option_traceurs.eq.20) then     
    2308         call isotrac_recolorise_extra(xt_seri,latitude_deg)
     2308        CALL isotrac_recolorise_extra(xt_seri,latitude_deg)
    23092309      endif !if (option_traceurs.eq.4) then
    23102310
    23112311        return
    2312         end subroutine isotrac_recolorise_general
     2312        END SUBROUTINE isotrac_recolorise_general
    23132313
    23142314
    23152315       
    23162316
    2317         subroutine iso_verif_traceur_jbid_vect(x,n,m)
     2317        SUBROUTINE iso_verif_traceur_jbid_vect(x,n,m)
    23182318        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    23192319        !use isotrac_mod, only: ntraceurs_zone=>nzone
    … …  
    23782378        endif !if (bidouille_anti_divergence) then   
    23792379
    2380         end subroutine iso_verif_traceur_jbid_vect
    2381 
    2382         subroutine iso_verif_traceur_jbidouille(x)
     2380        END SUBROUTINE iso_verif_traceur_jbid_vect
     2381
     2382        SUBROUTINE iso_verif_traceur_jbidouille(x)
    23832383        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    23842384        implicit none
    … …  
    24162416        endif !if (bidouille_anti_divergence) then     
    24172417
    2418         end subroutine iso_verif_traceur_jbidouille
    2419 
    2420 
    2421         subroutine iso_verif_traceur_jbid_pos(x)
     2418        END SUBROUTINE iso_verif_traceur_jbidouille
     2419
     2420
     2421        SUBROUTINE iso_verif_traceur_jbid_pos(x)
    24222422        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    24232423!#ifdef ISOVERIF
    … …  
    24842484        enddo !do iiso=1,ntraceurs_iso
    24852485#ifdef ISOVERIF       
    2486         call iso_verif_traceur_pbidouille(x,'iso_verif_trac 558')
     2486        CALL iso_verif_traceur_pbidouille(x,'iso_verif_trac 558')
    24872487#else
    2488         call iso_verif_traceur_jbidouille(x)
     2488        CALL iso_verif_traceur_jbidouille(x)
    24892489#endif
    24902490        endif !if (bidouille_anti_divergence) then     
    24912491
    2492         end subroutine iso_verif_traceur_jbid_pos
    2493 
    2494         subroutine iso_verif_traceur_jbid_pos_vect(n,m,x)
     2492        END SUBROUTINE iso_verif_traceur_jbid_pos
     2493
     2494        SUBROUTINE iso_verif_traceur_jbid_pos_vect(n,m,x)
    24952495        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    24962496#ifdef ISOVERIF
    … …  
    25642564        enddo !do iiso=1,ntraceurs_iso
    25652565#ifdef ISOVERIF       
    2566         call iso_verif_traceur_pbid_vect(x,n,m,'iso_verif_trac 558')
     2566        CALL iso_verif_traceur_pbid_vect(x,n,m,'iso_verif_trac 558')
    25672567#else
    2568         call iso_verif_traceur_jbid_vect(x,n,m)
     2568        CALL iso_verif_traceur_jbid_vect(x,n,m)
    25692569#endif
    25702570        endif !if (bidouille_anti_divergence) then     
    25712571
    2572         end subroutine iso_verif_traceur_jbid_pos_vect
    2573 
    2574         subroutine iso_verif_traceur_jbid_pos2(x,q)
     2572        END SUBROUTINE iso_verif_traceur_jbid_pos_vect
     2573
     2574        SUBROUTINE iso_verif_traceur_jbid_pos2(x,q)
    25752575        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    25762576#ifdef ISOVERIF
    … …  
    26312631              endif !if (dqtmp.gt.0.0) then   
    26322632#ifdef ISOVERIF
    2633               call iso_verif_traceur(x,'iso_verif_traceurs 612')
     2633              CALL iso_verif_traceur(x,'iso_verif_traceurs 612')
    26342634#endif
    26352635          endif !if (local_q(i,k).lt.0.0) then
    26362636
    26372637#ifdef ISOVERIF
    2638         call iso_verif_traceur_pbidouille(x,'iso_verif_trac 625')   
     2638        CALL iso_verif_traceur_pbidouille(x,'iso_verif_trac 625')
    26392639#endif
    26402640        endif ! if (bidouille_anti_divergence) then 
    26412641
    2642         end subroutine iso_verif_traceur_jbid_pos2
    2643 
    2644         subroutine iso_verif_traceur_jbid_vect1D(x,n)
     2642        END SUBROUTINE iso_verif_traceur_jbid_pos2
     2643
     2644        SUBROUTINE iso_verif_traceur_jbid_vect1D(x,n)
    26452645        USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule
    26462646        implicit none
    … …  
    26802680        endif !if (bidouille_anti_divergence) then             
    26812681
    2682         end subroutine iso_verif_traceur_jbid_vect1D
     2682        END SUBROUTINE iso_verif_traceur_jbid_vect1D
    26832683     
    26842684! on met ces routines ici pour éviter dépendances circulaires   
    26852685#ifdef ISOVERIF
    26862686
    2687         subroutine iso_verif_traceur_pbidouille(x,err_msg)
     2687        SUBROUTINE iso_verif_traceur_pbidouille(x,err_msg)
    26882688        use isotopes_verif_mod
    26892689        implicit none
    … …  
    27082708        endif
    27092709
    2710         end subroutine iso_verif_traceur_pbidouille
     2710        END SUBROUTINE iso_verif_traceur_pbidouille
    27112711
    27122712        function iso_verif_traceur_pbid_ns(x,err_msg)
    … …  
    27582758            ! on réajuste pour que les traceurs fasses bien la somme
    27592759            ! des traceurs
    2760             call iso_verif_traceur_jbidouille(x)
     2760            CALL iso_verif_traceur_jbidouille(x)
    27612761        endif !if (bidouille_anti_divergence) then   
    27622762
    … …  
    27722772        end function iso_verif_traceur_pbid_ns
    27732773
    2774         subroutine iso_verif_traceur_pbid_vect(x,n,m,err_msg)
     2774        SUBROUTINE iso_verif_traceur_pbid_vect(x,n,m,err_msg)
    27752775        use isotopes_mod, ONLY: iso_HDO,bidouille_anti_divergence
    27762776        use isotopes_verif_mod
    … …  
    27872787       
    27882788        ! verif noNaN
    2789         call iso_verif_traceur_noNaN_vect(x,n,m,err_msg)
     2789        CALL iso_verif_traceur_noNaN_vect(x,n,m,err_msg)
    27902790
    27912791        ! verif masse
    2792         call iso_verif_trac_masse_vect(x,n,m,err_msg,errmax*10, &
     2792        CALL iso_verif_trac_masse_vect(x,n,m,err_msg,errmax*10, &
    27932793                 errmaxrel*50)
    27942794
    … …  
    27962796            ! on réajuste pour que les traceurs fasses bien la somme
    27972797            ! des traceurs
    2798             call iso_verif_traceur_jbid_vect(x,n,m)
     2798            CALL iso_verif_traceur_jbid_vect(x,n,m)
    27992799        endif !if (bidouille_anti_divergence) then   
    28002800
    28012801       ! verif deltaD
    28022802       if (iso_HDO.gt.0) then
    2803         call iso_verif_tracdd_vect(x,n,m,err_msg) 
     2803        CALL iso_verif_tracdd_vect(x,n,m,err_msg)
    28042804       endif
    28052805
    2806         end subroutine iso_verif_traceur_pbid_vect
     2806        END SUBROUTINE iso_verif_traceur_pbid_vect
    28072807#endif
    28082808
Note: See TracChangeset for help on using the changeset viewer.