Changeset 5123 for LMDZ6/branches
- Timestamp:
- Jul 25, 2024, 8:45:50 AM (6 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf
- Files:
-
- 58 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d/addfi.F90
r5117 r5123 7 7 USE infotrac, ONLY: nqtot 8 8 USE control_mod, ONLY: planet_type 9 USE lmdz_ssum_scopy, ONLY: ssum 9 10 IMPLICIT NONE 10 11 ! … … 76 77 REAL, PARAMETER :: qtestw = 1.0e-15 77 78 REAL, PARAMETER :: qtestt = 1.0e-40 78 79 REAL :: SSUM80 79 ! 81 80 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F90
r5105 r5123 6 6 USE logic_mod, ONLY: conser 7 7 USE ener_mod, ONLY: gtot 8 USE lmdz_ssum_scopy, ONLY: ssum 8 9 9 10 IMPLICIT NONE … … 46 47 47 48 INTEGER :: ij, l 48 49 REAL :: SSUM50 49 51 50 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/dyn3d/dissip.F90
r5113 r5123 53 53 54 54 INTEGER :: l, ij 55 56 REAL :: SSUM57 55 58 56 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/dyn3d/integrd.F90
r5119 r5123 12 12 USE temps_mod, ONLY: dt 13 13 USE lmdz_iniprint, ONLY: lunout, prt_level 14 USE lmdz_ssum_scopy, ONLY: scopy 14 USE lmdz_ssum_scopy, ONLY: scopy, ssum 15 15 16 16 IMPLICIT NONE … … 73 73 74 74 INTEGER :: l, ij, iq, i, j 75 76 REAL :: SSUM77 75 78 76 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F90
r5119 r5123 27 27 USE lmdz_description, ONLY: descript 28 28 USE lmdz_iniprint, ONLY: lunout, prt_level 29 USE lmdz_ssum_scopy, ONLY: scopy 29 USE lmdz_ssum_scopy, ONLY: scopy, ssum 30 30 31 31 IMPLICIT NONE … … 117 117 REAL :: time 118 118 119 REAL :: SSUM120 119 ! REAL finvmaold(ip1jmp1,llm) 121 120 -
LMDZ6/branches/Amaury_dev/libf/dyn3d/qminimum.F90
r5117 r5123 6 6 USE lmdz_strings, ONLY: strIdx 7 7 USE lmdz_readTracFiles, ONLY: addPhase 8 USE lmdz_ssum_scopy, ONLY: ssum 9 8 10 IMPLICIT none 9 11 ! … … 35 37 REAL :: zx_defau_diag(ip1jmp1, llm, 2) 36 38 REAL :: q_follow(ip1jmp1, llm, 2) 37 !38 REAL :: SSUM39 39 ! 40 40 INTEGER :: imprim -
LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F90
r5119 r5123 421 421 USE infotrac, ONLY: nqtot, tracers, & ! CRisi 422 422 min_qParent, min_qMass, min_ratio ! MVals et CRisi 423 ! 423 USE lmdz_ssum_scopy, ONLY: ssum 424 425 ! 424 426 ! Auteurs: P.Le Van, F.Hourdin, F.Forget 425 427 ! … … 469 471 REAL :: masseq(ip1jmp1, llm, nqtot), Ratio(ip1jmp1, llm, nqtot) ! CRisi 470 472 INTEGER :: ifils, iq2 ! CRisi 471 472 !473 !474 REAL :: SSUM475 473 476 474 DATA first/.TRUE./ -
LMDZ6/branches/Amaury_dev/libf/dyn3d/vlspltqs.F90
r5119 r5123 219 219 SAVE first 220 220 221 REAL :: SSUM222 221 REAL :: temps0, temps1, temps2, temps3, temps4, temps5 223 222 SAVE temps0, temps1, temps2, temps3, temps4, temps5 … … 522 521 523 522 USE comconst_mod, ONLY: pi 523 USE lmdz_ssum_scopy, ONLY: ssum 524 524 525 525 IMPLICIT NONE … … 563 563 REAL :: masseq(ip1jmp1, llm, nqtot), Ratio(ip1jmp1, llm, nqtot) ! CRisi 564 564 INTEGER :: ifils, iq2 ! CRisi 565 !566 !567 REAL :: SSUM568 565 569 566 DATA first/.TRUE./ -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advn.f90
r5118 r5123 15 15 ! -------------------------------------------------------------------- 16 16 USE lmdz_iniprint, ONLY: lunout, prt_level 17 USE lmdz_ssum_scopy, ONLY: ssum 18 17 19 IMPLICIT NONE 18 20 ! … … 46 48 REAL :: zqh(ip1jmp1, llm), zqb(ip1jmp1, llm) 47 49 REAL :: temps0, temps1, temps2, temps3 48 REAL :: ztemps1, ztemps2 , ssum50 REAL :: ztemps1, ztemps2 49 51 save temps1, temps2, temps3 50 52 REAL :: zzpbar, zzw … … 726 728 ! -------------------------------------------------------------------- 727 729 USE lmdz_iniprint, ONLY: lunout, prt_level 730 USE lmdz_ssum_scopy, ONLY: ssum 731 728 732 IMPLICIT NONE 729 733 ! … … 749 753 REAL :: convpn, convps, convmpn, convmps, massen, masses 750 754 REAL :: zm, zq, zsigm, zsigp, zqm, zqp 751 REAL :: ssum752 755 REAL :: prec 753 756 save prec -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advy.f90
r5116 r5123 92 92 REAL :: z1(iim), zcos(iim), zsin(iim) 93 93 REAL :: smpn,smps,s0pn,s0ps 94 REAL :: SSUM95 EXTERNAL SSUM96 94 ! 97 95 REAL :: sqi,sqf -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advyp.f90
r5105 r5123 109 109 REAL :: cy1(llm,ntra), cyLAT(llm,ntra) 110 110 REAL :: z1(iim), zcos(iim), zsin(iim) 111 REAL :: SSUM112 EXTERNAL SSUM113 111 ! 114 112 REAL :: sqi,sqf -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.f90
r5105 r5123 18 18 ! nbniv est le nombre de niveaux vert. de xflu et de yflu . 19 19 ! 20 USE lmdz_ssum_scopy, ONLY: ssum 21 22 20 23 IMPLICIT NONE 21 24 ! … … 26 29 DIMENSION xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) , & 27 30 convfl( ip1jmp1,nbniv ) 28 ! 29 REAL :: SSUM 30 ! 31 ! 31 32 32 33 INCLUDE "comgeom.h" 33 34 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90
r5118 r5123 36 36 REAL :: x 37 37 38 REAL :: SSUM39 EXTERNAL SSUM40 38 REAL :: newsig 41 39 REAL :: dz0,dz1,nhaut,sig1,esig,csig,zz -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.f90
r5105 r5123 11 11 ! x et y etant des composantes covariantes ... 12 12 ! ********************************************************************* 13 USE lmdz_ssum_scopy, ONLY: ssum 14 13 15 IMPLICIT NONE 14 16 ! … … 37 39 REAL :: sumypn,sumyps 38 40 ! ................................................................... 39 ! 40 REAL :: SSUM 41 ! 42 ! 41 43 42 DO l = 1,klevel 44 43 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.f90
r5105 r5123 12 12 ! x et y etant des composantes covariantes ... 13 13 ! ********************************************************************* 14 USE lmdz_ssum_scopy, ONLY: ssum 15 14 16 IMPLICIT NONE 15 17 ! … … 40 42 INTEGER :: l,ij 41 43 ! ................................................................... 42 ! 43 REAL :: SSUM 44 ! 44 45 45 ! 46 46 DO l = 1,klevel -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.f90
r5106 r5123 12 12 ! ********************************************************************* 13 13 USE lmdz_filtreg, ONLY: filtreg 14 USE lmdz_ssum_scopy, ONLY: ssum 15 14 16 IMPLICIT NONE 15 17 ! … … 38 40 REAL :: sumypn,sumyps 39 41 ! ................................................................... 40 ! 41 REAL :: SSUM 42 ! 42 43 43 ! 44 44 DO l = 1,klevel -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.f90
r5105 r5123 1 2 1 ! $Header$ 3 2 4 SUBROUTINE divergst(klevel,x,y,div) 3 SUBROUTINE divergst(klevel, x, y, div) 4 USE lmdz_ssum_scopy, ONLY: ssum 5 5 6 IMPLICIT NONE 6 7 ! … … 17 18 ! ------------------------------------------------------------------- 18 19 ! 20 19 21 INCLUDE "dimensions.h" 20 22 INCLUDE "paramet.h" … … 22 24 23 25 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 30 30 ! 31 31 ! 32 DO l = 1, klevel33 !34 DO ij = iip2, ip1jm - 135 div( ij + 1, l ) = x(ij+1,l) - x(ij,l)+ y(ij-iim,l)-y(ij+1,l)36 END DO37 !38 ! .... correction pour div( 1,j,l) ......39 ! .... div(1,j,l)= div(iip1,j,l) ....40 !41 !DIR$ IVDEP42 DO ij = iip2,ip1jm,iip143 div( ij,l ) = div( ij + iim,l)44 END DO45 !46 ! .... calcul aux poles .....47 !48 !49 DO i = 1,iim50 aiy1(i)= y(i,l)51 aiy2(i)= y(i+ip1jmi1,l)52 END DO53 sumypn = SSUM ( iim,aiy1,1)54 sumyps = SSUM ( iim,aiy2,1)55 DO i = 1,iip156 div( i , l ) = - sumypn/iim57 div( i + ip1jm, l ) = sumyps/iim58 END DO59 !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 ! 60 60 END DO 61 61 RETURN -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.f90
r5117 r5123 22 22 USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, & 23 23 alphax,alphay,taux,tauy,transx,transy,pxo,pyo 24 USE lmdz_ssum_scopy, ONLY: ssum 24 25 IMPLICIT NONE 25 26 ! … … 50 51 SAVE rlatu1,yprimu1,rlatu2,yprimu2,yprimv,yprimu 51 52 SAVE rlonm025,xprimm025,rlonp025,xprimp025 52 53 REAL :: SSUM54 53 ! 55 54 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limx.f90
r5117 r5123 42 42 Logical :: extremum,first 43 43 save first 44 45 REAL :: SSUM46 INTEGER :: ismax,ismin47 EXTERNAL SSUM, ismin,ismax48 49 44 data first/.TRUE./ 50 45 -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limy.f90
r5117 r5123 16 16 USE comconst_mod, ONLY: pi 17 17 USE lmdz_libmath, ONLY: ismax, ismin 18 USE lmdz_ssum_scopy, ONLY: ssum 19 18 20 IMPLICIT NONE 19 21 ! … … 47 49 REAL :: coslon(iip1),coslondlon(iip1) 48 50 save sinlon,coslon,sinlondlon,coslondlon 49 !50 !51 REAL :: SSUM52 EXTERNAL SSUM53 51 54 52 data first/.TRUE./ -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limz.f90
r5117 r5123 42 42 Logical :: extremum,first 43 43 save first 44 45 REAL :: SSUM46 INTEGER :: ismax,ismin47 EXTERNAL SSUM, ismin,ismax48 44 49 45 data first/.TRUE./ -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.f90
r5105 r5123 31 31 REAL :: massemoyn, massemoys 32 32 33 REAL :: SSUM34 33 ! 35 34 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.f90
r5117 r5123 5 5 6 6 USE comconst_mod, ONLY: pi, dtvr 7 USE lmdz_ssum_scopy, ONLY: ssum 7 8 8 9 IMPLICIT NONE … … 55 56 REAL :: qmin,zq,pente_max 56 57 ! 57 REAL :: SSUM 58 INTEGER :: ismax,ismin,lati,latf 59 EXTERNAL SSUM, ismin,ismax 58 INTEGER :: lati,latf 60 59 LOGICAL :: first 61 60 save first -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.f90
r5117 r5123 5 5 6 6 USE comconst_mod, ONLY: pi 7 USE lmdz_ssum_scopy, ONLY: ssum 7 8 8 9 IMPLICIT NONE … … 58 59 REAL :: dyn1,dyn2,dys1,dys2,qpn,qps,dqzpn,dqzps 59 60 REAL :: masn,mass 60 ! 61 REAL :: SSUM 62 INTEGER :: ismax,ismin 63 EXTERNAL SSUM, ismin,ismax 61 64 62 LOGICAL :: first 65 63 save first -
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/sortvarc.f90
r5119 r5123 13 13 USE lmdz_filtreg, ONLY: filtreg 14 14 USE lmdz_iniprint, ONLY: lunout, prt_level 15 USE lmdz_ssum_scopy, ONLY: scopy 15 USE lmdz_ssum_scopy, ONLY: scopy, ssum 16 16 17 17 IMPLICIT NONE … … 64 64 INTEGER :: l, ij, imjmp1 65 65 66 REAL :: SSUM67 66 LOGICAL, SAVE :: firstcal = .TRUE. 68 67 CHARACTER(LEN = *), PARAMETER :: modname = "sortvarc" -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/addfi_loc.f90
r5117 r5123 1 2 1 ! $Id$ 3 2 4 3 SUBROUTINE 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) 7 6 USE parallel_lmdz 8 7 USE infotrac, ONLY: nqtot 9 8 USE control_mod, ONLY: planet_type 9 USE lmdz_ssum_scopy, ONLY: ssum 10 10 11 IMPLICIT NONE 11 12 ! … … 53 54 ! ----------- 54 55 ! 55 REAL, INTENT(IN) :: pdt ! time step for the integration (s)56 ! 57 REAL, INTENT(INOUT) :: pvcov(ijb_v:ije_v,llm) ! covariant meridional wind58 REAL, INTENT(INOUT) :: pucov(ijb_u:ije_u,llm) ! covariant zonal wind59 REAL, INTENT(INOUT) :: pteta(ijb_u:ije_u,llm) ! potential temperature60 REAL, INTENT(INOUT) :: pq(ijb_u:ije_u,llm,nqtot) ! tracers61 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) 62 63 ! 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 used64 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 70 71 ! 71 72 ! … … 73 74 ! ----------------- 74 75 ! 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 84 82 ! 85 83 !----------------------------------------------------------------------- 86 84 87 ijb =ij_begin88 ije =ij_end89 90 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)91 DO k = 1, llm92 DO j = ijb,ije93 pteta(j,k)= pteta(j,k) + pdhfi(j,k) * pdt94 95 ENDDO 96 !$OMP END DO NOWAIT85 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 97 95 98 96 IF (pole_nord) THEN 99 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)100 DO k 101 DO ij= 1, iim102 xpn(ij) = aire( ij ) * pteta( ij ,k)103 ENDDO104 tpn = SSUM(iim,xpn,1)/ apoln105 106 DO ij= 1, iip1107 pteta( ij ,k)= tpn108 ENDDO109 ENDDO110 !$OMP END DO NOWAIT97 !$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 111 109 ENDIF 112 110 113 111 IF (pole_sud) THEN 114 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)115 DO k 116 DO ij= 1, iim117 xps(ij) = aire(ij+ip1jm) * pteta(ij+ip1jm,k)118 ENDDO119 tps = SSUM(iim,xps,1)/ apols120 121 DO ij= 1, iip1122 pteta(ij+ip1jm,k)= tps123 ENDDO124 ENDDO125 !$OMP END DO NOWAIT126 ENDIF 127 ! 128 129 ijb =ij_begin130 ije =ij_end131 IF (pole_nord) ijb =ij_begin+iip1132 IF (pole_sud) ije =ij_end-iip1133 134 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)135 DO k = 1, llm136 DO j = ijb,ije137 pucov(j,k)= pucov(j,k) + pdufi(j,k) * pdt138 139 ENDDO 140 !$OMP END DO NOWAIT141 142 IF (pole_nord) ijb =ij_begin143 144 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)145 DO k = 1, llm146 DO j = ijb,ije147 pvcov(j,k)= pvcov(j,k) + pdvfi(j,k) * pdt148 149 ENDDO 150 !$OMP END DO NOWAIT151 152 ! 153 IF (pole_sud) ije =ij_end154 !$OMP MASTER155 DO j = ijb, ije156 157 ENDDO 158 !$OMP END MASTER112 !$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 159 157 160 158 IF (planet_type=="earth") THEN 161 ! earth case, special treatment for first 2 tracers (water)162 DO iq = 1, 2163 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)164 DO k = 1,llm165 DO j = ijb, ije166 pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt167 pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestw)168 ENDDO 169 ENDDO170 !$OMP END DO NOWAIT171 ENDDO172 173 DO iq = 3, nqtot174 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)175 DO k = 1,llm176 DO j = ijb, ije177 pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt178 pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestt)179 ENDDO 180 ENDDO181 !$OMP END DO NOWAIT182 ENDDO159 ! 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 183 181 else 184 ! general case, treat all tracers equally)185 DO iq = 1, nqtot186 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)187 DO k = 1,llm188 DO j = ijb, ije189 pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt190 pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestt)191 ENDDO 192 ENDDO193 !$OMP END DO NOWAIT194 ENDDO182 ! 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 195 193 ENDIF ! of if (planet_type=="earth") 196 194 197 !$OMP MASTER195 !$OMP MASTER 198 196 IF (pole_nord) THEN 199 DO ij 200 xpn(ij) = aire( ij ) * pps( ij)201 ENDDO 202 203 tpn = SSUM(iim,xpn,1)/apoln204 205 DO ij 206 pps ( ij )= tpn197 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 207 205 ENDDO 208 206 … … 210 208 211 209 IF (pole_sud) THEN 212 DO ij 213 xps(ij) = aire(ij +ip1jm) * pps(ij+ip1jm)214 ENDDO 215 216 tps = SSUM(iim,xps,1)/apols217 218 DO ij 219 pps ( ij+ip1jm )= tps220 ENDDO 221 222 ENDIF 223 !$OMP END MASTER210 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 224 222 225 223 IF (pole_nord) THEN 226 224 DO iq = 1, nqtot 227 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)228 DO k 229 DO ij 230 xpn(ij) = aire( ij ) * pq( ij ,k,iq)231 ENDDO 232 tpn = SSUM(iim,xpn,1)/apoln233 234 DO ij 235 pq ( ij ,k,iq)= tpn236 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 239 237 ENDDO 240 238 ENDIF … … 242 240 IF (pole_sud) THEN 243 241 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 260 256 261 257 END SUBROUTINE addfi_loc -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/advect_new_loc.f90
r5117 r5123 51 51 52 52 INTEGER :: ij, l, ijb, ije 53 EXTERNAL SSUM54 REAL :: SSUM55 56 57 53 58 54 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/call_dissip_mod.f90
r5117 r5123 85 85 USE write_field_loc 86 86 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 87 USE lmdz_ssum_scopy, ONLY: ssum 88 87 89 88 90 IMPLICIT NONE … … 97 99 REAL :: tpn, tps 98 100 99 REAL SSUM100 101 LOGICAL, PARAMETER :: dissip_conservative = .TRUE. 101 102 TYPE(Request), SAVE :: Request_dissip -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/convflu_loc.f90
r5105 r5123 16 16 ! 17 17 USE parallel_lmdz 18 USE lmdz_ssum_scopy, ONLY: ssum 19 18 20 IMPLICIT NONE 19 21 ! … … 26 28 ! 27 29 INTEGER :: ijb,ije 28 EXTERNAL SSUM29 REAL :: SSUM30 30 ! 31 31 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dissip_loc.f90
r5117 r5123 56 56 57 57 INTEGER :: l, ij 58 59 REAL :: SSUM60 58 INTEGER :: ijb, ije 61 59 -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.f90
r5117 r5123 10 10 ! ********************************************************************* 11 11 USE parallel_lmdz 12 USE lmdz_ssum_scopy, ONLY: ssum 13 12 14 IMPLICIT NONE 13 15 ! … … 39 41 INTEGER :: l,ij 40 42 ! ................................................................... 41 !42 EXTERNAL SSUM43 REAL :: SSUM44 43 INTEGER :: ijb,ije,jjb,jje 45 44 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.f90
r5117 r5123 9 9 ! ********************************************************************* 10 10 USE parallel_lmdz 11 USE lmdz_ssum_scopy, ONLY: ssum 12 11 13 IMPLICIT NONE 12 14 ! … … 36 38 INTEGER :: ijb,ije 37 39 ! ................................................................... 38 !39 EXTERNAL SSUM40 REAL :: SSUM41 40 ! 42 41 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/divergf_loc.f90
r5117 r5123 10 10 USE parallel_lmdz 11 11 USE lmdz_filtreg_p 12 USE lmdz_ssum_scopy, ONLY: ssum 13 12 14 IMPLICIT NONE 13 15 ! … … 37 39 REAL :: sumypn,sumyps 38 40 ! ................................................................... 39 !40 EXTERNAL SSUM41 REAL :: SSUM42 41 INTEGER :: ijb,ije,jjb,jje 43 42 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.f90
r5118 r5123 17 17 USE lmdz_strings, ONLY: int2str 18 18 USE lmdz_iniprint, ONLY: lunout, prt_level 19 USE lmdz_ssum_scopy, ONLY: ssum 19 20 20 21 IMPLICIT NONE … … 78 79 INTEGER :: l,ij,iq,i,j 79 80 80 REAL :: SSUM81 EXTERNAL SSUM82 81 INTEGER :: ijb,ije,jjb,jje 83 82 LOGICAL :: checksum -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/leapfrog_loc.F90
r5118 r5123 134 134 REAL :: time 135 135 136 REAL :: SSUM137 136 ! REAL,SAVE,ALLOCATABLE :: finvmaold(:,:) 138 137 -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/massdair_loc.f90
r5117 r5123 30 30 REAL :: massemoyn, massemoys 31 31 32 REAL :: SSUM33 EXTERNAL SSUM34 32 ! 35 33 ! -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/qminimum_loc.f90
r5118 r5123 38 38 REAL :: zx_defau_diag(ijb_u:ije_u, llm, 2) 39 39 REAL :: q_follow(ijb_u:ije_u, llm, 2) 40 !41 REAL :: SSUM42 EXTERNAL SSUM43 40 ! 44 41 INTEGER :: imprim -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90
r5118 r5123 46 46 47 47 Logical :: extremum 48 49 REAL :: SSUM50 48 51 49 REAL :: z1, z2, z3 … … 384 382 min_qParent, min_qMass, min_ratio ! MVals et CRisi 385 383 USE comconst_mod, ONLY: pi 384 USE lmdz_ssum_scopy, ONLY: ssum 385 386 386 IMPLICIT NONE 387 387 ! … … 428 428 REAL :: Ratio(ijb_u:ije_u, llm, nqtot) ! CRisi 429 429 INTEGER :: ifils, iq2 ! CRisi 430 !431 !432 REAL :: SSUM433 EXTERNAL SSUM434 430 435 431 DATA first/.TRUE./ … … 850 846 !$OMP THREADPRIVATE(temps0,temps1,temps2,temps3,temps4,temps5) 851 847 852 REAL :: SSUM853 854 848 DATA temps0, temps1, temps2, temps3, temps4, temps5/0., 0., 0., 0., 0., 0./ 855 849 INTEGER :: ijb, ije, ijb_x, ije_x -
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_loc.f90
r5118 r5123 40 40 REAL :: Ratio(ijb_u:ije_u, llm, nqtot) ! CRisi 41 41 INTEGER :: ifils, iq2 ! CRisi 42 43 REAL :: SSUM44 42 45 43 INTEGER :: ijb, ije, ijb_x, ije_x … … 392 390 USE comconst_mod, ONLY: pi 393 391 USE lmdz_iniprint, ONLY: lunout, prt_level 392 USE lmdz_ssum_scopy, ONLY: ssum 393 394 394 IMPLICIT NONE 395 395 ! … … 437 437 INTEGER :: ijb, ije 438 438 INTEGER :: ijbm, ijem 439 440 REAL :: ssum441 439 442 440 ijb = ij_begin - 2 * iip1 -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.f90
r5119 r5123 35 35 USE comvert_mod, ONLY: preff, presnivs 36 36 USE lmdz_iniprint, ONLY: lunout, prt_level 37 USE lmdz_ssum_scopy, ONLY: scopy 37 USE lmdz_ssum_scopy, ONLY: scopy, ssum 38 38 39 39 IMPLICIT NONE … … 164 164 REAL :: flxwfi(ngridmx, llm) ! Flux de masse verticale sur la grille physiq 165 165 ! 166 167 REAL :: SSUM168 169 166 LOGICAL, SAVE :: firstcal = .TRUE., debut = .TRUE. 170 167 ! REAL rdayvrai -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90
r5118 r5123 51 51 USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, kappa, pi 52 52 USE lmdz_iniprint, ONLY: lunout, prt_level 53 USE lmdz_ssum_scopy, ONLY: ssum 54 53 55 54 56 !======================================================================= … … 216 218 REAL :: zsinbis(iim),zcosbis(iim),z1bis(iim) 217 219 REAL :: unskap, pksurcp 218 219 REAL :: SSUM220 220 221 221 LOGICAL,SAVE :: firstcal=.TRUE., debut=.TRUE. -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90
r5118 r5123 460 460 var3d(iml, :, :) = var3d(1, :, :) 461 461 462 DO il = 1, lml463 CALL minmax(iml * jml, var3d(1, 1, il), chmin, chmax)464 WRITE(lunout, *)' ' // TRIM(var) // ' min max l ', il, chmin, chmax465 END DO466 467 462 END SUBROUTINE start_inter_3d 468 463 -
LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_filtreg.F90
r5120 r5123 330 330 ! 331 331 USE lmdz_coefils, ONLY: sddv, sddu, unsddu, unsddv, eignfnv, eignfnu 332 USE lmdz_ssum_scopy, ONLY: ssum 332 333 IMPLICIT NONE 333 334 ! … … 341 342 REAL :: pi 342 343 INTEGER :: i, j, k, imm1, nrot 343 EXTERNAL SSUM344 REAL :: SSUM345 344 ! 346 345 … … 557 556 558 557 SUBROUTINE acc(vec, d, im) 558 USE lmdz_ssum_scopy, ONLY: ssum 559 559 IMPLICIT NONE 560 560 INTEGER :: im … … 562 562 INTEGER :: i, j 563 563 REAL :: sum 564 REAL, external :: ssum565 564 do j = 1, im 566 565 do i = 1, im -
LMDZ6/branches/Amaury_dev/libf/grid/lmdz_fxy_prim.f90
r5122 r5123 1 ! This module replaces grid/fxy_ reg.h1 ! This module replaces grid/fxy_prim.h 2 2 3 3 MODULE lmdz_fxy_prim -
LMDZ6/branches/Amaury_dev/libf/grid/lmdz_fxy_sin.f90
r5122 r5123 1 ! This module replaces grid/fxy_ reg.h1 ! This module replaces grid/fxy_sin.h 2 2 3 3 MODULE lmdz_fxy_sin -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_libmath_pch.f90
r5117 r5123 915 915 REAL :: G, HALF, ONE, STEMP(3), THREE, TWO, XTEMP(4), ZERO 916 916 SAVE ZERO, HALF, ONE, TWO, THREE 917 REAL :: PCHDF918 917 ! 919 918 DATA ZERO /0./, HALF /0.5/, ONE /1./, TWO /2./, THREE /3./ -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_wxios.F90
r5119 r5123 443 443 USE lmdz_grid_phy, ONLY: nvertex, klon_glo 444 444 USE lmdz_physical_constants, ONLY: PI 445 USE lmdz_phys_para 445 446 USE lmdz_ioipsl_getin_p, ONLY: getin_p 446 447 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_xer.f90
r5117 r5123 201 201 CHARACTER(len = 20) :: LFIRST 202 202 INTEGER :: NERR, LEVEL, LKNTRL 203 INTEGER :: J4SAVE,MAXMES, KDUMMY, I, KOUNT, LERR, LLEVEL203 INTEGER :: MAXMES, KDUMMY, I, KOUNT, LERR, LLEVEL 204 204 INTEGER :: MKNTRL, LTEMP 205 205 !***FIRST EXECUTABLE STATEMENT XERMSG … … 549 549 !***END PROLOGUE XGETUA 550 550 DIMENSION IUNITA(5) 551 INTEGER :: IUNITA, N, J4SAVE,INDEX, I551 INTEGER :: IUNITA, N, INDEX, I 552 552 !***FIRST EXECUTABLE STATEMENT XGETUA 553 553 N = J4SAVE(5, 0, .FALSE.) … … 630 630 INTEGER :: NERR, LEVEL, KONTRL 631 631 INTEGER :: NERTAB, LEVTAB, KOUNT, KOUNTX, NMSG 632 INTEGER :: KFLAG, ICOUNT, NUNIT, KUNIT, IUNIT, I 1MACH, I632 INTEGER :: KFLAG, ICOUNT, NUNIT, KUNIT, IUNIT, I 633 633 !***FIRST EXECUTABLE STATEMENT XERSVE 634 634 ! … … 803 803 CHARACTER(len = 2) :: NEWLIN 804 804 PARAMETER (NEWLIN = '$$') 805 INTEGER :: N, I 1MACH, I, LPREF, LWRAP, LENMSG, NEXTC805 INTEGER :: N, I, LPREF, LWRAP, LENMSG, NEXTC 806 806 INTEGER :: LPIECE, IDELTA 807 807 !***FIRST EXECUTABLE STATEMENT XERPRN … … 1056 1056 1057 1057 !DECK J4SAVE 1058 FUNCTION J4SAVE(IWHICH, IVALUE, ISET)1058 INTEGER FUNCTION J4SAVE(IWHICH, IVALUE, ISET) 1059 1059 IMPLICIT NONE 1060 1060 !***BEGIN PROLOGUE J4SAVE … … 1116 1116 DATA IPARAM(5)/1/ 1117 1117 DATA IPARAM(6), IPARAM(7), IPARAM(8), IPARAM(9)/0, 0, 0, 0/ 1118 INTEGER :: J4SAVE,IWHICH, IVALUE1118 INTEGER :: IWHICH, IVALUE 1119 1119 !***FIRST EXECUTABLE STATEMENT J4SAVE 1120 1120 J4SAVE = IPARAM(IWHICH) -
LMDZ6/branches/Amaury_dev/libf/phy_common/lmdz_abort_physic.F90
r5117 r5123 7 7 8 8 USE IOIPSL 9 USE lmdz_phys_ para9 USE lmdz_phys_mpi_data 10 10 USE lmdz_print_control, ONLY: lunout 11 11 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phys_output_write_spl_mod.F90
r5118 r5123 1486 1486 DO i = 1, klon 1487 1487 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 1490 1489 freq_moyNMC(iff-6) =(mth_len*un_jour)/freq_calNMC(iff-6) 1491 1490 ELSE -
LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90
r5117 r5123 3 3 4 4 MODULE iophy 5 USE lmdz_abort_physic, ONLY: abort_physic 5 6 6 7 ! abd REAL,PRIVATE,ALLOCATABLE,DIMENSION(:),save :: io_lat … … 809 810 USE lmdz_print_control, ONLY: prt_level,lunout 810 811 USE lmdz_grid_phy, ONLY: nbp_lon, nbp_lat 811 USE lmdz_abort_physic, ONLY: abort_physic812 812 813 813 IMPLICIT NONE … … 955 955 USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios 956 956 USE lmdz_print_control, ONLY: lunout, prt_level 957 USE lmdz_abort_physic, ONLY: abort_physic958 957 959 958 IMPLICIT NONE … … 1168 1167 USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios 1169 1168 USE lmdz_print_control, ONLY: prt_level,lunout 1170 USE lmdz_abort_physic, ONLY: abort_physic1171 1169 1172 1170 IMPLICIT NONE … … 1365 1363 USE lmdz_xios, ONLY: xios_send_field 1366 1364 USE lmdz_print_control, ONLY: prt_level,lunout 1367 USE lmdz_abort_physic, ONLY: abort_physic1368 1365 1369 1366 IMPLICIT NONE … … 1449 1446 USE lmdz_grid_phy, ONLY: nbp_lon, nbp_lat, grid_type, unstructured 1450 1447 USE lmdz_print_control, ONLY: prt_level,lunout 1451 USE lmdz_abort_physic, ONLY: abort_physic1452 1448 1453 1449 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/limit_read_mod.F90
r5116 r5123 170 170 USE phys_cal_mod, ONLY: calend, year_len 171 171 USE lmdz_print_control, ONLY: lunout, prt_level 172 USE lmdz_xios, ONLY: xios_recv_field, using_xios w172 USE lmdz_xios, ONLY: xios_recv_field, using_xios 173 173 USE lmdz_abort_physic, ONLY: abort_physic 174 174 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_main.F90
r5117 r5123 7 7 CONTAINS 8 8 9 SUBROUTINE thermcell_main 10 (itap, ngrid, nlay, ptimestep & 9 SUBROUTINE thermcell_main(itap, ngrid, nlay, ptimestep & 11 10 , pplay,pplev, pphi, debut & 12 11 , puwind, pvwind,ptemp, p_o, ptemp_env, po_env & -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_old.F90
r5119 r5123 1 1 MODULE lmdz_thermcell_old 2 USE lmdz_abort_physic, ONLY: abort_physic 3 2 4 CONTAINS 3 5 … … 10 12 USE lmdz_thermcell_dv2, ONLY: thermcell_dv2 11 13 USE lmdz_thermcell_dq, ONLY: thermcell_dq 12 USE lmdz_abort_physic, ONLY: abort_physic13 14 IMPLICIT NONE 14 15 -
LMDZ6/branches/Amaury_dev/libf/phylmd/pbl_surface_mod.F90
r5117 r5123 256 256 !**************************************************************************************** 257 257 258 SUBROUTINE pbl_surface 259 ( & 258 SUBROUTINE pbl_surface( & 260 259 dtime, date0, itap, jour, & 261 260 debut, lafin, & … … 4186 4185 !**************************************************************************************** 4187 4186 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 & 4190 4188 #ifdef ISO 4191 4189 ,xtsnow_rst,Rland_ice_rst & … … 4259 4257 4260 4258 !albedo SB >>> 4261 SUBROUTINE pbl_surface_newfrac 4262 (itime, pctsrf_new, pctsrf_old, & 4259 SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, & 4263 4260 evap, z0m, z0h, agesno, & 4264 4261 tsurf, alb_dir, alb_dif, ustar, u10m, v10m, tke & -
LMDZ6/branches/Amaury_dev/libf/phylmd/slab_heat_transp_mod.F90
r5119 r5123 986 986 ! computes the divergence of a vector field of components 987 987 ! x,y. x and y being covariant components 988 USE lmdz_ssum_scopy, ONLY: ssum 989 988 990 IMPLICIT NONE 989 991 … … 998 1000 REAL :: aiy1(nbp_lon + 1), aiy2(nbp_lon + 1) 999 1001 REAL :: sumypn, sumyps 1000 REAL, EXTERNAL :: SSUM1001 1002 1002 1003 iim = nbp_lon -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90
r5117 r5123 7 7 !**************************************************************************************** 8 8 9 SUBROUTINE surf_landice 10 (itime, dtime, knon, knindex, & 9 SUBROUTINE surf_landice(itime, dtime, knon, knindex, & 11 10 rlon, rlat, debut, lafin, & 12 11 rmu0, lwdownm, albedo, pphi1, & -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90
r5117 r5123 9 9 !****************************************************************************** 10 10 11 SUBROUTINE surf_ocean 12 (rlon, rlat, swnet, lwnet, alb1, & 11 SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, & 13 12 windsp, rmu0, fder, tsurf_in, & 14 13 itime, dtime, jour, knon, knindex, &
Note: See TracChangeset
for help on using the changeset viewer.