- Timestamp:
- Aug 23, 2016, 5:21:06 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/physiq_mod.F90
r2611 r2613 578 578 LOGICAL,SAVE :: ok_adjwk=.FALSE. 579 579 !$OMP THREADPRIVATE(ok_adjwk) 580 REAL,SAVE :: oliqmax=999. 581 !$OMP THREADPRIVATE(oliqmax) 580 582 REAL, SAVE :: alp_offset 581 583 !$OMP THREADPRIVATE(alp_offset) … … 1197 1199 CALL getin_p('random_notrig_max',random_notrig_max) 1198 1200 CALL getin_p('ok_adjwk',ok_adjwk) 1201 CALL getin_p('oliqmax',oliqmax) 1199 1202 CALL getin_p('ratqsp0',ratqsp0) 1200 1203 CALL getin_p('ratqsdp',ratqsdp) … … 3056 3059 CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, & 3057 3060 'lsc',abortphy,flag_inhib_tend) 3061 rain_num(:)=0. 3062 DO k = 1, klev ! re-evaporation de l'eau liquide nuageuse 3063 DO i = 1, klon 3064 IF (ql_seri(i,k)>oliqmax) THEN 3065 rain_num(i)=rain_num(i)+(ql_seri(i,k)-oliqmax)*zmasse(i,k)/pdtphys 3066 ql_seri(i,k)=oliqmax 3067 ENDIF 3068 ENDDO 3069 ENDDO 3070 3058 3071 !--------------------------------------------------------------------------- 3059 3072 DO k = 1, klev
Note: See TracChangeset
for help on using the changeset viewer.