Changeset 3763 for LMDZ6/branches/Optimisation_LMDZ

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:

• cv3_routines.F90 (modified) (10 diffs)
• cv3a_driver.f90 (modified) (8 diffs)

LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3_routines.F90

@@ -1101 +1101 @@
 
 !JAM--------------------------------------------------------------------
-! Calcul de la quantité d'eau sous forme de glace
+! Calcul de la quantit� d'eau sous forme de glace
 ! --------------------------------------------------------------------
   INTEGER nl, len
@@ -3072 +3072 @@
 !!---end jyg---
 
-! --------retour à la formulation originale d'Emanuel.
+! --------retour � la formulation originale d'Emanuel.
         IF (cvflag_ice) THEN
 
@@ -3082 +3082 @@
 
 !JAM  Attention: evap=sigt*E
-!    Modification: evap devient l'évaporation en milieu de couche
-!    car nécessaire dans cv3_yield
-!    Du coup, il faut modifier pas mal d'équations...
+!    Modification: evap devient l'�vaporation en milieu de couche
+!    car n�cessaire dans cv3_yield
+!    Du coup, il faut modifier pas mal d'�quations...
 !    et l'expression de afac qui devient afac1
 !    revap=sqrt((prec(i+1)+prec(i))/2)
@@ -3103 +3103 @@
 !JYG    Dans sa formulation originale, Emanuel calcule l'evaporation par:
 ! c             evap(il,i)=sigt*afac*revap
-! ce qui n'est pas correct. Dans cv_routines, la formulation a été modifiee.
+! ce qui n'est pas correct. Dans cv_routines, la formulation a �t� modifiee.
 ! Ici,l'evaporation evap est simplement calculee par l'equation de
 ! conservation.
@@ -3550 +3550 @@
 ! print*,'cv3_yield initialisation delt', delt
 
-  DO il = 1, ncum
-    precip(il) = 0.0
-    wd(il) = 0.0 ! gust
-  END DO
+  precip(:) = 0
+  wd(:) = 0
+  ft(:,:) = 0
+  fr(:,:) = 0
+  fu(:,:) = 0
+  fv(:,:) = 0
+  ma(:,:) = 0
+  mip(:,:) = 0
+  upwd(:,:) = 0
+  dnwd(:,:) = 0
+  dnwd0(:,:) = 0
+  qcondc(:,:) = 0
+  ftd(:,:) = 0
+  fqd(:,:) = 0
+  qtc(:,:) = 0
+  sigt(:,:) = 0
 
 !   Fluxes are on a staggered grid : loops extend up to nl+1
@@ -3560 +3572 @@
       Vprecip(il, i) = 0.0
       Vprecipi(il, i) = 0.0                               ! jyg
-      upwd(il, i) = 0.0
-      dnwd(il, i) = 0.0
-      dnwd0(il, i) = 0.0
-      mip(il, i) = 0.0
     END DO
   END DO
   DO i = 1, nl
     DO il = 1, ncum
-      ft(il, i) = 0.0
-      fr(il, i) = 0.0
-      fu(il, i) = 0.0
-      fv(il, i) = 0.0
-      ftd(il, i) = 0.0
-      fqd(il, i) = 0.0
-      qcondc(il, i) = 0.0 ! cld
       qcond(il, i) = 0.0 ! cld
-      qtc(il, i) = 0.0 ! cld
       qtment(il, i) = 0.0 ! cld
       sigment(il, i) = 0.0 ! cld
-      sigt(il, i) = 0.0 ! cld
       nqcond(il, i) = 0.0 ! cld
     END DO
@@ -4888 +4887 @@
     ENDDO       ! k
 
-! 14/01/15 AJ delta n'a rien à faire là...                                                 
+! 14/01/15 AJ delta n'a rien � faire l�...                                                 
     DO il = 1, ncum                                                  ! cld
 !!      IF (wa(il,i)>0.0 .AND. iflag(il)<=1) &                         ! cld
@@ -4904 +4903 @@
 
 ! IM cf. FH 
-! 14/01/15 AJ ne correspond pas à ce qui a été codé par JYG et SB           
+! 14/01/15 AJ ne correspond pas � ce qui a �t� cod� par JYG et SB           
                                                          
       IF (iflag_clw==0) THEN                                         ! cld
@@ -5101 +5100 @@
 
         ! On fait varier epmax en fn de la cape
-        ! Il faut donc recalculer ep, et hp qui a déjà été calculé et
-        ! qui en dépend
-        ! Toutes les autres variables fn de ep sont calculées plus bas.
+        ! Il faut donc recalculer ep, et hp qui a d�j� �t� calcul� et
+        ! qui en d�pend
+        ! Toutes les autres variables fn de ep sont calcul�es plus bas.
 
   include "cvthermo.h"
@@ -5138 +5137 @@
 
         ! il faut calculer la cape: on fait un calcule simple car tant qu'on ne
