Ignore:
Timestamp:
Jul 23, 2016, 7:45:38 AM (8 years ago)
Author:
Ehouarn Millour
Message:

Cleanup in the dynamics: turn comvert.h into module comvert_mod.F90
EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dynphy_lonlat/calfis_p.F

    r2597 r2600  
    4747      USE callphysiq_mod, ONLY: call_physiq
    4848#endif
     49      USE comvert_mod, ONLY: preff, presnivs
    4950      USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi
    5051
    … …  
    111112      PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    112113
    113 #include "comvert.h"
    114114#include "comgeom2.h"
    115115#include "iniprint.h"
    … …  
    277277      ALLOCATE(flxwfi(klon,llm))
    278278c$OMP END MASTER
    279 c$OMP BARRIER    
     279c$OMP BARRIER         
    280280      ELSE
    281281          debut = .FALSE.
    … …  
    365365           enddo
    366366         ENDDO
    367 c$OMP END DO NOWAIT     
     367c$OMP END DO NOWAIT       
    368368      ENDDO
    369369
    … …  
    380380c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    381381      DO l=1,llm
    382         DO ig=1,klon
    383            zphi(ig,l)=zphi(ig,l)-zphis(ig)
    384         ENDDO
     382        DO ig=1,klon
     383           zphi(ig,l)=zphi(ig,l)-zphis(ig)
     384        ENDDO
    385385      ENDDO
    386386c$OMP END DO NOWAIT
    … …  
    510510           DO i=2,iim
    511511             z1(i)   =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm)
    512            ENDDO
     512           ENDDO
    513513 
    514514           DO i=1,iim
    515515              zcos(i)    = COS(rlonv(i))*z1(i)
    516516              zsin(i)    = SIN(rlonv(i))*z1(i)
    517            ENDDO
     517           ENDDO
    518518 
    519519           zufi(klon,l)  = SSUM(iim,zcos,1)/pi
    … …  
    556556        allocate(zdpsrf_omp(klon))
    557557        allocate(flxwfi_omp(klon,llm))
    558         first_omp=.false.
     558        first_omp=.false.
    559559      endif
    560560       
    561            
     561           
    562562      klon=klon_omp
    563563      offset=klon_omp_begin-1
    … …  
    566566        do i=1,klon
    567567          zplev_omp(i,l)=zplev(offset+i,l)
    568         enddo
    569       enddo
    570          
     568        enddo
     569      enddo
     570         
    571571       do l=1,llm
    572572        do i=1,klon 
    573           zplay_omp(i,l)=zplay(offset+i,l)
    574         enddo
    575       enddo
    576        
    577       do l=1,llm
    578         do i=1,klon
    579           zphi_omp(i,l)=zphi(offset+i,l)
    580         enddo
    581       enddo
    582        
     573          zplay_omp(i,l)=zplay(offset+i,l)
     574        enddo
     575      enddo
     576       
     577      do l=1,llm
     578        do i=1,klon
     579          zphi_omp(i,l)=zphi(offset+i,l)
     580        enddo
     581      enddo
     582       
    583583      do i=1,klon
    584         zphis_omp(i)=zphis(offset+i)
     584        zphis_omp(i)=zphis(offset+i)
    585585      enddo
    586586     
    587        
     587       
    588588      do l=1,llm
    589589        presnivs_omp(l)=presnivs(l)
    590590      enddo
    591        
    592       do l=1,llm
    593         do i=1,klon
    594           zufi_omp(i,l)=zufi(offset+i,l)
    595         enddo
    596       enddo
    597        
    598       do l=1,llm
    599         do i=1,klon
    600           zvfi_omp(i,l)=zvfi(offset+i,l)
    601         enddo
    602       enddo
    603        
    604       do l=1,llm
    605         do i=1,klon
    606           zrfi_omp(i,l)=zrfi(offset+i,l)
    607         enddo
    608       enddo
    609        
    610        
    611       do l=1,llm
    612         do i=1,klon
    613           ztfi_omp(i,l)=ztfi(offset+i,l)
    614         enddo
    615       enddo
    616        
     591       
     592      do l=1,llm
     593        do i=1,klon
     594          zufi_omp(i,l)=zufi(offset+i,l)
     595        enddo
     596      enddo
     597       
     598      do l=1,llm
     599        do i=1,klon
     600          zvfi_omp(i,l)=zvfi(offset+i,l)
     601        enddo
     602      enddo
     603       
     604      do l=1,llm
     605        do i=1,klon
     606          zrfi_omp(i,l)=zrfi(offset+i,l)
     607        enddo
     608      enddo
     609       
     610       
     611      do l=1,llm
     612        do i=1,klon
     613          ztfi_omp(i,l)=ztfi(offset+i,l)
     614        enddo
     615      enddo
     616       
    617617      do iq=1,nqtot
    618618        do l=1,llm
    619619          do i=1,klon
    620620            zqfi_omp(i,l,iq)=zqfi(offset+i,l,iq)
    621           enddo
    622         enddo
    623       enddo
    624        
    625       do l=1,llm
    626         do i=1,klon
    627           zdufi_omp(i,l)=zdufi(offset+i,l)
    628         enddo
    629       enddo
    630        
    631       do l=1,llm
    632         do i=1,klon
    633           zdvfi_omp(i,l)=zdvfi(offset+i,l)
    634         enddo
    635       enddo
    636        
     621          enddo
     622        enddo
     623      enddo
     624       
     625      do l=1,llm
     626        do i=1,klon
     627          zdufi_omp(i,l)=zdufi(offset+i,l)
     628        enddo
     629      enddo
     630       
     631      do l=1,llm
     632        do i=1,klon
     633          zdvfi_omp(i,l)=zdvfi(offset+i,l)
     634        enddo
     635      enddo
     636       
    637637      do l=1,llm
    638638        do i=1,klon
    639639          zdtfi_omp(i,l)=zdtfi(offset+i,l)
    640         enddo
    641       enddo
    642        
     640        enddo
     641      enddo
     642       
    643643      do iq=1,nqtot
    644644        do l=1,llm
    645645          do i=1,klon
    646             zdqfi_omp(i,l,iq)=zdqfi(offset+i,l,iq)
    647           enddo
     646            zdqfi_omp(i,l,iq)=zdqfi(offset+i,l,iq)
     647          enddo
    648648        enddo
    649649      enddo
    650        
     650             
    651651      do i=1,klon
    652         zdpsrf_omp(i)=zdpsrf(offset+i)
     652        zdpsrf_omp(i)=zdpsrf(offset+i)
    653653      enddo
    654654
    … …  
    656656        do i=1,klon
    657657          flxwfi_omp(i,l)=flxwfi(offset+i,l)
    658         enddo
     658        enddo
    659659      enddo
    660660     
    … …  
    713713        do i=1,klon
    714714          zplev(offset+i,l)=zplev_omp(i,l)
    715         enddo
    716       enddo
    717          
     715        enddo
     716      enddo
     717         
    718718       do l=1,llm
    719719        do i=1,klon 
    720           zplay(offset+i,l)=zplay_omp(i,l)
    721         enddo
    722       enddo
    723        
    724       do l=1,llm
    725         do i=1,klon
    726           zphi(offset+i,l)=zphi_omp(i,l)
    727         enddo
    728       enddo
    729        
     720          zplay(offset+i,l)=zplay_omp(i,l)
     721        enddo
     722      enddo
     723       
     724      do l=1,llm
     725        do i=1,klon
     726          zphi(offset+i,l)=zphi_omp(i,l)
     727        enddo
     728      enddo
     729       
    730730
    731731      do i=1,klon
    732         zphis(offset+i)=zphis_omp(i)
     732        zphis(offset+i)=zphis_omp(i)
    733733      enddo
    734734     
    735        
     735       
    736736      do l=1,llm
    737737        presnivs(l)=presnivs_omp(l)
    738738      enddo
    739        
    740       do l=1,llm
    741         do i=1,klon
    742           zufi(offset+i,l)=zufi_omp(i,l)
    743         enddo
    744       enddo
    745        
    746       do l=1,llm
    747         do i=1,klon
    748           zvfi(offset+i,l)=zvfi_omp(i,l)
    749         enddo
    750       enddo
    751        
    752       do l=1,llm
    753         do i=1,klon
    754           ztfi(offset+i,l)=ztfi_omp(i,l)
    755         enddo
    756       enddo
    757        
     739       
     740      do l=1,llm
     741        do i=1,klon
     742          zufi(offset+i,l)=zufi_omp(i,l)
     743        enddo
     744      enddo
     745       
     746      do l=1,llm
     747        do i=1,klon
     748          zvfi(offset+i,l)=zvfi_omp(i,l)
     749        enddo
     750      enddo
     751       
     752      do l=1,llm
     753        do i=1,klon
     754          ztfi(offset+i,l)=ztfi_omp(i,l)
     755        enddo
     756      enddo
     757       
    758758      do iq=1,nqtot
    759759        do l=1,llm
    760760          do i=1,klon
    761761            zqfi(offset+i,l,iq)=zqfi_omp(i,l,iq)
    762           enddo
    763         enddo
    764       enddo
    765        
    766       do l=1,llm
    767         do i=1,klon
    768           zdufi(offset+i,l)=zdufi_omp(i,l)
    769         enddo
    770       enddo
    771        
    772       do l=1,llm
    773         do i=1,klon
    774           zdvfi(offset+i,l)=zdvfi_omp(i,l)
    775         enddo
    776       enddo
    777        
     762          enddo
     763        enddo
     764      enddo
     765       
     766      do l=1,llm
     767        do i=1,klon
     768          zdufi(offset+i,l)=zdufi_omp(i,l)
     769        enddo
     770      enddo
     771       
     772      do l=1,llm
     773        do i=1,klon
     774          zdvfi(offset+i,l)=zdvfi_omp(i,l)
     775        enddo
     776      enddo
     777       
    778778      do l=1,llm
    779779        do i=1,klon
    780780          zdtfi(offset+i,l)=zdtfi_omp(i,l)
    781         enddo
    782       enddo
    783        
     781        enddo
     782      enddo
     783       
    784784      do iq=1,nqtot
    785785        do l=1,llm
    786786          do i=1,klon
    787             zdqfi(offset+i,l,iq)=zdqfi_omp(i,l,iq)
    788           enddo
     787            zdqfi(offset+i,l,iq)=zdqfi_omp(i,l,iq)
     788          enddo
    789789        enddo
    790790      enddo
    791        
     791             
    792792      do i=1,klon
    793         zdpsrf(offset+i)=zdpsrf_omp(i)
     793        zdpsrf(offset+i)=zdpsrf_omp(i)
    794794      enddo
    795795     
    … …  
    892892       pdpsfi(:,jj_begin)=0   
    893893       if (.not. is_south_pole_dyn) then
    894         pdpsfi(:,jj_end)=0
     894        pdpsfi(:,jj_end)=0
    895895       endif
    896896c$OMP END MASTER
    … …  
    976976        pdqfi(:,:,l,:)=0.
    977977      ENDDO
    978 c$OMP END DO NOWAIT     
     978c$OMP END DO NOWAIT       
    979979
    980980C
    … …  
    987987!CDIR ON_ADB(index_j)
    988988!cdir NODEP           
    989              DO ig0=kstart,kend
     989             DO ig0=kstart,kend
    990990              i=index_i(ig0)
    991991              j=index_j(ig0)
    … …  
    993993              if (i==1) pdqfi(iip1,j,l,iiq) = zdqfi(ig0,l,iq)
    994994            ENDDO
    995            
    996             IF (is_north_pole_dyn) then
    997               DO i=1,iip1
     995           
     996            IF (is_north_pole_dyn) then
     997              DO i=1,iip1
    998998                pdqfi(i,1,l,iiq)    = zdqfi(1,l,iq)
    999               ENDDO
    1000             ENDIF
    1001            
    1002             IF (is_south_pole_dyn) then
    1003               DO i=1,iip1
     999              ENDDO
     1000            ENDIF
     1001           
     1002            IF (is_south_pole_dyn) then
     1003              DO i=1,iip1
    10041004                pdqfi(i,jjp1,l,iiq) = zdqfi(klon,l,iq)
    1005               ENDDO
    1006             ENDIF
    1007            
     1005              ENDDO
     1006            ENDIF
     1007           
    10081008         ENDDO
    1009 c$OMP END DO NOWAIT     
     1009c$OMP END DO NOWAIT       
    10101010      ENDDO
    10111011     
    … …  
    10671067           pdvfi(i,j,l)=0.5*(zdvfi2(ig0,l)+zdvfi2(ig0+iim,l))*cv(i,j)
    10681068           if (i==1) pdvfi(iip1,j,l) = 0.5*(zdvfi2(ig0,l)+
    1069      $                                      zdvfi2(ig0+iim,l))
    1070      $                                    *cv(i,j)
     1069     $                                            zdvfi2(ig0+iim,l))
     1070     $                                          *cv(i,j)
    10711071        enddo
    10721072         
Note: See TracChangeset for help on using the changeset viewer.