Changeset 1299 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_old.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_old.F (modified) (19 diffs)

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

@@ -112 +112 @@
       character (len=10) :: str10
 
+      character (len=20) :: modname='thermcell2002'
+      character (len=80) :: abort_message
+
       LOGICAL vtest(klon),down
 
@@ -336 +339 @@
             if(w2di.eq.2) then
                entr(ig,k)=entr(ig,k)+
-     s         ptimestep*(zzz-entr(ig,k))/float(tho)
+     s         ptimestep*(zzz-entr(ig,k))/REAL(tho)
             else
                entr(ig,k)=zzz
@@ -379 +382 @@
 c     print*,'ig,l+1,ztv(ig,l+1)'
 c     print*, ig,l+1,ztv(ig,l+1)
-c        stop'dans thermiques'
 c     endif
                zw2(ig,l+1)=2.*RG*(ztv(ig,l)-ztv(ig,l+1))/ztv(ig,l+1)
@@ -395 +397 @@
 c     print*,'Tv ',(ztv(ig,ll),ll=1,klev)
 c     print*,'Entr ',(entr(ig,ll),ll=1,klev)
-c        stop'dans thermiques'
 c     endif
                ztva(ig,l)=(fmc(ig,l)*ztva(ig,l-1)+entr(ig,l)*ztv(ig,l))
@@ -517 +518 @@
          do ig=1,ngrid
             if(fracd(ig,l).lt.0.1) then
-               stop'fracd trop petit'
-            else
+              abort_message = 'fracd trop petit'
+              CALL abort_gcm (modname,abort_message,1) 
+           else
 c    vitesse descendante "diagnostique"
                wd(ig,l)=fm(ig,l)/(fracd(ig,l)*rhobarz(ig,l))
@@ -588 +590 @@
 
       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
@@ -1000 +1002 @@
       character*2 str2
       character*10 str10
+
+      character (len=20) :: modname='thermcell_cld'
+      character (len=80) :: abort_message
 
       LOGICAL vtest(klon),down
@@ -1855 +1860 @@
        if (l.eq.klev) then
           print*,'THERMCELL PB ig=',ig,'   l=',l
-          stop
+          abort_message = 'THERMCELL PB'
+          CALL abort_gcm (modname,abort_message,1)
        endif
 !       if ((zw2(ig,l+1).gt.1.e-10).and.(zw2(ig,l).gt.1.e-10).and.
@@ -2164 +2170 @@
          do ig=1,ngrid
             if(fracd(ig,l).lt.0.1.and.(test(ig).eq.1)) then
-               stop'fracd trop petit'
+              abort_message = 'fracd trop petit'
+              CALL abort_gcm (modname,abort_message,1)
             else
 c    vitesse descendante "diagnostique"
@@ -2262 +2269 @@
 
       if (w2di.eq.1) then
-         fm0=fm0+ptimestep*(fm-fm0)/float(tho)
-         entr0=entr0+ptimestep*(alim+entr-entr0)/float(tho)
+         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
+         entr0=entr0+ptimestep*(alim+entr-entr0)/REAL(tho)
       else
          fm0=fm
@@ -2747 +2754 @@
       character*10 str10
 
+      character (len=20) :: modname='thermcell_eau'
+      character (len=80) :: abort_message
+
       LOGICAL vtest(klon),down
       LOGICAL Zsat(klon)
@@ -3410 +3420 @@
          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"
@@ -3481 +3492 @@
 
       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
@@ -3848 +3859 @@
       character*10 str10
 
+      character (len=20) :: modname='thermcell'
+      character (len=80) :: abort_message
+
       LOGICAL vtest(klon),down
 
@@ -4394 +4408 @@
          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"
@@ -4477 +4492 @@
 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
@@ -5257 +5272 @@
       character*10 str10
 
+      character (len=20) :: modname='thermcell_sec'
+      character (len=80) :: abort_message
+
       LOGICAL vtest(klon),down
 
@@ -5822 +5840 @@
          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"
@@ -5905 +5924 @@
 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