Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d/integrd.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/integrd.F90
r5136 r5159 15 15 USE lmdz_comgeom 16 16 17 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 18 USE lmdz_paramet 17 19 IMPLICIT NONE 18 20 19 21 20 22 !======================================================================= 21 ! 23 22 24 ! Auteur: P. Le Van 23 25 ! ------- 24 ! 26 25 27 ! objet: 26 28 ! ------ 27 ! 29 28 30 ! Incrementation des tendances dynamiques 29 ! 31 30 32 !======================================================================= 31 33 !----------------------------------------------------------------------- … … 33 35 ! ------------- 34 36 35 INCLUDE "dimensions.h" 36 INCLUDE "paramet.h" 37 38 37 39 38 40 ! Arguments: … … 94 96 ps (ij) = psm1(ij) + dt * dp(ij) 95 97 ENDDO 96 ! 98 97 99 DO ij = 1, ip1jmp1 98 100 IF(ps(ij)<0.) THEN … … 107 109 ENDIF 108 110 ENDDO 109 ! 111 110 112 DO ij = 1, iim 111 113 tppn(ij) = aire(ij) * ps(ij) … … 118 120 ps(ij + ip1jm) = tps 119 121 ENDDO 120 ! 122 121 123 ! ... Calcul de la nouvelle masse d'air au dernier temps integre t+1 ... 122 ! 124 123 125 CALL pression (ip1jmp1, ap, bp, ps, p) 124 126 CALL massdair (p, masse) … … 151 153 152 154 ! .... Calcul de la valeur moyenne, unique aux poles pour teta ...... 153 ! 154 ! 155 156 155 157 DO ij = 1, iim 156 158 tppn(ij) = aire(ij) * teta(ij, l) … … 175 177 176 178 177 ! 179 178 180 ! ....... integration de q ...... 179 ! 181 180 182 !$$$ IF( iadv(1).NE.3.AND.iadv(2).NE.3 ) THEN 181 183 !$$$c … … 212 214 CALL qminimum(q, nq, deltap) 213 215 214 ! 216 215 217 ! ..... Calcul de la valeur moyenne, unique aux poles pour q ..... 216 218 ! … … 238 240 239 241 ENDIF ! of if (planet_type.EQ."earth") 240 ! 241 ! 242 243 242 244 ! ..... FIN de l'integration de q ....... 243 245
Note: See TracChangeset
for help on using the changeset viewer.