Changeset 2600


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
Files:
1 added
1 deleted
123 edited
1 moved

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3d/advect.F

    r2597 r2600  
    2828      include "dimensions.h"
    2929      include "paramet.h"
    30       include "comvert.h"
    3130      include "comgeom.h"
    3231      include "logic.h"

  • LMDZ5/trunk/libf/dyn3d/advtrac.F90

    r2597 r2600  
    1717  include "dimensions.h"
    1818  include "paramet.h"
    19   include "comvert.h"
    2019  include "comdissip.h"
    2120  include "comgeom2.h"

  • LMDZ5/trunk/libf/dyn3d/bilan_dyn.F

    r2597 r2600  
    1414#endif
    1515      USE comconst_mod, ONLY: pi, cpp
     16      USE comvert_mod, ONLY: presnivs
    1617
    1718      IMPLICIT NONE
    … …  
    1920      include "dimensions.h"
    2021      include "paramet.h"
    21       include "comvert.h"
    2222      include "comgeom2.h"
    2323      include "temps.h"

  • LMDZ5/trunk/libf/dyn3d/caldyn.F

    r2597 r2600  
    66     $  phi,conser,du,dv,dteta,dp,w,pbaru,pbarv,time )
    77
     8     
     9      USE comvert_mod, ONLY: ap, bp
     10     
    811      IMPLICIT NONE
    912
    … …  
    2629      include "dimensions.h"
    2730      include "paramet.h"
    28       include "comvert.h"
    2931      include "comgeom.h"
    3032

  • LMDZ5/trunk/libf/dyn3d/dudv2.F

    r1907 r2600  
    2626c=======================================================================
    2727c
    28 #include "dimensions.h"
    29 #include "paramet.h"
    30 #include "comvert.h"
     28      include "dimensions.h"
     29      include "paramet.h"
    3130
    3231      REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ),

  • LMDZ5/trunk/libf/dyn3d/dynetat0.f90

    r2598 r2600  
    1111  USE control_mod, ONLY: planet_type
    1212  USE assert_eq_m, ONLY: assert_eq
     13  USE comvert_mod, ONLY: pa,preff
    1314  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa, lllm, omeg, rad
    1415  USE serre_mod, ONLY: clon, clat, grossismx, grossismy
     16
    1517  IMPLICIT NONE
    1618  include "dimensions.h"
    1719  include "paramet.h"
    1820  include "temps.h"
    19   include "comvert.h"
    2021  include "comgeom2.h"
    2122  include "ener.h"

  • LMDZ5/trunk/libf/dyn3d/dynredem.F90

    r2598 r2600  
    1111                    NF90_CLOSE,  NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER
    1212  USE dynredem_mod, ONLY: cre_var, put_var1, put_var2, err, modname, fil
     13  USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff, &
     14                              nivsig,nivsigs
    1315  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    1416  USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, &
    15                         taux,tauy
     17                              taux,tauy
    1618  IMPLICIT NONE
    1719  include "dimensions.h"
    1820  include "paramet.h"
    19   include "comvert.h"
    2021  include "comgeom2.h"
    2122  include "temps.h"
    … …  
    173174  include "paramet.h"
    174175  include "description.h"
    175   include "comvert.h"
    176176  include "comgeom.h"
    177177  include "temps.h"

  • LMDZ5/trunk/libf/dyn3d/fluxstokenc.F

    r2597 r2600  
    1818      include "dimensions.h"
    1919      include "paramet.h"
    20       include "comvert.h"
    2120      include "comgeom.h"
    2221      include "tracstoke.h"
    … …  
    3534
    3635      REAL pbarvst(iip1,jjp1,llm),zistdyn
    37         real dtcum
     36        real dtcum
    3837
    3938      INTEGER iadvtr,ndex(1)
    … …  
    5958      if(first) then
    6059
    61         CALL initfluxsto( 'fluxstoke',
     60        CALL initfluxsto( 'fluxstoke',
    6261     .  time_step,istdyn* time_step,istdyn* time_step,
    6362     .  fluxid,fluxvid,fluxdid)
    64        
    65         ndex(1) = 0
     63       
     64        ndex(1) = 0
    6665        call histwrite(fluxid, 'phis', 1, phis, iip1*jjp1, ndex)
    6766        call histwrite(fluxid, 'aire', 1, aire, iip1*jjp1, ndex)
    68        
    69         ndex(1) = 0
     67       
     68        ndex(1) = 0
    7069        nscal = 1
    7170        tst(1) = time_step
    … …  
    7574        istp(1)= istphy
    7675        call histwrite(fluxdid, 'istphy', 1, istp, nscal, ndex)
    77        
    78         first = .false.
     76       
     77        first = .false.
    7978
    8079      endif
    … …  
    140139
    141140         iadvtr=0
    142         write(lunout,*)'ITAU auquel on stoke les fluxmasses',itau
    143        
    144         call histwrite(fluxid, 'masse', itau, massem,
     141        write(lunout,*)'ITAU auquel on stoke les fluxmasses',itau
     142       
     143        call histwrite(fluxid, 'masse', itau, massem,
    145144     .               iip1*jjp1*llm, ndex)
    146        
    147         call histwrite(fluxid, 'pbaru', itau, pbarug,
     145       
     146        call histwrite(fluxid, 'pbaru', itau, pbarug,
    148147     .               iip1*jjp1*llm, ndex)
    149        
    150         call histwrite(fluxvid, 'pbarv', itau, pbarvg,
     148       
     149        call histwrite(fluxvid, 'pbarv', itau, pbarvg,
    151150     .               iip1*jjm*llm, ndex)
    152        
     151       
    153152        call histwrite(fluxid, 'w' ,itau, wg,
    154153     .             iip1*jjp1*llm, ndex)
    155        
    156         call histwrite(fluxid, 'teta' ,itau, tetac,
     154       
     155        call histwrite(fluxid, 'teta' ,itau, tetac,
    157156     .             iip1*jjp1*llm, ndex)
    158        
    159         call histwrite(fluxid, 'phi' ,itau, phic,
     157       
     158        call histwrite(fluxid, 'phi' ,itau, phic,
    160159     .             iip1*jjp1*llm, ndex)
    161        
     160       
    162161C
    163162

  • LMDZ5/trunk/libf/dyn3d/gcm.F90

    r2598 r2600  
    6969  include "paramet.h"
    7070  include "comdissnew.h"
    71   include "comvert.h"
    7271  include "comgeom.h"
    7372  include "logic.h"

  • LMDZ5/trunk/libf/dyn3d/groupe.F

    r2597 r2600  
    2222      include "paramet.h"
    2323      include "comgeom2.h"
    24       include "comvert.h"
    2524
    2625!     integer ngroup

  • LMDZ5/trunk/libf/dyn3d/guide_mod.F90

    r2598 r2600  
    314314    USE control_mod, ONLY: day_step, iperiod
    315315    USE comconst_mod, ONLY: dtvr, daysec
     316    USE comvert_mod, ONLY: ap, bp, preff, presnivs
    316317 
    317318    IMPLICIT NONE
    … …  
    319320    INCLUDE "dimensions.h"
    320321    INCLUDE "paramet.h"
    321     INCLUDE "comvert.h"
    322322
    323323    ! Variables entree
    … …  
    607607  use exner_milieu_m, only: exner_milieu
    608608  use comconst_mod, only: kappa, cpp
     609  use comvert_mod, only: preff, pressure_exner, bp, ap
    609610  IMPLICIT NONE
    610611
    611612  include "dimensions.h"
    612613  include "paramet.h"
    613   include "comvert.h"
    614614  include "comgeom2.h"
    615615
    … …  
    15171517
    15181518    USE comconst_mod, ONLY: pi
     1519    USE comvert_mod, ONLY: presnivs
    15191520   
    15201521    IMPLICIT NONE
    … …  
    15241525    INCLUDE "netcdf.inc"
    15251526    INCLUDE "comgeom2.h"
    1526     INCLUDE "comvert.h"
    15271527   
    15281528    ! Variables entree

  • LMDZ5/trunk/libf/dyn3d/iniacademic.F90

    r2597 r2600  
    1717  use exner_milieu_m, only: exner_milieu
    1818  USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
     19  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
    1920
    2021  !   Author:    Frederic Hourdin      original: 15/01/93
    … …  
    2930  include "dimensions.h"
    3031  include "paramet.h"
    31   include "comvert.h"
    3232  include "comgeom.h"
    3333  include "academic.h"

  • LMDZ5/trunk/libf/dyn3d/integrd.F

    r2598 r2600  
    99      use control_mod, only : planet_type
    1010      use comconst_mod, only: pi
     11      use comvert_mod, only: ap, bp
    1112
    1213      IMPLICIT NONE
    … …  
    3132      include "paramet.h"
    3233      include "comgeom.h"
    33       include "comvert.h"
    3434      include "logic.h"
    3535      include "temps.h"

  • LMDZ5/trunk/libf/dyn3d/leapfrog.F

    r2598 r2600  
    2020      use exner_hyb_m, only: exner_hyb
    2121      use exner_milieu_m, only: exner_milieu
     22      USE comvert_mod, ONLY: ap,bp,pressure_exner,presnivs
    2223      USE comconst_mod, ONLY: cpp, dtphys, dtvr, pi, ihf
    2324
    … …  
    5859      include "paramet.h"
    5960      include "comdissnew.h"
    60       include "comvert.h"
    6161      include "comgeom.h"
    6262      include "logic.h"
    … …  
    712712            IF( MOD(itau,iecri).EQ.0) THEN
    713713             ! Ehouarn: output only during LF or Backward Matsuno
    714              if (leapf.or.(.not.leapf.and.(.not.forward))) then
     714             if (leapf.or.(.not.leapf.and.(.not.forward))) then
    715715              CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
    716716              unat=0.
    … …  
    722722              if (ok_dyn_ins) then
    723723!               write(lunout,*) "leapfrog: call writehist, itau=",itau
    724                CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
     724               CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
    725725!               call WriteField('ucov',reshape(ucov,(/iip1,jmp1,llm/)))
    726726!               call WriteField('vcov',reshape(vcov,(/iip1,jjm,llm/)))
    … …  
    852852!                write(lunout,*) "leapfrog: call writehist (b)",
    853853!     &                        itau,iecri
    854                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
     854                CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
    855855              endif ! of if (ok_dyn_ins)
    856856#endif

  • LMDZ5/trunk/libf/dyn3d/qminimum.F

    r2286 r2600  
    1010c             pour l'eau vapeur et l'eau liquide
    1111c
    12 #include "dimensions.h"
    13 #include "paramet.h"
    14 #include "comvert.h"
     12      include "dimensions.h"
     13      include "paramet.h"
    1514c
    1615      INTEGER nqtot

  • LMDZ5/trunk/libf/dyn3d/sw_case_williamson91_6.F

    r2597 r2600  
    2727c=======================================================================
    2828      USE comconst_mod, ONLY: cpp, omeg, rad
     29      USE comvert_mod, ONLY: ap, bp, preff
    2930     
    3031      IMPLICIT NONE
    … …  
    3536      include "dimensions.h"
    3637      include "paramet.h"
    37       include "comvert.h"
    3838      include "comgeom.h"
    3939      include "iniprint.h"

  • LMDZ5/trunk/libf/dyn3d/top_bound.F

    r2597 r2600  
    66      USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound,
    77     &                        tau_top_bound
     8      USE comvert_mod, ONLY: presnivs, preff, scaleheight
    89     
    910      IMPLICIT NONE
    … …  
    1112      include "dimensions.h"
    1213      include "paramet.h"
    13       include "comvert.h"
    1414      include "comgeom2.h"
    1515

  • LMDZ5/trunk/libf/dyn3d/vlsplt.F

    r2597 r2600  
    2424      include "paramet.h"
    2525      include "logic.h"
    26       include "comvert.h"
    2726
    2827c
    … …  
    157156      include "paramet.h"
    158157      include "logic.h"
    159       include "comvert.h"
    160158      include "iniprint.h"
    161159c
    … …  
    532530      include "paramet.h"
    533531      include "logic.h"
    534       include "comvert.h"
    535532      include "comgeom.h"
    536533c
    … …  
    894891      include "paramet.h"
    895892      include "logic.h"
    896       include "comvert.h"
    897893c
    898894c

  • LMDZ5/trunk/libf/dyn3d/vlspltqs.F

    r2597 r2600  
    3131      include "paramet.h"
    3232      include "logic.h"
    33       include "comvert.h"
    3433
    3534c
    … …  
    197196      include "paramet.h"
    198197      include "logic.h"
    199       include "comvert.h"
    200198c
    201199c
    … …  
    568566      include "paramet.h"
    569567      include "logic.h"
    570       include "comvert.h"
    571568      include "comgeom.h"
    572569c

  • LMDZ5/trunk/libf/dyn3d_common/adaptdt.F

    r2597 r2600  
    99      IMPLICIT NONE
    1010
    11 #include "dimensions.h"
    12 c#include "paramr2.h"
    13 #include "paramet.h"
    14 #include "comdissip.h"
    15 #include "comvert.h"
    16 #include "comgeom2.h"
    17 #include "logic.h"
    18 #include "temps.h"
    19 #include "ener.h"
    20 #include "description.h"
     11      include "dimensions.h"
     12      include "paramet.h"
     13      include "comdissip.h"
     14      include "comgeom2.h"
     15      include "logic.h"
     16      include "temps.h"
     17      include "ener.h"
     18      include "description.h"
    2119
    2220c----------------------------------------------------------

  • LMDZ5/trunk/libf/dyn3d_common/advn.F

    r2597 r2600  
    1717      IMPLICIT NONE
    1818c
    19 #include "dimensions.h"
    20 #include "paramet.h"
    21 #include "logic.h"
    22 #include "comvert.h"
    23 #include "comgeom.h"
    24 #include "iniprint.h"
     19      include "dimensions.h"
     20      include "paramet.h"
     21      include "logic.h"
     22      include "comgeom.h"
     23      include "iniprint.h"
    2524
    2625c
    … …  
    483482      IMPLICIT NONE
    484483c
    485 #include "dimensions.h"
    486 #include "paramet.h"
    487 #include "logic.h"
    488 #include "comvert.h"
    489 #include "iniprint.h"
     484      include "dimensions.h"
     485      include "paramet.h"
     486      include "logic.h"
     487      include "iniprint.h"
    490488c
    491489c

  • LMDZ5/trunk/libf/dyn3d_common/advx.F

    r2597 r2600  
    1818C  sm,s0,sx,sy,sz                                                C
    1919C  sont les arguments de sortie pour le s-pg                     C
    20 C                                                                C
     20C                                                                C
    2121CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    2222C
    2323C  parametres principaux du modele
    2424C
    25 #include "dimensions.h"
    26 #include "paramet.h"
    27 #include "comvert.h"
     25      include "dimensions.h"
     26      include "paramet.h"
    2827
    2928C  Arguments :

  • LMDZ5/trunk/libf/dyn3d_common/advxp.F

    r2597 r2600  
    1313C  parametres principaux du modele
    1414C
    15 #include "dimensions.h"
    16 #include "paramet.h"
    17 #include "comvert.h"
     15      include "dimensions.h"
     16      include "paramet.h"
    1817
    1918       INTEGER ntra
    … …  
    102101c         IF (S0(i,j,l,ntra) .lt. 0. ) THEN
    103102c         PRINT*,'S0(',i,j,l,')=',S0(i,j,l,ntra)
    104 c            print*, 'SSX(',i,j,l,')=',SSX(i,j,l,ntra)
     103c             print*, 'SSX(',i,j,l,')=',SSX(i,j,l,ntra)
    105104c         print*, 'SY(',i,j,l,')=',SY(i,j,l,ntra)
    106105c         print*, 'SZ(',i,j,l,')=',SZ(i,j,l,ntra)
    … …  
    119118      DO j = 1, jjp1
    120119      DO i = 1, iim
    121         sqi = sqi + S0(i,j,l,ntra)
     120        sqi = sqi + S0(i,j,l,ntra)
    122121      END DO
    123122      END DO
    … …  
    611610c      DO 9999 j = 1, jjp1
    612611c      DO 9999 i = 1, iip1
    613 c          IF (S0(i,j,l,ntra).lt.0..and.LIMIT) THEN
     612c           IF (S0(i,j,l,ntra).lt.0..and.LIMIT) THEN
    614613c           PRINT*, '-------------------'
    615 c               PRINT*, 'En fin de ADVXP'
     614c                PRINT*, 'En fin de ADVXP'
    616615c           PRINT*,'S0(',i,j,l,')=',S0(i,j,l,ntra)
    617 c               print*, 'SSX(',i,j,l,')=',SSX(i,j,l,ntra)
     616c                print*, 'SSX(',i,j,l,')=',SSX(i,j,l,ntra)
    618617c           print*, 'SY(',i,j,l,')=',SY(i,j,l,ntra)
    619 c               print*, 'SZ(',i,j,l,')=',SZ(i,j,l,ntra)
     618c               print*, 'SZ(',i,j,l,')=',SZ(i,j,l,ntra)
    620619c            WRITE (*,*) 'On arrete !! - pbl en fin de ADVXP'
    621620c            STOP
    … …  
    628627         SM(iip1,j,l) = SM(1,j,l)
    629628         S0(iip1,j,l,ntra) = S0(1,j,l,ntra)
    630          SSX(iip1,j,l,ntra) = SSX(1,j,l,ntra)
    631         SY(iip1,j,l,ntra) = SY(1,j,l,ntra)
    632         SZ(iip1,j,l,ntra) = SZ(1,j,l,ntra)
     629             SSX(iip1,j,l,ntra) = SSX(1,j,l,ntra)
     630            SY(iip1,j,l,ntra) = SY(1,j,l,ntra)
     631            SZ(iip1,j,l,ntra) = SZ(1,j,l,ntra)
    633632      END DO
    634633      END DO

  • LMDZ5/trunk/libf/dyn3d_common/advy.F

    r2597 r2600  
    99C  first-order moments (SOM) advection of tracer in Y direction  C
    1010C                                                                C
    11 C  Source : Pascal Simon ( Meteo, CNRM )                        C
    12 C  Adaptation : A.A. (LGGE)                                     C
     11C  Source : Pascal Simon ( Meteo, CNRM )                        C
     12C  Adaptation : A.A. (LGGE)                                     C
    1313C  Derniere Modif : 15/12/94 LAST
    14 C                                                                C
    15 C  sont les arguments d'entree pour le s-pg                      C
    16 C                                                                C
    17 C  argument de sortie du s-pg                                    C
    18 C                                                                C
     14C                                                                C
     15C  sont les arguments d'entree pour le s-pg                      C
     16C                                                                C
     17C  argument de sortie du s-pg                                    C
     18C                                                                C
    1919CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    2020CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    … …  
    2626C
    2727C
    28 #include "dimensions.h"
    29 #include "paramet.h"
    30 #include "comvert.h"
    31 #include "comgeom2.h"
     28      include "dimensions.h"
     29      include "paramet.h"
     30      include "comgeom2.h"
    3231 
    3332C  Arguments :

  • LMDZ5/trunk/libf/dyn3d_common/advyp.F

    r2597 r2600  
    1111CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    1212C                                                                C
    13 C  Source : Pascal Simon ( Meteo, CNRM )                        C
    14 C  Adaptation : A.A. (LGGE)                                     C
     13C  Source : Pascal Simon ( Meteo, CNRM )                        C
     14C  Adaptation : A.A. (LGGE)                                     C
    1515C  Derniere Modif : 19/10/95 LAST
    16 C                                                                C
    17 C  sont les arguments d'entree pour le s-pg                      C
    18 C                                                                C
    19 C  argument de sortie du s-pg                                    C
    20 C                                                                C
     16C                                                                C
     17C  sont les arguments d'entree pour le s-pg                      C
     18C                                                                C
     19C  argument de sortie du s-pg                                    C
     20C                                                                C
    2121CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    2222CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    … …  
    2828C
    2929C
    30 #include "dimensions.h"
    31 #include "paramet.h"
    32 #include "comvert.h"
    33 #include "comgeom.h"
     30      include "dimensions.h"
     31      include "paramet.h"
     32      include "comgeom.h"
    3433 
    3534C  Arguments :

  • LMDZ5/trunk/libf/dyn3d_common/advz.F

    r2597 r2600  
    1616C                                                                C
    1717C  dq est l'argument de sortie pour le s-pg                      C
    18 C                                                                C
     18C                                                                C
    1919CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    2020C
    2121C  parametres principaux du modele
    2222C
    23 #include "dimensions.h"
    24 #include "paramet.h"
    25 #include "comvert.h"
     23      include "dimensions.h"
     24      include "paramet.h"
    2625
    2726C    #include "traceur.h"

  • LMDZ5/trunk/libf/dyn3d_common/advzp.F

    r2597 r2600  
    3131C  parametres principaux du modele
    3232C
    33 #include "dimensions.h"
    34 #include "paramet.h"
    35 #include "comvert.h"
    36 #include "comgeom.h"
     33      include "dimensions.h"
     34      include "paramet.h"
     35      include "comgeom.h"
    3736C
    3837C  Arguments :
    … …  
    355354      DO j = 1,jjp1
    356355          SM(iip1,j,l) = SM(1,j,l)
    357           S0(iip1,j,l,ntra) = S0(1,j,l,ntra)
     356          S0(iip1,j,l,ntra) = S0(1,j,l,ntra)
    358357          SSX(iip1,j,l,ntra) = SSX(1,j,l,ntra)
    359           SY(iip1,j,l,ntra) = SY(1,j,l,ntra)
     358          SY(iip1,j,l,ntra) = SY(1,j,l,ntra)
    360359          SZ(iip1,j,l,ntra) = SZ(1,j,l,ntra)
    361360      ENDDO
    362361      ENDDO
    363 c                                                                               C-------------------------------------------------------------
     362c                                                                                C-------------------------------------------------------------
    364363C *** Test : diag de la qqtite totale de tarceur
    365364C            dans l'atmosphere avant l'advection en z

  • LMDZ5/trunk/libf/dyn3d_common/caldyn0.F90

    r2597 r2600  
    77!-------------------------------------------------------------------------------
    88  USE control_mod, ONLY: resetvarc
     9  USE comvert_mod, ONLY: ap, bp
    910  IMPLICIT NONE
    1011  include "dimensions.h"
    1112  include "paramet.h"
    12   include "comvert.h"
    1313  include "comgeom.h"
    1414!===============================================================================

  • LMDZ5/trunk/libf/dyn3d_common/comvert_mod.F90

    r2599 r2600  
    22! $Id$
    33!
    4 !-----------------------------------------------------------------------
    5 !   INCLUDE 'comvert.h'
     4MODULE comvert_mod
    65
    7       COMMON/comvertr/ap(llm+1),bp(llm+1),presnivs(llm),dpres(llm),     &
    8      &               pa,preff,nivsigs(llm),nivsig(llm+1),               &
    9      &               aps(llm),bps(llm),scaleheight,pseudoalt(llm)
     6IMPLICIT NONE 
    107
    11       common/comverti/disvert_type, pressure_exner
     8INCLUDE "dimensions.h"
    129
    13       real ap     ! hybrid pressure contribution at interlayers
    14       real bp     ! hybrid sigma contribution at interlayer
    15       real presnivs ! (reference) pressure at mid-layers
    16       real dpres
    17       real pa     ! reference pressure (Pa) at which hybrid coordinates
    18                   ! become purely pressure
    19       real preff  ! reference surface pressure (Pa)
    20       real nivsigs
    21       real nivsig
    22       real aps    ! hybrid pressure contribution at mid-layers
    23       real bps    ! hybrid sigma contribution at mid-layers
    24       real scaleheight ! atmospheric (reference) scale height (km)
    25       real pseudoalt ! pseudo-altitude of model levels (km), based on presnivs(),
     10REAL ap(llm+1) ! hybrid pressure contribution at interlayers
     11REAL bp (llm+1) ! hybrid sigma contribution at interlayer
     12REAL presnivs(llm) ! (reference) pressure at mid-layers
     13REAL dpres(llm)
     14REAL sig(llm+1)
     15REAL ds(llm)
     16REAL pa ! reference pressure (Pa) at which hybrid coordinates
     17        ! become purely pressure (more or less)
     18REAL preff  ! reference surface pressure (Pa)
     19REAL nivsigs(llm)
     20REAL nivsig(llm+1)
     21REAL aps(llm) ! hybrid pressure contribution at mid-layers
     22REAL bps(llm) ! hybrid sigma contribution at mid-layers
     23REAL scaleheight ! atmospheric (reference) scale height (km)
     24REAL pseudoalt(llm) ! pseudo-altitude of model levels (km), based on presnivs(),
    2625                     ! preff and scaleheight
    2726
    28       integer disvert_type ! type of vertical discretization:
    29                            ! 1: Earth (default for planet_type==earth),
    30                            !     automatic generation
    31                            ! 2: Planets (default for planet_type!=earth),
    32                            !     using 'z2sig.def' (or 'esasig.def) file
     27INTEGER disvert_type ! type of vertical discretization:
     28                     ! 1: Earth (default for planet_type==earth),
     29                     !     automatic generation
     30                     ! 2: Planets (default for planet_type!=earth),
     31                     !     using 'z2sig.def' (or 'esasig.def) file
    3332
    34       logical pressure_exner
     33LOGICAL pressure_exner
    3534!     compute pressure inside layers using Exner function, else use mean
    3635!     of pressure values at interfaces
    3736
    38  !-----------------------------------------------------------------------
     37END MODULE comvert_mod

  • LMDZ5/trunk/libf/dyn3d_common/conf_planete.F90

    r2597 r2600  
    1212USE comconst_mod, ONLY: pi, g, molmass, kappa, cpp, omeg, rad, &
    1313                        year_day, daylen, daysec, ihf
     14USE comvert_mod, ONLY: preff, pa
    1415IMPLICIT NONE
    1516!
    … …  
    1718!   Declarations :
    1819!   --------------
    19 #include "dimensions.h"
    20 #include "comvert.h"
     20
    2121!
    2222!   local:

  • LMDZ5/trunk/libf/dyn3d_common/disvert.F90

    r2597 r2600  
    1010  use new_unit_m, only: new_unit
    1111  use assert_m, only: assert
     12  USE comvert_mod, ONLY: ap, bp, nivsigs, nivsig, dpres, presnivs, &
     13                         pa, preff, scaleheight
    1214
    1315  IMPLICIT NONE
    … …  
    1517  include "dimensions.h"
    1618  include "paramet.h"
    17   include "comvert.h"
    1819  include "iniprint.h"
    1920  include "logic.h"
    … …  
    2324!          Triggered by the levels number llm.
    2425!-------------------------------------------------------------------------------
    25 ! Read    in "comvert.h":
     26! Read    in "comvert_mod":
    2627
    2728! pa !--- vertical coordinate is close to a PRESSURE COORDINATE FOR P
    … …  
    2930
    3031! preff                      !--- REFERENCE PRESSURE                 (101325 Pa)
    31 ! Written in "comvert.h":
     32! Written in "comvert_mod":
    3233! ap(llm+1), bp(llm+1)       !--- Ap, Bp HYBRID COEFFICIENTS AT INTERFACES
    3334! aps(llm),  bps(llm)        !--- Ap, Bp HYBRID COEFFICIENTS AT MID-LAYERS

  • LMDZ5/trunk/libf/dyn3d_common/disvert_noterre.F

    r2597 r2600  
    1212      use ioipsl_getincom
    1313#endif
     14      USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,
     15     &                       nivsig,nivsigs,pa,preff,scaleheight
    1416      USE comconst_mod, ONLY: kappa
    1517
    1618      IMPLICIT NONE
    1719
    18 #include "dimensions.h"
    19 #include "paramet.h"
    20 #include "comvert.h"
    21 #include "logic.h"
    22 #include "iniprint.h"
     20      include "dimensions.h"
     21      include "paramet.h"
     22      include "logic.h"
     23      include "iniprint.h"
    2324c
    2425c=======================================================================

  • LMDZ5/trunk/libf/dyn3d_common/exner_hyb_m.F90

    r2598 r2600  
    3434    !
    3535    USE comconst_mod, ONLY: jmp1, cpp, kappa, r
     36    USE comvert_mod, ONLY: preff
     37   
    3638    IMPLICIT NONE
    3739   
    … …  
    3941    include "paramet.h"
    4042    include "comgeom.h"
    41     include "comvert.h"
    4243
    4344    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3d_common/exner_milieu_m.F90

    r2598 r2600  
    3131    !
    3232    USE comconst_mod, ONLY: jmp1, cpp, kappa, r
     33    USE comvert_mod, ONLY: preff
     34   
    3335    IMPLICIT NONE
    3436   
    … …  
    3638    include "paramet.h"
    3739    include "comgeom.h"
    38     include "comvert.h"
    3940
    4041    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3d_common/geopot.F

    r1945 r2600  
    2727c   -------------
    2828
    29 #include "dimensions.h"
    30 #include "paramet.h"
    31 #include "comvert.h"
     29      include "dimensions.h"
     30      include "paramet.h"
    3231
    3332c   Arguments:

  • LMDZ5/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90

    r2598 r2600  
    1010
    1111  USE comconst_mod, ONLY: cpp, kappa, g, omeg, daysec, rad, pi
     12  USE comvert_mod, ONLY: presnivs, preff, pa
     13 
    1214  IMPLICIT NONE
    1315
    … …  
    1517  INCLUDE "paramet.h"
    1618  INCLUDE "comgeom.h"
    17   INCLUDE "comvert.h"
    1819  INCLUDE "netcdf.inc"
    1920

  • LMDZ5/trunk/libf/dyn3d_common/iniconst.F90

    r2597 r2600  
    1313  USE comconst_mod, ONLY: im, imp1, jm, jmp1, lllm, lllmm1, lllmp1, &
    1414                          unsim, pi, r, kappa, cpp, dtvr, dtphys
     15  USE comvert_mod, ONLY: disvert_type, pressure_exner
    1516 
    1617  IMPLICIT NONE
    … …  
    2425  include "paramet.h"
    2526  include "temps.h"
    26   include "comvert.h"
    2727  include "iniprint.h"
    2828

  • LMDZ5/trunk/libf/dyn3d_common/inidissip.F90

    r2597 r2600  
    1414  USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
    1515                          dtdiss, dtvr, rad
     16  USE comvert_mod, ONLY: preff, presnivs
    1617
    1718  IMPLICIT NONE
    … …  
    1920  include "paramet.h"
    2021  include "comdissipn.h"
    21   include "comvert.h"
    2222  include "logic.h"
    2323  include "iniprint.h"

  • LMDZ5/trunk/libf/dyn3d_common/initdynav.F90

    r2598 r2600  
    1010       dynhistave_file,dynhistvave_file,dynhistuave_file
    1111  USE comconst_mod, ONLY: pi
     12  USE comvert_mod, ONLY: presnivs
    1213  implicit none
    1314
    … …  
    3536  include "dimensions.h"
    3637  include "paramet.h"
    37   include "comvert.h"
    3838  include "comgeom.h"
    3939  include "temps.h"

  • LMDZ5/trunk/libf/dyn3d_common/initfluxsto.F

    r2598 r2600  
    1010#endif
    1111      USE comconst_mod, ONLY: pi
     12      USE comvert_mod, ONLY: nivsigs
    1213      implicit none
    1314
    … …  
    4142      include "dimensions.h"
    4243      include "paramet.h"
    43       include "comvert.h"
    4444      include "comgeom.h"
    4545      include "temps.h"
    … …  
    8585      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    8686      tau0 = itau_dyn
    87        
    88         do jj = 1, jjp1
     87       
     88        do jj = 1, jjp1
    8989        do ii = 1, iip1
    9090          rlong(ii,jj) = rlonu(ii) * 180. / pi
    … …  
    112112     .             1, iip1, 1, jjm,
    113113     .             tau0, zjulian, tstep, vhoriid, filevid)
    114        
    115         rl(1,1) = 1.
     114       
     115        rl(1,1) = 1.
    116116      call histbeg('defstoke.nc', 1, rl, 1, rl,
    117117     .             1, 1, 1, 1,
    … …  
    130130      call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar',
    131131     .              'Grille points scalaires', thoriid)
    132        
     132       
    133133C
    134134C  Appel a histvert pour la grille verticale
    … …  
    149149C
    150150C  Appels a histdef pour la definition des variables a sauvegarder
    151        
    152         CALL histdef(fileid, "phis", "Surface geop. height", "-",
     151       
     152        CALL histdef(fileid, "phis", "Surface geop. height", "-",
    153153     .                iip1,jjp1,thoriid, 1,1,1, -99, 32,
    154154     .                "once", t_ops, t_wrt)
    … …  
    157157     .                iip1,jjp1,thoriid, 1,1,1, -99, 32,
    158158     .                "once", t_ops, t_wrt)
    159        
    160         CALL histdef(filedid, "dtvr", "tps dyn", "s",
     159       
     160        CALL histdef(filedid, "dtvr", "tps dyn", "s",
    161161     .                1,1,dhoriid, 1,1,1, -99, 32,
    162162     .                "once", t_ops, t_wrt)
    … …  
    222222        call histsync(filedid)
    223223      endif
    224        
     224       
    225225#else
    226226! tell the user this routine should be run with ioipsl

  • LMDZ5/trunk/libf/dyn3d_common/inithist.F

    r2598 r2600  
    1111     &                        dynhist_file,dynhistv_file,dynhistu_file
    1212       USE comconst_mod, ONLY: pi
     13       USE comvert_mod, ONLY: presnivs
    1314      implicit none
    1415
    … …  
    4041      include "dimensions.h"
    4142      include "paramet.h"
    42       include "comvert.h"
    4343      include "comgeom.h"
    4444      include "temps.h"

  • LMDZ5/trunk/libf/dyn3d_common/interpost.F

    r2597 r2600  
    99      include "dimensions.h"
    1010      include "paramet.h"
    11       include "comvert.h"
    1211      include "comgeom2.h"
    1312

  • LMDZ5/trunk/libf/dyn3d_common/interpre.F

    r2597 r2600  
    77
    88      USE comconst_mod, ONLY: g
     9      USE comvert_mod, ONLY: ap, bp
    910
    1011       implicit none
    … …  
    1314      include "paramet.h"
    1415      include "comdissip.h"
    15       include "comvert.h"
    1616      include "comgeom2.h"
    1717      include "logic.h"

  • LMDZ5/trunk/libf/dyn3d_common/limx.F

    r2597 r2600  
    1818      include "paramet.h"
    1919      include "logic.h"
    20       include "comvert.h"
    2120      include "comgeom.h"
    2221c

  • LMDZ5/trunk/libf/dyn3d_common/limy.F

    r2597 r2600  
    2020      include "paramet.h"
    2121      include "logic.h"
    22       include "comvert.h"
    2322      include "comgeom.h"
    2423c

  • LMDZ5/trunk/libf/dyn3d_common/limz.F

    r2597 r2600  
    1818      include "paramet.h"
    1919      include "logic.h"
    20       include "comvert.h"
    2120      include "comgeom.h"
    2221c

  • LMDZ5/trunk/libf/dyn3d_common/pentes_ini.F

    r2597 r2600  
    2727      include "dimensions.h"
    2828      include "paramet.h"
    29       include "comvert.h"
    3029      include "comgeom2.h"
    3130
    … …  
    439438
    440439        DO l = 1,llm
    441         DO j = 1,jjp1
    442           DO i = 1,iip1
     440            DO j = 1,jjp1
     441              DO i = 1,iip1
    443442                IF (q(i,j,l,0).lt.0.)  THEN
    444443c                    PRINT*,'------------ BIP-----------'
    … …  
    447446c                    PRINT*,'QY(',i,j,l,')=',q(i,j,l,2)
    448447c                    PRINT*,'QZ(',i,j,l,')=',q(i,j,l,3)
    449 c                            PRINT*,' PBL EN SORTIE DE PENTES'
     448c                            PRINT*,' PBL EN SORTIE DE PENTES'
    450449                     q(i,j,l,0)=0.
    451450c                    STOP

  • LMDZ5/trunk/libf/dyn3d_common/prather.F

    r2597 r2600  
    2424      include "dimensions.h"
    2525      include "paramet.h"
    26       include "comvert.h"
    2726      include "comgeom2.h"
    2827

  • LMDZ5/trunk/libf/dyn3d_common/sortvarc.F

    r2597 r2600  
    2828      INCLUDE "dimensions.h"
    2929      INCLUDE "paramet.h"
    30       INCLUDE "comvert.h"
    3130      INCLUDE "comgeom.h"
    3231      INCLUDE "ener.h"

  • LMDZ5/trunk/libf/dyn3d_common/traceurpole.F

    r2597 r2600  
    99      include "paramet.h"
    1010      include "comdissip.h"
    11       include "comvert.h"
    1211      include "comgeom2.h"
    1312      include "logic.h"

  • LMDZ5/trunk/libf/dyn3d_common/vitvert.F90

    r2336 r2600  
    55!-------------------------------------------------------------------------------
    66! Purpose: Compute vertical speed at sigma levels.
     7  USE comvert_mod, ONLY: bp
    78  IMPLICIT NONE
    89  include "dimensions.h"
    910  include "paramet.h"
    10   include "comvert.h"
    1111!===============================================================================
    1212! Arguments:

  • LMDZ5/trunk/libf/dyn3d_common/writedynav.F90

    r2598 r2600  
    3232  include "dimensions.h"
    3333  include "paramet.h"
    34   include "comvert.h"
    3534  include "comgeom.h"
    3635  include "temps.h"

  • LMDZ5/trunk/libf/dyn3d_common/writehist.F

    r2598 r2600  
    3535      include "dimensions.h"
    3636      include "paramet.h"
    37       include "comvert.h"
    3837      include "comgeom.h"
    3938      include "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/advect_new_loc.F

    r2597 r2600  
    3030      include "dimensions.h"
    3131      include "paramet.h"
    32       include "comvert.h"
    3332      include "comgeom.h"
    3433      include "logic.h"

  • LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F

    r2597 r2600  
    3232      include "dimensions.h"
    3333      include "paramet.h"
    34       include "comvert.h"
    3534      include "comdissip.h"
    3635      include "comgeom2.h"
    … …  
    172171#endif
    173172         
    174           GOTO 1234     
     173          GOTO 1234     
    175174c-----------------------------------------------------------
    176175c     Appel des sous programmes d'advection
    … …  
    184183        if(iadv(iq).eq.10) THEN
    185184     
    186 !LF         call vlsplt_p(q(1,1,iq),2.,massem,wg,pbarug,pbarvg,dtvr)
     185!LF            call vlsplt_p(q(1,1,iq),2.,massem,wg,pbarug,pbarvg,dtvr)
    187186
    188187c   ----------------------------------------------------------------

  • LMDZ5/trunk/libf/dyn3dmem/bilan_dyn_loc.F

    r2597 r2600  
    1818      USE write_field_loc
    1919      USE comconst_mod, ONLY: cpp, pi
     20      USE comvert_mod, ONLY: presnivs
     21     
    2022      IMPLICIT NONE
    2123
    2224      include "dimensions.h"
    2325      include "paramet.h"
    24       include "comvert.h"
    2526      include "comgeom2.h"
    2627      include "temps.h"
    … …  
    598599        do l=1,llm
    599600          Q_cum(:,jjb:jje,l,iQ)=Q_cum(:,jjb:jje,l,iQ)
    600      .                          /masse_cum(:,jjb:jje,l)
     601     .                                /masse_cum(:,jjb:jje,l)
    601602        enddo
    602603!$OMP ENDDO NOWAIT

  • LMDZ5/trunk/libf/dyn3dmem/caldyn_loc.F

    r2597 r2600  
    1212      USE caldyn_mod, ONLY: vcont, ucont, ang, p, massebx, masseby,
    1313     &                      vorpot, ecin, bern, massebxy, convm
     14      USE comvert_mod, ONLY: ap, bp
    1415     
    1516      IMPLICIT NONE
    … …  
    3334      include "dimensions.h"
    3435      include "paramet.h"
    35       include "comvert.h"
    3636      include "comgeom.h"
    3737

  • LMDZ5/trunk/libf/dyn3dmem/call_calfis_mod.F90

    r2597 r2600  
    8383  USE write_field
    8484  USE comconst_mod, ONLY: dtphys
     85  USE comvert_mod, ONLY: ap, bp, pressure_exner
     86 
    8587  IMPLICIT NONE
    86     INCLUDE "comvert.h"
    8788    INCLUDE "logic.h"
    8889    INCLUDE "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/dudv2_loc.F

    r1907 r2600  
    2323c=======================================================================
    2424c
    25 #include "dimensions.h"
    26 #include "paramet.h"
    27 #include "comvert.h"
     25      include "dimensions.h"
     26      include "paramet.h"
    2827
    2928      REAL teta( ijb_u:ije_u,llm ),pkf( ijb_u:ije_u,llm )

  • LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.f90

    r2598 r2600  
    1212  USE control_mod, ONLY: planet_type
    1313  USE assert_eq_m, ONLY: assert_eq
     14  USE comvert_mod, ONLY: pa,preff
    1415  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa, lllm, &
    1516                          omeg, rad
    … …  
    2021  include "paramet.h"
    2122  include "temps.h"
    22   include "comvert.h"
    2323  include "comgeom.h"
    2424  include "ener.h"

  • LMDZ5/trunk/libf/dyn3dmem/dynredem_loc.F90

    r2598 r2600  
    1313                    NF90_CLOSE,  NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER
    1414  USE dynredem_mod, ONLY: cre_var, put_var, err, modname, fil
     15  USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff, &
     16                         nivsig,nivsigs
    1517  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    1618  USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, &
    … …  
    2022  include "dimensions.h"
    2123  include "paramet.h"
    22   include "comvert.h"
    2324  include "comgeom.h"
    2425  include "temps.h"
    … …  
    182183  include "paramet.h"
    183184  include "description.h"
    184   include "comvert.h"
    185185  include "comgeom.h"
    186186  include "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/exner_hyb_loc_m.F90

    r2598 r2600  
    3636    USE write_field_loc
    3737    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     38    USE comvert_mod, ONLY: preff
     39   
    3840    IMPLICIT NONE
    3941    !
    … …  
    4143    include "paramet.h"
    4244    include "comgeom.h"
    43     include "comvert.h"
    4445
    4546    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3dmem/exner_milieu_loc_m.F90

    r2598 r2600  
    3232    USE mod_filtreg_p
    3333    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     34    USE comvert_mod, ONLY: preff
     35   
    3436    IMPLICIT NONE
    3537    !
    … …  
    3739    include "paramet.h"
    3840    include "comgeom.h"
    39     include "comvert.h"
    4041
    4142    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3dmem/gcm.F90

    r2598 r2600  
    5656  include "paramet.h"
    5757  include "comdissnew.h"
    58   include "comvert.h"
    5958  include "comgeom.h"
    6059  include "logic.h"

  • LMDZ5/trunk/libf/dyn3dmem/geopot_loc.F

    r1907 r2600  
    2626c   -------------
    2727
    28 #include "dimensions.h"
    29 #include "paramet.h"
    30 #include "comvert.h"
     28      include "dimensions.h"
     29      include "paramet.h"
    3130
    3231c   Arguments:

  • LMDZ5/trunk/libf/dyn3dmem/groupe_loc.F

    r2597 r2600  
    2020      include "paramet.h"
    2121      include "comgeom2.h"
    22       include "comvert.h"
    2322
    2423!     integer ngroup

  • LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

    r2598 r2600  
    345345    USE write_field_loc
    346346    USE comconst_mod, ONLY: cpp, daysec, dtvr, kappa
     347    USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
    347348   
    348349    IMPLICIT NONE
    … …  
    350351    INCLUDE "dimensions.h"
    351352    INCLUDE "paramet.h"
    352     INCLUDE "comvert.h"
    353353
    354354    ! Variables entree
    … …  
    891891  USE Bands
    892892  USE comconst_mod, ONLY: cpp, kappa
     893  USE comvert_mod, ONLY: preff, pressure_exner, bp, ap, disvert_type
    893894  IMPLICIT NONE
    894895
    895896  include "dimensions.h"
    896897  include "paramet.h"
    897   include "comvert.h"
    898898  include "comgeom2.h"
    899899
    … …  
    21292129    USE mod_hallo, ONLY : gather_field_u, gather_field_v
    21302130    USE comconst_mod, ONLY: pi
     2131    USE comvert_mod, ONLY: presnivs
    21312132    IMPLICIT NONE
    21322133
    … …  
    21352136    INCLUDE "netcdf.inc"
    21362137    INCLUDE "comgeom2.h"
    2137     INCLUDE "comvert.h"
    21382138   
    21392139    ! Variables entree

  • LMDZ5/trunk/libf/dyn3dmem/iniacademic_loc.F90

    r2597 r2600  
    1919  USE Write_Field
    2020  USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
     21  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
    2122
    2223  !   Author:    Frederic Hourdin      original: 15/01/93
    … …  
    3132  include "dimensions.h"
    3233  include "paramet.h"
    33   include "comvert.h"
    3434  include "comgeom.h"
    3535  include "academic.h"

  • LMDZ5/trunk/libf/dyn3dmem/initdynav_loc.F

    r2598 r2600  
    1515     &        dynhistave_file,dynhistvave_file,dynhistuave_file
    1616       USE comconst_mod, ONLY: pi
     17       USE comvert_mod, ONLY: presnivs
     18       
    1719       implicit none
    1820
    … …  
    4446      include "dimensions.h"
    4547      include "paramet.h"
    46       include "comvert.h"
    4748      include "comgeom.h"
    4849      include "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/initfluxsto_p.F

    r2598 r2600  
    1414       use misc_mod
    1515       USE comconst_mod, ONLY: pi
     16       USE comvert_mod, ONLY: nivsigs
    1617       
    1718      implicit none
    … …  
    4647      include "dimensions.h"
    4748      include "paramet.h"
    48       include "comvert.h"
    4949      include "comgeom.h"
    5050      include "temps.h"
    … …  
    105105      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    106106      tau0 = itau_dyn
    107        
    108         do jj = 1, jjp1
     107       
     108        do jj = 1, jjp1
    109109        do ii = 1, iip1
    110110          rlong(ii,jj) = rlonu(ii) * 180. / pi
    … …  
    164164     .             1, iip1, 1, jjn,tau0, zjulian, tstep, vhoriid,
    165165     .             filevid,dynv_domain_id)
    166        
     166       
    167167      rl(1,1) = 1.
    168168     
    … …  
    190190      call histhori(fileid, iip1, rlong(:,jjb:jje),jjn,rlat(:,jjb:jje),
    191191     .             'scalar','Grille points scalaires', thoriid)
    192        
     192       
    193193C
    194194C  Appel a histvert pour la grille verticale
    … …  
    210210C
    211211C  Appels a histdef pour la definition des variables a sauvegarder
    212        
    213         CALL histdef(fileid, "phis", "Surface geop. height", "-",
     212       
     213        CALL histdef(fileid, "phis", "Surface geop. height", "-",
    214214     .                iip1,jjn,thoriid, 1,1,1, -99, 32,
    215215     .                "once", t_ops, t_wrt)
    … …  
    218218     .                iip1,jjn,thoriid, 1,1,1, -99, 32,
    219219     .                "once", t_ops, t_wrt)
    220        
     220       
    221221        if (mpi_rank==0) then
    222        
    223         CALL histdef(filedid, "dtvr", "tps dyn", "s",
     222       
     223        CALL histdef(filedid, "dtvr", "tps dyn", "s",
    224224     .                1,1,dhoriid, 1,1,1, -99, 32,
    225225     .                "once", t_ops, t_wrt)
    … …  
    288288        if (mpi_rank==0) call histsync(filedid)
    289289      endif
    290        
     290       
    291291#else
    292292      write(lunout,*)'initfluxsto_p: Needs IOIPSL to function'

  • LMDZ5/trunk/libf/dyn3dmem/inithist_loc.F

    r2598 r2600  
    1515     &                        dynhist_file,dynhistv_file,dynhistu_file
    1616       USE comconst_mod, ONLY: pi
     17       USE comvert_mod, ONLY: presnivs
     18       
    1719       implicit none
    1820
    … …  
    4345      include "dimensions.h"
    4446      include "paramet.h"
    45       include "comvert.h"
    4647      include "comgeom.h"
    4748      include "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F

    r2598 r2600  
    1313      USE infotrac, ONLY: ok_iso_verif ! ajout CRisi
    1414      USE comconst_mod, ONLY: pi
     15      USE comvert_mod, ONLY: ap, bp
     16     
    1517      IMPLICIT NONE
    1618
    … …  
    3436      include "paramet.h"
    3537      include "comgeom.h"
    36       include "comvert.h"
    3738      include "logic.h"
    3839      include "temps.h"

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

    r2598 r2600  
    3333       use exner_milieu_loc_m, only: exner_milieu_loc
    3434       USE comconst_mod, ONLY: cpp, dtvr, ihf
     35       USE comvert_mod, ONLY: ap, bp, pressure_exner
     36       
    3537      IMPLICIT NONE
    3638
    … …  
    6971      include "paramet.h"
    7072      include "comdissnew.h"
    71       include "comvert.h"
    7273      include "comgeom.h"
    7374      include "logic.h"
    … …  
    16321633            IF( MOD(itau,iecri).EQ.0) THEN
    16331634             ! Ehouarn: output only during LF or Backward Matsuno
    1634              if (leapf.or.(.not.leapf.and.(.not.forward))) then
     1635             if (leapf.or.(.not.leapf.and.(.not.forward))) then
    16351636
    16361637c$OMP BARRIER
    … …  
    16421643#ifdef CPP_IOIPSL
    16431644             if (ok_dyn_ins) then
    1644                  CALL writehist_loc(itau,vcov,ucov,teta,pk,phi,q,
     1645                 CALL writehist_loc(itau,vcov,ucov,teta,pk,phi,q,
    16451646     &                              masse,ps,phis)
    16461647             endif

  • LMDZ5/trunk/libf/dyn3dmem/qminimum_loc.F

    r2286 r2600  
    77c             pour l'eau vapeur et l'eau liquide
    88c
    9 #include "dimensions.h"
    10 #include "paramet.h"
    11 #include "comvert.h"
     9      include "dimensions.h"
     10      include "paramet.h"
    1211c
    1312      INTEGER nqtot ! CRisi: on remplace nq par nqtot

  • LMDZ5/trunk/libf/dyn3dmem/sw_case_williamson91_6_loc.F

    r2597 r2600  
    2828      USE parallel_lmdz
    2929      USE comconst_mod, ONLY: cpp, omeg, rad
     30      USE comvert_mod, ONLY: ap, bp, preff
    3031
    3132      IMPLICIT NONE
    … …  
    3637      include "dimensions.h"
    3738      include "paramet.h"
    38       include "comvert.h"
    3939      include "comgeom.h"
    4040      include "iniprint.h"

  • LMDZ5/trunk/libf/dyn3dmem/top_bound_loc.F

    r2597 r2600  
    66      USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound,
    77     &                        tau_top_bound
     8      USE comvert_mod, ONLY: presnivs, preff, scaleheight
     9
    810      IMPLICIT NONE
    911c
    1012      include "dimensions.h"
    1113      include "paramet.h"
    12       include "comvert.h"
    1314      include "comgeom2.h"
    1415

  • LMDZ5/trunk/libf/dyn3dmem/vitvert_loc.F90

    r2336 r2600  
    66! Purpose: Compute vertical speed at sigma levels.
    77  USE parallel_lmdz
     8  USE comvert_mod, ONLY: bp
     9 
    810  IMPLICIT NONE
    911  include "dimensions.h"
    1012  include "paramet.h"
    11   include "comvert.h"
    1213!===============================================================================
    1314! Arguments:

  • LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F

    r2597 r2600  
    2020      include "paramet.h"
    2121      include "logic.h"
    22       include "comvert.h"
    2322c
    2423c
    … …  
    422421      include "paramet.h"
    423422      include "logic.h"
    424       include "comvert.h"
    425423      include "comgeom.h"
    426424c
    … …  
    877875      include "paramet.h"
    878876      include "logic.h"
    879       include "comvert.h"
    880877c
    881878c

  • LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F

    r2597 r2600  
    3838      include "paramet.h"
    3939      include "logic.h"
    40       include "comvert.h"
    4140
    4241c

  • LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F

    r2597 r2600  
    1515      include "paramet.h"
    1616      include "logic.h"
    17       include "comvert.h"
    1817c
    1918c
    … …  
    173172c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    174173        DO l=1,llm
    175           iadvplus(ip1jm+1:ip1jmp1,l)=0
     174          iadvplus(ip1jm+1:ip1jmp1,l)=0
    176175        ENDDO
    177176c$OMP END DO NOWAIT
    178177      endif
    179        
     178             
    180179c   calcul des flux a gauche et a droite
    181180
    … …  
    413412c    ********************************************************************
    414413c     q,masse_adv_v,w sont des arguments d'entree  pour le s-pg ....
    415 c     qsat             est   un argument de sortie pour le s-pg ....
     414c     qsat                est   un argument de sortie pour le s-pg ....
    416415c
    417416c
    … …  
    425424      include "paramet.h"
    426425      include "logic.h"
    427       include "comvert.h"
    428426      include "comgeom.h"
    429427c
    … …  
    581579         dyq(ij,l)=fn*dyq(ij,l)
    582580        ENDDO
    583          
     581         
    584582      ENDIF
    585583     
    … …  
    602600        ENDDO
    603601       
    604 c   calcul des pentes limites aux poles 
     602c   calcul des pentes limites aux poles       
    605603        fs=1.
    606604        DO ij=1,iim
    … …  
    613611         dyq(ip1jm+ij,l)=fs*dyq(ip1jm+ij,l)
    614612        ENDDO
    615        
     613       
    616614      ENDIF
    617615
    … …  
    778776              masse(ij,l,iq)=newmasse
    779777           ENDDO
    780          
    781          ENDIF
    782778         
    783          IF (pole_sud) THEN
    784          
    785            convps  = -SSUM(iim,qbyv(ip1jm-iim,l,iq),iq,1)/apols
     779         ENDIF
     780         
     781         IF (pole_sud) THEN
     782         
     783           convps  = -SSUM(iim,qbyv(ip1jm-iim,l,iq),iq,1)/apols
    786784           convmps = -SSUM(iim,masse_adv_v(ip1jm-iim,l),1)/apols
    787785           DO ij = ip1jm+1,ip1jmp1
    … …  
    791789              masse(ij,l,iq)=newmasse
    792790           ENDDO
    793        
    794         ENDIF
     791       
     792        ENDIF
    795793c.-. fin ancienne version
    796794

  • LMDZ5/trunk/libf/dyn3dmem/writedynav_loc.F

    r2598 r2600  
    4444      include "dimensions.h"
    4545      include "paramet.h"
    46       include "comvert.h"
    4746      include "comgeom.h"
    4847      include "temps.h"

  • LMDZ5/trunk/libf/dyn3dmem/writehist_loc.F

    r2598 r2600  
    4444      include "dimensions.h"
    4545      include "paramet.h"
    46       include "comvert.h"
    4746      include "comgeom.h"
    4847      include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/advect_new_p.F

    r2597 r2600  
    2929#include "dimensions.h"
    3030#include "paramet.h"
    31 #include "comvert.h"
    3231#include "comgeom.h"
    3332#include "logic.h"

  • LMDZ5/trunk/libf/dyn3dpar/advect_p.F

    r2597 r2600  
    2828#include "dimensions.h"
    2929#include "paramet.h"
    30 #include "comvert.h"
    3130#include "comgeom.h"
    3231#include "logic.h"
    … …  
    8685                 
    8786           uav(ij,l)=0.25*(ucov(ij,l)+ucov(ij-iip1,l))
    88      .               +0.25*(ucov(ij+iip1,l)+ucov(ij,l))
     87     .                     +0.25*(ucov(ij+iip1,l)+ucov(ij,l))
    8988         ENDDO
    9089         

  • LMDZ5/trunk/libf/dyn3dpar/advtrac_p.F90

    r2597 r2600  
    2424  include "dimensions.h"
    2525  include "paramet.h"
    26   include "comvert.h"
    2726  include "comdissip.h"
    2827  include "comgeom2.h"

  • LMDZ5/trunk/libf/dyn3dpar/bilan_dyn_p.F

    r2597 r2600  
    1818      use write_field_p
    1919      USE comconst_mod, ONLY: cpp, pi
     20      USE comvert_mod, ONLY: presnivs
     21     
    2022      IMPLICIT NONE
    2123
    2224#include "dimensions.h"
    2325#include "paramet.h"
    24 #include "comvert.h"
    2526#include "comgeom2.h"
    2627#include "temps.h"
    … …  
    202203
    203204!        if (i_sortie.eq.1) then
    204 !        file='dynzon'
     205!        file='dynzon'
    205206!         if (mpi_rank==0) then
    206 !        call inigrads(ifile,1
     207!        call inigrads(ifile,1
    207208!     s  ,0.,180./pi,0.,0.,jjm,rlatv,-90.,90.,180./pi
    208209!     s  ,llm,presnivs,1.
    … …  
    568569        do l=1,llm
    569570          Q_cum(:,jjb:jje,l,iQ)=Q_cum(:,jjb:jje,l,iQ)
    570      .                          /masse_cum(:,jjb:jje,l)
     571     .                                /masse_cum(:,jjb:jje,l)
    571572        enddo
    572573!$OMP ENDDO NOWAIT

  • LMDZ5/trunk/libf/dyn3dpar/caldyn_p.F

    r2597 r2600  
    1010      USE parallel_lmdz
    1111      USE Write_Field_p
     12      USE comvert_mod, ONLY: ap, bp
    1213     
    1314      IMPLICIT NONE
    … …  
    3132#include "dimensions.h"
    3233#include "paramet.h"
    33 #include "comvert.h"
    3434#include "comgeom.h"
    3535

  • LMDZ5/trunk/libf/dyn3dpar/convmas1_p.F

    r1907 r2600  
    3232#include "dimensions.h"
    3333#include "paramet.h"
    34 #include "comvert.h"
    3534#include "logic.h"
    3635

  • LMDZ5/trunk/libf/dyn3dpar/convmas2_p.F

    r1907 r2600  
    3232#include "dimensions.h"
    3333#include "paramet.h"
    34 #include "comvert.h"
    3534#include "logic.h"
    3635

  • LMDZ5/trunk/libf/dyn3dpar/convmas_p.F

    r1907 r2600  
    3232#include "dimensions.h"
    3333#include "paramet.h"
    34 #include "comvert.h"
    3534#include "logic.h"
    3635

  • LMDZ5/trunk/libf/dyn3dpar/dudv2_p.F

    r1907 r2600  
    2525#include "dimensions.h"
    2626#include "paramet.h"
    27 #include "comvert.h"
    2827
    2928      REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ),

  • LMDZ5/trunk/libf/dyn3dpar/dynetat0.F

    r2598 r2600  
    99
    1010      use control_mod, only : planet_type
     11      USE comvert_mod, ONLY: pa,preff
    1112      USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa,
    1213     &                        lllm, omeg, rad
     14
    1315      USE serre_mod, ONLY: clon,clat,grossismx,grossismy
    1416
    … …  
    3335#include "paramet.h"
    3436#include "temps.h"
    35 #include "comvert.h"
    3637#include "comgeom2.h"
    3738#include "ener.h"

  • LMDZ5/trunk/libf/dyn3dpar/dynredem.F

    r2598 r2600  
    99      USE infotrac
    1010      use netcdf95, only: NF95_PUT_VAR
     11      USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff,
     12     &                        nivsig,nivsigs
    1113      USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    1214      USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy,
    … …  
    2123#include "dimensions.h"
    2224#include "paramet.h"
    23 #include "comvert.h"
    2425#include "comgeom2.h"
    2526#include "temps.h"
    … …  
    481482#include "description.h"
    482483#include "netcdf.inc"
    483 #include "comvert.h"
    484484#include "comgeom.h"
    485485#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/dynredem_p.F

    r2598 r2600  
    1010      USE infotrac
    1111      use netcdf95, only: NF95_PUT_VAR
     12      USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff,
     13     &                        nivsig,nivsigs
    1214      USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    1315      USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy,
    … …  
    2224#include "dimensions.h"
    2325#include "paramet.h"
    24 #include "comvert.h"
    2526#include "comgeom2.h"
    2627#include "temps.h"
    … …  
    482483#include "description.h"
    483484#include "netcdf.inc"
    484 #include "comvert.h"
    485485#include "comgeom.h"
    486486#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/exner_hyb_p_m.F90

    r2598 r2600  
    3434    USE parallel_lmdz
    3535    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     36    USE comvert_mod, ONLY: preff
    3637    !
    3738    include "dimensions.h"
    3839    include "paramet.h"
    3940    include "comgeom.h"
    40     include "comvert.h"
    4141
    4242    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3dpar/exner_milieu_p_m.F90

    r2598 r2600  
    3131    USE parallel_lmdz
    3232    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     33    USE comvert_mod, ONLY: preff
    3334    !
    3435    include "dimensions.h"
    3536    include "paramet.h"
    3637    include "comgeom.h"
    37     include "comvert.h"
    3838
    3939    INTEGER  ngrid

  • LMDZ5/trunk/libf/dyn3dpar/fluxstokenc_p.F

    r2597 r2600  
    2121#include "dimensions.h"
    2222#include "paramet.h"
    23 #include "comvert.h"
    2423#include "comgeom.h"
    2524#include "tracstoke.h"

  • LMDZ5/trunk/libf/dyn3dpar/gcm.F

    r2598 r2600  
    6868#include "paramet.h"
    6969#include "comdissnew.h"
    70 #include "comvert.h"
    7170#include "comgeom.h"
    7271#include "logic.h"

  • LMDZ5/trunk/libf/dyn3dpar/geopot_p.F

    r1907 r2600  
    2828#include "dimensions.h"
    2929#include "paramet.h"
    30 #include "comvert.h"
    3130
    3231c   Arguments:

  • LMDZ5/trunk/libf/dyn3dpar/groupe_p.F

    r2597 r2600  
    1818#include "paramet.h"
    1919#include "comgeom2.h"
    20 #include "comvert.h"
    2120
    2221!     integer ngroup

  • LMDZ5/trunk/libf/dyn3dpar/guide_p_mod.F90

    r2598 r2600  
    340340    USE control_mod
    341341    USE comconst_mod, ONLY: daysec, dtvr, cpp, kappa
     342    USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
    342343   
    343344    IMPLICIT NONE
    … …  
    345346    INCLUDE "dimensions.h"
    346347    INCLUDE "paramet.h"
    347     INCLUDE "comvert.h"
    348348
    349349    ! Variables entree
    … …  
    707707  USE Bands
    708708  USE comconst_mod, ONLY: cpp, kappa
     709  USE comvert_mod, ONLY: preff, pressure_exner, bp, ap
    709710  IMPLICIT NONE
    710711
    711712  include "dimensions.h"
    712713  include "paramet.h"
    713   include "comvert.h"
    714714  include "comgeom2.h"
    715715
    … …  
    18141814    USE parallel_lmdz
    18151815    USE comconst_mod, ONLY: pi
     1816    USE comvert_mod, ONLY: presnivs
    18161817    IMPLICIT NONE
    18171818
    … …  
    18201821    INCLUDE "netcdf.inc"
    18211822    INCLUDE "comgeom2.h"
    1822     INCLUDE "comvert.h"
    18231823   
    18241824    ! Variables entree

  • LMDZ5/trunk/libf/dyn3dpar/iniacademic.F90

    r2597 r2600  
    1717  use exner_milieu_m, only: exner_milieu
    1818  USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
     19  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
    1920
    2021  !   Author:    Frederic Hourdin      original: 15/01/93
    … …  
    2930  include "dimensions.h"
    3031  include "paramet.h"
    31   include "comvert.h"
    3232  include "comgeom.h"
    3333  include "academic.h"

  • LMDZ5/trunk/libf/dyn3dpar/initdynav_p.F

    r2598 r2600  
    1313       USE infotrac
    1414       USE comconst_mod, ONLY: pi
     15       USE comvert_mod, ONLY: nivsigs
    1516
    1617      implicit none
    … …  
    4445#include "dimensions.h"
    4546#include "paramet.h"
    46 #include "comvert.h"
    4747#include "comgeom.h"
    4848#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/initfluxsto_p.F

    r2598 r2600  
    1414       use misc_mod
    1515       USE comconst_mod, ONLY: pi
     16       USE comvert_mod, ONLY: nivsigs
    1617       
    1718      implicit none
    … …  
    4647#include "dimensions.h"
    4748#include "paramet.h"
    48 #include "comvert.h"
    4949#include "comgeom.h"
    5050#include "temps.h"
    … …  
    105105      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    106106      tau0 = itau_dyn
    107        
    108         do jj = 1, jjp1
     107       
     108        do jj = 1, jjp1
    109109        do ii = 1, iip1
    110110          rlong(ii,jj) = rlonu(ii) * 180. / pi
    … …  
    164164     .             1, iip1, 1, jjn,tau0, zjulian, tstep, vhoriid,
    165165     .             filevid,dynv_domain_id)
    166        
    167       rl(1,1) = 1.     
     166       
     167      rl(1,1) = 1.       
    168168     
    169169      if (mpi_rank==0) then
    … …  
    190190      call histhori(fileid, iip1, rlong(:,jjb:jje),jjn,rlat(:,jjb:jje),
    191191     .             'scalar','Grille points scalaires', thoriid)
    192        
     192       
    193193C
    194194C  Appel a histvert pour la grille verticale
    … …  
    210210C
    211211C  Appels a histdef pour la definition des variables a sauvegarder
    212        
    213         CALL histdef(fileid, "phis", "Surface geop. height", "-",
     212       
     213        CALL histdef(fileid, "phis", "Surface geop. height", "-",
    214214     .                iip1,jjn,thoriid, 1,1,1, -99, 32,
    215215     .                "once", t_ops, t_wrt)
    … …  
    218218     .                iip1,jjn,thoriid, 1,1,1, -99, 32,
    219219     .                "once", t_ops, t_wrt)
    220        
     220       
    221221        if (mpi_rank==0) then
    222        
    223         CALL histdef(filedid, "dtvr", "tps dyn", "s",
     222       
     223        CALL histdef(filedid, "dtvr", "tps dyn", "s",
    224224     .                1,1,dhoriid, 1,1,1, -99, 32,
    225225     .                "once", t_ops, t_wrt)
    … …  
    288288        if (mpi_rank==0) call histsync(filedid)
    289289      endif
    290        
     290       
    291291#else
    292292      write(lunout,*)'initfluxsto_p: Needs IOIPSL to function'

  • LMDZ5/trunk/libf/dyn3dpar/inithist_p.F

    r2598 r2600  
    1414       USE infotrac
    1515       USE comconst_mod, ONLY: pi
     16       USE comvert_mod, ONLY: nivsigs
    1617
    1718      implicit none
    … …  
    4647#include "dimensions.h"
    4748#include "paramet.h"
    48 #include "comvert.h"
    4949#include "comgeom.h"
    5050#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/integrd_p.F

    r2598 r2600  
    88      USE control_mod, only : planet_type
    99      USE comconst_mod, ONLY: pi
     10      USE comvert_mod, ONLY: ap, bp
    1011      IMPLICIT NONE
    1112
    … …  
    2930#include "paramet.h"
    3031#include "comgeom.h"
    31 #include "comvert.h"
    3232#include "logic.h"
    3333#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F

    r2598 r2600  
    2828     &                       periodav, ok_dyn_ave, output_grads_dyn,
    2929     &                       iapp_tracvl
     30       USE comvert_mod, ONLY: ap,bp,pressure_exner,presnivs
    3031       USE comconst_mod, ONLY: cpp, dtvr, ihf, dtphys, pi, jmp1
    3132      IMPLICIT NONE
    … …  
    6566#include "paramet.h"
    6667#include "comdissnew.h"
    67 #include "comvert.h"
    6868#include "comgeom.h"
    6969#include "logic.h"
    … …  
    14541454                call Gather_Field(teta,ip1jmp1,llm,0)
    14551455                call Gather_Field(pk,ip1jmp1,llm,0)
    1456                 call Gather_Field(phi,ip1jmp1,llm,0)
     1456                call Gather_Field(phi,ip1jmp1,llm,0)
    14571457                do iq=1,nqtot
    14581458                  call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
    14591459                enddo
    1460                 call Gather_Field(masse,ip1jmp1,llm,0)
     1460                call Gather_Field(masse,ip1jmp1,llm,0)
    14611461                call Gather_Field(ps,ip1jmp1,1,0)
    1462                 call Gather_Field(phis,ip1jmp1,1,0)
     1462                call Gather_Field(phis,ip1jmp1,1,0)
    14631463                if (mpi_rank==0) then
    14641464                 CALL writedynav(itau,vcov,
    14651465     &                 ucov,teta,pk,phi,q,masse,ps,phis)
    1466                 endif
     1466                endif
    14671467#endif
    14681468!$OMP END MASTER
    … …  
    14761476            IF( MOD(itau,iecri).EQ.0) THEN
    14771477             ! Ehouarn: output only during LF or Backward Matsuno
    1478              if (leapf.or.(.not.leapf.and.(.not.forward))) then
     1478             if (leapf.or.(.not.leapf.and.(.not.forward))) then
    14791479c$OMP BARRIER
    14801480c$OMP MASTER
    … …  
    15141514                call Gather_Field(ucov,ip1jmp1,llm,0)
    15151515                call Gather_Field(teta,ip1jmp1,llm,0)
    1516                 call Gather_Field(phi,ip1jmp1,llm,0)
     1516                call Gather_Field(phi,ip1jmp1,llm,0)
    15171517                do iq=1,nqtot
    15181518                  call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
    15191519                enddo
    1520                 call Gather_Field(masse,ip1jmp1,llm,0)
     1520                call Gather_Field(masse,ip1jmp1,llm,0)
    15211521                call Gather_Field(ps,ip1jmp1,1,0)
    1522                 call Gather_Field(phis,ip1jmp1,1,0)
     1522                call Gather_Field(phis,ip1jmp1,1,0)
    15231523                if (mpi_rank==0) then
    1524                  CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
    1525                 endif
     1524                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
     1525                endif
    15261526!              CALL writehist_p(histid,histvid, itau,vcov,
    15271527!     &                         ucov,teta,phi,q,masse,ps,phis)
    … …  
    16491649                call Gather_Field(teta,ip1jmp1,llm,0)
    16501650                call Gather_Field(pk,ip1jmp1,llm,0)
    1651                 call Gather_Field(phi,ip1jmp1,llm,0)
     1651                call Gather_Field(phi,ip1jmp1,llm,0)
    16521652                do iq=1,nqtot
    16531653                  call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
    16541654                enddo
    1655                 call Gather_Field(masse,ip1jmp1,llm,0)
     1655                call Gather_Field(masse,ip1jmp1,llm,0)
    16561656                call Gather_Field(ps,ip1jmp1,1,0)
    1657                 call Gather_Field(phis,ip1jmp1,1,0)
     1657                call Gather_Field(phis,ip1jmp1,1,0)
    16581658                if (mpi_rank==0) then
    16591659                 CALL writedynav(itau,vcov,
    16601660     &                 ucov,teta,pk,phi,q,masse,ps,phis)
    1661                 endif
     1661                endif
    16621662#endif
    16631663!$OMP END MASTER
    … …  
    17051705                call Gather_Field(ucov,ip1jmp1,llm,0)
    17061706                call Gather_Field(teta,ip1jmp1,llm,0)
    1707                 call Gather_Field(phi,ip1jmp1,llm,0)
     1707                call Gather_Field(phi,ip1jmp1,llm,0)
    17081708                do iq=1,nqtot
    17091709                  call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
    17101710                enddo
    1711                 call Gather_Field(masse,ip1jmp1,llm,0)
     1711                call Gather_Field(masse,ip1jmp1,llm,0)
    17121712                call Gather_Field(ps,ip1jmp1,1,0)
    1713                 call Gather_Field(phis,ip1jmp1,1,0)
     1713                call Gather_Field(phis,ip1jmp1,1,0)
    17141714                if (mpi_rank==0) then
    17151715                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
    1716                 endif
     1716                endif
    17171717!                CALL writehist_p(histid, histvid, itau,vcov ,
    17181718!     &                           ucov,teta,phi,q,masse,ps,phis)

  • LMDZ5/trunk/libf/dyn3dpar/qminimum_p.F

    r1907 r2600  
    88#include "dimensions.h"
    99#include "paramet.h"
    10 #include "comvert.h"
    1110c
    1211      INTEGER nq

  • LMDZ5/trunk/libf/dyn3dpar/sw_case_williamson91_6.F

    r2597 r2600  
    2727c=======================================================================
    2828      USE comconst_mod, ONLY: cpp, omeg, rad
     29      USE comvert_mod, ONLY: ap, bp, preff
     30     
    2931      IMPLICIT NONE
    3032c-----------------------------------------------------------------------
    … …  
    3436#include "dimensions.h"
    3537#include "paramet.h"
    36 #include "comvert.h"
    3738#include "comgeom.h"
    3839#include "iniprint.h"

  • LMDZ5/trunk/libf/dyn3dpar/top_bound_p.F

    r2597 r2600  
    66      USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound,
    77     &                        tau_top_bound
     8      USE comvert_mod, ONLY: presnivs, preff, scaleheight
     9     
    810      IMPLICIT NONE
    911c
    1012#include "dimensions.h"
    1113#include "paramet.h"
    12 #include "comvert.h"
    1314#include "comgeom2.h"
    1415
    … …  
    7071      REAL uzon(jjp1,llm),vzon(jjm,llm),tzon(jjp1,llm)
    7172     
    72       integer i 
     73      integer i
    7374      REAL,SAVE :: rdamp(llm) ! quenching coefficient
    7475      real,save :: lambda(llm) ! inverse or quenching time scale (Hz)

  • LMDZ5/trunk/libf/dyn3dpar/vitvert_p.F

    r1907 r2600  
    22c
    33      USE parallel_lmdz
     4      USE comvert_mod, ONLY: bp
     5     
    46      IMPLICIT NONE
    57
    … …  
    2830#include "dimensions.h"
    2931#include "paramet.h"
    30 #include "comvert.h"
    3132
    3233      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)

  • LMDZ5/trunk/libf/dyn3dpar/vlsplt_p.F

    r2597 r2600  
    2626#include "paramet.h"
    2727#include "logic.h"
    28 #include "comvert.h"
    2928
    3029c
    … …  
    208207#include "paramet.h"
    209208#include "logic.h"
    210 #include "comvert.h"
    211209c
    212210c
    … …  
    543541#include "paramet.h"
    544542#include "logic.h"
    545 #include "comvert.h"
    546543#include "comgeom.h"
    547544c
    … …  
    931928#include "paramet.h"
    932929#include "logic.h"
    933 #include "comvert.h"
    934930c
    935931c

  • LMDZ5/trunk/libf/dyn3dpar/vlspltgen_p.F

    r2597 r2600  
    3434#include "paramet.h"
    3535#include "logic.h"
    36 #include "comvert.h"
    3736
    3837c
    … …  
    105104
    106105       
    107         ijb=ij_begin-iip1
    108         ije=ij_end+iip1
    109         if (pole_nord) ijb=ij_begin
    110         if (pole_sud) ije=ij_end
    111        
     106        ijb=ij_begin-iip1
     107        ije=ij_end+iip1
     108        if (pole_nord) ijb=ij_begin
     109        if (pole_sud) ije=ij_end
     110       
    112111c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    113         DO l = 1, llm
     112        DO l = 1, llm
    114113         DO ij = ijb, ije
    115114          tempe(ij) = teta(ij,l) * pk(ij,l) /cpp
    … …  
    191190
    192191        if(iadv(iq) == 0) then
    193        
    194           cycle
    195        
    196         else if (iadv(iq)==10) then
     192       
     193          cycle
     194       
     195        else if (iadv(iq)==10) then
    197196
    198197#ifdef _ADV_HALO       
    199           call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
    200      &               ij_begin,ij_begin+2*iip1-1)
     198          call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
     199     &                     ij_begin,ij_begin+2*iip1-1)
    201200          call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
    202201     &               ij_end-2*iip1+1,ij_end)
    203202#else
    204           call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
    205      &               ij_begin,ij_end)
     203          call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
     204     &                     ij_begin,ij_end)
    206205#endif
    207206
    … …  
    215214          call VTe(VTHallo)
    216215c$OMP END MASTER
    217         else if (iadv(iq)==14) then
     216        else if (iadv(iq)==14) then
    218217
    219218#ifdef _ADV_HALO           
    … …  
    239238c$OMP END MASTER
    240239        else
    241        
    242           stop 'vlspltgen_p : schema non parallelise'
     240       
     241          stop 'vlspltgen_p : schema non parallelise'
    243242     
    244243        endif
    … …  
    261260
    262261        if(iadv(iq) == 0) then
    263        
    264           cycle
    265        
    266         else if (iadv(iq)==10) then
     262       
     263          cycle
     264       
     265        else if (iadv(iq)==10) then
    267266
    268267#ifdef _ADV_HALLO
    … …  
    270269     &               ij_begin+2*iip1,ij_end-2*iip1)
    271270#endif       
    272         else if (iadv(iq)==14) then
     271        else if (iadv(iq)==14) then
    273272#ifdef _ADV_HALLO
    274273          call vlxqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,qsat,
    … …  
    276275#endif   
    277276        else
    278        
    279           stop 'vlspltgen_p : schema non parallelise'
     277       
     278          stop 'vlspltgen_p : schema non parallelise'
    280279     
    281280        endif
    … …  
    301300
    302301        if(iadv(iq) == 0) then
    303        
    304           cycle
    305        
    306         else if (iadv(iq)==10) then
     302       
     303          cycle
     304       
     305        else if (iadv(iq)==10) then
    307306       
    308307          call vly_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv)
    309308 
    310         else if (iadv(iq)==14) then
     309        else if (iadv(iq)==14) then
    311310     
    312311          call vlyqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv,qsat)
    313312 
    314313        else
    315        
    316           stop 'vlspltgen_p : schema non parallelise'
     314       
     315          stop 'vlspltgen_p : schema non parallelise'
    317316     
    318317        endif
    … …  
    324323
    325324        if(iadv(iq) == 0) then
    326          
    327           cycle
    328        
    329         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     325         
     326          cycle
     327       
     328        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    330329
    331330c$OMP BARRIER       
    … …  
    350349c$OMP MASTER
    351350          call VTe(VTHallo)
    352 c$OMP END MASTER       
    353 c$OMP BARRIER
    354         else
    355        
    356           stop 'vlspltgen_p : schema non parallelise'
     351c$OMP END MASTER       
     352c$OMP BARRIER
     353        else
     354       
     355          stop 'vlspltgen_p : schema non parallelise'
    357356     
    358357        endif
    … …  
    369368c$OMP MASTER
    370369      call VTe(VTHallo)
    371 c$OMP END MASTER       
     370c$OMP END MASTER       
    372371
    373372c$OMP BARRIER
    … …  
    375374
    376375        if(iadv(iq) == 0) then
    377          
    378           cycle
    379        
    380         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     376         
     377          cycle
     378       
     379        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    381380c$OMP BARRIER       
    382381
    … …  
    388387c$OMP BARRIER       
    389388        else
    390        
    391           stop 'vlspltgen_p : schema non parallelise'
     389       
     390          stop 'vlspltgen_p : schema non parallelise'
    392391     
    393392        endif
    … …  
    414413
    415414        if(iadv(iq) == 0) then
    416        
    417           cycle
    418        
    419         else if (iadv(iq)==10) then
     415       
     416          cycle
     417       
     418        else if (iadv(iq)==10) then
    420419       
    421420          call vly_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv)
    422421 
    423         else if (iadv(iq)==14) then
     422        else if (iadv(iq)==14) then
    424423     
    425424          call vlyqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mv,qsat)
    426425 
    427426        else
    428        
    429           stop 'vlspltgen_p : schema non parallelise'
     427       
     428          stop 'vlspltgen_p : schema non parallelise'
    430429     
    431430        endif
    … …  
    436435
    437436        if(iadv(iq) == 0) then
    438          
    439           cycle
    440        
    441         else if (iadv(iq)==10) then
     437         
     438          cycle
     439       
     440        else if (iadv(iq)==10) then
    442441       
    443442          call vlx_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,
    444443     &               ij_begin,ij_end)
    445444 
    446         else if (iadv(iq)==14) then
     445        else if (iadv(iq)==14) then
    447446     
    448447          call vlxqs_p(zq(1,1,iq),pente_max,zm(1,1,iq),mu,qsat,
    … …  
    450449 
    451450        else
    452        
     451       
    453452          stop 'vlspltgen_p : schema non parallelise'
    454453     
    … …  
    469468           DO ij=ijb,ije
    470469c             print *,'zq-->',ij,l,iq,zq(ij,l,iq)
    471 c            print *,'q-->',ij,l,iq,q(ij,l,iq)
    472              q(ij,l,iq)=zq(ij,l,iq)
     470c             print *,'q-->',ij,l,iq,q(ij,l,iq)
     471             q(ij,l,iq)=zq(ij,l,iq)
    473472           ENDDO
    474473        ENDDO

  • LMDZ5/trunk/libf/dyn3dpar/vlspltqs_p.F

    r2597 r2600  
    3232#include "paramet.h"
    3333#include "logic.h"
    34 #include "comvert.h"
    3534
    3635c
    … …  
    237236#include "paramet.h"
    238237#include "logic.h"
    239 #include "comvert.h"
    240238c
    241239c
    … …  
    589587#include "paramet.h"
    590588#include "logic.h"
    591 #include "comvert.h"
    592589#include "comgeom.h"
    593590c

  • LMDZ5/trunk/libf/dyn3dpar/writedynav_p.F

    r2598 r2600  
    4343#include "dimensions.h"
    4444#include "paramet.h"
    45 #include "comvert.h"
    4645#include "comgeom.h"
    4746#include "temps.h"

  • LMDZ5/trunk/libf/dyn3dpar/writehist_p.F

    r2598 r2600  
    4343#include "dimensions.h"
    4444#include "paramet.h"
    45 #include "comvert.h"
    4645#include "comgeom.h"
    4746#include "temps.h"

  • 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

  • LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

    r2565 r2600  
    950950!---------------------------------------------------------------------
    951951       call physiq(ngrid,llm, &
    952                     firstcall,lastcall,timestep, &
    953                     plev,play,phi,phis,presnivs, &
    954                     u,v, rot, temp,q,omega2, &
    955                     du_phys,dv_phys,dt_phys,dq,dpsrf)
    956                 firstcall=.false.
     952                    firstcall,lastcall,timestep, &
     953                    plev,play,phi,phis,presnivs, &
     954                    u,v, rot, temp,q,omega2, &
     955                    du_phys,dv_phys,dt_phys,dq,dpsrf)
     956                firstcall=.false.
    957957
    958958!---------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.