- Timestamp:
- Jul 24, 2024, 4:23:34 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inidissip.F90
r5116 r5117 16 16 USE lmdz_filtreg, ONLY: filtreg 17 17 USE lmdz_libmath, ONLY: minmax 18 USE lmdz_ran1, ONLY: ran1 18 19 19 20 IMPLICIT NONE … … 23 24 include "iniprint.h" 24 25 25 LOGICAL, INTENT( in) :: lstardis26 INTEGER, INTENT( in) :: nitergdiv, nitergrot, niterh27 REAL, INTENT( in) :: tetagdiv, tetagrot, tetatemp28 29 integer, INTENT(in) :: vert_prof_dissip26 LOGICAL, INTENT(IN) :: lstardis 27 INTEGER, INTENT(IN) :: nitergdiv, nitergrot, niterh 28 REAL, INTENT(IN) :: tetagdiv, tetagrot, tetatemp 29 30 INTEGER, INTENT(IN) :: vert_prof_dissip 30 31 ! Vertical profile of horizontal dissipation 31 32 ! Allowed values: … … 36 37 REAL fact, zvert(llm), zz 37 38 REAL zh(ip1jmp1), zu(ip1jmp1), gx(ip1jmp1), divgra(ip1jmp1) 38 realzv(ip1jm), gy(ip1jm), deltap(ip1jmp1, llm)39 REAL zv(ip1jm), gy(ip1jm), deltap(ip1jmp1, llm) 39 40 REAL ullm, vllm, umin, vmin, zhmin, zhmax 40 41 REAL zllm … … 43 44 REAL tetamin 44 45 REAL pseudoz 45 character (len = 80) :: abort_message 46 47 REAL ran1 48 46 CHARACTER (LEN = 80) :: abort_message 49 47 50 48 !----------------------------------------------------------------------- … … 174 172 ! -------------------------------------------------- 175 173 176 if(vert_prof_dissip == 1) THEN174 IF (vert_prof_dissip == 1) THEN 177 175 do l = 1, llm 178 176 pseudoz = 8. * log(preff / presnivs(l)) … … 190 188 zvert(l) = fact - (fact - 1.) / (1. + zz * zz) 191 189 ENDDO 192 endif190 ENDIF 193 191 194 192 WRITE(lunout, *)'inidissip: Constantes de temps de la diffusion horizontale'
Note: See TracChangeset
for help on using the changeset viewer.