Changeset 5123 for LMDZ6/branches


Ignore:
Timestamp:
Jul 25, 2024, 8:45:50 AM (6 months ago)
Author:
abarral
Message:

Correct various minor mistakes from previous commits

Location:
LMDZ6/branches/Amaury_dev/libf
Files:
58 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/addfi.F90

    r5117 r5123  
    77  USE infotrac, ONLY: nqtot
    88  USE control_mod, ONLY: planet_type
     9  USE lmdz_ssum_scopy, ONLY: ssum
    910  IMPLICIT NONE
    1011  !
     
    7677  REAL, PARAMETER :: qtestw = 1.0e-15
    7778  REAL, PARAMETER :: qtestt = 1.0e-40
    78 
    79   REAL :: SSUM
    8079  !
    8180  !-----------------------------------------------------------------------
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F90

    r5105 r5123  
    66  USE logic_mod, ONLY: conser
    77  USE ener_mod, ONLY: gtot
     8  USE lmdz_ssum_scopy, ONLY: ssum
    89
    910  IMPLICIT NONE
     
    4647
    4748  INTEGER :: ij, l
    48 
    49   REAL :: SSUM
    5049
    5150  !-----------------------------------------------------------------------
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dissip.F90

    r5113 r5123  
    5353
    5454  INTEGER :: l, ij
    55 
    56   REAL :: SSUM
    5755
    5856  !-----------------------------------------------------------------------
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/integrd.F90

    r5119 r5123  
    1212  USE temps_mod, ONLY: dt
    1313  USE lmdz_iniprint, ONLY: lunout, prt_level
    14   USE lmdz_ssum_scopy, ONLY: scopy
     14  USE lmdz_ssum_scopy, ONLY: scopy, ssum
    1515
    1616  IMPLICIT NONE
     
    7373
    7474  INTEGER :: l, ij, iq, i, j
    75 
    76   REAL :: SSUM
    7775
    7876  !-----------------------------------------------------------------------
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F90

    r5119 r5123  
    2727  USE lmdz_description, ONLY: descript
    2828  USE lmdz_iniprint, ONLY: lunout, prt_level
    29   USE lmdz_ssum_scopy, ONLY: scopy
     29  USE lmdz_ssum_scopy, ONLY: scopy, ssum
    3030
    3131  IMPLICIT NONE
     
    117117  REAL :: time
    118118
    119   REAL :: SSUM
    120119  ! REAL finvmaold(ip1jmp1,llm)
    121120
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/qminimum.F90

    r5117 r5123  
    66  USE lmdz_strings, ONLY: strIdx
    77  USE lmdz_readTracFiles, ONLY: addPhase
     8  USE lmdz_ssum_scopy, ONLY: ssum
     9
    810  IMPLICIT none
    911  !
     
    3537  REAL :: zx_defau_diag(ip1jmp1, llm, 2)
    3638  REAL :: q_follow(ip1jmp1, llm, 2)
    37   !
    38   REAL :: SSUM
    3939  !
    4040  INTEGER :: imprim
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F90

    r5119 r5123  
    421421  USE infotrac, ONLY: nqtot, tracers, & ! CRisi
    422422          min_qParent, min_qMass, min_ratio ! MVals et CRisi
    423   !
     423  USE lmdz_ssum_scopy, ONLY: ssum
     424
     425 !
    424426  ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
    425427  !
     
    469471  REAL :: masseq(ip1jmp1, llm, nqtot), Ratio(ip1jmp1, llm, nqtot) ! CRisi
    470472  INTEGER :: ifils, iq2 ! CRisi
    471 
    472   !
    473   !
    474   REAL :: SSUM
    475473
    476474  DATA first/.TRUE./
  • LMDZ6/branches/Amaury_dev/libf/dyn3d/vlspltqs.F90

    r5119 r5123  
    219219  SAVE first
    220220
    221   REAL :: SSUM
    222221  REAL :: temps0, temps1, temps2, temps3, temps4, temps5
    223222  SAVE temps0, temps1, temps2, temps3, temps4, temps5
     
    522521
    523522  USE comconst_mod, ONLY: pi
     523  USE lmdz_ssum_scopy, ONLY: ssum
    524524
    525525  IMPLICIT NONE
     
    563563  REAL :: masseq(ip1jmp1, llm, nqtot), Ratio(ip1jmp1, llm, nqtot) ! CRisi
    564564  INTEGER :: ifils, iq2 ! CRisi
    565   !
    566   !
    567   REAL :: SSUM
    568565
    569566  DATA first/.TRUE./
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advn.f90

    r5118 r5123  
    1515  !   --------------------------------------------------------------------
    1616  USE lmdz_iniprint, ONLY: lunout, prt_level
     17  USE lmdz_ssum_scopy, ONLY: ssum
     18
    1719  IMPLICIT NONE
    1820  !
     
    4648  REAL :: zqh(ip1jmp1, llm), zqb(ip1jmp1, llm)
    4749  REAL :: temps0, temps1, temps2, temps3
    48   REAL :: ztemps1, ztemps2, ssum
     50  REAL :: ztemps1, ztemps2
    4951  save temps1, temps2, temps3
    5052  REAL :: zzpbar, zzw
     
    726728  !   --------------------------------------------------------------------
    727729  USE lmdz_iniprint, ONLY: lunout, prt_level
     730  USE lmdz_ssum_scopy, ONLY: ssum
     731
    728732  IMPLICIT NONE
    729733  !
     
    749753  REAL :: convpn, convps, convmpn, convmps, massen, masses
    750754  REAL :: zm, zq, zsigm, zsigp, zqm, zqp
    751   REAL :: ssum
    752755  REAL :: prec
    753756  save prec
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advy.f90

    r5116 r5123  
    9292  REAL :: z1(iim), zcos(iim), zsin(iim)
    9393  REAL :: smpn,smps,s0pn,s0ps
    94   REAL :: SSUM
    95   EXTERNAL SSUM
    9694  !
    9795  REAL :: sqi,sqf
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advyp.f90

    r5105 r5123  
    109109  REAL :: cy1(llm,ntra), cyLAT(llm,ntra)
    110110  REAL :: z1(iim), zcos(iim), zsin(iim)
    111   REAL :: SSUM
    112   EXTERNAL SSUM
    113111  !
    114112  REAL :: sqi,sqf
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.f90

    r5105 r5123  
    1818  ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
    1919  !
     20  USE lmdz_ssum_scopy, ONLY: ssum
     21
     22
    2023  IMPLICIT NONE
    2124  !
     
    2629  DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) , &
    2730        convfl( ip1jmp1,nbniv )
    28   !
    29   REAL :: SSUM
    30   !
    31   !
     31
     32
    3233  INCLUDE "comgeom.h"
    3334  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90

    r5118 r5123  
    3636  REAL :: x
    3737
    38   REAL :: SSUM
    39   EXTERNAL SSUM
    4038  REAL :: newsig
    4139  REAL :: dz0,dz1,nhaut,sig1,esig,csig,zz
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.f90

    r5105 r5123  
    1111  !          x et y  etant des composantes covariantes   ...
    1212  !  *********************************************************************
     13  USE lmdz_ssum_scopy, ONLY: ssum
     14
    1315  IMPLICIT NONE
    1416  !
     
    3739  REAL :: sumypn,sumyps
    3840  !    ...................................................................
    39   !
    40   REAL :: SSUM
    41   !
    42   !
     41
    4342  DO l = 1,klevel
    4443  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.f90

    r5105 r5123  
    1212  !          x et y  etant des composantes covariantes   ...
    1313  !  *********************************************************************
     14  USE lmdz_ssum_scopy, ONLY: ssum
     15
    1416  IMPLICIT NONE
    1517  !
     
    4042  INTEGER :: l,ij
    4143  !    ...................................................................
    42   !
    43   REAL :: SSUM
    44   !
     44
    4545  !
    4646  DO l = 1,klevel
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.f90

    r5106 r5123  
    1212  !  *********************************************************************
    1313  USE lmdz_filtreg, ONLY: filtreg
     14  USE lmdz_ssum_scopy, ONLY: ssum
     15
    1416  IMPLICIT NONE
    1517  !
     
    3840  REAL :: sumypn,sumyps
    3941  !    ...................................................................
    40   !
    41   REAL :: SSUM
    42   !
     42
    4343  !
    4444  DO l = 1,klevel
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.f90

    r5105 r5123  
    1 
    21! $Header$
    32
    4 SUBROUTINE divergst(klevel,x,y,div)
     3SUBROUTINE divergst(klevel, x, y, div)
     4  USE lmdz_ssum_scopy, ONLY: ssum
     5
    56  IMPLICIT NONE
    67  !
     
    1718  !   -------------------------------------------------------------------
    1819  !
     20
    1921  INCLUDE "dimensions.h"
    2022  INCLUDE "paramet.h"
     
    2224
    2325  INTEGER :: klevel
    24   REAL :: x( ip1jmp1,klevel ),y( ip1jm,klevel ),div( ip1jmp1,klevel )
    25   INTEGER :: ij,l,i
    26   REAL :: aiy1( iip1 ) , aiy2( iip1 )
    27   REAL :: sumypn,sumyps
    28 
    29   REAL :: SSUM
     26  REAL :: x(ip1jmp1, klevel), y(ip1jm, klevel), div(ip1jmp1, klevel)
     27  INTEGER :: ij, l, i
     28  REAL :: aiy1(iip1), aiy2(iip1)
     29  REAL :: sumypn, sumyps
    3030  !
    3131  !
    32   DO l = 1,klevel
    33   !
    34   DO ij = iip2, ip1jm - 1
    35   div( ij + 1, l ) = x(ij+1,l) - x(ij,l)+ y(ij-iim,l)-y(ij+1,l)
    36   END DO
    37   !
    38   ! ....  correction pour  div( 1,j,l)  ......
    39   ! ....   div(1,j,l)= div(iip1,j,l) ....
    40   !
    41   !DIR$ IVDEP
    42   DO ij = iip2,ip1jm,iip1
    43   div( ij,l ) = div( ij + iim,l )
    44   END DO
    45   !
    46   ! ....  calcul  aux poles  .....
    47   !
    48   !
    49   DO i  = 1,iim
    50   aiy1(i)= y(i,l)
    51   aiy2(i)= y(i+ip1jmi1,l)
    52   END DO
    53   sumypn = SSUM ( iim,aiy1,1 )
    54   sumyps = SSUM ( iim,aiy2,1 )
    55   DO i = 1,iip1
    56   div(     i    , l ) = - sumypn/iim
    57   div( i + ip1jm, l ) =   sumyps/iim
    58   END DO
    59   !
     32  DO l = 1, klevel
     33    !
     34    DO ij = iip2, ip1jm - 1
     35      div(ij + 1, l) = x(ij + 1, l) - x(ij, l) + y(ij - iim, l) - y(ij + 1, l)
     36    END DO
     37    !
     38    ! ....  correction pour  div( 1,j,l)  ......
     39    ! ....   div(1,j,l)= div(iip1,j,l) ....
     40    !
     41    !DIR$ IVDEP
     42    DO ij = iip2, ip1jm, iip1
     43      div(ij, l) = div(ij + iim, l)
     44    END DO
     45    !
     46    ! ....  calcul  aux poles  .....
     47    !
     48    !
     49    DO i = 1, iim
     50      aiy1(i) = y(i, l)
     51      aiy2(i) = y(i + ip1jmi1, l)
     52    END DO
     53    sumypn = SSUM (iim, aiy1, 1)
     54    sumyps = SSUM (iim, aiy2, 1)
     55    DO i = 1, iip1
     56      div(i, l) = - sumypn / iim
     57      div(i + ip1jm, l) = sumyps / iim
     58    END DO
     59    !
    6060  END DO
    6161  RETURN
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.f90

    r5117 r5123  
    2222  USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, &
    2323        alphax,alphay,taux,tauy,transx,transy,pxo,pyo
     24  USE lmdz_ssum_scopy, ONLY: ssum
    2425  IMPLICIT NONE
    2526  !
     
    5051  SAVE rlatu1,yprimu1,rlatu2,yprimu2,yprimv,yprimu
    5152  SAVE rlonm025,xprimm025,rlonp025,xprimp025
    52 
    53   REAL :: SSUM
    5453  !
    5554  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limx.f90

    r5117 r5123  
    4242  Logical :: extremum,first
    4343  save first
    44 
    45   REAL :: SSUM
    46   INTEGER :: ismax,ismin
    47   EXTERNAL  SSUM, ismin,ismax
    48 
    4944  data first/.TRUE./
    5045
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limy.f90

    r5117 r5123  
    1616  USE comconst_mod, ONLY: pi
    1717  USE lmdz_libmath, ONLY: ismax, ismin
     18  USE lmdz_ssum_scopy, ONLY: ssum
     19
    1820  IMPLICIT NONE
    1921  !
     
    4749  REAL :: coslon(iip1),coslondlon(iip1)
    4850  save sinlon,coslon,sinlondlon,coslondlon
    49   !
    50   !
    51   REAL :: SSUM
    52   EXTERNAL  SSUM
    5351
    5452  data first/.TRUE./
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limz.f90

    r5117 r5123  
    4242  Logical :: extremum,first
    4343  save first
    44 
    45   REAL :: SSUM
    46   INTEGER :: ismax,ismin
    47   EXTERNAL  SSUM, ismin,ismax
    4844
    4945  data first/.TRUE./
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.f90

    r5105 r5123  
    3131  REAL :: massemoyn, massemoys
    3232
    33   REAL :: SSUM
    3433  !
    3534  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.f90

    r5117 r5123  
    55
    66  USE comconst_mod, ONLY: pi, dtvr
     7  USE lmdz_ssum_scopy, ONLY: ssum
    78
    89  IMPLICIT NONE
     
    5556  REAL :: qmin,zq,pente_max
    5657  !
    57   REAL :: SSUM
    58   INTEGER :: ismax,ismin,lati,latf
    59   EXTERNAL  SSUM, ismin,ismax
     58  INTEGER :: lati,latf
    6059  LOGICAL :: first
    6160  save first
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.f90

    r5117 r5123  
    55
    66  USE comconst_mod, ONLY: pi
     7  USE lmdz_ssum_scopy, ONLY: ssum
    78
    89  IMPLICIT NONE
     
    5859  REAL :: dyn1,dyn2,dys1,dys2,qpn,qps,dqzpn,dqzps
    5960  REAL :: masn,mass
    60   !
    61   REAL :: SSUM
    62   INTEGER :: ismax,ismin
    63   EXTERNAL  SSUM, ismin,ismax
     61
    6462  LOGICAL :: first
    6563  save first
  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/sortvarc.f90

    r5119 r5123  
    1313  USE lmdz_filtreg, ONLY: filtreg
    1414  USE lmdz_iniprint, ONLY: lunout, prt_level
    15   USE lmdz_ssum_scopy, ONLY: scopy
     15  USE lmdz_ssum_scopy, ONLY: scopy, ssum
    1616
    1717  IMPLICIT NONE
     
    6464  INTEGER :: l, ij, imjmp1
    6565
    66   REAL :: SSUM
    6766  LOGICAL, SAVE :: firstcal = .TRUE.
    6867  CHARACTER(LEN = *), PARAMETER :: modname = "sortvarc"
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/addfi_loc.f90

    r5117 r5123  
    1 
    21! $Id$
    32
    43SUBROUTINE addfi_loc(pdt, leapf, forward, &
    5         pucov, pvcov, pteta, pq   , pps , &
    6         pdufi, pdvfi, pdhfi,pdqfi, pdpfi  )
     4        pucov, pvcov, pteta, pq, pps, &
     5        pdufi, pdvfi, pdhfi, pdqfi, pdpfi)
    76  USE parallel_lmdz
    87  USE infotrac, ONLY: nqtot
    98  USE control_mod, ONLY: planet_type
     9  USE lmdz_ssum_scopy, ONLY: ssum
     10
    1011  IMPLICIT NONE
    1112  !
     
    5354  !    -----------
    5455  !
    55   REAL,INTENT(IN) :: pdt ! time step for the integration (s)
    56   !
    57   REAL,INTENT(INOUT) :: pvcov(ijb_v:ije_v,llm) ! covariant meridional wind
    58   REAL,INTENT(INOUT) :: pucov(ijb_u:ije_u,llm) ! covariant zonal wind
    59   REAL,INTENT(INOUT) :: pteta(ijb_u:ije_u,llm) ! potential temperature
    60   REAL,INTENT(INOUT) :: pq(ijb_u:ije_u,llm,nqtot) ! tracers
    61   REAL,INTENT(INOUT) :: pps(ijb_u:ije_u) ! surface pressure (Pa)
     56  REAL, INTENT(IN) :: pdt ! time step for the integration (s)
     57  !
     58  REAL, INTENT(INOUT) :: pvcov(ijb_v:ije_v, llm) ! covariant meridional wind
     59  REAL, INTENT(INOUT) :: pucov(ijb_u:ije_u, llm) ! covariant zonal wind
     60  REAL, INTENT(INOUT) :: pteta(ijb_u:ije_u, llm) ! potential temperature
     61  REAL, INTENT(INOUT) :: pq(ijb_u:ije_u, llm, nqtot) ! tracers
     62  REAL, INTENT(INOUT) :: pps(ijb_u:ije_u) ! surface pressure (Pa)
    6263  ! respective tendencies (.../s) to add
    63   REAL,INTENT(IN) :: pdvfi(ijb_v:ije_v,llm)
    64   REAL,INTENT(IN) :: pdufi(ijb_u:ije_u,llm)
    65   REAL,INTENT(IN) :: pdqfi(ijb_u:ije_u,llm,nqtot)
    66   REAL,INTENT(IN) :: pdhfi(ijb_u:ije_u,llm)
    67   REAL,INTENT(IN) :: pdpfi(ijb_u:ije_u)
    68   !
    69   LOGICAL,INTENT(IN) :: leapf,forward ! not used
     64  REAL, INTENT(IN) :: pdvfi(ijb_v:ije_v, llm)
     65  REAL, INTENT(IN) :: pdufi(ijb_u:ije_u, llm)
     66  REAL, INTENT(IN) :: pdqfi(ijb_u:ije_u, llm, nqtot)
     67  REAL, INTENT(IN) :: pdhfi(ijb_u:ije_u, llm)
     68  REAL, INTENT(IN) :: pdpfi(ijb_u:ije_u)
     69  !
     70  LOGICAL, INTENT(IN) :: leapf, forward ! not used
    7071  !
    7172  !
     
    7374  !    -----------------
    7475  !
    75   REAL :: xpn(iim),xps(iim),tpn,tps
    76   INTEGER :: j,k,iq,ij
    77   REAL,PARAMETER :: qtestw = 1.0e-15
    78   REAL,PARAMETER :: qtestt = 1.0e-40
    79 
    80   REAL :: SSUM
    81   EXTERNAL SSUM
    82 
    83   INTEGER :: ijb,ije
     76  REAL :: xpn(iim), xps(iim), tpn, tps
     77  INTEGER :: j, k, iq, ij
     78  REAL, PARAMETER :: qtestw = 1.0e-15
     79  REAL, PARAMETER :: qtestt = 1.0e-40
     80
     81  INTEGER :: ijb, ije
    8482  !
    8583  !-----------------------------------------------------------------------
    8684
    87   ijb=ij_begin
    88   ije=ij_end
    89 
    90 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    91   DO k = 1,llm
    92      DO j = ijb,ije
    93         pteta(j,k)= pteta(j,k) + pdhfi(j,k) * pdt
    94      ENDDO
    95   ENDDO
    96 !$OMP END DO NOWAIT
     85  ijb = ij_begin
     86  ije = ij_end
     87
     88  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     89  DO k = 1, llm
     90    DO j = ijb, ije
     91      pteta(j, k) = pteta(j, k) + pdhfi(j, k) * pdt
     92    ENDDO
     93  ENDDO
     94  !$OMP END DO NOWAIT
    9795
    9896  IF (pole_nord) THEN
    99 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    100     DO  k    = 1, llm
    101      DO  ij  = 1, iim
    102        xpn(ij) = aire(   ij   ) * pteta(  ij    ,k)
    103      ENDDO
    104      tpn      = SSUM(iim,xpn,1)/ apoln
    105 
    106      DO ij  = 1, iip1
    107        pteta(   ij   ,k) = tpn
    108      ENDDO
    109    ENDDO
    110 !$OMP END DO NOWAIT
     97    !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     98    DO  k = 1, llm
     99      DO  ij = 1, iim
     100        xpn(ij) = aire(ij) * pteta(ij, k)
     101      ENDDO
     102      tpn = SSUM(iim, xpn, 1) / apoln
     103
     104      DO ij = 1, iip1
     105        pteta(ij, k) = tpn
     106      ENDDO
     107    ENDDO
     108    !$OMP END DO NOWAIT
    111109  ENDIF
    112110
    113111  IF (pole_sud) THEN
    114 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    115     DO  k    = 1, llm
    116      DO  ij  = 1, iim
    117        xps(ij) = aire(ij+ip1jm) * pteta(ij+ip1jm,k)
    118      ENDDO
    119      tps      = SSUM(iim,xps,1)/ apols
    120 
    121      DO ij  = 1, iip1
    122        pteta(ij+ip1jm,k) = tps
    123      ENDDO
    124    ENDDO
    125 !$OMP END DO NOWAIT
    126   ENDIF
    127   !
    128 
    129   ijb=ij_begin
    130   ije=ij_end
    131   IF (pole_nord) ijb=ij_begin+iip1
    132   IF (pole_sud)  ije=ij_end-iip1
    133 
    134 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    135   DO k = 1,llm
    136      DO j = ijb,ije
    137         pucov(j,k)= pucov(j,k) + pdufi(j,k) * pdt
    138      ENDDO
    139   ENDDO
    140 !$OMP END DO NOWAIT
    141 
    142   IF (pole_nord) ijb=ij_begin
    143 
    144 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    145   DO k = 1,llm
    146      DO j = ijb,ije
    147         pvcov(j,k)= pvcov(j,k) + pdvfi(j,k) * pdt
    148      ENDDO
    149   ENDDO
    150 !$OMP END DO NOWAIT
    151 
    152   !
    153   IF (pole_sud)  ije=ij_end
    154 !$OMP MASTER
    155   DO j = ijb,ije
    156      pps(j) = pps(j) + pdpfi(j) * pdt
    157   ENDDO
    158 !$OMP END MASTER
     112    !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     113    DO  k = 1, llm
     114      DO  ij = 1, iim
     115        xps(ij) = aire(ij + ip1jm) * pteta(ij + ip1jm, k)
     116      ENDDO
     117      tps = SSUM(iim, xps, 1) / apols
     118
     119      DO ij = 1, iip1
     120        pteta(ij + ip1jm, k) = tps
     121      ENDDO
     122    ENDDO
     123    !$OMP END DO NOWAIT
     124  ENDIF
     125  !
     126
     127  ijb = ij_begin
     128  ije = ij_end
     129  IF (pole_nord) ijb = ij_begin + iip1
     130  IF (pole_sud)  ije = ij_end - iip1
     131
     132  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     133  DO k = 1, llm
     134    DO j = ijb, ije
     135      pucov(j, k) = pucov(j, k) + pdufi(j, k) * pdt
     136    ENDDO
     137  ENDDO
     138  !$OMP END DO NOWAIT
     139
     140  IF (pole_nord) ijb = ij_begin
     141
     142  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     143  DO k = 1, llm
     144    DO j = ijb, ije
     145      pvcov(j, k) = pvcov(j, k) + pdvfi(j, k) * pdt
     146    ENDDO
     147  ENDDO
     148  !$OMP END DO NOWAIT
     149
     150  !
     151  IF (pole_sud)  ije = ij_end
     152  !$OMP MASTER
     153  DO j = ijb, ije
     154    pps(j) = pps(j) + pdpfi(j) * pdt
     155  ENDDO
     156  !$OMP END MASTER
    159157
    160158  IF (planet_type=="earth") THEN
    161   ! earth case, special treatment for first 2 tracers (water)
    162   DO iq = 1, 2
    163 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    164      DO k = 1,llm
    165         DO j = ijb,ije
    166            pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt
    167            pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestw )
    168         ENDDO
    169      ENDDO
    170 !$OMP END DO NOWAIT
    171   ENDDO
    172 
    173   DO iq = 3, nqtot
    174 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    175      DO k = 1,llm
    176         DO j = ijb,ije
    177            pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt
    178            pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestt )
    179         ENDDO
    180      ENDDO
    181 !$OMP END DO NOWAIT
    182   ENDDO
     159    ! earth case, special treatment for first 2 tracers (water)
     160    DO iq = 1, 2
     161      !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     162      DO k = 1, llm
     163        DO j = ijb, ije
     164          pq(j, k, iq) = pq(j, k, iq) + pdqfi(j, k, iq) * pdt
     165          pq(j, k, iq) = AMAX1(pq(j, k, iq), qtestw)
     166        ENDDO
     167      ENDDO
     168      !$OMP END DO NOWAIT
     169    ENDDO
     170
     171    DO iq = 3, nqtot
     172      !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     173      DO k = 1, llm
     174        DO j = ijb, ije
     175          pq(j, k, iq) = pq(j, k, iq) + pdqfi(j, k, iq) * pdt
     176          pq(j, k, iq) = AMAX1(pq(j, k, iq), qtestt)
     177        ENDDO
     178      ENDDO
     179      !$OMP END DO NOWAIT
     180    ENDDO
    183181  else
    184   ! general case, treat all tracers equally)
    185    DO iq = 1, nqtot
    186 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    187      DO k = 1,llm
    188         DO j = ijb,ije
    189            pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt
    190            pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestt )
    191         ENDDO
    192      ENDDO
    193 !$OMP END DO NOWAIT
    194    ENDDO
     182    ! general case, treat all tracers equally)
     183    DO iq = 1, nqtot
     184      !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     185      DO k = 1, llm
     186        DO j = ijb, ije
     187          pq(j, k, iq) = pq(j, k, iq) + pdqfi(j, k, iq) * pdt
     188          pq(j, k, iq) = AMAX1(pq(j, k, iq), qtestt)
     189        ENDDO
     190      ENDDO
     191      !$OMP END DO NOWAIT
     192    ENDDO
    195193  ENDIF ! of if (planet_type=="earth")
    196194
    197 !$OMP MASTER
     195  !$OMP MASTER
    198196  IF (pole_nord) THEN
    199     DO  ij   = 1, iim
    200       xpn(ij) = aire(   ij   ) * pps(  ij     )
    201     ENDDO
    202 
    203     tpn      = SSUM(iim,xpn,1)/apoln
    204 
    205     DO ij   = 1, iip1
    206       pps (   ij     ) = tpn
     197    DO  ij = 1, iim
     198      xpn(ij) = aire(ij) * pps(ij)
     199    ENDDO
     200
     201    tpn = SSUM(iim, xpn, 1) / apoln
     202
     203    DO ij = 1, iip1
     204      pps (ij) = tpn
    207205    ENDDO
    208206
     
    210208
    211209  IF (pole_sud) THEN
    212     DO  ij   = 1, iim
    213       xps(ij) = aire(ij+ip1jm) * pps(ij+ip1jm )
    214     ENDDO
    215 
    216     tps      = SSUM(iim,xps,1)/apols
    217 
    218     DO ij   = 1, iip1
    219       pps ( ij+ip1jm ) = tps
    220     ENDDO
    221 
    222   ENDIF
    223 !$OMP END MASTER
     210    DO  ij = 1, iim
     211      xps(ij) = aire(ij + ip1jm) * pps(ij + ip1jm)
     212    ENDDO
     213
     214    tps = SSUM(iim, xps, 1) / apols
     215
     216    DO ij = 1, iip1
     217      pps (ij + ip1jm) = tps
     218    ENDDO
     219
     220  ENDIF
     221  !$OMP END MASTER
    224222
    225223  IF (pole_nord) THEN
    226224    DO iq = 1, nqtot
    227 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    228       DO  k    = 1, llm
    229         DO  ij   = 1, iim
    230           xpn(ij) = aire(   ij   ) * pq(  ij    ,k,iq)
    231         ENDDO
    232         tpn      = SSUM(iim,xpn,1)/apoln
    233 
    234         DO ij   = 1, iip1
    235           pq (   ij   ,k,iq) = tpn
    236         ENDDO
    237       ENDDO
    238 !$OMP END DO NOWAIT     
     225      !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     226      DO  k = 1, llm
     227        DO  ij = 1, iim
     228          xpn(ij) = aire(ij) * pq(ij, k, iq)
     229        ENDDO
     230        tpn = SSUM(iim, xpn, 1) / apoln
     231
     232        DO ij = 1, iip1
     233          pq (ij, k, iq) = tpn
     234        ENDDO
     235      ENDDO
     236      !$OMP END DO NOWAIT
    239237    ENDDO
    240238  ENDIF
     
    242240  IF (pole_sud) THEN
    243241    DO iq = 1, nqtot
    244 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    245       DO  k    = 1, llm
    246         DO  ij   = 1, iim
    247           xps(ij) = aire(ij+ip1jm) * pq(ij+ip1jm,k,iq)
    248         ENDDO
    249         tps      = SSUM(iim,xps,1)/apols
    250 
    251         DO ij   = 1, iip1
    252           pq (ij+ip1jm,k,iq)  = tps
    253         ENDDO
    254       ENDDO
    255 !$OMP END DO NOWAIT     
    256     ENDDO
    257   ENDIF
    258 
    259 
     242      !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     243      DO  k = 1, llm
     244        DO  ij = 1, iim
     245          xps(ij) = aire(ij + ip1jm) * pq(ij + ip1jm, k, iq)
     246        ENDDO
     247        tps = SSUM(iim, xps, 1) / apols
     248
     249        DO ij = 1, iip1
     250          pq (ij + ip1jm, k, iq) = tps
     251        ENDDO
     252      ENDDO
     253      !$OMP END DO NOWAIT
     254    ENDDO
     255  ENDIF
    260256
    261257END SUBROUTINE addfi_loc
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/advect_new_loc.f90

    r5117 r5123  
    5151
    5252  INTEGER :: ij, l, ijb, ije
    53   EXTERNAL  SSUM
    54   REAL :: SSUM
    55 
    56 
    5753
    5854  !-----------------------------------------------------------------------
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/call_dissip_mod.f90

    r5117 r5123  
    8585    USE write_field_loc
    8686    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO
     87    USE lmdz_ssum_scopy, ONLY: ssum
     88
    8789
    8890    IMPLICIT NONE
     
    9799    REAL :: tpn, tps
    98100
    99     REAL  SSUM
    100101    LOGICAL, PARAMETER :: dissip_conservative = .TRUE.
    101102    TYPE(Request), SAVE :: Request_dissip
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/convflu_loc.f90

    r5105 r5123  
    1616  !
    1717  USE parallel_lmdz
     18  USE lmdz_ssum_scopy, ONLY: ssum
     19
    1820  IMPLICIT NONE
    1921  !
     
    2628  !
    2729  INTEGER :: ijb,ije
    28   EXTERNAL   SSUM
    29   REAL :: SSUM
    3030  !
    3131  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dissip_loc.f90

    r5117 r5123  
    5656
    5757  INTEGER :: l, ij
    58 
    59   REAL :: SSUM
    6058  INTEGER :: ijb, ije
    6159
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.f90

    r5117 r5123  
    1010  !  *********************************************************************
    1111  USE parallel_lmdz
     12  USE lmdz_ssum_scopy, ONLY: ssum
     13
    1214  IMPLICIT NONE
    1315  !
     
    3941  INTEGER :: l,ij
    4042  !    ...................................................................
    41   !
    42   EXTERNAL  SSUM
    43   REAL :: SSUM
    4443  INTEGER :: ijb,ije,jjb,jje
    4544  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.f90

    r5117 r5123  
    99  !  *********************************************************************
    1010  USE parallel_lmdz
     11  USE lmdz_ssum_scopy, ONLY: ssum
     12
    1113  IMPLICIT NONE
    1214  !
     
    3638  INTEGER :: ijb,ije
    3739  !    ...................................................................
    38   !
    39   EXTERNAL  SSUM
    40   REAL :: SSUM
    4140  !
    4241  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/divergf_loc.f90

    r5117 r5123  
    1010  USE parallel_lmdz
    1111  USE lmdz_filtreg_p
     12  USE lmdz_ssum_scopy, ONLY: ssum
     13
    1214  IMPLICIT NONE
    1315  !
     
    3739  REAL :: sumypn,sumyps
    3840  !    ...................................................................
    39   !
    40   EXTERNAL  SSUM
    41   REAL :: SSUM
    4241  INTEGER :: ijb,ije,jjb,jje
    4342  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.f90

    r5118 r5123  
    1717  USE lmdz_strings, ONLY: int2str
    1818  USE lmdz_iniprint, ONLY: lunout, prt_level
     19  USE lmdz_ssum_scopy, ONLY: ssum
    1920
    2021  IMPLICIT NONE
     
    7879  INTEGER :: l,ij,iq,i,j
    7980
    80   REAL :: SSUM
    81   EXTERNAL SSUM
    8281  INTEGER :: ijb,ije,jjb,jje
    8382  LOGICAL :: checksum
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/leapfrog_loc.F90

    r5118 r5123  
    134134  REAL :: time
    135135
    136   REAL :: SSUM
    137136  ! REAL,SAVE,ALLOCATABLE :: finvmaold(:,:)
    138137
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/massdair_loc.f90

    r5117 r5123  
    3030  REAL :: massemoyn, massemoys
    3131
    32   REAL :: SSUM
    33   EXTERNAL SSUM
    3432  !
    3533  !
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/qminimum_loc.f90

    r5118 r5123  
    3838  REAL :: zx_defau_diag(ijb_u:ije_u, llm, 2)
    3939  REAL :: q_follow(ijb_u:ije_u, llm, 2)
    40   !
    41   REAL :: SSUM
    42   EXTERNAL SSUM
    4340  !
    4441  INTEGER :: imprim
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90

    r5118 r5123  
    4646
    4747  Logical :: extremum
    48 
    49   REAL :: SSUM
    5048
    5149  REAL :: z1, z2, z3
     
    384382          min_qParent, min_qMass, min_ratio ! MVals et CRisi
    385383  USE comconst_mod, ONLY: pi
     384  USE lmdz_ssum_scopy, ONLY: ssum
     385
    386386  IMPLICIT NONE
    387387  !
     
    428428  REAL :: Ratio(ijb_u:ije_u, llm, nqtot) ! CRisi
    429429  INTEGER :: ifils, iq2 ! CRisi
    430   !
    431   !
    432   REAL :: SSUM
    433   EXTERNAL  SSUM
    434430
    435431  DATA first/.TRUE./
     
    850846  !$OMP THREADPRIVATE(temps0,temps1,temps2,temps3,temps4,temps5)
    851847
    852   REAL :: SSUM
    853 
    854848  DATA temps0, temps1, temps2, temps3, temps4, temps5/0., 0., 0., 0., 0., 0./
    855849  INTEGER :: ijb, ije, ijb_x, ije_x
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_loc.f90

    r5118 r5123  
    4040  REAL :: Ratio(ijb_u:ije_u, llm, nqtot) ! CRisi
    4141  INTEGER :: ifils, iq2 ! CRisi
    42 
    43   REAL :: SSUM
    4442
    4543  INTEGER :: ijb, ije, ijb_x, ije_x
     
    392390  USE comconst_mod, ONLY: pi
    393391  USE lmdz_iniprint, ONLY: lunout, prt_level
     392  USE lmdz_ssum_scopy, ONLY: ssum
     393
    394394  IMPLICIT NONE
    395395  !
     
    437437  INTEGER :: ijb, ije
    438438  INTEGER :: ijbm, ijem
    439 
    440   REAL :: ssum
    441439
    442440  ijb = ij_begin - 2 * iip1
  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.f90

    r5119 r5123  
    3535  USE comvert_mod, ONLY: preff, presnivs
    3636  USE lmdz_iniprint, ONLY: lunout, prt_level
    37   USE lmdz_ssum_scopy, ONLY: scopy
     37  USE lmdz_ssum_scopy, ONLY: scopy, ssum
    3838
    3939  IMPLICIT NONE
     
    164164  REAL :: flxwfi(ngridmx, llm)  ! Flux de masse verticale sur la grille physiq
    165165  !
    166 
    167   REAL :: SSUM
    168 
    169166  LOGICAL, SAVE :: firstcal = .TRUE., debut = .TRUE.
    170167  ! REAL rdayvrai
  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90

    r5118 r5123  
    5151    USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi
    5252    USE lmdz_iniprint, ONLY: lunout, prt_level
     53    USE lmdz_ssum_scopy, ONLY: ssum
     54
    5355
    5456    !=======================================================================
     
    216218  REAL :: zsinbis(iim),zcosbis(iim),z1bis(iim)
    217219  REAL :: unskap, pksurcp
    218 
    219   REAL :: SSUM
    220220
    221221  LOGICAL,SAVE :: firstcal=.TRUE., debut=.TRUE.
  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

    r5118 r5123  
    460460    var3d(iml, :, :) = var3d(1, :, :)
    461461
    462     DO il = 1, lml
    463       CALL minmax(iml * jml, var3d(1, 1, il), chmin, chmax)
    464       WRITE(lunout, *)' ' // TRIM(var) // '  min max l ', il, chmin, chmax
    465     END DO
    466 
    467462  END SUBROUTINE start_inter_3d
    468463
  • LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_filtreg.F90

    r5120 r5123  
    330330    !
    331331    USE lmdz_coefils, ONLY: sddv, sddu, unsddu, unsddv, eignfnv, eignfnu
     332    USE lmdz_ssum_scopy, ONLY: ssum
    332333    IMPLICIT NONE
    333334    !
     
    341342    REAL :: pi
    342343    INTEGER :: i, j, k, imm1, nrot
    343     EXTERNAL SSUM
    344     REAL :: SSUM
    345344    !
    346345
     
    557556
    558557  SUBROUTINE acc(vec, d, im)
     558    USE lmdz_ssum_scopy, ONLY: ssum
    559559    IMPLICIT NONE
    560560    INTEGER :: im
     
    562562    INTEGER :: i, j
    563563    REAL :: sum
    564     REAL, external :: ssum
    565564    do j = 1, im
    566565      do i = 1, im
  • LMDZ6/branches/Amaury_dev/libf/grid/lmdz_fxy_prim.f90

    r5122 r5123  
    1 ! This module replaces grid/fxy_reg.h
     1! This module replaces grid/fxy_prim.h
    22
    33MODULE lmdz_fxy_prim
  • LMDZ6/branches/Amaury_dev/libf/grid/lmdz_fxy_sin.f90

    r5122 r5123  
    1 ! This module replaces grid/fxy_reg.h
     1! This module replaces grid/fxy_sin.h
    22
    33MODULE lmdz_fxy_sin
  • LMDZ6/branches/Amaury_dev/libf/misc/lmdz_libmath_pch.f90

    r5117 r5123  
    915915    REAL :: G, HALF, ONE, STEMP(3), THREE, TWO, XTEMP(4), ZERO
    916916    SAVE ZERO, HALF, ONE, TWO, THREE
    917     REAL :: PCHDF
    918917    !
    919918    DATA  ZERO /0./, HALF /0.5/, ONE /1./, TWO /2./, THREE /3./
  • LMDZ6/branches/Amaury_dev/libf/misc/lmdz_wxios.F90

    r5119 r5123  
    443443    USE lmdz_grid_phy, ONLY: nvertex, klon_glo
    444444    USE lmdz_physical_constants, ONLY: PI
     445    USE lmdz_phys_para
    445446    USE lmdz_ioipsl_getin_p, ONLY: getin_p
    446447    IMPLICIT NONE
  • LMDZ6/branches/Amaury_dev/libf/misc/lmdz_xer.f90

    r5117 r5123  
    201201    CHARACTER(len = 20) :: LFIRST
    202202    INTEGER :: NERR, LEVEL, LKNTRL
    203     INTEGER :: J4SAVE, MAXMES, KDUMMY, I, KOUNT, LERR, LLEVEL
     203    INTEGER :: MAXMES, KDUMMY, I, KOUNT, LERR, LLEVEL
    204204    INTEGER :: MKNTRL, LTEMP
    205205    !***FIRST EXECUTABLE STATEMENT  XERMSG
     
    549549    !***END PROLOGUE  XGETUA
    550550    DIMENSION IUNITA(5)
    551     INTEGER :: IUNITA, N, J4SAVE, INDEX, I
     551    INTEGER :: IUNITA, N, INDEX, I
    552552    !***FIRST EXECUTABLE STATEMENT  XGETUA
    553553    N = J4SAVE(5, 0, .FALSE.)
     
    630630    INTEGER :: NERR, LEVEL, KONTRL
    631631    INTEGER :: NERTAB, LEVTAB, KOUNT, KOUNTX, NMSG
    632     INTEGER :: KFLAG, ICOUNT, NUNIT, KUNIT, IUNIT, I1MACH, I
     632    INTEGER :: KFLAG, ICOUNT, NUNIT, KUNIT, IUNIT, I
    633633    !***FIRST EXECUTABLE STATEMENT  XERSVE
    634634    !
     
    803803    CHARACTER(len = 2) :: NEWLIN
    804804    PARAMETER (NEWLIN = '$$')
    805     INTEGER :: N, I1MACH, I, LPREF, LWRAP, LENMSG, NEXTC
     805    INTEGER :: N, I, LPREF, LWRAP, LENMSG, NEXTC
    806806    INTEGER :: LPIECE, IDELTA
    807807    !***FIRST EXECUTABLE STATEMENT  XERPRN
     
    10561056
    10571057  !DECK J4SAVE
    1058   FUNCTION J4SAVE(IWHICH, IVALUE, ISET)
     1058  INTEGER FUNCTION J4SAVE(IWHICH, IVALUE, ISET)
    10591059    IMPLICIT NONE
    10601060    !***BEGIN PROLOGUE  J4SAVE
     
    11161116    DATA IPARAM(5)/1/
    11171117    DATA IPARAM(6), IPARAM(7), IPARAM(8), IPARAM(9)/0, 0, 0, 0/
    1118     INTEGER :: J4SAVE, IWHICH, IVALUE
     1118    INTEGER :: IWHICH, IVALUE
    11191119    !***FIRST EXECUTABLE STATEMENT  J4SAVE
    11201120    J4SAVE = IPARAM(IWHICH)
  • LMDZ6/branches/Amaury_dev/libf/phy_common/lmdz_abort_physic.F90

    r5117 r5123  
    77
    88    USE IOIPSL
    9     USE lmdz_phys_para
     9    USE lmdz_phys_mpi_data
    1010    USE lmdz_print_control, ONLY: lunout
    1111    IMPLICIT NONE
  • LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phys_output_write_spl_mod.F90

    r5118 r5123  
    14861486                                                                                                                                                                                                                                                                                                                                    DO i = 1, klon
    14871487                                                                                                                                                                                                                                                                                                                                      IF (tnondef(i, k, iff - 6).NE.missing_val) THEN
    1488                                                                                                                                                                                                                                                                                                                                         IF (
    1489                                                                                                                                                                                                                                                                                                                                           freq_outNMC(iff-6).LT.0) THEN
     1488                                                                                                                                                                                                                                                                                                                                        IF (                                                                                                                                                                                                                                                                                                                                          freq_outNMC(iff-6).LT.0) THEN
    14901489                                                                                                                                                                                                                                                                                                                                          freq_moyNMC(iff-6) =(mth_len*un_jour)/freq_calNMC(iff-6)
    14911490                                                                                                                                                                                                                                                                                                                                          ELSE
  • LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90

    r5117 r5123  
    33
    44MODULE iophy
     5    USE lmdz_abort_physic, ONLY: abort_physic
    56
    67! abd  REAL,PRIVATE,ALLOCATABLE,DIMENSION(:),save :: io_lat
     
    809810    USE lmdz_print_control, ONLY: prt_level,lunout
    810811    USE lmdz_grid_phy, ONLY: nbp_lon, nbp_lat
    811     USE lmdz_abort_physic, ONLY: abort_physic
    812812
    813813    IMPLICIT NONE
     
    955955  USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios
    956956  USE lmdz_print_control, ONLY: lunout, prt_level
    957   USE lmdz_abort_physic, ONLY: abort_physic
    958957
    959958  IMPLICIT NONE
     
    11681167  USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios
    11691168  USE lmdz_print_control, ONLY: prt_level,lunout
    1170   USE lmdz_abort_physic, ONLY: abort_physic
    11711169
    11721170  IMPLICIT NONE
     
    13651363  USE lmdz_xios, ONLY: xios_send_field
    13661364  USE lmdz_print_control, ONLY: prt_level,lunout
    1367   USE lmdz_abort_physic, ONLY: abort_physic
    13681365
    13691366  IMPLICIT NONE
     
    14491446  USE lmdz_grid_phy, ONLY: nbp_lon, nbp_lat, grid_type, unstructured
    14501447  USE lmdz_print_control, ONLY: prt_level,lunout
    1451   USE lmdz_abort_physic, ONLY: abort_physic
    14521448
    14531449  IMPLICIT NONE
  • LMDZ6/branches/Amaury_dev/libf/phylmd/limit_read_mod.F90

    r5116 r5123  
    170170    USE phys_cal_mod, ONLY: calend, year_len
    171171    USE lmdz_print_control, ONLY: lunout, prt_level
    172     USE lmdz_xios, ONLY: xios_recv_field, using_xiosw
     172    USE lmdz_xios, ONLY: xios_recv_field, using_xios
    173173    USE lmdz_abort_physic, ONLY: abort_physic
    174174   
  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_main.F90

    r5117 r5123  
    77CONTAINS
    88
    9   SUBROUTINE thermcell_main
    10     (itap, ngrid, nlay, ptimestep  &
     9  SUBROUTINE thermcell_main(itap, ngrid, nlay, ptimestep  &
    1110            , pplay,pplev, pphi, debut  &
    1211            , puwind, pvwind,ptemp, p_o, ptemp_env, po_env  &
  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_old.F90

    r5119 r5123  
    11MODULE lmdz_thermcell_old
     2      USE lmdz_abort_physic, ONLY: abort_physic
     3
    24CONTAINS
    35
     
    1012    USE lmdz_thermcell_dv2, ONLY: thermcell_dv2
    1113    USE lmdz_thermcell_dq, ONLY: thermcell_dq
    12     USE lmdz_abort_physic, ONLY: abort_physic
    1314    IMPLICIT NONE
    1415
  • LMDZ6/branches/Amaury_dev/libf/phylmd/pbl_surface_mod.F90

    r5117 r5123  
    256256  !****************************************************************************************
    257257
    258   SUBROUTINE pbl_surface
    259     ( &
     258  SUBROUTINE pbl_surface( &
    260259            dtime, date0, itap, jour, &
    261260            debut,     lafin, &
     
    41864185  !****************************************************************************************
    41874186
    4188   SUBROUTINE pbl_surface_final
    4189     (fder_rst, snow_rst, qsurf_rst, ftsoil_rst &
     4187  SUBROUTINE pbl_surface_final(fder_rst, snow_rst, qsurf_rst, ftsoil_rst &
    41904188#ifdef ISO
    41914189       ,xtsnow_rst,Rland_ice_rst &
     
    42594257
    42604258  !albedo SB >>>
    4261   SUBROUTINE pbl_surface_newfrac
    4262     (itime, pctsrf_new, pctsrf_old, &
     4259  SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &
    42634260            evap, z0m, z0h, agesno, &
    42644261            tsurf, alb_dir, alb_dif, ustar, u10m, v10m, tke &
  • LMDZ6/branches/Amaury_dev/libf/phylmd/slab_heat_transp_mod.F90

    r5119 r5123  
    986986    ! computes the divergence of a vector field of components
    987987    ! x,y. x and y being covariant components
     988    USE lmdz_ssum_scopy, ONLY: ssum
     989
    988990    IMPLICIT NONE
    989991
     
    9981000    REAL :: aiy1(nbp_lon + 1), aiy2(nbp_lon + 1)
    9991001    REAL :: sumypn, sumyps
    1000     REAL, EXTERNAL :: SSUM
    10011002
    10021003    iim = nbp_lon
  • LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90

    r5117 r5123  
    77  !****************************************************************************************
    88
    9   SUBROUTINE surf_landice
    10     (itime, dtime, knon, knindex, &
     9  SUBROUTINE surf_landice(itime, dtime, knon, knindex, &
    1110            rlon, rlat, debut, lafin, &
    1211            rmu0, lwdownm, albedo, pphi1, &
  • LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90

    r5117 r5123  
    99  !******************************************************************************
    1010
    11   SUBROUTINE surf_ocean
    12     (rlon, rlat, swnet, lwnet, alb1, &
     11  SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, &
    1312            windsp, rmu0, fder, tsurf_in, &
    1413            itime, dtime, jour, knon, knindex, &
Note: See TracChangeset for help on using the changeset viewer.