Changeset 1176 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/inidissip.F
- Timestamp:
- Jun 11, 2009, 10:54:10 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/inidissip.F
r774 r1176 1 1 ! 2 ! $ Header$2 ! $Id$ 3 3 ! 4 4 SUBROUTINE inidissip ( lstardis,nitergdiv,nitergrot,niterh , … … 18 18 #include "comvert.h" 19 19 #include "control.h" 20 #include "logic.h" 20 21 21 22 LOGICAL lstardis … … 29 30 INTEGER l,ij,idum,ii 30 31 REAL tetamin 32 REAL pseudoz 31 33 32 34 REAL ran1 … … 174 176 c -------------------------------------------------- 175 177 176 DO l=1,llm 177 zvert(l)=1. 178 ENDDO 179 180 fact=2. 181 c 182 DO l = 1, llm 183 zz = 1. - preff/presnivs(l) 184 zvert(l)= fact -( fact-1.)/( 1.+zz*zz ) 185 ENDDO 178 if (ok_strato .and. llm==39) then 179 do l=1,llm 180 pseudoz=8.*log(preff/presnivs(l)) 181 zvert(l)=1+ 182 s (tanh((pseudoz-dissip_zref)/dissip_deltaz)+1.)/2. 183 s *(dissip_factz-1.) 184 enddo 185 else 186 DO l=1,llm 187 zvert(l)=1. 188 ENDDO 189 fact=2. 190 DO l = 1, llm 191 zz = 1. - preff/presnivs(l) 192 zvert(l)= fact -( fact-1.)/( 1.+zz*zz ) 193 ENDDO 194 endif 186 195 187 196
Note: See TracChangeset
for help on using the changeset viewer.