Ignore:
Timestamp:
Jul 19, 2024, 4:15:44 PM (4 months ago)
Author:
abarral
Message:

Commit linked to correcting relevant warnings during gfortran compilation

[minor]
Correct remnants of fixed-form "+" in massbarxy.F90
Correct COMMON alignment in clesphys.h, flux_arp.h, cv*param.h, YOECUMF.h, alpale.h
Correct obsolete logical operators
Restrict use of iso_fortran_env in nf95_abort.f90
Remove redundant save in module declarations
Remove <continue> without labels
Fix nonstandard kind selectors

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/quickbeam.F90

    r3491 r5081  
    179179         
    180180          ! Attenuation due to gaseous absorption between radar and volume
    181           if ((rcfg%use_gas_abs == 1) .or. (rcfg%use_gas_abs == 2 .and. pr .eq. 1)) then
     181          if ((rcfg%use_gas_abs == 1) .or. (rcfg%use_gas_abs == 2 .and. pr == 1)) then
    182182             if (d_gate==1) then
    183183                if (k>1) then
     
    272272
    273273    ! Which platforms to create diagnostics for?
    274     if (platform .eq. 'cloudsat') lcloudsat=.true.
     274    if (platform == 'cloudsat') lcloudsat=.true.
    275275
    276276    ! Create Cloudsat diagnostics.
     
    289289             enddo
    290290          enddo
    291           where(cfad_ze .ne. R_UNDEF) cfad_ze = cfad_ze/Ncolumns
     291          where(cfad_ze /= R_UNDEF) cfad_ze = cfad_ze/Ncolumns
    292292
    293293          ! Compute cloudsat near-surface precipitation diagnostics
     
    306306             enddo
    307307          enddo
    308           where(cfad_ze .ne. R_UNDEF) cfad_ze = cfad_ze/Ncolumns
     308          where(cfad_ze /= R_UNDEF) cfad_ze = cfad_ze/Ncolumns
    309309       endif
    310310    endif
     
    402402       do pr=1,Ncolumns
    403403          ! 1) Compute the PIA in all profiles containing hydrometeors
    404           if ( (Ze_non_out(i,pr,cloudsat_preclvl).gt.-100) .and. (Ze_out(i,pr,cloudsat_preclvl).gt.-100) ) then
    405              if ( (Ze_non_out(i,pr,cloudsat_preclvl).lt.100) .and. (Ze_out(i,pr,cloudsat_preclvl).lt.100) ) then
     404          if ( (Ze_non_out(i,pr,cloudsat_preclvl)>-100) .and. (Ze_out(i,pr,cloudsat_preclvl)>-100) ) then
     405             if ( (Ze_non_out(i,pr,cloudsat_preclvl)<100) .and. (Ze_out(i,pr,cloudsat_preclvl)<100) ) then
    406406                cloudsat_precip_pia(i,pr) = Ze_non_out(i,pr,cloudsat_preclvl) - Ze_out(i,pr,cloudsat_preclvl)
    407407             endif
     
    412412          ! 2a) Oceanic points.
    413413          ! ################################################################################
    414           if (land(i) .eq. 0) then
     414          if (land(i) == 0) then
    415415!             print*, 'aaa i, pr, fracPrecipIce(i,pr) : ', i, pr, fracPrecipIce(i,pr) !Artem
    416416             ! Snow
    417              if(fracPrecipIce(i,pr).gt.0.9) then
    418                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(2)) then
     417             if(fracPrecipIce(i,pr)>0.9) then
     418                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(2)) then
    419419                   cloudsat_pflag(i,pr) = pClass_Snow2                   ! TSL: Snow certain
    420420                endif
    421                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(4).and. &
    422                      Ze_non_out(i,pr,cloudsat_preclvl).le.Zenonbinval(2)) then
     421                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(4).and. &
     422                     Ze_non_out(i,pr,cloudsat_preclvl)<=Zenonbinval(2)) then
    423423                   cloudsat_pflag(i,pr) = pClass_Snow1                   ! TSL: Snow possible
    424424                endif
     
    426426             
    427427             ! Mixed
    428              if(fracPrecipIce(i,pr).gt.0.1.and.fracPrecipIce(i,pr).le.0.9) then
    429                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(2)) then
     428             if(fracPrecipIce(i,pr)>0.1.and.fracPrecipIce(i,pr)<=0.9) then
     429                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(2)) then
    430430                   cloudsat_pflag(i,pr) = pClass_Mixed2                  ! TSL: Mixed certain
    431431                endif
    432                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(4).and. &
    433                      Ze_non_out(i,pr,cloudsat_preclvl).le.Zenonbinval(2)) then
     432                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(4).and. &
     433                     Ze_non_out(i,pr,cloudsat_preclvl)<=Zenonbinval(2)) then
    434434                   cloudsat_pflag(i,pr) = pClass_Mixed1                  ! TSL: Mixed possible
    435435                endif
     
    437437             
    438438             ! Rain
    439              if(fracPrecipIce(i,pr).le.0.1) then
    440                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(1)) then
     439             if(fracPrecipIce(i,pr)<=0.1) then
     440                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(1)) then
    441441                   cloudsat_pflag(i,pr) = pClass_Rain3                   ! TSL: Rain certain
    442442                endif
    443                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(3).and. &
    444                      Ze_non_out(i,pr,cloudsat_preclvl).le.Zenonbinval(1)) then
     443                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(3).and. &
     444                     Ze_non_out(i,pr,cloudsat_preclvl)<=Zenonbinval(1)) then
    445445                   cloudsat_pflag(i,pr) = pClass_Rain2                   ! TSL: Rain probable
    446446                endif
    447                 if(Ze_non_out(i,pr,cloudsat_preclvl).gt.Zenonbinval(4).and. &
    448                      Ze_non_out(i,pr,cloudsat_preclvl).le.Zenonbinval(3)) then
     447                if(Ze_non_out(i,pr,cloudsat_preclvl)>Zenonbinval(4).and. &
     448                     Ze_non_out(i,pr,cloudsat_preclvl)<=Zenonbinval(3)) then
    449449                   cloudsat_pflag(i,pr) = pClass_Rain1                   ! TSL: Rain possible
    450450                endif
    451                 if(cloudsat_precip_pia(i,pr).gt.40) then
     451                if(cloudsat_precip_pia(i,pr)>40) then
    452452                   cloudsat_pflag(i,pr) = pClass_Rain4                   ! TSL: Heavy Rain
    453453                endif
     
    455455             
    456456             ! No precipitation
    457              if(Ze_non_out(i,pr,cloudsat_preclvl).le.-15) then
     457             if(Ze_non_out(i,pr,cloudsat_preclvl)<=-15) then
    458458                cloudsat_pflag(i,pr) = pClass_noPrecip                   ! TSL: Not Raining
    459459             endif
     
    463463          ! 2b) Land points.
    464464          ! ################################################################################
    465           if (land(i) .eq. 1) then
     465          if (land(i) == 1) then
    466466             ! Find Zmax, the maximum reflectivity value in the attenuated profile (Ze_out);
    467467             Zmax=maxval(Ze_out(i,pr,:))
    468468
    469469             ! Snow (T<273)
    470              if(t2m(i) .lt. 273._wp) then
    471                 if(Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(5)) then
     470             if(t2m(i) < 273._wp) then
     471                if(Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(5)) then
    472472                   cloudsat_pflag(i,pr) = pClass_Snow2                      ! JEK: Snow certain
    473473                endif
    474                 if(Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(6) .and. &
    475                      Ze_out(i,pr,cloudsat_preclvl).le.Zbinvallnd(5)) then
     474                if(Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(6) .and. &
     475                     Ze_out(i,pr,cloudsat_preclvl)<=Zbinvallnd(5)) then
    476476                   cloudsat_pflag(i,pr) = pClass_Snow1                      ! JEK: Snow possible
    477477                endif
     
    479479             
    480480             ! Mized phase (273<T<275)
    481              if(t2m(i) .ge. 273._wp .and. t2m(i) .le. 275._wp) then
    482                 if ((Zmax .gt. Zbinvallnd(1) .and. cloudsat_precip_pia(i,pr).gt.30) .or. &
    483                      (Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(4))) then
     481             if(t2m(i) >= 273._wp .and. t2m(i) <= 275._wp) then
     482                if ((Zmax > Zbinvallnd(1) .and. cloudsat_precip_pia(i,pr)>30) .or. &
     483                     (Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(4))) then
    484484                   cloudsat_pflag(i,pr) = pClass_Mixed2                     ! JEK: Mixed certain
    485485                endif
    486                 if ((Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(6)  .and. &
    487                      Ze_out(i,pr,cloudsat_preclvl) .le. Zbinvallnd(4)) .and. &
    488                      (Zmax .gt. Zbinvallnd(5)) ) then
     486                if ((Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(6)  .and. &
     487                     Ze_out(i,pr,cloudsat_preclvl) <= Zbinvallnd(4)) .and. &
     488                     (Zmax > Zbinvallnd(5)) ) then
    489489                   cloudsat_pflag(i,pr) = pClass_Mixed1                     ! JEK: Mixed possible
    490490                endif
     
    492492
    493493             ! Rain (T>275)
    494              if(t2m(i) .gt. 275) then
    495                 if ((Zmax .gt. Zbinvallnd(1) .and. cloudsat_precip_pia(i,pr).gt.30) .or. &
    496                      (Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(2))) then
     494             if(t2m(i) > 275) then
     495                if ((Zmax > Zbinvallnd(1) .and. cloudsat_precip_pia(i,pr)>30) .or. &
     496                     (Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(2))) then
    497497                   cloudsat_pflag(i,pr) = pClass_Rain3                      ! JEK: Rain certain
    498498                endif
    499                 if((Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(6)) .and. &
    500                      (Zmax .gt. Zbinvallnd(3))) then
     499                if((Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(6)) .and. &
     500                     (Zmax > Zbinvallnd(3))) then
    501501                   cloudsat_pflag(i,pr) = pClass_Rain2                      ! JEK: Rain probable
    502502                endif
    503                 if((Ze_out(i,pr,cloudsat_preclvl) .gt. Zbinvallnd(6)) .and. &
    504                      (Zmax.lt.Zbinvallnd(3))) then
     503                if((Ze_out(i,pr,cloudsat_preclvl) > Zbinvallnd(6)) .and. &
     504                     (Zmax<Zbinvallnd(3))) then
    505505                   cloudsat_pflag(i,pr) = pClass_Rain1                      ! JEK: Rain possible
    506506                endif
    507                 if(cloudsat_precip_pia(i,pr).gt.40) then
     507                if(cloudsat_precip_pia(i,pr)>40) then
    508508                   cloudsat_pflag(i,pr) = pClass_Rain4                      ! JEK: Heavy Rain
    509509                endif
     
    511511             
    512512             ! No precipitation
    513              if(Ze_out(i,pr,cloudsat_preclvl).le.-15) then
     513             if(Ze_out(i,pr,cloudsat_preclvl)<=-15) then
    514514                cloudsat_pflag(i,pr) =  pClass_noPrecip                     ! JEK: Not Precipitating
    515515             endif         
     
    526526       ! Gridmean precipitation fraction for each precipitation type
    527527       do k=1,nCloudsatPrecipClass
    528           if (any(cloudsat_pflag(i,:) .eq. k-1)) then
    529              cloudsat_precip_cover(i,k) = count(cloudsat_pflag(i,:) .eq. k-1)
     528          if (any(cloudsat_pflag(i,:) == k-1)) then
     529             cloudsat_precip_cover(i,k) = count(cloudsat_pflag(i,:) == k-1)
    530530          endif
    531531       enddo
Note: See TracChangeset for help on using the changeset viewer.