Changeset 2922 for LMDZ5/trunk


Ignore:
Timestamp:
Jun 29, 2017, 5:45:27 PM (7 years ago)
Author:
fhourdin
Message:

Nouvelle option iflag_wk_check_trgl pour tester la bonne sante des poches.
L'ancien flag logique flag_wk_check_trgl est toujours utilisable.
Il vaut y/n si iflag_wk_check_trgl=1/0.
L'option 1 avait tendance à faire disparaitre les poches dans des cas 1D.
L'option 2 est donc précaunisée pour les réglages à venir.
Fredho pour Jean-Yves

File:
1 edited

Legend:

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

    r2761 r2922  
    178178  LOGICAL, SAVE                                         :: flag_wk_check_trgl
    179179  !$OMP THREADPRIVATE(flag_wk_check_trgl)
     180  INTEGER, SAVE                                         :: iflag_wk_check_trgl
     181  !$OMP THREADPRIVATE(iflag_wk_check_trgl)
    180182
    181183  REAL                                                  :: delta_t_min
     
    347349  CALL getin_p('flag_wk_check_trgl ', flag_wk_check_trgl)
    348350  WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
     351  WRITE(*,*) 'flag_wk_check_trgl OBSOLETE. Utilisr iflag_wk_check_trgl plutot'
     352  iflag_wk_check_trgl=0 ; IF (flag_wk_check_trgl) iflag_wk_check_trgl=1
     353  CALL getin_p('iflag_wk_check_trgl ', iflag_wk_check_trgl)
     354  WRITE(*,*) 'iflag_wk_check_trgl=', iflag_wk_check_trgl
    349355
    350356  first=.false.
     
    17931799  ! Filter out bad wakes
    17941800
    1795   IF (flag_wk_check_trgl) THEN
     1801  IF (iflag_wk_check_trgl>=1) THEN
    17961802    ! Check triangular shape of dth profile
    17971803    DO i = 1, klon
     
    18051811          wape2(i) = -1.
    18061812          !! print *,'wake, rej 1'
    1807         ELSE IF (abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN
     1813        ELSE IF (iflag_wk_check_trgl==1.AND.abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN
    18081814          wape2(i) = -1.
    18091815          !! print *,'wake, rej 2'
Note: See TracChangeset for help on using the changeset viewer.