Ignore:
Timestamp:
Sep 11, 2024, 4:27:07 PM (3 months ago)
Author:
abarral
Message:

Replace REPROBUS CPP KEY by logical using handmade wonky wrapper

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dates.F90

    r5159 r5185  
    136136!-------------------------------------------------
    137137
    138 if(imi == 0 .and. imiv == 0) then
     138if(imi == 0 .AND. imiv == 0) then
    139139
    140140!-------------------------------------------------
    … …  
    885885ia100=100*(iaaaa/100)
    886886ia4=4*(iaaaa/4)
    887 if((iaaaa == ia400).or.((iaaaa == ia4).and.(iaaaa /= ia100)))then
     887if((iaaaa == ia400).or.((iaaaa == ia4).AND.(iaaaa /= ia100)))then
    888888  ibissext=1
    889889else
    890890  ibissext=0
    891891endif
    892 if ((ibissext == 1).and.(imm > 2)) then
     892if ((ibissext == 1).AND.(imm > 2)) then
    893893  ijourp=1
    894894else
    … …  
    924924ia100=100*(iaaaa/100)
    925925ia4=4*(iaaaa/4)
    926 if((iaaaa == ia400).or.((iaaaa == ia4).and.(iaaaa /= ia100)))then
     926if((iaaaa == ia400).or.((iaaaa == ia4).AND.(iaaaa /= ia100)))then
    927927  ibissext=1
    928928else
    929929  ibissext=0
    930930endif
    931 if ((ibissext == 1).and.(imm > 2)) then
     931if ((ibissext == 1).AND.(imm > 2)) then
    932932  ijourp=1
    933933else
    … …  
    10501050ia100=100*(iaaaa/100)
    10511051ia4=4*(iaaaa/4)
    1052 if((iaaaa == ia400).or.((iaaaa == ia4).and.(iaaaa /= ia100)))then
     1052if((iaaaa == ia400).or.((iaaaa == ia4).AND.(iaaaa /= ia100)))then
    10531053  ibissext=1
    10541054else

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F

    r5133 r5185  
    88c NAN et INF ajoute aux plots                           Pat fin 2006
    99c ==================================================================
    10 c Comme dump2d sauf que le signe est préservé, la valeur zero
     10c Comme dump2d sauf que le signe est pr�serv�, la valeur zero
    1111c identifiee par un blanc.
    1212c detection des Infty (= ou -) et NaN (?)
    … …  
    7979      kchar(i)=16
    8080      GOTO 10022
    81 10021 IF(.NOT.((az.ne.0.and.icheck(1).eq.0.and.icheck(2).eq.2146435072))
     8110021 IF(.NOT.((az.ne.0.AND.icheck(1).eq.0.AND.icheck(2).eq.2146435072))
    8282     *)GOTO 10023
    8383      kchar(i)=31
    8484      zinf=.true.
    8585      GOTO 10022
    86 10023 IF(.NOT.((az.ne.0.and.icheck(1).eq.0.and.icheck(2).eq.2146959360))
     8610023 IF(.NOT.((az.ne.0.AND.icheck(1).eq.0.AND.icheck(2).eq.2146959360))
    8787     *)GOTO 10024
    8888      kchar(i)=32

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90

    r5159 r5185  
    170170integer(kind=jpim) :: num_c
    171171logical enough
    172 enough = (N > 2) .and. (a_ideal > 0)
     172enough = (N > 2) .AND. (a_ideal > 0)
    173173if( enough )then
    174174  num_c = max(1,nint((pi-2.*c_polar)/a_ideal))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwu.F90

    r5160 r5185  
    1 
    21! $Id$
    32
    … …  
    7473  !USE YOERDI   , ONLY : RCH4     ,RN2O     ,RCFC11   ,RCFC12
    7574  USE YOERDU, ONLY: R10E, REPSCO, REPSCQ
    76 #ifdef REPROBUS
    77 USE chem_rep, ONLY: rch42d, rn2o2d, rcfc112d, rcfc122d, ok_rtime2d
    78 USE infotrac_phy, ONLY : type_trac
    79 #endif
     75  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_REPROBUS
     76  USE lmdz_reprobus_wrappers, ONLY: rch42d, rn2o2d, rcfc112d, rcfc122d, ok_rtime2d
     77  USE infotrac_phy, ONLY: type_trac
    8078  USE lmdz_clesphys
    8179
    … …  
    319317        PABCU(JL, 17, IC) = PABCU(JL, 17, ICP1) + ZUAER(JL, 4) * ZDUC(JL, IC) * ZDIFF
    320318        PABCU(JL, 18, IC) = PABCU(JL, 18, ICP1) + ZUAER(JL, 5) * ZDUC(JL, IC) * ZDIFF
    321 #ifdef REPROBUS
    322         IF (type_trac=='repr'.and. ok_rtime2d) THEN
    323 !- CH4
    324       PABCU(JL,19,IC)=PABCU(JL,19,ICP1)&
    325        & + ZABLY(JL,2,IC)*RCH42D(JL, IC)/PCCO2*ZPHM6(JL)*ZDIFF
    326       PABCU(JL,20,IC)=PABCU(JL,20,ICP1)&
    327        & + ZABLY(JL,3,IC)*RCH42D(JL, IC)/PCCO2*ZPSM6(JL)*ZDIFF
    328 !- N2O
    329       PABCU(JL,21,IC)=PABCU(JL,21,ICP1)&
    330        & + ZABLY(JL,2,IC)*RN2O2D(JL, IC)/PCCO2*ZPHN6(JL)*ZDIFF
    331       PABCU(JL,22,IC)=PABCU(JL,22,ICP1)&
    332        & + ZABLY(JL,3,IC)*RN2O2D(JL, IC)/PCCO2*ZPSN6(JL)*ZDIFF
    333 !- CFC11
    334       PABCU(JL,23,IC)=PABCU(JL,23,ICP1)&
    335        & + ZABLY(JL,2,IC)*RCFC112D(JL, IC)/PCCO2        *ZDIFF
    336 !- CFC12
    337       PABCU(JL,24,IC)=PABCU(JL,24,ICP1)&
    338        & + ZABLY(JL,2,IC)*RCFC122D(JL, IC)/PCCO2        *ZDIFF
    339 
    340          ELSE
    341 #endif
    342         !- CH4
    343         PABCU(JL, 19, IC) = PABCU(JL, 19, ICP1)&
    344                 & + ZABLY(JL, 2, IC) * RCH4 / PCCO2 * ZPHM6(JL) * ZDIFF
    345         PABCU(JL, 20, IC) = PABCU(JL, 20, ICP1)&
    346                 & + ZABLY(JL, 3, IC) * RCH4 / PCCO2 * ZPSM6(JL) * ZDIFF
    347         !- N2O
    348         PABCU(JL, 21, IC) = PABCU(JL, 21, ICP1)&
    349                 & + ZABLY(JL, 2, IC) * RN2O / PCCO2 * ZPHN6(JL) * ZDIFF
    350         PABCU(JL, 22, IC) = PABCU(JL, 22, ICP1)&
    351                 & + ZABLY(JL, 3, IC) * RN2O / PCCO2 * ZPSN6(JL) * ZDIFF
    352         !- CFC11
    353         PABCU(JL, 23, IC) = PABCU(JL, 23, ICP1)&
    354                 & + ZABLY(JL, 2, IC) * RCFC11 / PCCO2 * ZDIFF
    355         !- CFC12
    356         PABCU(JL, 24, IC) = PABCU(JL, 24, ICP1)&
    357                 & + ZABLY(JL, 2, IC) * RCFC12 / PCCO2 * ZDIFF
    358 #ifdef REPROBUS
     319        IF (CPPKEY_REPROBUS .AND. type_trac=='repr'.AND. ok_rtime2d) THEN
     320          !- CH4
     321          PABCU(JL, 19, IC) = PABCU(JL, 19, ICP1)&
     322                  & + ZABLY(JL, 2, IC) * RCH42D(JL, IC) / PCCO2 * ZPHM6(JL) * ZDIFF
     323          PABCU(JL, 20, IC) = PABCU(JL, 20, ICP1)&
     324                  & + ZABLY(JL, 3, IC) * RCH42D(JL, IC) / PCCO2 * ZPSM6(JL) * ZDIFF
     325          !- N2O
     326          PABCU(JL, 21, IC) = PABCU(JL, 21, ICP1)&
     327                  & + ZABLY(JL, 2, IC) * RN2O2D(JL, IC) / PCCO2 * ZPHN6(JL) * ZDIFF
     328          PABCU(JL, 22, IC) = PABCU(JL, 22, ICP1)&
     329                  & + ZABLY(JL, 3, IC) * RN2O2D(JL, IC) / PCCO2 * ZPSN6(JL) * ZDIFF
     330          !- CFC11
     331          PABCU(JL, 23, IC) = PABCU(JL, 23, ICP1)&
     332                  & + ZABLY(JL, 2, IC) * RCFC112D(JL, IC) / PCCO2 * ZDIFF
     333          !- CFC12
     334          PABCU(JL, 24, IC) = PABCU(JL, 24, ICP1)&
     335                  & + ZABLY(JL, 2, IC) * RCFC122D(JL, IC) / PCCO2 * ZDIFF
     336
     337        ELSE
     338          !- CH4
     339          PABCU(JL, 19, IC) = PABCU(JL, 19, ICP1)&
     340                  & + ZABLY(JL, 2, IC) * RCH4 / PCCO2 * ZPHM6(JL) * ZDIFF
     341          PABCU(JL, 20, IC) = PABCU(JL, 20, ICP1)&
     342                  & + ZABLY(JL, 3, IC) * RCH4 / PCCO2 * ZPSM6(JL) * ZDIFF
     343          !- N2O
     344          PABCU(JL, 21, IC) = PABCU(JL, 21, ICP1)&
     345                  & + ZABLY(JL, 2, IC) * RN2O / PCCO2 * ZPHN6(JL) * ZDIFF
     346          PABCU(JL, 22, IC) = PABCU(JL, 22, ICP1)&
     347                  & + ZABLY(JL, 3, IC) * RN2O / PCCO2 * ZPSN6(JL) * ZDIFF
     348          !- CFC11
     349          PABCU(JL, 23, IC) = PABCU(JL, 23, ICP1)&
     350                  & + ZABLY(JL, 2, IC) * RCFC11 / PCCO2 * ZDIFF
     351          !- CFC12
     352          PABCU(JL, 24, IC) = PABCU(JL, 24, ICP1)&
     353                  & + ZABLY(JL, 2, IC) * RCFC12 / PCCO2 * ZDIFF
    359354        END IF
    360 #endif
    361355      ENDDO
    362356    ENDDO

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/recmwf_aero.F90

    r5159 r5185  
    592592
    593593     !--Case 4
    594      IF (ok_ade .and. ok_aie) THEN
     594     IF (ok_ade .AND. ok_aie) THEN
    595595
    596596        ! total aerosols for direct indirect effect
    … …  
    628628        LWDN_AERO(:,:,4) = PFLUX(:,2,:)
    629629
    630      ENDIF ! ok_ade .and. ok_aie
     630     ENDIF ! ok_ade .AND. ok_aie
    631631
    632632  ENDIF !--if flag_aerosol GT 0 OR flag_aerosol_strat
    … …  
    734734  IF  ( AEROSOLFEEDBACK_ACTIVE .AND. (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) ) THEN
    735735
    736      IF ( ok_ade .and. ok_aie  ) THEN
     736     IF ( ok_ade .AND. ok_aie  ) THEN
    737737        PFSUP(:,:) =    ZFSUP_AERO(:,:,4)
    738738        PFSDN(:,:) =    ZFSDN_AERO(:,:,4)
    … …  
    746746     ENDIF
    747747
    748      IF ( ok_ade .and. (.not. ok_aie) )  THEN
     748     IF ( ok_ade .AND. (.not. ok_aie) )  THEN
    749749        PFSUP(:,:) =    ZFSUP_AERO(:,:,3)
    750750        PFSDN(:,:) =    ZFSDN_AERO(:,:,3)
    … …  
    758758     ENDIF
    759759
    760      IF ( (.not. ok_ade) .and. ok_aie  )  THEN
     760     IF ( (.not. ok_ade) .AND. ok_aie  )  THEN
    761761        PFSUP(:,:) =    ZFSUP_AERO(:,:,2)
    762762        PFSDN(:,:) =    ZFSDN_AERO(:,:,2)
    … …  
    770770     ENDiF
    771771
    772      IF ((.not. ok_ade) .and. (.not. ok_aie)) THEN
     772     IF ((.not. ok_ade) .AND. (.not. ok_aie)) THEN
    773773        PFSUP(:,:) =    ZFSUP_AERO(:,:,1)
    774774        PFSDN(:,:) =    ZFSDN_AERO(:,:,1)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/rrtm_rtrn1a_140gp.F90

    r5160 r5185  
    383383!       & (1.0_JPRB - Z_CLDFRAC(I_LEV-1))   
    384384!    ENDIF
    385      if(istcld(i_lev).ne.1.and.i_lev.ne.1) then
     385     if(istcld(i_lev).ne.1.AND.i_lev.ne.1) then
    386386        z_faccmb1(i_lev+1) = max(0.,min(z_cldfrac(i_lev+1)-z_cldfrac(i_lev), &
    387387               z_cldfrac(i_lev-1)-z_cldfrac(i_lev)))
    … …  
    496496!    Z_FACCMB2D(I_LEV-1) = Z_FACCLD1D(I_LEV-1) * Z_FACCLR2D(I_LEV) *&
    497497!     & (1.0_JPRB - Z_CLDFRAC(I_LEV+1)) 
    498     if (istcldd(i_lev).ne.1.and.i_lev.ne.1) then
     498    if (istcldd(i_lev).ne.1.AND.i_lev.ne.1) then
    499499       z_faccmb1d(i_lev-1) = max(0.,min(z_cldfrac(i_lev+1)-z_cldfrac(i_lev), &
    500500                            z_cldfrac(i_lev-1)-z_cldfrac(i_lev)))
Note: See TracChangeset for help on using the changeset viewer.