Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/limy.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_common/limy.f90
r5158 r5159 1 ! 1 2 2 ! $Id$ 3 ! 3 4 4 SUBROUTINE limy(s0,sy,sm,pente_max) 5 ! 5 6 6 ! Auteurs: P.Le Van, F.Hourdin, F.Forget 7 ! 7 8 8 ! ******************************************************************** 9 9 ! Shema d'advection " pseudo amont " . … … 11 11 ! q,w sont des arguments d'entree pour le s-pg .... 12 12 ! dq sont des arguments de sortie pour le s-pg .... 13 ! 14 ! 13 14 15 15 ! -------------------------------------------------------------------- 16 16 USE comconst_mod, ONLY: pi … … 19 19 USE lmdz_comgeom 20 20 21 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 22 USE lmdz_paramet 21 23 IMPLICIT NONE 22 24 ! 23 INCLUDE "dimensions.h" 24 INCLUDE "paramet.h" 25 ! 26 ! 25 26 27 28 27 29 ! Arguments: 28 30 ! ---------- 29 31 REAL :: pente_max 30 32 REAL :: s0(ip1jmp1,llm),sy(ip1jmp1,llm),sm(ip1jmp1,llm) 31 ! 33 32 34 ! Local 33 35 ! --------- 34 ! 36 35 37 INTEGER :: i,ij,l 36 ! 38 37 39 REAL :: q(ip1jmp1,llm) 38 40 REAL :: airej2,airejjm,airescb(iim),airesch(iim) … … 70 72 71 73 DO l = 1, llm 72 ! 74 73 75 DO ij=1,ip1jmp1 74 76 q(ij,l) = s0(ij,l) / sm ( ij,l ) 75 77 dyq(ij) = sy(ij,l) / sm ( ij,l ) 76 78 ENDDO 77 ! 79 78 80 ! -------------------------------- 79 81 ! CALCUL EN LATITUDE
Note: See TracChangeset
for help on using the changeset viewer.