Changeset 2613 for LMDZ5


Ignore:
Timestamp:
Aug 23, 2016, 5:21:06 PM (8 years ago)
Author:
fhourdin
Message:

Intronduction d'une limite sur l'eau liquide nuageuse, oliqmax
Typiquement oliqmax>1e-3 (1g/kg)
Peut permettre d'éviter certains plantages.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2611 r2613  
    578578    LOGICAL,SAVE :: ok_adjwk=.FALSE.
    579579    !$OMP THREADPRIVATE(ok_adjwk)
     580    REAL,SAVE :: oliqmax=999.
     581    !$OMP THREADPRIVATE(oliqmax)
    580582    REAL, SAVE :: alp_offset
    581583    !$OMP THREADPRIVATE(alp_offset)
     
    11971199       CALL getin_p('random_notrig_max',random_notrig_max)
    11981200       CALL getin_p('ok_adjwk',ok_adjwk)
     1201       CALL getin_p('oliqmax',oliqmax)
    11991202       CALL getin_p('ratqsp0',ratqsp0)
    12001203       CALL getin_p('ratqsdp',ratqsdp)
     
    30563059    CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, &
    30573060         '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
    30583071    !---------------------------------------------------------------------------
    30593072    DO k = 1, klev
Note: See TracChangeset for help on using the changeset viewer.