Changeset 4085 for trunk/LMDZ.VENUS/libf


Ignore:
Timestamp:
Feb 26, 2026, 10:20:22 AM (5 weeks ago)
Author:
emillour
Message:

Venus PCM:
Follow-up of commit r4080: some adaptations for the gfortran compiler which
is more strict than ifort: logicals should not be compared with ".eq." and
array sizes must be integers.
EM

Location:
trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/evap.f90

    r4080 r4085  
    66  IMPLICIT NONE
    77
    8   real, intent(in) :: nq
     8  integer, intent(in) :: nq
    99  real, intent(in) :: RHOSA
    1010  real, intent(inout), dimension(nq) :: ztrac   ! Zeroth moment of mode 1
    … …  
    4343  IMPLICIT NONE
    4444
    45   real, intent(in) :: nq
     45  integer, intent(in) :: nq
    4646  real, intent(in) :: RHOSA
    4747  real, intent(inout), dimension(nq) :: ztrac   ! Zeroth moment of mode 1

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/flux.f90

    r4080 r4085  
    155155    ENDIF
    156156
    157     IF (checkm1 .eq. .false.) then 
     157    IF (checkm1 .eqv. .false.) then 
    158158                k1 =  dt_eff*dM3_dt(m3sa_m1,              m3w_m1, m3sa_m2, M3_m1(3), M0_m1(1), mmean, MSA, RHOAIR, PAIR, MRSATOTtemp, sigma1, RHOSA, RHOsasat, ST, RGAS, TAIR, RDSA+RCSA)
    159159               
    … …  
    168168    ENDIF
    169169
    170     IF (checkm2 .eq. .false.) then
     170    IF (checkm2 .eqv. .false.) then
    171171                k1 =  dt_eff*dM3_dt(m3sa_m2,              m3w_m2, m3sa_m1, M3_m2(3), M0_m2(1), mmean, MSA, RHOAIR, PAIR, MRSATOTtemp, sigma2, RHOSA, RHOsasat, ST, RGAS, TAIR, RDSA+RCSA)
    172172               
    … …  
    181181    ENDIF
    182182
    183     IF (checkm1 .eq. .false.) then 
     183    IF (checkm1 .eqv. .false.) then 
    184184               
    185185               
    … …  
    227227                s1pred1 = s1pred2
    228228    ENDIF
    229     IF (checkm2 .eq. .false.) then
     229    IF (checkm2 .eqv. .false.) then
    230230
    231231               

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/init_moment.f90

    r4080 r4085  
    2020          real    :: M3_m1sa(nlev)  , M3_m2sa(nlev)
    2121          real    :: M3_m1wv(nlev)  , M3_m2wv(nlev)
    22           real    :: r1i,r1j,r2i,r2j,J,ROSA,dummy,WSAM
     22          real    :: r1i,r1j,r2i,r2j,ROSA,dummy,WSAM
     23          integer :: J
    2324          real    :: KMASS_mode1(nlev), KMASS_mode2(nlev)
    2425          real    :: RHOAIR(nlev)

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/mad_muphy.f90

    r4080 r4085  
    8585
    8686  !-------------
    87   IF (debug_cloud .eq. .false.) then
     87  IF (debug_cloud .eqv. .false.) then
    8888          dbg_cons = .false.
    8989          dbg_coag = .false.
    … …  
    201201  ! obtenues a l'equilibre
    202202 
    203   IF (dbg_glob .eq. .true.) then
     203  IF (dbg_glob) then
    204204  write(*,*)'RHOpi et RHOeqi', RHOSAp, RHOSA   
    205205  ENDIF
    206206  call change_wsa(TAIR,PAIR,MRSAv,MRSAss,WSAv)
    207207 
    208   IF (dbg_glob .eq. .true.) then
     208  IF (dbg_glob) then
    209209  write(*,*)'WSAp, WSAVin, WSAeq, WSAVeq', WSAp,WSAVin, WSA, WSAV
    210210  write(*,*)'RHOp et RHOeq', RHOSAp, RHOSA   
    … …  
    226226  IF((NUCLEA) .AND. (SH2SO4 .gt. 1)) then ! HOmogeneous NUcleation with MOment
    227227 
    228       IF (dbg_nucl .eq. .true.) write(*,*)'[NUCLEA] Here in NUCLEA (hom)', NUCLEA
     228      IF (dbg_nucl) write(*,*)'[NUCLEA] Here in NUCLEA (hom)', NUCLEA
    229229      IF ((TAIR .GT. 165.e0) .and. (TAIR .LT. 400.e0) .and. (h2so4_m3 .GT. 1.e10) .and. (h2so4_m3 .LT. 1.e19)  .and.  (RH .gt. 1e-5) .AND. (RH .LT. 1.D0)) then
    230230         CALL HONUMO(TAIR,PAIR,dt,M0_m1,M3_m1,dM0_hom,dM3_hom) ! #.m(air)-3 and m3.m(air)-3
    … …  
    238238  IF((HETNUCLEA) .AND. (SH2SO4 .gt. 1)) then ! HETerogeneous NUcleation with MOment
    239239   
    240      IF (dbg_nucl .eq. .true.) write(*,*)'[NUCLEA] Here in NUCLEA (het)', HETNUCLEA
     240     IF (dbg_nucl) write(*,*)'[NUCLEA] Here in NUCLEA (het)', HETNUCLEA
    241241   
    242242   ! Calculation of ra (mean radius of aerosols) before heterogeneous nucleation
    … …  
    392392     dM3_m1_ccn = 0.D0
    393393
    394      IF (dbg_tend .eq. .true.) THEN
     394     IF (dbg_tend) THEN
    395395        write(*,*)'[TEND  ] Tendencies calculation'
    396396        write(*,*)'[TEND  ] SUM tendance M3(m1) SA/WV ',dM3_m1_SA,dM3_m1_WV
    … …  
    406406   
    407407
    408      IF (dbg_tend .eq. .true.) THEN
     408     IF (dbg_tend) THEN
    409409        write(*,*)'[TEND  ] Tendencies calculation'
    410410        write(*,*)'[TEND  ] SUM tendance M3(m2) SA/WV ',dM3_m2_SA,dM3_m2_WV
    … …  
    514514     dvmr_WV_m2 = dvmr_WV_m2 * dvmr_WV
    515515
    516      IF (dbg_tend .eq. .true.) THEN
     516     IF (dbg_tend) THEN
    517517        write(*,*)'[TEND  ] M3_m1  (SA, W, vmr)', dvmr_SA_m1, dvmr_WV_m1
    518518        write(*,*)'[TEND  ] M3_m2  (SA, W, vmr)', dvmr_SA_m2, dvmr_WV_m2
    … …  
    608608     ! Partition vapeur
    609609     
    610      IF (dbg_cons .eq. .true.) THEN
     610     IF (dbg_cons) THEN
    611611        write(*,*)'[VapUpdate prev] MRSAv, MRWVv',MRSAv, MRWVv
    612612        write(*,*)'[VapUpdate prev] MRSAl, MRWVl',MRSAl, MRWVl
    … …  
    620620     MRSAl = MRSAl + dvmr_SA ! acid
    621621     MRWVl = MRWVl + dvmr_WV ! water
    622      IF (dbg_cons .eq. .true.) THEN
     622     IF (dbg_cons) THEN
    623623        write(*,*)'[VapUpdate aft] MRSAv, MRWVv',MRSAv, MRWVv
    624624        write(*,*)'[VapUpdate aft] MRSAl, MRWVl',MRSAl, MRWVl
    … …  
    677677     diff = abs(check - MRTOT)  ! any differences ?
    678678
    679      IF (dbg_tend .eq. .true.) THEN
     679     IF (dbg_tend) THEN
    680680        write(*,*)'[TEND  ] vap: PI, RHOSA',PI, RHOSA
    681681        write(*,*)'[TEND  ] vap: CSA, CWV, RHOAIR', RHOAIR
    … …  
    684684     ENDIF
    685685
    686      IF (dbg_cons .eq. .true.) THEN
     686     IF (dbg_cons) THEN
    687687        write(*,*)'[CONS  ] WSA AFTER',WSAV, WSA
    688688        write(*,*)'[CONS  ] MRSA V AFTER', MRSAv
    … …  
    715715
    716716     ! Changements des moments effectues directement dans la subroutines ?
    717      IF (dbg_coag .eq. .true.) THEN
     717     IF (dbg_coag) THEN
    718718        write(*,*)'[COAG  ] Here in COAG', COAG
    719719        write(*,*)'[COAG  ] r1,r2 tot BEF COAG',r1,r2
    … …  
    732732          &  dM3_m1_coag,dM0_m2_coag,dM3_m2_coag)
    733733     
    734      IF (dbg_coag .eq. .true.) THEN
     734     IF (dbg_coag) THEN
    735735        write(*,*)'[COAG  ] M0 tot BEF UP', M0_m1(1),M0_m2(1), M0_m1(1)+M0_m2(1)
    736736        write(*,*)'[COAG  ] M3 tot BEF UP', M3_m1(1)+M3_m1(2),M3_m2(1)+M3_m2(2), &
    … …  
    763763!     somme = dM3_m1_coag(1)+dM3_m1_coag(2)+dM3_m1_coag(3)+dM3_m2_coag(1)+dM3_m2_coag(2)+dM3_m2_coag(3)
    764764!     write(*,*),'coag',somme
    765      IF (dbg_coag .eq. .true.) THEN
     765     IF (dbg_coag) THEN
    766766        write(*,*)'[COAG  ] M0 tot AFT UP', M0_m1(1),M0_m2(1), M0_m1(1)+M0_m2(1)
    767767        write(*,*)'[COAG  ] M3 tot AFT UP', M3_m1(1)+M3_m1(2),M3_m2(1)+M3_m2(2), &
    … …  
    808808     IF ((M0_m1(1) .gt. 1e-15) .and. (somme .gt. 1e-30) .and. (r_merge .gt. (ri1*ri2)**0.5)) THEN
    809809!        write(*,*)'[COAG, MERGING  ] r1,r2 BEF MERGING',r1, r2
    810         IF (dbg_coag .eq. .true.) THEN
     810        IF (dbg_coag) THEN
    811811           write(*,*)'[COAG, MERGING  ] Mode merging 1 vers 2'
    812812        ENDIF
    … …  
    859859      IF ((M0_m2(1) .gt. 1e-15) .and. (somme .gt. 1e-30) .and. (r_merge .lt. (ri1*ri2)**0.5)) THEN
    860860!        write(*,*)'[COAG, MERGING  ] r1,r2 BEF MERGING',r1, r2
    861          IF (dbg_coag .eq. .true.) THEN
     861         IF (dbg_coag) THEN
    862862            write(*,*)'[COAG, MERGING  ] Mode merging 2 vers 1',abs(sum(dM3_merge_m1)+sum(dM3_merge_m2))
    863863         ENDIF
Note: See TracChangeset for help on using the changeset viewer.