Changeset 1299 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/calltherm.F90

Timestamp:

Jan 20, 2010, 3:27:21 PM (14 years ago)

Author:

Laurent Fairhead

Message:

Nettoyage general pour se rapprocher des normes et éviter des erreurs a la
compilation:

• tous les FLOAT() sont remplacés par des REAL()
• tous les STOP dans phylmd sont remplacés par des appels à abort_gcm
• le common control défini dans le fichier control.h est remplacé par le module control_mod pour éviter des messages sur l'alignement des variables dans les déclarations
• des $Header$ remplacés par des $Id$ pour svn

Quelques remplacements à faire ont pu m'échapper

---

General cleanup of the code to try and adhere to norms and to prevent some
compilation errors:

• all FLOAT() instructions have been replaced by REAL() instructions
• all STOP instructions in phylmd have been replaced by calls to abort_gcm
• the common block control defined in the control.h file has been replaced by the control_mod to prevent compilation warnings on the alignement of declared variables
• $Header$ replaced by $Id$ for svn

Some changes which should have been made might have escaped me

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/calltherm.F90 (modified) (8 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/calltherm.F90

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
       subroutine calltherm(dtime  &
@@ -83 +83 @@
 !      save zentr_therm,zfm_therm
 
+      character (len=20) :: modname='calltherm'
+      character (len=80) :: abort_message
+
       integer i,k
       logical, save :: first=.true.
@@ -137 +140 @@
          if(nbptspb.GT.0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb   
 
-         zdt=dtime/float(nsplit_thermals)
+         zdt=dtime/REAL(nsplit_thermals)
          do isplit=1,nsplit_thermals
 
@@ -173 +176 @@
      &      ,tau_thermals,3)
           else if (iflag_thermals.eq.11) then
-            stop 'cas non prevu dans calltherm'
+              abort_message = 'cas non prevu dans calltherm'
+              CALL abort_gcm (modname,abort_message,1)
+
 !           CALL thermcell_pluie(klon,klev,zdt  &
 !   &      ,pplay,paprs,pphi,zlev  &
@@ -210 +215 @@
            if (prt_level.gt.10) write(lunout,*)'Apres thermcell_main OK'
          else
-            STOP'Cas des thermiques non prevu'
+           abort_message = 'Cas des thermiques non prevu'
+           CALL abort_gcm (modname,abort_message,1)
          endif
 
@@ -218 +224 @@
       DO i=1,klon
        logexpr1(i)=flag_bidouille_stratocu.or.weak_inversion(i).gt.0.5
-       IF(logexpr1(i)) fact(i)=1./float(nsplit_thermals)
+       IF(logexpr1(i)) fact(i)=1./REAL(nsplit_thermals)
       ENDDO
 
@@ -255 +261 @@
         if(prt_level.GE.10) print*,'calltherm i Alp_bl Alp Ale_bl Ale',i,Alp_bl(i),Alp(i),Ale_bl(i),Ale(i)
             fm_therm(i,klev+1)=0.
-            Ale_bl(i)=Ale_bl(i)+Ale(i)/float(nsplit_thermals)
+            Ale_bl(i)=Ale_bl(i)+Ale(i)/REAL(nsplit_thermals)
 !            write(22,*)'ALE CALLTHERM',Ale_bl(i),Ale(i)
-            Alp_bl(i)=Alp_bl(i)+Alp(i)/float(nsplit_thermals)
+            Alp_bl(i)=Alp_bl(i)+Alp(i)/REAL(nsplit_thermals)
 !            write(23,*)'ALP CALLTHERM',Alp_bl(i),Alp(i)
        ENDDO
@@ -276 +282 @@
 !    &         'fm=',zfm_therm(i,k),'entr=',entr_therm(i,k)
                  endif
-!       stop
             ENDDO
             ENDDO