-        ! connait pas ep, on ne connait pas les mélanges, ddfts etc... qui sont
+        ! connait pas ep, on ne connait pas les m�langes, ddfts etc... qui sont
         ! necessaires au calcul de la cape dans la nouvelle physique
         

LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_driver.f90

@@ -634 +634 @@
     else !compress
       nloc = len
-      ktop1 = 0
-      precip1 = 0
-      cbmf1 = 0
-      plfc1 = 0
-      wbeff1 = 0
-      ptop21 = 0
-      sigd1 = 0
-      wd1 = 0
-      cape1 = Cape_noconv
-      cin1 = Cin_noconv
-      Plim11 = 0
-      plim21 = 0
-      supmax01 = 0
-      asupmaxmin1 = 0
-      epmax_diag1 = 0
-      ft1 = 0
-      fq1 = 0
-      fu1 = 0
-      fv1 = 0
-      ma1 = 0
-      mip1 = 0
-      upwd1 = 0
-      dnwd1 = 0
-      dnwd01 = 0
-      qcondc1 = 0
-      ftd1 = 0
-      fqd1 = 0
-      asupmax1 = 0
-      da1 = 0
-      mp1 = 0
-      d1a1 = 0
-      dam1 = 0
-      qta1 = 0
-      evap1 = 0
-      ep1 = 0
-      eplaMm1 = 0
-      wdtrainA1 = 0
-      wdtrainS1 = 0
-      wdtrainM1 = 0
-      qtc1 = 0
-      sigt1 = 0
-      vprecip1 = 0
-      vprecipi1 = 0
-      phi1 = 0
-      phi21 = 0
-      sigij1 = 0
-      elij1 = 0
-      epmlmMm1 = 0
       coef_convective(:) = merge(1, 0, iflag1(:) == 0)
       call cv3a_driver_compressed(nloc, nd, ntra, &
@@ -810 +762 @@
 
       call driver_log('cv3_undilute2')
+      qta(:,:) = 0
+      ep(:,:) = 0
       CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &
                          tnk, qnk, gznk, hnk, t, q, qs, gz, &
                          p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
                          inb, tp, tvp, clw, hp, ep, sigp, buoy, &
-                         frac_a, frac_s, qpreca, qta)
-
+                         frac_a, frac_s, qpreca, qta)      
+      
 ! epmax_cape
 ! on recalcule ep et hp
       call driver_log('cv3_epmax_cape')
+      epmax_diag(:) = 0
       call cv3_epmax_fn_cape(nloc, ncum, nd &
                              , ep, hp, icb, inb, clw, nk, t, h, hnk, lv, lf, frac_s &
@@ -828 +783 @@
 ! -------------------------------------------------------------------
       call enter_profile("cv3p_mixing")
+      elij(:,:,:) = 0
+      supmax(:,:) = 0
       IF (iflag_mix >= 1) THEN
-        CALL zilch(supmax, nloc*nd)
         call driver_log('cv3p_mixing')
         CALL cv3p_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &
@@ -836 +792 @@
                          ment, qent, hent, uent, vent, nent, &
                          sigij, elij, supmax, ments, qents, traent)
-      ELSE
-        CALL zilch(supmax, nloc*nd)
       END IF
       call exit_profile("cv3p_mixing")
@@ -845 +799 @@
 ! -------------------------------------------------------------------
 
+      cape(:) = -1
+      cin(:) = -100000.
       ptop2(:) = 0
       coef_clos(:) = 1.
+      Plim1(:) = 0
+      plim2(:) = 0
+      supmax0(:) = 0
+      asupmaxmin(:) = 0
+      cbmf(:) = 0
+      plfc(:) = 0
+      wbeff(:) = 0
+      asupmax(:,:) = 0
+
+      ok_inhib = (iflag_mix == 2)
+      
       IF (iflag_clos == 0) THEN
         call driver_log('cv3_closure')
@@ -852 +819 @@
                          pbase, p, ph, tv, buoy, &
                          sig, w0, cape, m, iflag)
-      END IF   ! iflag_clos==0
-
-      ok_inhib = iflag_mix == 2
+      END IF   ! iflag_clos==0     
 
       IF (iflag_clos == 1) PRINT *, ' pas d appel cv3p_closure'
@@ -958 +923 @@
       ENDIF
 
+    
+
 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 !--- passive tracers
@@ -964 +931 @@
       call driver_log('cv3_tracer')
       call enter_profile("cv3_tracer")
+ 
+      sigij(:,:,:) = 0    
+      vprecip(:,:) = 0
+      ! GLITCHY : vprecip is unused in cv3_tracer and sigij is intent in
       CALL cv3_tracer(nloc, -1, ncum, nd, nd, &
                       ment, sigij, da, phi, phi2, d1a, dam, &
                       ep, vprecip, elij, clw, epmlmMm, eplaMm, &
                       icb, inb)
+
+      vprecipi = 0
       call exit_profile("cv3_tracer")
     END IF ! ncum>0