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

Location:
LMDZ5/trunk/libf/dynphy_lonlat
Files:
7 edited

Legend:

Unmodified
Added
Removed

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

    r2597 r2600  
    3535#endif
    3636      USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi
     37      USE comvert_mod, ONLY: preff, presnivs
    3738     
    3839      IMPLICIT NONE
    … …  
    9091c    ------------------
    9192
    92 #include "dimensions.h"
    93 #include "paramet.h"
    94 #include "temps.h"
     93      include "dimensions.h"
     94      include "paramet.h"
     95      include "temps.h"
    9596
    9697      INTEGER ngridmx
    9798      PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    9899
    99 #include "comvert.h"
    100 #include "comgeom2.h"
    101 #include "iniprint.h"
     100      include "comgeom2.h"
     101      include "iniprint.h"
    102102
    103103c    Arguments :

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

    r2597 r2600  
    5050      USE callphysiq_mod, ONLY: call_physiq
    5151#endif
     52      USE comvert_mod, ONLY: preff, presnivs
    5253      USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi
    5354
    … …  
    115116      PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    116117
    117       include "comvert.h"
    118118      include "comgeom2.h"
    119119      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
    … …  
    382382           enddo
    383383         ENDDO
    384 c$OMP END DO NOWAIT     
     384c$OMP END DO NOWAIT       
    385385
    386386c      CALL gr_dyn_fi_p(llm,iip1,jjp1,klon,pphi,zphi)
    … …  
    403403c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    404404      DO l=1,llm
    405         DO ig=1,klon
    406            zphi(ig,l)=zphi(ig,l)-zphis(ig)
    407         ENDDO
     405        DO ig=1,klon
     406           zphi(ig,l)=zphi(ig,l)-zphis(ig)
     407        ENDDO
    408408      ENDDO
    409409c$OMP END DO NOWAIT
    … …  
    532532           DO i=2,iim
    533533             z1(i)   =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm)
    534            ENDDO
     534           ENDDO
    535535 
    536536           DO i=1,iim
    537537              zcos(i)    = COS(rlonv(i))*z1(i)
    538538              zsin(i)    = SIN(rlonv(i))*z1(i)
    539            ENDDO
     539           ENDDO
    540540 
    541541           zufi(klon,l)  = SSUM(iim,zcos,1)/pi
    … …  
    590590        allocate(zdpsrf_omp(klon))
    591591        allocate(flxwfi_omp(klon,llm))
    592         first_omp=.false.
     592        first_omp=.false.
    593593      endif
    594594       
    595            
     595           
    596596      klon=klon_omp
    597597      offset=klon_omp_begin-1
    … …  
    600600        do i=1,klon
    601601          zplev_omp(i,l)=zplev(offset+i,l)
    602         enddo
    603       enddo
    604          
     602        enddo
     603      enddo
     604         
    605605       do l=1,llm
    606606        do i=1,klon 
    607           zplay_omp(i,l)=zplay(offset+i,l)
    608         enddo
    609       enddo
    610        
    611       do l=1,llm
    612         do i=1,klon
    613           zphi_omp(i,l)=zphi(offset+i,l)
    614         enddo
    615       enddo
    616        
     607          zplay_omp(i,l)=zplay(offset+i,l)
     608        enddo
     609      enddo
     610       
     611      do l=1,llm
     612        do i=1,klon
     613          zphi_omp(i,l)=zphi(offset+i,l)
     614        enddo
     615      enddo
     616       
    617617      do i=1,klon
    618         zphis_omp(i)=zphis(offset+i)
     618        zphis_omp(i)=zphis(offset+i)
    619619      enddo
    620620     
    621        
     621       
    622622      do l=1,llm
    623623        presnivs_omp(l)=presnivs(l)
    624624      enddo
    625        
    626       do l=1,llm
    627         do i=1,klon
    628           zufi_omp(i,l)=zufi(offset+i,l)
    629         enddo
    630       enddo
    631        
    632       do l=1,llm
    633         do i=1,klon
    634           zvfi_omp(i,l)=zvfi(offset+i,l)
    635         enddo
    636       enddo
    637        
    638       do l=1,llm
    639         do i=1,klon
    640           zrfi_omp(i,l)=zrfi(offset+i,l)
    641         enddo
    642       enddo
    643        
    644       do l=1,llm
    645         do i=1,klon
    646           ztfi_omp(i,l)=ztfi(offset+i,l)
    647         enddo
    648       enddo
    649        
     625       
     626      do l=1,llm
     627        do i=1,klon
     628          zufi_omp(i,l)=zufi(offset+i,l)
     629        enddo
     630      enddo
     631       
     632      do l=1,llm
     633        do i=1,klon
     634          zvfi_omp(i,l)=zvfi(offset+i,l)
     635        enddo
     636      enddo
     637       
     638      do l=1,llm
     639        do i=1,klon
     640          zrfi_omp(i,l)=zrfi(offset+i,l)
     641        enddo
     642      enddo
     643       
     644      do l=1,llm
     645        do i=1,klon
     646          ztfi_omp(i,l)=ztfi(offset+i,l)
     647        enddo
     648      enddo
     649       
    650650      do iq=1,nqtot
    651651        do l=1,llm
    652652          do i=1,klon
    653653            zqfi_omp(i,l,iq)=zqfi(offset+i,l,iq)
    654           enddo
    655         enddo
    656       enddo
    657        
    658       do l=1,llm
    659         do i=1,klon
    660           zdufi_omp(i,l)=zdufi(offset+i,l)
    661         enddo
    662       enddo
    663        
    664       do l=1,llm
    665         do i=1,klon
    666           zdvfi_omp(i,l)=zdvfi(offset+i,l)
    667         enddo
    668       enddo
    669        
     654          enddo
     655        enddo
     656      enddo
     657       
     658      do l=1,llm
     659        do i=1,klon
     660          zdufi_omp(i,l)=zdufi(offset+i,l)
     661        enddo
     662      enddo
     663       
     664      do l=1,llm
     665        do i=1,klon
     666          zdvfi_omp(i,l)=zdvfi(offset+i,l)
     667        enddo
     668      enddo
     669       
    670670      do l=1,llm
    671671        do i=1,klon
    672672          zdtfi_omp(i,l)=zdtfi(offset+i,l)
    673         enddo
    674       enddo
    675        
     673        enddo
     674      enddo
     675       
    676676      do iq=1,nqtot
    677677        do l=1,llm
    678678          do i=1,klon
    679             zdqfi_omp(i,l,iq)=zdqfi(offset+i,l,iq)
    680           enddo
     679            zdqfi_omp(i,l,iq)=zdqfi(offset+i,l,iq)
     680          enddo
    681681        enddo
    682682      enddo
    683        
     683             
    684684      do i=1,klon
    685         zdpsrf_omp(i)=zdpsrf(offset+i)
     685        zdpsrf_omp(i)=zdpsrf(offset+i)
    686686      enddo
    687687
    … …  
    689689        do i=1,klon
    690690          flxwfi_omp(i,l)=flxwfi(offset+i,l)
    691         enddo
     691        enddo
    692692      enddo
    693693     
    … …  
    749749        do i=1,klon
    750750          zplev(offset+i,l)=zplev_omp(i,l)
    751         enddo
    752       enddo
    753          
     751        enddo
     752      enddo
     753         
    754754       do l=1,llm
    755755        do i=1,klon 
    756           zplay(offset+i,l)=zplay_omp(i,l)
    757         enddo
    758       enddo
    759        
    760       do l=1,llm
    761         do i=1,klon
    762           zphi(offset+i,l)=zphi_omp(i,l)
    763         enddo
    764       enddo
    765        
     756          zplay(offset+i,l)=zplay_omp(i,l)
     757        enddo
     758      enddo
     759       
     760      do l=1,llm
     761        do i=1,klon
     762          zphi(offset+i,l)=zphi_omp(i,l)
     763        enddo
     764      enddo
     765       
    766766
    767767      do i=1,klon
    768         zphis(offset+i)=zphis_omp(i)
     768        zphis(offset+i)=zphis_omp(i)
    769769      enddo
    770770     
    771        
     771       
    772772      do l=1,llm
    773773        presnivs(l)=presnivs_omp(l)
    774774      enddo
    775        
    776       do l=1,llm
    777         do i=1,klon
    778           zufi(offset+i,l)=zufi_omp(i,l)
    779         enddo
    780       enddo
    781        
    782       do l=1,llm
    783         do i=1,klon
    784           zvfi(offset+i,l)=zvfi_omp(i,l)
    785         enddo
    786       enddo
    787        
    788       do l=1,llm
    789         do i=1,klon
    790           ztfi(offset+i,l)=ztfi_omp(i,l)
    791         enddo
    792       enddo
    793        
     775       
     776      do l=1,llm
     777        do i=1,klon
     778          zufi(offset+i,l)=zufi_omp(i,l)
     779        enddo
     780      enddo
     781       
     782      do l=1,llm
     783        do i=1,klon
     784          zvfi(offset+i,l)=zvfi_omp(i,l)
     785        enddo
     786      enddo
     787       
     788      do l=1,llm
     789        do i=1,klon
     790          ztfi(offset+i,l)=ztfi_omp(i,l)
     791        enddo
     792      enddo
     793       
    794794      do iq=1,nqtot
    795795        do l=1,llm
    796796          do i=1,klon
    797797            zqfi(offset+i,l,iq)=zqfi_omp(i,l,iq)
    798           enddo
    799         enddo
    800       enddo
    801        
    802       do l=1,llm
    803         do i=1,klon
    804           zdufi(offset+i,l)=zdufi_omp(i,l)
    805         enddo
    806       enddo
    807        
    808       do l=1,llm
    809         do i=1,klon
    810           zdvfi(offset+i,l)=zdvfi_omp(i,l)
    811         enddo
    812       enddo
    813        
     798          enddo
     799        enddo
     800      enddo
     801       
     802      do l=1,llm
     803        do i=1,klon
     804          zdufi(offset+i,l)=zdufi_omp(i,l)
     805        enddo
     806      enddo
     807       
     808      do l=1,llm
     809        do i=1,klon
     810          zdvfi(offset+i,l)=zdvfi_omp(i,l)
     811        enddo
     812      enddo
     813       
    814814      do l=1,llm
    815815        do i=1,klon
    816816          zdtfi(offset+i,l)=zdtfi_omp(i,l)
    817         enddo
    818       enddo
    819        
     817        enddo
     818      enddo
     819       
    820820      do iq=1,nqtot
    821821        do l=1,llm
    822822          do i=1,klon
    823             zdqfi(offset+i,l,iq)=zdqfi_omp(i,l,iq)
    824           enddo
     823            zdqfi(offset+i,l,iq)=zdqfi_omp(i,l,iq)
     824          enddo
    825825        enddo
    826826      enddo
    827        
     827             
    828828      do i=1,klon
    829         zdpsrf(offset+i)=zdpsrf_omp(i)
     829        zdpsrf(offset+i)=zdpsrf_omp(i)
    830830      enddo
    831831     
    … …  
    929929       
    930930       if (.not. is_south_pole_dyn) then
    931         pdpsfi(:,jj_end:jj_end+1)=0
     931        pdpsfi(:,jj_end:jj_end+1)=0
    932932       endif
    933933c$OMP END MASTER
    … …  
    10451045        pdqfi(:,jj_begin:jj_end,l,:)=0.
    10461046      ENDDO
    1047 c$OMP END DO NOWAIT     
     1047c$OMP END DO NOWAIT       
    10481048
    10491049C
    … …  
    10561056!CDIR ON_ADB(index_j)
    10571057!cdir NODEP           
    1058              DO ig0=kstart,kend
     1058             DO ig0=kstart,kend
    10591059              i=index_i(ig0)
    10601060              j=index_j(ig0)
    … …  
    10621062              if (i==1) pdqfi(iip1,j,l,iiq) = zdqfi(ig0,l,iq)
    10631063            ENDDO
    1064            
    1065             IF (is_north_pole_dyn) then
    1066               DO i=1,iip1
     1064           
     1065            IF (is_north_pole_dyn) then
     1066              DO i=1,iip1
    10671067                pdqfi(i,1,l,iiq)    = zdqfi(1,l,iq)
    1068               ENDDO
    1069             ENDIF
    1070            
    1071             IF (is_south_pole_dyn) then
    1072               DO i=1,iip1
     1068              ENDDO
     1069            ENDIF
     1070           
     1071            IF (is_south_pole_dyn) then
     1072              DO i=1,iip1
    10731073                pdqfi(i,jjp1,l,iiq) = zdqfi(klon,l,iq)
    1074               ENDDO
    1075             ENDIF
    1076            
     1074              ENDDO
     1075            ENDIF
     1076           
    10771077         ENDDO
    1078 c$OMP END DO NOWAIT     
     1078c$OMP END DO NOWAIT       
    10791079      ENDDO
    10801080     
    … …  
    11361136           pdvfi(i,j,l)=0.5*(zdvfi2(ig0,l)+zdvfi2(ig0+iim,l))*cv(i,j)
    11371137           if (i==1) pdvfi(iip1,j,l) = 0.5*(zdvfi2(ig0,l)+
    1138      $                                      zdvfi2(ig0+iim,l))
    1139      $                                    *cv(i,j)
     1138     $                                            zdvfi2(ig0+iim,l))
     1139     $                                          *cv(i,j)
    11401140        enddo
    11411141         

  • 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         

  • LMDZ5/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90

    r2597 r2600  
    3636  USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, kappa, omeg, r, rad, &
    3737                          pi, jmp1
     38  USE comvert_mod, ONLY: pa, preff, pressure_exner
    3839
    3940  IMPLICIT NONE
    … …  
    4445  include "paramet.h"
    4546  include "comgeom2.h"
    46   include "comvert.h"
    4747  include "iniprint.h"
    4848  include "temps.h"

  • LMDZ5/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

    r2598 r2600  
    3737  USE assert_eq_m,    ONLY: assert_eq
    3838  USE comconst_mod, ONLY: pi, cpp, kappa
     39  USE comvert_mod, ONLY: ap, bp, preff, pressure_exner
    3940  IMPLICIT NONE
    4041
    … …  
    4647  include "paramet.h"
    4748  include "comgeom2.h"
    48   include "comvert.h"
    4949  include "temps.h"
    5050  include "comdissnew.h"

  • LMDZ5/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90

    r2597 r2600  
    1212                     prad,pg,pr,pcpp,iflag_phys)
    1313  USE dimphy, ONLY: init_dimphy
     14  USE comvert_mod, ONLY: preff, ap, bp, presnivs, scaleheight, pseudoalt
    1415  USE inigeomphy_mod, ONLY: inigeomphy
    1516  USE mod_grid_phy_lmdz, ONLY: klon_glo ! number of atmospheric columns (on full grid)
    … …  
    5051
    5152  include "dimensions.h"
    52   include "comvert.h"
    5353  include "iniprint.h"
    5454  include "temps.h"

  • LMDZ5/trunk/libf/dynphy_lonlat/phymar/iniphysiq_mod.F90

    r2599 r2600  
    1212                     prad,pg,pr,pcpp,iflag_phys)
    1313  USE dimphy, ONLY: init_dimphy
     14  USE comvert_mod, ONLY: preff, ap, bp, presnivs, scaleheight, pseudoalt
    1415  USE inigeomphy_mod, ONLY: inigeomphy
    1516  USE vertical_layers_mod, ONLY : init_vertical_layers
    … …  
    2930 
    3031 
    31   include "dimensions.h"
    32   include "comvert.h"
    3332  include "iniprint.h"
    3433
Note: See TracChangeset for help on using the changeset viewer.