Ignore:
Timestamp:
Aug 2, 2024, 2:12:03 PM (7 weeks ago)
Author:
abarral
Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bands.F90

    r5134 r5158  
    6262     
    6363      IF (ierr==0) THEN
    64          do i=0,mpi_size-1
     64         DO i=0,mpi_size-1
    6565          read (unit_number,*) j,jj_nb_caldyn(i)
    6666        enddo
    6767     
    68         do i=0,mpi_size-1
     68        DO i=0,mpi_size-1
    6969          read (unit_number,*) j,jj_nb_vanleer(i)
    7070        enddo
    7171     
    72         do i=0,mpi_size-1
     72        DO i=0,mpi_size-1
    7373          read (unit_number,*) j,jj_nb_dissip(i)
    7474        enddo
    7575     
    76         do i=0,mpi_size-1
     76        DO i=0,mpi_size-1
    7777          read (unit_number,*) j,distrib_phys(i)
    7878        enddo
    79        
    80         CLOSE(unit_number) 
     79
     80    CLOSE(unit_number)
    8181 
    8282      else
    83         do i=0,mpi_size-1
     83        DO i=0,mpi_size-1
    8484          jj_nb_caldyn(i)=(jjm+1)/mpi_size
    85           IF (i<MOD(jjm+1,mpi_size)) jj_nb_caldyn(i)=jj_nb_caldyn(i)+1
     85      IF (i<MOD(jjm+1,mpi_size)) jj_nb_caldyn(i)=jj_nb_caldyn(i)+1
    8686        enddo
    8787     
    … …  
    8989        jj_nb_dissip(:)=jj_nb_caldyn(:)
    9090       
    91         do i=0,mpi_size-1
    92           distrib_phys(i)=(iim*(jjm-1)+2)/mpi_size
    93           IF (i<MOD(iim*(jjm-1)+2,mpi_size)) distrib_phys(i)=distrib_phys(i)+1
    94         enddo
     91    DO i=0,mpi_size-1
     92      distrib_phys(i)=(iim*(jjm-1)+2)/mpi_size
     93      IF (i<MOD(iim*(jjm-1)+2,mpi_size)) distrib_phys(i)=distrib_phys(i)+1
     94    enddo
    9595      endif
    9696     
    … …  
    110110     INTEGER :: jj_para_end(0:mpi_size-1)
    111111       
    112       do i=0,mpi_size-1
     112      DO i=0,mpi_size-1
    113113         jj_nb_vanleer2(i)=(jjm+1)/mpi_size
    114         IF (i<MOD(jjm+1,mpi_size)) jj_nb_vanleer2(i)=jj_nb_vanleer2(i)+1
     114    IF (i<MOD(jjm+1,mpi_size)) jj_nb_vanleer2(i)=jj_nb_vanleer2(i)+1
    115115      enddo
    116116         
    … …  
    126126      ENDDO
    127127 
    128        do i=0,MPI_Size-1
     128       DO i=0,MPI_Size-1
    129129        jj_Nb_physic(i)=jj_para_end(i)-jj_para_begin(i)+1
    130130        IF (i/=0) THEN
    … …  
    135135      enddo
    136136     
    137       do i=0,MPI_Size-1
     137      DO i=0,MPI_Size-1
    138138        jj_Nb_physic_bis(i)=jj_para_end(i)-jj_para_begin(i)+1
    139139        IF (i/=0) THEN
    … …  
    141141            jj_Nb_physic_bis(i)=jj_Nb_physic_bis(i)-1
    142142          else
    143             jj_Nb_physic_bis(i-1)=jj_Nb_physic_bis(i-1)+1
    144             jj_Nb_physic_bis(i)=jj_Nb_physic_bis(i)-1
    145           ENDIF
     143        jj_Nb_physic_bis(i-1)=jj_Nb_physic_bis(i-1)+1
     144        jj_Nb_physic_bis(i)=jj_Nb_physic_bis(i)-1
     145      ENDIF
    146146        endif
    147147      enddo
    … …  
    193193      CALL allgather_timer_average
    194194
    195       do i=0,mpi_size-1
     195      DO i=0,mpi_size-1
    196196        value(i)=timer_average(jj_nb_caldyn(i),timer_caldyn,i)
    197         index(i)=i
    198       enddo
    199      
    200       do i=0,mpi_size-2
    201         do j=i+1,mpi_size-1
    202           IF (value(i)>value(j)) THEN
    203             tmpvalue=value(i)
    204             value(i)=value(j)
    205             value(j)=tmpvalue
    206            
    207             tmpindex=index(i)
    208             index(i)=index(j)
    209             index(j)=tmpindex
    210            endif
    211         enddo
     197    index(i)=i
     198      enddo
     199     
     200      DO i=0,mpi_size-2
     201        DO j=i+1,mpi_size-1
     202      IF (value(i)>value(j)) THEN
     203        tmpvalue=value(i)
     204        value(i)=value(j)
     205        value(j)=tmpvalue
     206
     207        tmpindex=index(i)
     208        index(i)=index(j)
     209        index(j)=tmpindex
     210       endif
     211    enddo
    212212      enddo
    213213     
    … …  
    215215      max_proc=index(mpi_size-1)           
    216216           
    217       do i=0,mpi_size-2
     217      DO i=0,mpi_size-2
    218218        minvalue=value(i)
    219219        min_proc=index(i)
    … …  
    222222             jj_nb_caldyn(min_proc)=jj_nb_caldyn(min_proc)+1
    223223             jj_nb_caldyn(max_proc)=jj_nb_caldyn(max_proc)-1
    224              exit
     224         exit
    225225           else
    226226             IF (timer_average(jj_nb_caldyn(min_proc)+1,timer_caldyn,min_proc)                 &
    227                 -timer_delta(jj_nb_caldyn(min_proc)+1,timer_caldyn,min_proc) < maxvalue) THEN
     227            -timer_delta(jj_nb_caldyn(min_proc)+1,timer_caldyn,min_proc) < maxvalue) THEN
    228228               jj_nb_caldyn(min_proc)=jj_nb_caldyn(min_proc)+1
    229229               jj_nb_caldyn(max_proc)=jj_nb_caldyn(max_proc)-1
    230230               exit
    231              endif
     231         endif
    232232           endif
    233233         endif
    … …  
    260260      CALL allgather_timer_average
    261261
    262       do i=0,mpi_size-1
     262      DO i=0,mpi_size-1
    263263        value(i)=timer_average(jj_nb_vanleer(i),timer_vanleer,i)
    264         index(i)=i
    265       enddo
    266      
    267       do i=0,mpi_size-2
    268         do j=i+1,mpi_size-1
    269           IF (value(i)>value(j)) THEN
    270             tmpvalue=value(i)
    271             value(i)=value(j)
    272             value(j)=tmpvalue
    273            
    274             tmpindex=index(i)
    275             index(i)=index(j)
    276             index(j)=tmpindex
    277            endif
    278         enddo
     264    index(i)=i
     265      enddo
     266     
     267      DO i=0,mpi_size-2
     268        DO j=i+1,mpi_size-1
     269      IF (value(i)>value(j)) THEN
     270        tmpvalue=value(i)
     271        value(i)=value(j)
     272        value(j)=tmpvalue
     273
     274        tmpindex=index(i)
     275        index(i)=index(j)
     276        index(j)=tmpindex
     277       endif
     278    enddo
    279279      enddo
    280280     
    … …  
    282282      max_proc=index(mpi_size-1)           
    283283           
    284       do i=0,mpi_size-2
     284      DO i=0,mpi_size-2
    285285        minvalue=value(i)
    286286        min_proc=index(i)
    … …  
    291291             jj_nb_vanleer(min_proc)=jj_nb_vanleer(min_proc)+1
    292292             jj_nb_vanleer(max_proc)=jj_nb_vanleer(max_proc)-1
    293              exit
     293         exit
    294294           else
    295295             IF (timer_average(jj_nb_vanleer(min_proc)+1,timer_vanleer,min_proc) < maxvalue) THEN
    … …  
    297297               jj_nb_vanleer(max_proc)=jj_nb_vanleer(max_proc)-1
    298298               exit
    299              endif
     299         endif
    300300           endif
    301301         endif
    … …  
    329329      CALL allgather_timer_average
    330330
    331       do i=0,mpi_size-1
     331      DO i=0,mpi_size-1
    332332        value(i)=timer_average(jj_nb_dissip(i),timer_dissip,i)
    333         index(i)=i
    334       enddo
    335      
    336       do i=0,mpi_size-2
    337         do j=i+1,mpi_size-1
    338           IF (value(i)>value(j)) THEN
    339             tmpvalue=value(i)
    340             value(i)=value(j)
    341             value(j)=tmpvalue
    342            
    343             tmpindex=index(i)
    344             index(i)=index(j)
    345             index(j)=tmpindex
    346            endif
    347         enddo
     333    index(i)=i
     334      enddo
     335     
     336      DO i=0,mpi_size-2
     337        DO j=i+1,mpi_size-1
     338      IF (value(i)>value(j)) THEN
     339        tmpvalue=value(i)
     340        value(i)=value(j)
     341        value(j)=tmpvalue
     342
     343        tmpindex=index(i)
     344        index(i)=index(j)
     345        index(j)=tmpindex
     346       endif
     347    enddo
    348348      enddo
    349349     
    … …  
    351351      max_proc=index(mpi_size-1)           
    352352           
    353       do i=0,mpi_size-2
     353      DO i=0,mpi_size-2
    354354        minvalue=value(i)
    355355        min_proc=index(i)
    … …  
    359359             jj_nb_dissip(min_proc)=jj_nb_dissip(min_proc)+1
    360360             jj_nb_dissip(max_proc)=jj_nb_dissip(max_proc)-1
    361              exit
     361         exit
    362362           else
    363363             IF (timer_average(jj_nb_dissip(min_proc)+1,timer_dissip,min_proc)         &
    364                 - timer_delta(jj_nb_dissip(min_proc)+1,timer_dissip,min_proc) < maxvalue) THEN
     364            - timer_delta(jj_nb_dissip(min_proc)+1,timer_dissip,min_proc) < maxvalue) THEN
    365365               jj_nb_dissip(min_proc)=jj_nb_dissip(min_proc)+1
    366366               jj_nb_dissip(max_proc)=jj_nb_dissip(max_proc)-1
    367367               exit
    368              endif
     368         endif
    369369           endif
    370370         endif
    … …  
    401401     
    402402      medium=0
    403       do i=0,mpi_size-1
     403      DO i=0,mpi_size-1
    404404        value(i)=timer_average(jj_nb_physic(i),timer_physic,i)
    405         medium=medium+value(i)
     405    medium=medium+value(i)
    406406      enddo   
    407407     
    … …  
    409409      NbTot=0
    410410      IF (CPPKEY_PHYS) THEN
    411       do i=0,mpi_size-1
     411      DO i=0,mpi_size-1
    412412        Inc(i)=nint(klon_mpi_para_nb(i)*(medium-value(i))/value(i))
    413413        NbTot=NbTot+Inc(i) 
    … …  
    418418      else
    419419        Sgn=-1
    420         NbTot=-NbTot
     420    NbTot=-NbTot
    421421      endif
    422422     
    423423      Index=0
    424       do i=1,NbTot
     424      DO i=1,NbTot
    425425        Inc(Index)=Inc(Index)-Sgn
    426         Index=Index+1
    427         IF (Index>mpi_size-1) Index=0
    428       enddo
    429      
    430       do i=0,mpi_size-1
     426    Index=Index+1
     427    IF (Index>mpi_size-1) Index=0
     428      enddo
     429     
     430      DO i=0,mpi_size-1
    431431        distrib_phys(i)=klon_mpi_para_nb(i)+inc(i)
    432432      enddo
    … …  
    457457     
    458458      IF (ierr==0) THEN
    459 !       write (unit_number,*) '*** Bandes caldyn ***'
    460         do i=0,mpi_size-1
     459!    write (unit_number,*) '*** Bandes caldyn ***'
     460    DO i=0,mpi_size-1
    461461          write (unit_number,*) i,jj_nb_caldyn(i)
    462462        enddo
    463463       
    464 !       write (unit_number,*) '*** Bandes vanleer ***'
    465         do i=0,mpi_size-1
     464!    write (unit_number,*) '*** Bandes vanleer ***'
     465        DO i=0,mpi_size-1
    466466          write (unit_number,*) i,jj_nb_vanleer(i)
    467467        enddo
    468468       
    469469!        write (unit_number,*) '*** Bandes dissip ***'
    470         do i=0,mpi_size-1
     470        DO i=0,mpi_size-1
    471471          write (unit_number,*) i,jj_nb_dissip(i)
    472472        enddo
    473473       
    474         do i=0,mpi_size-1
     474    DO i=0,mpi_size-1
    475475          write (unit_number,*) i,distrib_phys(i)
    476476        enddo
    477        
     477
    478478        CLOSE(unit_number)   
    479479      else
Note: See TracChangeset for help on using the changeset viewer.