Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d/addfi.F90
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d/addfi.F90
r5136 r5159 1 1 ! $Id$ 2 2 3 SUBROUTINE addfi(pdt, leapf, forward, & 4 pucov, pvcov, pteta, pq, pps, & 5 pdufi, pdvfi, pdhfi, pdqfi, pdpfi) 3 SUBROUTINE addfi(pdt, leapf, forward, pucov, pvcov, pteta, pq, pps, pdufi, pdvfi, pdhfi, pdqfi, pdpfi) 6 4 7 5 USE infotrac, ONLY: nqtot … … 9 7 USE lmdz_ssum_scopy, ONLY: ssum 10 8 USE lmdz_comgeom 9 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 10 USE lmdz_paramet 11 11 12 IMPLICIT NONE 12 ! 13 13 14 !======================================================================= 14 ! 15 15 16 ! Addition of the physical tendencies 16 ! 17 17 18 ! Interface : 18 19 ! ----------- 19 ! 20 20 21 ! Input : 21 22 ! ------- … … 31 32 ! pdhfi(ip1jmp1) | tendencies 32 33 ! pdtsfi(ip1jmp1) | 33 ! 34 34 35 ! Output : 35 36 ! -------- … … 38 39 ! ph 39 40 ! pts 40 ! 41 ! 41 42 42 43 !======================================================================= 43 ! 44 !----------------------------------------------------------------------- 45 ! 46 ! 0. Declarations : 47 ! ------------------ 48 ! 49 INCLUDE "dimensions.h" 50 INCLUDE "paramet.h" 51 ! 44 ! ! 52 45 ! Arguments : 53 46 ! ----------- 54 ! 47 55 48 REAL, INTENT(IN) :: pdt ! time step for the integration (s) 56 ! 49 57 50 REAL, INTENT(INOUT) :: pvcov(ip1jm, llm) ! covariant meridional wind 58 51 REAL, INTENT(INOUT) :: pucov(ip1jmp1, llm) ! covariant zonal wind … … 66 59 REAL, INTENT(IN) :: pdhfi(ip1jmp1, llm) 67 60 REAL, INTENT(IN) :: pdpfi(ip1jmp1) 68 ! 61 69 62 LOGICAL, INTENT(IN) :: leapf, forward ! not used 70 ! 71 ! 63 64 72 65 ! Local variables : 73 66 ! ----------------- 74 ! 67 75 68 REAL :: xpn(iim), xps(iim), tpn, tps 76 69 INTEGER :: j, k, iq, ij 77 70 REAL, PARAMETER :: qtestw = 1.0e-15 78 71 REAL, PARAMETER :: qtestt = 1.0e-40 79 ! 72 80 73 !----------------------------------------------------------------------- 81 74 … … 113 106 ENDDO 114 107 115 ! 108 116 109 DO j = 1, ip1jmp1 117 110 pps(j) = pps(j) + pdpfi(j) * pdt … … 177 170 ENDDO 178 171 179 180 172 END SUBROUTINE addfi
Note: See TracChangeset
for help on using the changeset viewer.