Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diagedyn.f90

Timestamp:

Aug 2, 2024, 9:58:25 PM (7 weeks ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diagedyn.f90 (modified) (21 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diagedyn.f90

@@ -6 +6 @@
         , ucov    , vcov , ps, p ,pk , teta , q, ql)
   !======================================================================
-  !
+
   ! Purpose:
   !    Calcul la difference d'enthalpie et de masse d'eau entre 2 appels,
@@ -14 +14 @@
   !    Outil pour diagnostiquer la conservation de l'energie
   !    et de la masse dans la dynamique.
-  !
-  !
+
+
   !======================================================================
   ! Arguments:
@@ -36 +36 @@
   ! ql-------input-R- liquid watter (kg/kg)
   ! aire-----input-R- mesh surafce (m2)
-  !
+
   ! the following total value are computed by UNIT of earth surface
-  !
+
   ! d_h_vcol--output-R- Heat flux (W/m2) define as the Enthalpy
   !        change (J/m2) during one time step (dtime) for the whole
@@ -47 +47 @@
   ! d_ql------output-R- same, for the liquid watter only (kg/m2/s)
   ! d_ec------output-R- Cinetic Energy Budget (W/m2) for vertical air column
-  !
-  !
+
+
   ! J.L. Dufresne, July 2002
   !======================================================================
@@ -56 +56 @@
   USE lmdz_comgeom
 
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
+
+
 
   ! Ehouarn: for now set these parameters to what is in Earth physics...
@@ -73 +75 @@
   REAL,PARAMETER :: RLVTT=2.5008E+6
 
-  !
+
   INTEGER :: imjmp1
   PARAMETER( imjmp1=iim*jjp1)
@@ -91 +93 @@
   ! Output variables
   REAL :: d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec
-  !
+
   ! Local variables
-  !
+
   REAL :: h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot &
         , h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot
@@ -106 +108 @@
   ! qs_tot------  total mass of solid watter (kg/m2)
   ! ec_tot------  total cinetic energy (kg/m2)
-  !
+
   REAL :: masse(ip1jmp1,llm)                ! masse d'air
   REAL :: vcont(ip1jm,llm),ucont(ip1jmp1,llm)
@@ -129 +131 @@
   REAL :: zh_dair_col(imjmp1)
   REAL :: zh_qw_col(imjmp1), zh_ql_col(imjmp1), zh_qs_col(imjmp1)
-  !
+
   REAL :: d_h_dair, d_h_qw, d_h_ql, d_h_qs
-  !
+
   REAL :: airetot, zcpvap, zcwat, zcice
-  !
+
   INTEGER :: i, k, jj, ij , l ,ip1jjm1
-  !
+
   INTEGER :: ndiag     ! max number of diagnostic in parallel
   PARAMETER (ndiag=10)
@@ -141 +143 @@
   save pas
   data pas/ndiag*0/
-  !
+
   REAL :: h_vcol_pre(ndiag), h_dair_pre(ndiag), h_qw_pre(ndiag) &
         , h_ql_pre(ndiag), h_qs_pre(ndiag), qw_pre(ndiag) &
@@ -157 +159 @@
   CALL massdair( p, masse )
   !======================================================================
-  !
-  !
+
+
   PRINT*,'MAIS POURQUOI DONC DIAGEDYN NE MARCHE PAS ?'
   RETURN
@@ -187 +189 @@
     ENDDO
   ENDDO
-  !
+
   ! Reset variables
   DO i = 1, imjmp1
@@ -199 +201 @@
     zh_qs_col(i) = 0.
   ENDDO
-  !
+
   zcpvap=RCPV
   zcwat=RCW
   zcice=RCS
-  !
+
   ! Compute vertical sum for each atmospheric column
   ! ================================================
@@ -230 +232 @@
     END DO
   ENDDO
-  !
+
   ! Mean over the planete surface
   ! =============================
@@ -243 +245 @@
   h_qs_tot = 0.
   airetot=0.
-  !
+
   DO i=1,imjmp1
     qw_tot = qw_tot + zqw_col(i)
@@ -255 +257 @@
     airetot=airetot+zaire(i)
   END DO
-  !
+
   qw_tot = qw_tot/airetot
   ql_tot = ql_tot/airetot
@@ -264 +266 @@
   h_ql_tot = h_ql_tot/airetot
   h_qs_tot = h_qs_tot/airetot
-  !
+
   h_vcol_tot = h_dair_tot+h_qw_tot+h_ql_tot+h_qs_tot
-  !
+
   ! Compute the change of the atmospheric state compare to the one
   ! stored in "idiag2", and convert it in flux. THis computation
@@ -272 +274 @@
   ! for "idiag"
   ! ===================================
-  !
+
   IF ( (idiag2>0) .AND. (pas(idiag2) /= 0) ) THEN
     d_h_vcol  = (h_vcol_tot - h_vcol_pre(idiag2) )/dtime
@@ -296 +298 @@
     d_qt     = 0.
   ENDIF
-  !
+
   IF (iprt>=2) THEN
     WRITE(6,9000) tit,pas(idiag),d_qt,d_qw,d_ql,d_qs
@@ -310 +312 @@
  9004   format('Dyn3d. Total Energy Budget (W/m2) ',A15,1i6,10(F8.2))
   END IF
-  !
+
   ! Store the new atmospheric state in "idiag"
-  !
+
   pas(idiag)=pas(idiag)+1
   h_vcol_pre(idiag)  = h_vcol_tot
@@ -323 +325 @@
   qs_pre(idiag)     = qs_tot
   ec_pre (idiag)    = ec_tot
-  !
+
   ELSE
     WRITE(lunout,*)'diagedyn: set to function with Earth parameters'