Changeset 3763 for LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3_routines.F90

Timestamp:

Jul 15, 2020, 10:14:36 PM (4 years ago)

Author:

adurocher

Message:

Remove useless zero-initializations

File:

• LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3_routines.F90 (modified) (10 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