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

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

@@ -39 +39 @@
       use conf_phys_m, only: conf_phys
       use radlwsw_m, only: radlwsw
+      USE control_mod
+
 
       IMPLICIT none
@@ -97 +99 @@
 #include "dimsoil.h"
 #include "clesphys.h"
-#include "control.h"
 #include "temps.h"
 #include "iniprint.h"
@@ -1348 +1349 @@
          ENDIF
 c
-         IF (dtime*FLOAT(radpas).GT.21600..AND.cycle_diurne) THEN 
+         IF (dtime*REAL(radpas).GT.21600..AND.cycle_diurne) THEN 
            WRITE(lunout,*)'Nbre d appels au rayonnement insuffisant'
            WRITE(lunout,*)"Au minimum 4 appels par jour si cycle diurne"
@@ -1518 +1519 @@
          CALL VTb(VTinca)
 !         iii = MOD(NINT(xjour),360)
-!         calday = FLOAT(iii) + jH_cur
-         calday = FLOAT(days_elapsed) + jH_cur
+!         calday = REAL(iii) + jH_cur
+         calday = REAL(days_elapsed) + jH_cur
          WRITE(lunout,*) 'initial time chemini', days_elapsed, calday
 
@@ -1801 +1802 @@
 !   solarlong0
         if (solarlong0<-999.) then
-           CALL orbite(FLOAT(days_elapsed+1),zlongi,dist)
+           CALL orbite(REAL(days_elapsed+1),zlongi,dist)
         else
            zlongi=solarlong0  ! longitude solaire vraie
@@ -1812 +1813 @@
 !  Avec ou sans cycle diurne
       IF (cycle_diurne) THEN
-        zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
+        zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
         CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract)
       ELSE
@@ -1961 +1962 @@
 
       IF (iflag_con.EQ.1) THEN
-          stop'reactiver le call conlmd dans physiq.F'
+        abort_message ='reactiver le call conlmd dans physiq.F'
+        CALL abort_gcm (modname,abort_message,1)
 c     CALL conlmd (dtime, paprs, pplay, t_seri, q_seri, conv_q,
 c    .             d_t_con, d_q_con,
@@ -2205 +2207 @@
           za = 0.0
           DO i = 1, klon
-            za = za + airephy(i)/FLOAT(klon)
+            za = za + airephy(i)/REAL(klon)
             zx_t = zx_t + (rain_con(i)+
-     .                   snow_con(i))*airephy(i)/FLOAT(klon)
+     .                   snow_con(i))*airephy(i)/REAL(klon)
           ENDDO
           zx_t = zx_t/za*dtime
@@ -2599 +2601 @@
          za = 0.0
          DO i = 1, klon
-            za = za + airephy(i)/FLOAT(klon)
+            za = za + airephy(i)/REAL(klon)
             zx_t = zx_t + (rain_lsc(i)
-     .                  + snow_lsc(i))*airephy(i)/FLOAT(klon)
+     .                  + snow_lsc(i))*airephy(i)/REAL(klon)
         ENDDO
          zx_t = zx_t/za*dtime
@@ -2823 +2825 @@
          CALL VTe(VTphysiq)
          CALL VTb(VTinca)
-         calday = FLOAT(days_elapsed + 1) + jH_cur
+         calday = REAL(days_elapsed + 1) + jH_cur
 
          call chemtime(itap+itau_phy-1, date0, dtime)