Changeset 3763


Ignore:
Timestamp:
Jul 15, 2020, 10:14:36 PM (4 years ago)
Author:
adurocher
Message:

Remove useless zero-initializations

Location:
LMDZ6/branches/Optimisation_LMDZ/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3_routines.F90

    r3759 r3763  
    11011101
    11021102!JAM--------------------------------------------------------------------
    1103 ! Calcul de la quantité d'eau sous forme de glace
     1103! Calcul de la quantit d'eau sous forme de glace
    11041104! --------------------------------------------------------------------
    11051105  INTEGER nl, len
     
    30723072!!---end jyg---
    30733073
    3074 ! --------retour à la formulation originale d'Emanuel.
     3074! --------retour la formulation originale d'Emanuel.
    30753075        IF (cvflag_ice) THEN
    30763076
     
    30823082
    30833083!JAM  Attention: evap=sigt*E
    3084 !    Modification: evap devient l'évaporation en milieu de couche
    3085 !    car nécessaire dans cv3_yield
    3086 !    Du coup, il faut modifier pas mal d'équations...
     3084!    Modification: evap devient l'vaporation en milieu de couche
     3085!    car ncessaire dans cv3_yield
     3086!    Du coup, il faut modifier pas mal d'quations...
    30873087!    et l'expression de afac qui devient afac1
    30883088!    revap=sqrt((prec(i+1)+prec(i))/2)
     
    31033103!JYG    Dans sa formulation originale, Emanuel calcule l'evaporation par:
    31043104! c             evap(il,i)=sigt*afac*revap
    3105 ! ce qui n'est pas correct. Dans cv_routines, la formulation a été modifiee.
     3105! ce qui n'est pas correct. Dans cv_routines, la formulation a �t� modifiee.
    31063106! Ici,l'evaporation evap est simplement calculee par l'equation de
    31073107! conservation.
     
    35503550! print*,'cv3_yield initialisation delt', delt
    35513551
    3552   DO il = 1, ncum
    3553     precip(il) = 0.0
    3554     wd(il) = 0.0 ! gust
    3555   END DO
     3552  precip(:) = 0
     3553  wd(:) = 0
     3554  ft(:,:) = 0
     3555  fr(:,:) = 0
     3556  fu(:,:) = 0
     3557  fv(:,:) = 0
     3558  ma(:,:) = 0
     3559  mip(:,:) = 0
     3560  upwd(:,:) = 0
     3561  dnwd(:,:) = 0
     3562  dnwd0(:,:) = 0
     3563  qcondc(:,:) = 0
     3564  ftd(:,:) = 0
     3565  fqd(:,:) = 0
     3566  qtc(:,:) = 0
     3567  sigt(:,:) = 0
    35563568
    35573569!   Fluxes are on a staggered grid : loops extend up to nl+1
     
    35603572      Vprecip(il, i) = 0.0
    35613573      Vprecipi(il, i) = 0.0                               ! jyg
    3562       upwd(il, i) = 0.0
    3563       dnwd(il, i) = 0.0
    3564       dnwd0(il, i) = 0.0
    3565       mip(il, i) = 0.0
    35663574    END DO
    35673575  END DO
    35683576  DO i = 1, nl
    35693577    DO il = 1, ncum
    3570       ft(il, i) = 0.0
    3571       fr(il, i) = 0.0
    3572       fu(il, i) = 0.0
    3573       fv(il, i) = 0.0
    3574       ftd(il, i) = 0.0
    3575       fqd(il, i) = 0.0
    3576       qcondc(il, i) = 0.0 ! cld
    35773578      qcond(il, i) = 0.0 ! cld
    3578       qtc(il, i) = 0.0 ! cld
    35793579      qtment(il, i) = 0.0 ! cld
    35803580      sigment(il, i) = 0.0 ! cld
    3581       sigt(il, i) = 0.0 ! cld
    35823581      nqcond(il, i) = 0.0 ! cld
    35833582    END DO
     
    48884887    ENDDO       ! k
    48894888
    4890 ! 14/01/15 AJ delta n'a rien à faire là...                                                 
     4889! 14/01/15 AJ delta n'a rien � faire l�...                                                 
    48914890    DO il = 1, ncum                                                  ! cld
    48924891!!      IF (wa(il,i)>0.0 .AND. iflag(il)<=1) &                         ! cld
     
    49044903
    49054904! IM cf. FH
    4906 ! 14/01/15 AJ ne correspond pas à ce qui a été codé par JYG et SB           
     4905! 14/01/15 AJ ne correspond pas � ce qui a �t� cod� par JYG et SB           
    49074906                                                         
    49084907      IF (iflag_clw==0) THEN                                         ! cld
     
    51015100
    51025101        ! On fait varier epmax en fn de la cape
    5103         ! Il faut donc recalculer ep, et hp qui a déjà été calculé et
    5104         ! qui en dépend
    5105         ! Toutes les autres variables fn de ep sont calculées plus bas.
     5102        ! Il faut donc recalculer ep, et hp qui a d�j� �t� calcul� et
     5103        ! qui en dpend
     5104        ! Toutes les autres variables fn de ep sont calcules plus bas.
    51065105
    51075106  include "cvthermo.h"
     
    51385137
    51395138        ! il faut calculer la cape: on fait un calcule simple car tant qu'on ne
    5140         ! connait pas ep, on ne connait pas les mélanges, ddfts etc... qui sont
     5139        ! connait pas ep, on ne connait pas les mlanges, ddfts etc... qui sont
    51415140        ! necessaires au calcul de la cape dans la nouvelle physique
    51425141       
  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_driver.f90

    r3762 r3763  
    634634    else !compress
    635635      nloc = len
    636       ktop1 = 0
    637       precip1 = 0
    638       cbmf1 = 0
    639       plfc1 = 0
    640       wbeff1 = 0
    641       ptop21 = 0
    642       sigd1 = 0
    643       wd1 = 0
    644       cape1 = Cape_noconv
    645       cin1 = Cin_noconv
    646       Plim11 = 0
    647       plim21 = 0
    648       supmax01 = 0
    649       asupmaxmin1 = 0
    650       epmax_diag1 = 0
    651       ft1 = 0
    652       fq1 = 0
    653       fu1 = 0
    654       fv1 = 0
    655       ma1 = 0
    656       mip1 = 0
    657       upwd1 = 0
    658       dnwd1 = 0
    659       dnwd01 = 0
    660       qcondc1 = 0
    661       ftd1 = 0
    662       fqd1 = 0
    663       asupmax1 = 0
    664       da1 = 0
    665       mp1 = 0
    666       d1a1 = 0
    667       dam1 = 0
    668       qta1 = 0
    669       evap1 = 0
    670       ep1 = 0
    671       eplaMm1 = 0
    672       wdtrainA1 = 0
    673       wdtrainS1 = 0
    674       wdtrainM1 = 0
    675       qtc1 = 0
    676       sigt1 = 0
    677       vprecip1 = 0
    678       vprecipi1 = 0
    679       phi1 = 0
    680       phi21 = 0
    681       sigij1 = 0
    682       elij1 = 0
    683       epmlmMm1 = 0
    684636      coef_convective(:) = merge(1, 0, iflag1(:) == 0)
    685637      call cv3a_driver_compressed(nloc, nd, ntra, &
     
    810762
    811763      call driver_log('cv3_undilute2')
     764      qta(:,:) = 0
     765      ep(:,:) = 0
    812766      CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &
    813767                         tnk, qnk, gznk, hnk, t, q, qs, gz, &
    814768                         p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
    815769                         inb, tp, tvp, clw, hp, ep, sigp, buoy, &
    816                          frac_a, frac_s, qpreca, qta)
    817 
     770                         frac_a, frac_s, qpreca, qta)     
     771     
    818772! epmax_cape
    819773! on recalcule ep et hp
    820774      call driver_log('cv3_epmax_cape')
     775      epmax_diag(:) = 0
    821776      call cv3_epmax_fn_cape(nloc, ncum, nd &
    822777                             , ep, hp, icb, inb, clw, nk, t, h, hnk, lv, lf, frac_s &
     
    828783! -------------------------------------------------------------------
    829784      call enter_profile("cv3p_mixing")
     785      elij(:,:,:) = 0
     786      supmax(:,:) = 0
    830787      IF (iflag_mix >= 1) THEN
    831         CALL zilch(supmax, nloc*nd)
    832788        call driver_log('cv3p_mixing')
    833789        CALL cv3p_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &
     
    836792                         ment, qent, hent, uent, vent, nent, &
    837793                         sigij, elij, supmax, ments, qents, traent)
    838       ELSE
    839         CALL zilch(supmax, nloc*nd)
    840794      END IF
    841795      call exit_profile("cv3p_mixing")
     
    845799! -------------------------------------------------------------------
    846800
     801      cape(:) = -1
     802      cin(:) = -100000.
    847803      ptop2(:) = 0
    848804      coef_clos(:) = 1.
     805      Plim1(:) = 0
     806      plim2(:) = 0
     807      supmax0(:) = 0
     808      asupmaxmin(:) = 0
     809      cbmf(:) = 0
     810      plfc(:) = 0
     811      wbeff(:) = 0
     812      asupmax(:,:) = 0
     813
     814      ok_inhib = (iflag_mix == 2)
     815     
    849816      IF (iflag_clos == 0) THEN
    850817        call driver_log('cv3_closure')
     
    852819                         pbase, p, ph, tv, buoy, &
    853820                         sig, w0, cape, m, iflag)
    854       END IF   ! iflag_clos==0
    855 
    856       ok_inhib = iflag_mix == 2
     821      END IF   ! iflag_clos==0     
    857822
    858823      IF (iflag_clos == 1) PRINT *, ' pas d appel cv3p_closure'
     
    958923      ENDIF
    959924
     925   
     926
    960927!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    961928!--- passive tracers
     
    964931      call driver_log('cv3_tracer')
    965932      call enter_profile("cv3_tracer")
     933 
     934      sigij(:,:,:) = 0   
     935      vprecip(:,:) = 0
     936      ! GLITCHY : vprecip is unused in cv3_tracer and sigij is intent in
    966937      CALL cv3_tracer(nloc, -1, ncum, nd, nd, &
    967938                      ment, sigij, da, phi, phi2, d1a, dam, &
    968939                      ep, vprecip, elij, clw, epmlmMm, eplaMm, &
    969940                      icb, inb)
     941
     942      vprecipi = 0
    970943      call exit_profile("cv3_tracer")
    971944    END IF ! ncum>0
Note: See TracChangeset for help on using the changeset viewer.