Changeset 1299 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell.F

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/thermcell.F (modified) (5 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell.F

@@ +1 @@
+!
+! $Id$
+!
       SUBROUTINE calcul_sec(ngrid,nlay,ptimestep
      s                  ,pplay,pplev,pphi,zlev
@@ -132 +135 @@
       character*10 str10
 
+      character (len=20) :: modname='calcul_sec'
+      character (len=80) :: abort_message
+
+
 !      LOGICAL vtest(klon),down
 
@@ -530 +537 @@
 c      write(8,*)zmax(ig),lmax(ig),lentr(ig),lmin(ig)
       enddo
-con stoppe après les calculs de zmax et wmax
+con stope après les calculs de zmax et wmax
       RETURN
 
@@ -776 +783 @@
          do ig=1,ngrid
             if(fracd(ig,l).lt.0.1) then
-               stop'fracd trop petit'
+              abort_message = 'fracd trop petit'
+              CALL abort_gcm (modname,abort_message,1)
+
             else
 c    vitesse descendante "diagnostique"
@@ -860 +869 @@
 cRC
       if (w2di.eq.1) then
-         fm0=fm0+ptimestep*(fm-fm0)/float(tho)
-         entr0=entr0+ptimestep*(entr-entr0)/float(tho)
+         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
+         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
       else
          fm0=fm