Changeset 2100 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jul 22, 2014, 6:33:56 PM (10 years ago)

Author:

lguez

Message:

Removed "on rentre dans guide_main" from guide_main in dyn3dpar, was
already commented out in the dyn3dmem version.

Keeping length of lines under 80 characters in physiq (for
readability). Removed wrong comments "ajout des tendances de la
diffusion turbulente". Replaced "con" by "convection" as an argument
of add_phys_tend.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• hgardfou.F90 (modified) (2 diffs)
• physiq.F90 (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/hgardfou.F90

@@ -13 +13 @@
   include "iniprint.h"
   REAL t(klon, klev), tsol(klon, nbsrf)
-  CHARACTER *(*) text
+  CHARACTER(len=*), intent(in):: text
   CHARACTER (LEN=20) :: modname = 'hgardfou'
-  CHARACTER (LEN=80) :: abort_message
 
   INTEGER i, k, nsrf
@@ -129 +128 @@
   END DO
 
-  IF (.NOT. ok) THEN
-    abort_message = 'hgardfou s arrete ' // text
-    CALL abort_gcm(modname, abort_message, 1)
-  END IF
+  IF (.NOT. ok) CALL abort_gcm(modname, text, 1)
 
-  RETURN
 END SUBROUTINE hgardfou

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1929 +1929 @@
      ntra = 1
 
-     !=====================================================================================
-     !ajout pour la parametrisation des poches froides: 
-     !calcul de t_wake et t_undi: si pas de poches froides, t_wake=t_undi=t_seri 
+     !=========================================================================
+     !ajout pour la parametrisation des poches froides: calcul de
+     !t_wake et t_undi: si pas de poches froides, t_wake=t_undi=t_seri
      do k=1,klev
         do i=1,klon
@@ -1952 +1952 @@
      enddo
 
-     !c--   Calcul de l'energie disponible ALE (J/kg) et de la puissance disponible ALP (W/m2)
-     !c--    pour le soulevement des particules dans le modele convectif
+     ! Calcul de l'energie disponible ALE (J/kg) et de la puissance
+     ! disponible ALP (W/m2) pour le soulevement des particules dans
+     ! le modele convectif
      !
      do i = 1,klon
@@ -1981 +1982 @@
         enddo
      endif
-     !combinaison avec ale et alp de couche limite: constantes si pas de couplage, valeurs calculees
-     !dans le thermique sinon
+     !combinaison avec ale et alp de couche limite: constantes si pas
+     !de couplage, valeurs calculees dans le thermique sinon
      if (iflag_coupl.eq.0) then
         if (debut.and.prt_level.gt.9) &
@@ -2004 +2005 @@
         !         enddo
 
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
         ! Modif FH 2010/04/27. Sans doute temporaire.
-        ! Deux options pour le alp_offset : constant si >?? 0 ou proportionnel ??a
-        ! w si <0
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+        ! Deux options pour le alp_offset : constant si >?? 0 ou
+        ! proportionnel ??a w si <0
+        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
         do i = 1,klon
            ALE(i) = max(ale_wake(i),Ale_bl(i))
@@ -2026 +2027 @@
            endif
         enddo
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
 
@@ -2046 +2047 @@
 
      !fin calcul ale et alp
-     !=================================================================================================
+     !=======================================================================
 
 
@@ -2174 +2175 @@
   !    .              d_u_con, d_v_con)
 
-  !-----------------------------------------------------------------------------------------
-  ! ajout des tendances de la diffusion turbulente
-  CALL add_phys_tend(d_u_con,d_v_con,d_t_con,d_q_con,dql0,dqi0,paprs,'con')
-  !-----------------------------------------------------------------------------------------
+  CALL add_phys_tend(d_u_con, d_v_con, d_t_con, d_q_con, dql0, dqi0, paprs, &
+       'convection')
+  !----------------------------------------------------------------------------
 
   if (mydebug) then
@@ -2287 +2287 @@
           ,wake_ddeltat,wake_ddeltaq)
      !
-     !-----------------------------------------------------------------------------------------
+     !-------------------------------------------------------------------------
      ! ajout des tendances des poches froides
      ! Faire rapidement disparaitre l'ancien dt_wake pour garder un d_t_wake
@@ -2294 +2294 @@
      d_q_wake(:,:)=dq_wake(:,:)*dtime
      CALL add_phys_tend(du0,dv0,d_t_wake,d_q_wake,dql0,dqi0,paprs,'wake')
-     !-----------------------------------------------------------------------------------------
+     !------------------------------------------------------------------------
 
   endif
@@ -2551 +2551 @@
         endif
 
-        !-----------------------------------------------------------------------------------------
+        !-----------------------------------------------------------------------
         ! ajout des tendances de l'ajustement sec ou des thermiques
         CALL add_phys_tend(du0,dv0,d_t_ajsb,d_q_ajsb,dql0,dqi0,paprs,'ajsb')
@@ -2557 +2557 @@
         d_q_ajs(:,:)=d_q_ajs(:,:)+d_q_ajsb(:,:)
 
-        !-----------------------------------------------------------------------------------------
+        !---------------------------------------------------------------------
 
      endif
@@ -2609 +2609 @@
   WHERE (rain_lsc < 0) rain_lsc = 0.
   WHERE (snow_lsc < 0) snow_lsc = 0.
-  !-----------------------------------------------------------------------------------------
-  ! ajout des tendances de la diffusion turbulente
+
   CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs,'lsc')
-  !-----------------------------------------------------------------------------------------
+  !---------------------------------------------------------------------------
   DO k = 1, klev
      DO i = 1, klon
@@ -2740 +2739 @@
 #else
 
-              abort_message='You should compile with -rrtm if running with iflag_rrtm=1'
+              abort_message = 'You should compile with -rrtm if running ' &
+                   // 'with iflag_rrtm=1'
               call abort_gcm(modname,abort_message,1)
 #endif