Ignore:
Timestamp:
Jul 22, 2016, 8:44:47 AM (9 years ago)
Author:
Ehouarn Millour
Message:

Cleanup in the dynamics: get rid of comconst.h, make it a module comconst_mod.
EM

Location:
LMDZ5/trunk/libf/dyn3dmem
Files:
33 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dmem/addfi_loc.F

    r1987 r2597  
    4646c    ------------------
    4747c
    48 #include "dimensions.h"
    49 #include "paramet.h"
    50 #include "comconst.h"
    51 #include "comgeom.h"
    52 #include "serre.h"
     48      include "dimensions.h"
     49      include "paramet.h"
     50      include "comgeom.h"
     51      include "serre.h"
    5352c
    5453c    Arguments :
  • LMDZ5/trunk/libf/dyn3dmem/advect_new_loc.F

    r1907 r2597  
    77      USE write_field_loc
    88      USE advect_new_mod
     9      USE comconst_mod, ONLY: daysec
    910      IMPLICIT NONE
    1011c=======================================================================
     
    2728c   -------------
    2829
    29 #include "dimensions.h"
    30 #include "paramet.h"
    31 #include "comconst.h"
    32 #include "comvert.h"
    33 #include "comgeom.h"
    34 #include "logic.h"
    35 #include "ener.h"
     30      include "dimensions.h"
     31      include "paramet.h"
     32      include "comvert.h"
     33      include "comgeom.h"
     34      include "logic.h"
     35      include "ener.h"
    3636
    3737c   Arguments:
     
    122122                 
    123123           uav(ij,l)=0.25*(ucov(ij,l)+ucov(ij-iip1,l))
    124      .               +0.25*(ucov(ij+iip1,l)+ucov(ij,l))
     124     .               +0.25*(ucov(ij+iip1,l)+ucov(ij,l))
    125125         ENDDO
    126126         
  • LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F

    r2286 r2597  
    2727      USE control_mod, ONLY: iapp_tracvl, day_step, planet_type
    2828      USE advtrac_mod, ONLY: finmasse
     29      USE comconst_mod, ONLY: dtvr
    2930      IMPLICIT NONE
    3031c
    31 #include "dimensions.h"
    32 #include "paramet.h"
    33 #include "comconst.h"
    34 #include "comvert.h"
    35 #include "comdissip.h"
    36 #include "comgeom2.h"
    37 #include "logic.h"
    38 #include "temps.h"
    39 #include "ener.h"
    40 #include "description.h"
     32      include "dimensions.h"
     33      include "paramet.h"
     34      include "comvert.h"
     35      include "comdissip.h"
     36      include "comgeom2.h"
     37      include "logic.h"
     38      include "temps.h"
     39      include "ener.h"
     40      include "description.h"
    4141
    4242c-------------------------------------------------------------------
  • LMDZ5/trunk/libf/dyn3dmem/bilan_dyn_loc.F

    r2475 r2597  
    1717      use misc_mod
    1818      USE write_field_loc
     19      USE comconst_mod, ONLY: cpp, pi
    1920      IMPLICIT NONE
    2021
    21 #include "dimensions.h"
    22 #include "paramet.h"
    23 #include "comconst.h"
    24 #include "comvert.h"
    25 #include "comgeom2.h"
    26 #include "temps.h"
    27 #include "iniprint.h"
     22      include "dimensions.h"
     23      include "paramet.h"
     24      include "comvert.h"
     25      include "comgeom2.h"
     26      include "temps.h"
     27      include "iniprint.h"
    2828
    2929c====================================================================
  • LMDZ5/trunk/libf/dyn3dmem/caladvtrac_loc.F

    r2286 r2597  
    2929
    3030
    31 #include "dimensions.h"
    32 #include "paramet.h"
    33 #include "comconst.h"
     31      include "dimensions.h"
     32      include "paramet.h"
    3433
    3534c   Arguments:
     
    210209        call resume_timer(timer_caldyn)
    211210c$OMP END MASTER
    212 c$OMP BARRIER   
     211c$OMP BARRIER
    213212          iadvtr=0
    214213       ENDIF ! if iadvtr.EQ.iapp_tracvl
  • LMDZ5/trunk/libf/dyn3dmem/caldyn_loc.F

    r1987 r2597  
    3131!   ----------------
    3232
    33 #include "dimensions.h"
    34 #include "paramet.h"
    35 #include "comconst.h"
    36 #include "comvert.h"
    37 #include "comgeom.h"
     33      include "dimensions.h"
     34      include "paramet.h"
     35      include "comvert.h"
     36      include "comgeom.h"
    3837
    3938!   Arguments:
  • LMDZ5/trunk/libf/dyn3dmem/call_calfis_mod.F90

    r2375 r2597  
    8282  USE write_field_loc
    8383  USE write_field
     84  USE comconst_mod, ONLY: dtphys
    8485  IMPLICIT NONE
    85     INCLUDE "comconst.h"
    8686    INCLUDE "comvert.h"
    8787    INCLUDE "logic.h"
  • LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F90

    r2444 r2597  
    1717  USE infotrac, ONLY : type_trac
    1818  use assert_m, only: assert
     19  USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
     20                          iflag_top_bound, mode_top_bound, tau_top_bound, &
     21                          ngroup
    1922
    2023  IMPLICIT NONE
     
    3942  include "comdissnew.h"
    4043  include "temps.h"
    41   include "comconst.h"
    4244  include "iniprint.h"
    4345
  • LMDZ5/trunk/libf/dyn3dmem/dissip_loc.F

    r1987 r2597  
    77      USE write_field_loc
    88      USE dissip_mod, ONLY: dissip_allocate
     9      USE comconst_mod, ONLY: dtdiss
    910      IMPLICIT NONE
    1011
     
    2829c   -------------
    2930
    30 #include "dimensions.h"
    31 #include "paramet.h"
    32 #include "comconst.h"
    33 #include "comgeom.h"
    34 #include "comdissnew.h"
    35 #include "comdissipn.h"
     31      include "dimensions.h"
     32      include "paramet.h"
     33      include "comgeom.h"
     34      include "comdissnew.h"
     35      include "comdissipn.h"
    3636
    3737c   Arguments:
  • LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.f90

    r2299 r2597  
    1212  USE control_mod, ONLY: planet_type
    1313  USE assert_eq_m, ONLY: assert_eq
     14  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa, lllm, &
     15                          omeg, rad
    1416  IMPLICIT NONE
    1517  include "dimensions.h"
    1618  include "paramet.h"
    1719  include "temps.h"
    18   include "comconst.h"
    1920  include "comvert.h"
    2021  include "comgeom.h"
  • LMDZ5/trunk/libf/dyn3dmem/dynredem_loc.F90

    r2584 r2597  
    1313                    NF90_CLOSE,  NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER
    1414  USE dynredem_mod, ONLY: cre_var, put_var, err, modname, fil
     15  USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    1516  IMPLICIT NONE
    1617  include "dimensions.h"
    1718  include "paramet.h"
    18   include "comconst.h"
    1919  include "comvert.h"
    2020  include "comgeom.h"
  • LMDZ5/trunk/libf/dyn3dmem/exner_hyb_loc_m.F90

    r2021 r2597  
    3535    USE mod_filtreg_p
    3636    USE write_field_loc
     37    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     38    IMPLICIT NONE
    3739    !
    3840    include "dimensions.h"
    3941    include "paramet.h"
    40     include "comconst.h"
    4142    include "comgeom.h"
    4243    include "comvert.h"
  • LMDZ5/trunk/libf/dyn3dmem/exner_milieu_loc_m.F90

    r2021 r2597  
    3131    USE parallel_lmdz
    3232    USE mod_filtreg_p
     33    USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     34    IMPLICIT NONE
    3335    !
    3436    include "dimensions.h"
    3537    include "paramet.h"
    36     include "comconst.h"
    3738    include "comgeom.h"
    3839    include "comvert.h"
  • LMDZ5/trunk/libf/dyn3dmem/friction_loc.F

    r1907 r2597  
    1212      USE ioipsl_getincom
    1313#endif
     14      USE comconst_mod, ONLY: pi
    1415      IMPLICIT NONE
    1516
     
    2627!=======================================================================
    2728
    28 #include "dimensions.h"
    29 #include "paramet.h"
    30 #include "comgeom2.h"
    31 #include "comconst.h"
    32 #include "iniprint.h"
    33 #include "academic.h"
     29      include "dimensions.h"
     30      include "paramet.h"
     31      include "comgeom2.h"
     32      include "iniprint.h"
     33      include "academic.h"
    3434
    3535! arguments:
  • LMDZ5/trunk/libf/dyn3dmem/gcm.F90

    r2475 r2597  
    2121  USE iniphysiq_mod, ONLY: iniphysiq
    2222#endif
     23  USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, r, rad
    2324  IMPLICIT NONE
    2425
     
    5455  include "dimensions.h"
    5556  include "paramet.h"
    56   include "comconst.h"
    5757  include "comdissnew.h"
    5858  include "comvert.h"
  • LMDZ5/trunk/libf/dyn3dmem/groupe_loc.F

    r2442 r2597  
    33      USE Write_field_loc
    44      USE groupe_mod
     5      USE comconst_mod, ONLY: ngroup
    56      implicit none
    67
     
    1617c   pas besoin de w en entree.
    1718
    18 #include "dimensions.h"
    19 #include "paramet.h"
    20 #include "comconst.h"
    21 #include "comgeom2.h"
    22 #include "comvert.h"
     19      include "dimensions.h"
     20      include "paramet.h"
     21      include "comgeom2.h"
     22      include "comvert.h"
    2323
    2424!     integer ngroup
  • LMDZ5/trunk/libf/dyn3dmem/groupeun_loc.F

    r2442 r2597  
    22      USE parallel_lmdz
    33      USE Write_Field_p
     4      USE comconst_mod, ONLY: ngroup
    45      IMPLICIT NONE
    56
    6 #include "dimensions.h"
    7 #include "paramet.h"
    8 #include "comconst.h"
    9 #include "comgeom2.h"
     7      include "dimensions.h"
     8      include "paramet.h"
     9      include "comgeom2.h"
    1010
    1111      INTEGER jjmax,llmax,sb,se,jjb,jje
     
    136136
    137137      USE parallel_lmdz
     138      USE comconst_mod, ONLY: ngroup
    138139      IMPLICIT NONE
    139140
    140 #include "dimensions.h"
    141 #include "paramet.h"
    142 #include "comconst.h"
    143 #include "comgeom2.h"
     141      include "dimensions.h"
     142      include "paramet.h"
     143      include "comgeom2.h"
    144144
    145145!     INTEGER ngroup
  • LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

    r2263 r2597  
    346346    USE control_mod
    347347    USE write_field_loc
     348    USE comconst_mod, ONLY: cpp, daysec, dtvr, kappa
    348349   
    349350    IMPLICIT NONE
     
    351352    INCLUDE "dimensions.h"
    352353    INCLUDE "paramet.h"
    353     INCLUDE "comconst.h"
    354354    INCLUDE "comvert.h"
    355355
     
    748748  SUBROUTINE guide_zonave_u(typ,vsize,field)
    749749
     750    USE comconst_mod, ONLY: pi
     751   
    750752    IMPLICIT NONE
    751753
     
    753755    INCLUDE "paramet.h"
    754756    INCLUDE "comgeom.h"
    755     INCLUDE "comconst.h"
    756757   
    757758    ! input/output variables
     
    819820  SUBROUTINE guide_zonave_v(typ,hsize,vsize,field)
    820821
     822    USE comconst_mod, ONLY: pi
     823   
    821824    IMPLICIT NONE
    822825
     
    824827    INCLUDE "paramet.h"
    825828    INCLUDE "comgeom.h"
    826     INCLUDE "comconst.h"
    827829   
    828830    ! input/output variables
     
    890892  USE mod_hallo
    891893  USE Bands
     894  USE comconst_mod, ONLY: cpp, kappa
    892895  IMPLICIT NONE
    893896
     
    896899  include "comvert.h"
    897900  include "comgeom2.h"
    898   include "comconst.h"
    899901
    900902  REAL, DIMENSION (iip1,jjb_u:jje_u),     INTENT(IN) :: psi ! Psol gcm
     
    10121014!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    10131015        DO l=1,llm
    1014             DO j=jjbu,jjeu
    1015                 DO i =1, iip1
     1016            DO j=jjbu,jjeu
     1017                DO i =1, iip1
    10161018                    pls(i,j,l)=(ap(l)+ap(l+1))/2.+psi(i,j)*(bp(l)+bp(l+1))/2.
    10171019                ENDDO
     
    10301032   DO l = 1, llm
    10311033       DO j=jjbu,jjeu
    1032         DO i =1, iip1
    1033             pls(i,j,l) = preff * ( pk(i,j,l)/cpp) ** unskap
    1034         ENDDO
     1034           DO i =1, iip1
     1035               pls(i,j,l) = preff * ( pk(i,j,l)/cpp) ** unskap
     1036           ENDDO
    10351037       ENDDO
    10361038   ENDDO
     
    13021304        ! Calcul des nouvelles valeurs des niveaux de pression du guidage
    13031305        IF (guide_plevs.EQ.1) THEN
    1304         CALL Register_Hallo_u(psnat1,1,1,2,2,1,Req)
    1305         CALL Register_Hallo_u(psnat2,1,1,2,2,1,Req)
    1306         CALL SendRequest(Req)
    1307 !$OMP BARRIER
    1308         CALL WaitRequest(Req)
     1306        CALL Register_Hallo_u(psnat1,1,1,2,2,1,Req)
     1307        CALL Register_Hallo_u(psnat2,1,1,2,2,1,Req)
     1308        CALL SendRequest(Req)
     1309!$OMP BARRIER
     1310        CALL WaitRequest(Req)
    13091311!$OMP BARRIER
    13101312!$OMP DO
     
    13201322            ENDDO
    13211323        ELSE IF (guide_plevs.EQ.2) THEN
    1322         CALL Register_Hallo_u(pnat1,llm,1,2,2,1,Req)
    1323         CALL Register_Hallo_u(pnat2,llm,1,2,2,1,Req)
    1324         CALL SendRequest(Req)
    1325 !$OMP BARRIER
    1326         CALL WaitRequest(Req)
     1324        CALL Register_Hallo_u(pnat1,llm,1,2,2,1,Req)
     1325        CALL Register_Hallo_u(pnat2,llm,1,2,2,1,Req)
     1326        CALL SendRequest(Req)
     1327!$OMP BARRIER
     1328        CALL WaitRequest(Req)
    13271329!$OMP BARRIER
    13281330!$OMP DO
     
    13701372! Calcul des constantes de rappel alpha (=1/tau)
    13711373
     1374    use comconst_mod, only: pi
     1375   
    13721376    implicit none
    13731377
    13741378    include "dimensions.h"
    13751379    include "paramet.h"
    1376     include "comconst.h"
    13771380    include "comgeom2.h"
    13781381    include "serre.h"
     
    21272130    USE parallel_lmdz
    21282131    USE mod_hallo, ONLY : gather_field_u, gather_field_v
     2132    USE comconst_mod, ONLY: pi
    21292133    IMPLICIT NONE
    21302134
     
    21332137    INCLUDE "netcdf.inc"
    21342138    INCLUDE "comgeom2.h"
    2135     INCLUDE "comconst.h"
    21362139    INCLUDE "comvert.h"
    21372140   
  • LMDZ5/trunk/libf/dyn3dmem/iniacademic_loc.F90

    r2270 r2597  
    1818#endif
    1919  USE Write_Field
     20  USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
    2021
    2122  !   Author:    Frederic Hourdin      original: 15/01/93
     
    3132  include "paramet.h"
    3233  include "comvert.h"
    33   include "comconst.h"
    3434  include "comgeom.h"
    3535  include "academic.h"
  • LMDZ5/trunk/libf/dyn3dmem/initdynav_loc.F

    r2475 r2597  
    1414       use com_io_dyn_mod, only : histaveid,histvaveid,histuaveid,       &
    1515     &        dynhistave_file,dynhistvave_file,dynhistuave_file
     16       USE comconst_mod, ONLY: pi
    1617       implicit none
    1718
     
    4142C
    4243C   Declarations
    43 #include "dimensions.h"
    44 #include "paramet.h"
    45 #include "comconst.h"
    46 #include "comvert.h"
    47 #include "comgeom.h"
    48 #include "temps.h"
    49 #include "ener.h"
    50 #include "logic.h"
    51 #include "description.h"
    52 #include "serre.h"
    53 #include "iniprint.h"
     44      include "dimensions.h"
     45      include "paramet.h"
     46      include "comvert.h"
     47      include "comgeom.h"
     48      include "temps.h"
     49      include "ener.h"
     50      include "logic.h"
     51      include "description.h"
     52      include "serre.h"
     53      include "iniprint.h"
    5454
    5555C   Arguments
  • LMDZ5/trunk/libf/dyn3dmem/initfluxsto_p.F

    r1907 r2597  
    1313       use Write_field
    1414       use misc_mod
     15       USE comconst_mod, ONLY: pi
    1516       
    1617      implicit none
     
    4344C
    4445C   Declarations
    45 #include "dimensions.h"
    46 #include "paramet.h"
    47 #include "comconst.h"
    48 #include "comvert.h"
    49 #include "comgeom.h"
    50 #include "temps.h"
    51 #include "ener.h"
    52 #include "logic.h"
    53 #include "description.h"
    54 #include "serre.h"
    55 #include "iniprint.h"
     46      include "dimensions.h"
     47      include "paramet.h"
     48      include "comvert.h"
     49      include "comgeom.h"
     50      include "temps.h"
     51      include "ener.h"
     52      include "logic.h"
     53      include "description.h"
     54      include "serre.h"
     55      include "iniprint.h"
    5656
    5757C   Arguments
     
    166166     .             filevid,dynv_domain_id)
    167167       
    168       rl(1,1) = 1.     
     168      rl(1,1) = 1.
    169169     
    170170      if (mpi_rank==0) then
  • LMDZ5/trunk/libf/dyn3dmem/inithist_loc.F

    r2475 r2597  
    1414       use com_io_dyn_mod, only : histid,histvid,histuid,               &
    1515     &                        dynhist_file,dynhistv_file,dynhistu_file
     16       USE comconst_mod, ONLY: pi
    1617       implicit none
    1718
     
    4041C
    4142C   Declarations
    42 #include "dimensions.h"
    43 #include "paramet.h"
    44 #include "comconst.h"
    45 #include "comvert.h"
    46 #include "comgeom.h"
    47 #include "temps.h"
    48 #include "ener.h"
    49 #include "logic.h"
    50 #include "description.h"
    51 #include "serre.h"
    52 #include "iniprint.h"
     43      include "dimensions.h"
     44      include "paramet.h"
     45      include "comvert.h"
     46      include "comgeom.h"
     47      include "temps.h"
     48      include "ener.h"
     49      include "logic.h"
     50      include "description.h"
     51      include "serre.h"
     52      include "iniprint.h"
    5353
    5454C   Arguments
  • LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F

    r2461 r2597  
    1212      USE integrd_mod
    1313      USE infotrac, ONLY: ok_iso_verif ! ajout CRisi
     14      USE comconst_mod, ONLY: pi
    1415      IMPLICIT NONE
    1516
     
    3031c   -------------
    3132
    32 #include "dimensions.h"
    33 #include "paramet.h"
    34 #include "comconst.h"
    35 #include "comgeom.h"
    36 #include "comvert.h"
    37 #include "logic.h"
    38 #include "temps.h"
    39 #include "serre.h"
    40 #include "iniprint.h"
    41 !      include 'mpif.h'
     33      include "dimensions.h"
     34      include "paramet.h"
     35      include "comgeom.h"
     36      include "comvert.h"
     37      include "logic.h"
     38      include "temps.h"
     39      include "serre.h"
     40      include "iniprint.h"
    4241
    4342c   Arguments:
     
    335334      ije=ij_end
    336335
    337         if (planet_type.eq."earth") then
     336        if (planet_type.eq."earth") then
    338337! Earth-specific treatment of first 2 tracers (water)
    339338c$OMP BARRIER
     
    452451          DO l = 1, llm
    453452            massem1(ijb:ije,l)=massescr(ijb:ije,l)
    454           ENDDO
    455 c$OMP END DO NOWAIT      
     453          ENDDO
     454c$OMP END DO NOWAIT         
    456455      END IF
    457456c$OMP BARRIER
  • LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F

    r2475 r2597  
    3232       use exner_hyb_loc_m, only: exner_hyb_loc
    3333       use exner_milieu_loc_m, only: exner_milieu_loc
     34       USE comconst_mod, ONLY: cpp, dtvr, ihf
    3435      IMPLICIT NONE
    3536
     
    6566c   -------------
    6667
    67 #include "dimensions.h"
    68 #include "paramet.h"
    69 #include "comconst.h"
    70 #include "comdissnew.h"
    71 #include "comvert.h"
    72 #include "comgeom.h"
    73 #include "logic.h"
    74 #include "temps.h"
    75 #include "ener.h"
    76 #include "description.h"
    77 #include "serre.h"
    78 !#include "com_io_dyn.h"
    79 #include "iniprint.h"
    80 #include "academic.h"
    81 !      include "mpif.h"
     68      include "dimensions.h"
     69      include "paramet.h"
     70      include "comdissnew.h"
     71      include "comvert.h"
     72      include "comgeom.h"
     73      include "logic.h"
     74      include "temps.h"
     75      include "ener.h"
     76      include "description.h"
     77      include "serre.h"
     78      include "iniprint.h"
     79      include "academic.h"
    8280     
    8381      REAL,INTENT(IN) :: time_0 ! not used
  • LMDZ5/trunk/libf/dyn3dmem/massbarxy_loc.F90

    r2336 r2597  
    1010  include "dimensions.h"
    1111  include "paramet.h"
    12   include "comconst.h"
    1312  include "comgeom.h"
    1413!===============================================================================
  • LMDZ5/trunk/libf/dyn3dmem/massdair_loc.F

    r1907 r2597  
    1616      IMPLICIT NONE
    1717c
    18 #include "dimensions.h"
    19 #include "paramet.h"
    20 #include "comconst.h"
    21 #include "comgeom.h"
     18      include "dimensions.h"
     19      include "paramet.h"
     20      include "comgeom.h"
    2221c
    2322c  .....   arguments  ....
  • LMDZ5/trunk/libf/dyn3dmem/sw_case_williamson91_6_loc.F

    r1907 r2597  
    2727c=======================================================================
    2828      USE parallel_lmdz
     29      USE comconst_mod, ONLY: cpp, omeg, rad
    2930
    3031      IMPLICIT NONE
     
    3334c   ---------------
    3435
    35 #include "dimensions.h"
    36 #include "paramet.h"
    37 #include "comvert.h"
    38 #include "comconst.h"
    39 #include "comgeom.h"
    40 #include "iniprint.h"
     36      include "dimensions.h"
     37      include "paramet.h"
     38      include "comvert.h"
     39      include "comgeom.h"
     40      include "iniprint.h"
    4141
    4242c   Arguments:
  • LMDZ5/trunk/libf/dyn3dmem/top_bound_loc.F

    r1907 r2597  
    44      SUBROUTINE top_bound_loc(vcov,ucov,teta,masse,dt)
    55      USE parallel_lmdz
     6      USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound,
     7     &                        tau_top_bound
    68      IMPLICIT NONE
    79c
    8 #include "dimensions.h"
    9 #include "paramet.h"
    10 #include "comconst.h"
    11 #include "comvert.h"
    12 #include "comgeom2.h"
     10      include "dimensions.h"
     11      include "paramet.h"
     12      include "comvert.h"
     13      include "comgeom2.h"
    1314
    1415
     
    4041! NB: top_bound sponge is only called from leapfrog if ok_strato=.true.
    4142
    42 ! sponge parameters: (loaded/set in conf_gcm.F ; stored in comconst.h)
     43! sponge parameters: (loaded/set in conf_gcm.F ; stored in comconst_mod)
    4344!    iflag_top_bound=0 for no sponge
    4445!    iflag_top_bound=1 for sponge over 4 topmost layers
     
    7475      REAL tzon(jjb_u:jje_u,llm)
    7576     
    76       integer i 
     77      integer i
    7778      REAL,SAVE :: rdamp(llm)
    7879      real,save :: lambda(llm) ! inverse or quenching time scale (Hz)
  • LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F

    r2286 r2597  
    1717      IMPLICIT NONE
    1818c
    19 #include "dimensions.h"
    20 #include "paramet.h"
    21 #include "logic.h"
    22 #include "comvert.h"
    23 #include "comconst.h"
     19      include "dimensions.h"
     20      include "paramet.h"
     21      include "logic.h"
     22      include "comvert.h"
    2423c
    2524c
     
    417416      USE parallel_lmdz
    418417      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
     418      USE comconst_mod, ONLY: pi
    419419      IMPLICIT NONE
    420420c
    421 #include "dimensions.h"
    422 #include "paramet.h"
    423 #include "logic.h"
    424 #include "comvert.h"
    425 #include "comconst.h"
    426 #include "comgeom.h"
     421      include "dimensions.h"
     422      include "paramet.h"
     423      include "logic.h"
     424      include "comvert.h"
     425      include "comgeom.h"
    427426c
    428427c
     
    875874      IMPLICIT NONE
    876875c
    877 #include "dimensions.h"
    878 #include "paramet.h"
    879 #include "logic.h"
    880 #include "comvert.h"
    881 #include "comconst.h"
     876      include "dimensions.h"
     877      include "paramet.h"
     878      include "logic.h"
     879      include "comvert.h"
    882880c
    883881c
  • LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F

    r2286 r2597  
    3131     &    ok_iso_verif
    3232      USE vlspltgen_mod
     33      USE comconst_mod, ONLY: cpp
    3334      IMPLICIT NONE
    3435
    3536c
    36 #include "dimensions.h"
    37 #include "paramet.h"
    38 #include "logic.h"
    39 #include "comvert.h"
    40 #include "comconst.h"
     37      include "dimensions.h"
     38      include "paramet.h"
     39      include "logic.h"
     40      include "comvert.h"
    4141
    4242c
     
    100100
    101101       
    102         ijb=ij_begin-iip1
    103         ije=ij_end+iip1
    104         if (pole_nord) ijb=ij_begin
    105         if (pole_sud) ije=ij_end
    106        
     102        ijb=ij_begin-iip1
     103        ije=ij_end+iip1
     104        if (pole_nord) ijb=ij_begin
     105        if (pole_sud) ije=ij_end
     106       
    107107c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    108         DO l = 1, llm
     108        DO l = 1, llm
    109109         DO ij = ijb, ije
    110110          tempe(ij) = teta(ij,l) * pk(ij,l) /cpp
     
    208208#endif
    209209        if(iadv(iq) == 0) then
    210        
    211           cycle
    212        
    213         else if (iadv(iq)==10) then
     210       
     211          cycle
     212       
     213        else if (iadv(iq)==10) then
    214214
    215215#ifdef _ADV_HALO       
    216216! CRisi: on ajoute les nombres de fils et tableaux des fils
    217217! On suppose qu'on ne peut advecter les fils que par le schéma 10. 
    218           call vlx_loc(zq,pente_max,zm,mu,
    219      &               ij_begin,ij_begin+2*iip1-1,iq)
     218          call vlx_loc(zq,pente_max,zm,mu,
     219     &                     ij_begin,ij_begin+2*iip1-1,iq)
    220220          call vlx_loc(zq,pente_max,zm,mu,
    221221     &               ij_end-2*iip1+1,ij_end,iq)
    222222#else
    223           call vlx_loc(zq,pente_max,zm,mu,
    224      &               ij_begin,ij_end,iq)
     223          call vlx_loc(zq,pente_max,zm,mu,
     224     &                     ij_begin,ij_end,iq)
    225225#endif
    226226
     
    240240          call VTe(VTHallo)
    241241c$OMP END MASTER
    242         else if (iadv(iq)==14) then
     242        else if (iadv(iq)==14) then
    243243
    244244#ifdef _ADV_HALO           
     
    268268c$OMP END MASTER
    269269        else
    270        
    271           stop 'vlspltgen_p : schema non parallelise'
     270       
     271          stop 'vlspltgen_p : schema non parallelise'
    272272     
    273273        endif
     
    301301
    302302        if(iadv(iq) == 0) then
    303        
    304           cycle
    305        
    306         else if (iadv(iq)==10) then
     303       
     304          cycle
     305       
     306        else if (iadv(iq)==10) then
    307307
    308308#ifdef _ADV_HALLO
     
    310310     &                 ij_begin+2*iip1,ij_end-2*iip1,iq)
    311311#endif       
    312         else if (iadv(iq)==14) then
     312        else if (iadv(iq)==14) then
    313313#ifdef _ADV_HALLO
    314314          call vlxqs_loc(zq,pente_max,zm,mu,
     
    316316#endif   
    317317        else
    318        
    319           stop 'vlspltgen_p : schema non parallelise'
     318       
     319          stop 'vlspltgen_p : schema non parallelise'
    320320     
    321321        endif
     
    358358
    359359        if(iadv(iq) == 0) then
    360        
    361           cycle
    362        
    363         else if (iadv(iq)==10) then
     360       
     361          cycle
     362       
     363        else if (iadv(iq)==10) then
    364364       
    365365          call vly_loc(zq,pente_max,zm,mv,iq)
    366366 
    367         else if (iadv(iq)==14) then
     367        else if (iadv(iq)==14) then
    368368     
    369369          call vlyqs_loc(zq,pente_max,zm,mv,
     
    371371 
    372372        else
    373        
    374           stop 'vlspltgen_p : schema non parallelise'
     373       
     374          stop 'vlspltgen_p : schema non parallelise'
    375375     
    376376        endif
     
    389389#endif
    390390        if(iadv(iq) == 0) then
    391          
    392           cycle
    393        
    394         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     391         
     392          cycle
     393       
     394        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    395395
    396396c$OMP BARRIER       
     
    420420c$OMP MASTER
    421421          call VTe(VTHallo)
    422 c$OMP END MASTER       
     422c$OMP END MASTER       
    423423c$OMP BARRIER
    424424        else
    425        
    426           stop 'vlspltgen_p : schema non parallelise'
     425       
     426          stop 'vlspltgen_p : schema non parallelise'
    427427     
    428428        endif
     
    439439c$OMP MASTER
    440440      call VTe(VTHallo)
    441 c$OMP END MASTER       
     441c$OMP END MASTER       
    442442
    443443
     
    451451
    452452        if(iadv(iq) == 0) then
    453          
    454           cycle
    455        
    456         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     453         
     454          cycle
     455       
     456        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    457457c$OMP BARRIER       
    458458
     
    464464c$OMP BARRIER       
    465465        else
    466        
    467           stop 'vlspltgen_p : schema non parallelise'
     466       
     467          stop 'vlspltgen_p : schema non parallelise'
    468468     
    469469        endif
     
    501501#endif
    502502        if(iadv(iq) == 0) then
    503        
    504           cycle
    505        
    506         else if (iadv(iq)==10) then
     503       
     504          cycle
     505       
     506        else if (iadv(iq)==10) then
    507507       
    508508          call vly_loc(zq,pente_max,zm,mv,iq)
    509509 
    510         else if (iadv(iq)==14) then
     510        else if (iadv(iq)==14) then
    511511     
    512512          call vlyqs_loc(zq,pente_max,zm,mv,
     
    514514 
    515515        else
    516        
    517           stop 'vlspltgen_p : schema non parallelise'
     516       
     517          stop 'vlspltgen_p : schema non parallelise'
    518518     
    519519        endif
     
    532532#endif
    533533        if(iadv(iq) == 0) then
    534          
    535           cycle
    536        
    537         else if (iadv(iq)==10) then
     534         
     535          cycle
     536       
     537        else if (iadv(iq)==10) then
    538538       
    539539          call vlx_loc(zq,pente_max,zm,mu,
    540540     &               ij_begin,ij_end,iq)
    541541 
    542         else if (iadv(iq)==14) then
     542        else if (iadv(iq)==14) then
    543543     
    544544          call vlxqs_loc(zq,pente_max,zm,mu,
     
    546546 
    547547        else
    548        
     548       
    549549          stop 'vlspltgen_p : schema non parallelise'
    550550     
     
    574574           DO ij=ijb,ije
    575575c             print *,'zq-->',ij,l,iq,zq(ij,l,iq)
    576 c            print *,'q-->',ij,l,iq,q(ij,l,iq)
    577              q(ij,l,iq)=zq(ij,l,iq)
     576c             print *,'q-->',ij,l,iq,q(ij,l,iq)
     577             q(ij,l,iq)=zq(ij,l,iq)
    578578           ENDDO
    579579        ENDDO
  • LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F

    r2286 r2597  
    1212      IMPLICIT NONE
    1313c
    14 #include "dimensions.h"
    15 #include "paramet.h"
    16 #include "logic.h"
    17 #include "comvert.h"
    18 #include "comconst.h"
     14      include "dimensions.h"
     15      include "paramet.h"
     16      include "logic.h"
     17      include "comvert.h"
    1918c
    2019c
     
    420419      USE parallel_lmdz
    421420      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi 
     421      USE comconst_mod, ONLY: pi
    422422      IMPLICIT NONE
    423423c
    424 #include "dimensions.h"
    425 #include "paramet.h"
    426 #include "logic.h"
    427 #include "comvert.h"
    428 #include "comconst.h"
    429 #include "comgeom.h"
     424      include "dimensions.h"
     425      include "paramet.h"
     426      include "logic.h"
     427      include "comvert.h"
     428      include "comgeom.h"
    430429c
    431430c
  • LMDZ5/trunk/libf/dyn3dmem/writedynav_loc.F

    r2475 r2597  
    1313      USE infotrac, ONLY : nqtot, ttext
    1414      use com_io_dyn_mod, only : histaveid,histvaveid,histuaveid
     15      USE comconst_mod, ONLY: cpp
    1516      implicit none
    1617
     
    4142C
    4243C   Declarations
    43 #include "dimensions.h"
    44 #include "paramet.h"
    45 #include "comconst.h"
    46 #include "comvert.h"
    47 #include "comgeom.h"
    48 #include "temps.h"
    49 #include "ener.h"
    50 #include "logic.h"
    51 #include "description.h"
    52 #include "serre.h"
    53 #include "iniprint.h"
     44      include "dimensions.h"
     45      include "paramet.h"
     46      include "comvert.h"
     47      include "comgeom.h"
     48      include "temps.h"
     49      include "ener.h"
     50      include "logic.h"
     51      include "description.h"
     52      include "serre.h"
     53      include "iniprint.h"
    5454
    5555C
  • LMDZ5/trunk/libf/dyn3dmem/writehist_loc.F

    r2475 r2597  
    1313      USE infotrac, ONLY : nqtot, ttext
    1414      use com_io_dyn_mod, only : histid,histvid,histuid
     15      USE comconst_mod, ONLY: cpp
    1516      implicit none
    1617
     
    4142C
    4243C   Declarations
    43 #include "dimensions.h"
    44 #include "paramet.h"
    45 #include "comconst.h"
    46 #include "comvert.h"
    47 #include "comgeom.h"
    48 #include "temps.h"
    49 #include "ener.h"
    50 #include "logic.h"
    51 #include "description.h"
    52 #include "serre.h"
    53 #include "iniprint.h"
     44      include "dimensions.h"
     45      include "paramet.h"
     46      include "comvert.h"
     47      include "comgeom.h"
     48      include "temps.h"
     49      include "ener.h"
     50      include "logic.h"
     51      include "description.h"
     52      include "serre.h"
     53      include "iniprint.h"
    5454
    5555C
Note: See TracChangeset for help on using the changeset viewer.