- Timestamp:
- Oct 9, 2023, 2:23:04 PM (15 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/cdrag_mod.F90
r4662 r4721 174 174 INTEGER, SAVE :: iflag_corr_insta 175 175 !$OMP THREADPRIVATE(iflag_corr_insta) 176 LOGICAL, SAVE :: ok_cdrag_iter 177 !$OMP THREADPRIVATE(ok_cdrag_iter) 176 178 177 179 !===================================================================c … … 247 249 ! On choisit les fonctions de stabilite utilisees au premier appel 248 250 !************************************************************************** 249 251 IF (firstcall) THEN 250 252 iflag_corr_sta=2 251 253 iflag_corr_insta=2 254 ok_cdrag_iter = .FALSE. 252 255 253 256 CALL getin_p('iflag_corr_sta',iflag_corr_sta) 254 257 CALL getin_p('iflag_corr_insta',iflag_corr_insta) 258 CALL getin_p('ok_cdrag_iter',ok_cdrag_iter) 255 259 256 260 firstcall = .FALSE. … … 489 493 cdhh(i)=MIN(cdh(i),cdhmax) 490 494 ENDIF 491 rugos_itm(i,1) = rugos_itm(i,2) 492 rugos_ith(i,1) = rugos_ith(i,2) 493 rugos_itm(i,2) = 0.018*cdmm(i) * (speed(i))/RG & 495 IF (ok_cdrag_iter) THEN 496 rugos_itm(i,1) = rugos_itm(i,2) 497 rugos_ith(i,1) = rugos_ith(i,2) 498 rugos_itm(i,2) = 0.018*cdmm(i) * (speed(i))/RG & 494 499 + 0.11*14e-6 / SQRT(cdmm(i) * zdu2) 495 500 496 501 !---------- Version SEPARATION DES Z0 ---------------------- 497 IF (iflag_z0_oce==0) THEN 498 rugos_ith(i,2) = rugos_itm(i,2) 499 ELSE IF (iflag_z0_oce==1) THEN 500 rugos_ith(i,2) = 0.40*14e-6 / SQRT(cdmm(i) * zdu2) 502 IF (iflag_z0_oce==0) THEN 503 rugos_ith(i,2) = rugos_itm(i,2) 504 ELSE IF (iflag_z0_oce==1) THEN 505 rugos_ith(i,2) = 0.40*14e-6 / SQRT(cdmm(i) * zdu2) 506 ENDIF 501 507 ENDIF 502 508 ENDIF 503 rugos_itm(i,2) = MAX(1.5e-05,rugos_itm(i,2)) 504 rugos_ith(i,2) = MAX(1.5e-05,rugos_ith(i,2)) 509 IF (ok_cdrag_iter) THEN 510 rugos_itm(i,2) = MAX(1.5e-05,rugos_itm(i,2)) 511 rugos_ith(i,2) = MAX(1.5e-05,rugos_ith(i,2)) 512 ENDIF 505 513 ENDDO 506 514 IF (nsrf.EQ.is_oce) THEN
Note: See TracChangeset
for help on using the changeset viewer.