Ignore:
Timestamp:
Apr 14, 2017, 12:38:26 PM (7 years ago)
Author:
jyg
Message:

Changing from a binary code to a decimal code:
iflag_pbl_split= 0-> no splitting;
1-> vdf splitting;
10-> thermals splitting;
11-> full splitting

File:
1 edited

Legend:

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

    r2851 r2852  
    21372137
    21382138       !jyg+nrlmd<
    2139        IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,2) .eq. 1) THEN
     2139!!jyg       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,2) .eq. 1) THEN
     2140       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,10) .ge. 1) THEN
    21402141          print *,'debut du splitting de la PBL'
    21412142       ENDIF
     
    22172218       !
    22182219       !  Add turbulent diffusion tendency to the wake difference variables
    2219        IF (mod(iflag_pbl_split,2) .NE. 0) THEN
     2220!!jyg       IF (mod(iflag_pbl_split,2) .NE. 0) THEN
     2221       IF (mod(iflag_pbl_split,10) .NE. 0) THEN
    22202222!jyg<
    22212223          d_deltat_vdf(:,:) = d_t_vdf_w(:,:)-d_t_vdf_x(:,:)
     
    28292831       IF (iflag_thermals>=1) THEN
    28302832          !jyg<
    2831           IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
     2833!!       IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
     2834       IF (mod(iflag_pbl_split/10,10) .GE. 1) THEN
    28322835             !  Appel des thermiques avec les profils exterieurs aux poches
    28332836             DO k=1,klev
     
    28742877          !
    28752878          !jyg<
    2876           IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
     2879!!jyg          IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
     2880          IF (mod(iflag_pbl_split/10,10) .GE. 1) THEN
    28772881             !  Si les thermiques ne sont presents que hors des
    28782882             !  poches, la tendance moyenne associ\'ee doit etre
     
    28962900             call prt_enerbil('the',itap)
    28972901          !
    2898           ENDIF  ! (mod(iflag_pbl_split/2,2) .EQ. 1)
     2902          ENDIF  ! (mod(iflag_pbl_split/10,10) .GE. 1)
    28992903          !
    29002904          CALL add_phys_tend(d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs,  &
Note: See TracChangeset for help on using the changeset viewer.