Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d/caladvtrac.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/caladvtrac.F90
r5134 r5159 1 1 ! $Id$ 2 2 3 ! 4 ! 3 4 5 5 SUBROUTINE caladvtrac(q, pbaru, pbarv, & 6 6 p, masse, dq, teta, & 7 7 flxw, pk) 8 ! 8 9 9 USE infotrac, ONLY: nqtot 10 10 USE control_mod, ONLY: iapp_tracvl, planet_type … … 13 13 USE lmdz_ssum_scopy, ONLY: scopy 14 14 15 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 16 USE lmdz_paramet 15 17 IMPLICIT NONE 16 ! 18 17 19 ! Auteurs: F.Hourdin , P.Le Van, F.Forget, F.Codron 18 ! 20 19 21 ! F.Codron (10/99) : ajout humidite specifique pour eau vapeur 20 22 !======================================================================= 21 ! 23 22 24 ! Shema de Van Leer 23 ! 25 24 26 !======================================================================= 25 27 26 INCLUDE "dimensions.h" 27 INCLUDE "paramet.h" 28 29 28 30 29 31 ! Arguments: … … 36 38 37 39 ! .................................................................. 38 ! 40 39 41 ! .. dq n'est utilise et dimensionne que pour l'eau vapeur et liqu. 40 ! 42 41 43 ! .................................................................. 42 ! 44 43 45 ! Local: 44 46 ! ------ … … 49 51 50 52 !c 51 ! 53 52 54 ! Earth-specific stuff for the first 2 tracers (water) 53 55 IF (planet_type=="earth") THEN … … 73 75 IF (planet_type=="earth") THEN 74 76 ! Earth-specific treatment for the first 2 tracers (water) 75 ! 77 76 78 !c CALL minmaxq(q(1,1,1),1.e33,-1.e33,'Eau vapeur ') 77 79 !c CALL minmaxq(q(1,1,2),1.e33,-1.e33,'Eau liquide ') 78 80 79 81 !c .... Calcul de deltap qu'on stocke dans finmasse ... 80 ! 82 81 83 DO l = 1, llm 82 84 DO ij = 1, ip1jmp1 … … 91 93 CALL SCOPY (ip1jmp1 * llm, masse, 1, finmasse, 1) 92 94 CALL filtreg (finmasse, jjp1, llm, -2, 2, .TRUE., 1) 93 ! 95 94 96 ! ***** Calcul de dq pour l'eau , pour le passer a la physique ****** 95 97 ! ******************************************************************** 96 ! 98 97 99 dtvrtrac = iapp_tracvl * dtvr 98 ! 100 99 101 DO iq = 1, nqtot 100 102 DO l = 1, llm … … 105 107 ENDDO 106 108 ENDDO 107 ! 109 108 110 endif ! of if (planet_type.EQ."earth") 109 111 ELSE
Note: See TracChangeset
for help on using the changeset viewer.