Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecumev6_flux.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/ecumev6_flux.F90 (modified) (9 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/ecumev6_flux.F90

@@ -253 +253 @@
 !----------------------- Modif Olive calcul de PRHOA ---------------------------
 
-!write(*,*) "PZ0SEA ",PZ0SEA
-!write(*,*) "PTA ",PTA
-!write(*,*) "PSST ",PSST
-!write(*,*) "PQA ",PQA
-!write(*,*) "PVMOD ",PVMOD
-!write(*,*) "PZREF ",PZREF
-!write(*,*) "PUREF ",PUREF
-!write(*,*) "PPS ",PPS
-!write(*,*) "PPA ",PPA
-!write(*,*) "OPRECIP ",OPRECIP
-!write(*,*) "PZ0HSEA ",PZ0HSEA
-!write(*,*) "PRAIN ",PRAIN
+!WRITE(*,*) "PZ0SEA ",PZ0SEA
+!WRITE(*,*) "PTA ",PTA
+!WRITE(*,*) "PSST ",PSST
+!WRITE(*,*) "PQA ",PQA
+!WRITE(*,*) "PVMOD ",PVMOD
+!WRITE(*,*) "PZREF ",PZREF
+!WRITE(*,*) "PUREF ",PUREF
+!WRITE(*,*) "PPS ",PPS
+!WRITE(*,*) "PPA ",PPA
+!WRITE(*,*) "OPRECIP ",OPRECIP
+!WRITE(*,*) "PZ0HSEA ",PZ0HSEA
+!WRITE(*,*) "PRAIN ",PRAIN
 
 
 PRHOA(:) = PPS(:) / (287.1 * PTA(:) * (1.+.61*PQA(:)))
-!write(*,*) "klon klon ",klon,PTA
-!write(*,*) "PRHOA ",SIZE(PRHOA),PRHOA
+!WRITE(*,*) "klon klon ",klon,PTA
+!WRITE(*,*) "PRHOA ",SIZE(PRHOA),PRHOA
 
 PREF = 100900.                    ! = 1000 hPa
@@ -358 +358 @@
 ZVMOD = MAX(PVMOD , 0.1 * MIN(10.,PUREF) )    !set a minimum value to wind
 
-write(*,*) "ZVMOD ",SIZE(ZVMOD)
+WRITE(*,*) "ZVMOD ",SIZE(ZVMOD)
 
 !       2.0. Radiative fluxes - For warm layer & cool skin
@@ -386 +386 @@
 ZWORK2A    = XRD/XRV
 
-!write(*,*) "ZFOES ",ZFOES
-!write(*,*) "PPS ",PPS
-!write(*,*) "ZWORK1 ",ZWORK1
-!write(*,*) "XRD ",XRD
-!write(*,*) "XRV ",XRV
-!write(*,*) "PPA ",PPA
-!write(*,*) "ZWORK1A ",ZWORK1A
-
-write(*,*) "PQSAT : ",PQSAT
-write(*,*) "PQSATA : ",PQSATA
+!WRITE(*,*) "ZFOES ",ZFOES
+!WRITE(*,*) "PPS ",PPS
+!WRITE(*,*) "ZWORK1 ",ZWORK1
+!WRITE(*,*) "XRD ",XRD
+!WRITE(*,*) "XRV ",XRV
+!WRITE(*,*) "PPA ",PPA
+!WRITE(*,*) "ZWORK1A ",ZWORK1A
+
+WRITE(*,*) "PQSAT : ",PQSAT
+WRITE(*,*) "PQSATA : ",PQSATA
 
 !*       2.    COMPUTE SATURATION HUMIDITY
@@ -410 +410 @@
 ZDQ(:) = PQA(:)-PQSATA(:)
 
-write(*,*) "PQA ",PQA(:)
-write(*,*) "PQSAT",PQSAT(:)
-write(*,*) "ZDQ",ZDQ(:)
+WRITE(*,*) "PQA ",PQA(:)
+WRITE(*,*) "PQSAT",PQSAT(:)
+WRITE(*,*) "ZDQ",ZDQ(:)
 
 !       2.3. Latent heat of vaporisation
@@ -421 +421 @@
 
 
-write(*,*) "ZLVA ",ZLVA
-write(*,*) "ZLVS ",ZLVS
+WRITE(*,*) "ZLVA ",ZLVA
+WRITE(*,*) "ZLVS ",ZLVS
 
 
@@ -450 +450 @@
 ZDDQ(:) = SIGN(MAX(ABS(ZDDQ(:)),10.0*ZDQSR0),ZDDQ(:))
 
-write(*,*) "ZDDU ",ZDDU
-write(*,*) "ZDDQ ",ZDDQ
-write(*,*) "ZDDT ",ZDDT
+WRITE(*,*) "ZDDU ",ZDDU
+WRITE(*,*) "ZDDQ ",ZDDQ
+WRITE(*,*) "ZDDT ",ZDDT
 
 JCV (:) = -1
@@ -624 +624 @@
   ZEF(JLON)  = -PRHOA(JLON)*ZLVS(JLON)*ZUSR(JLON)*ZQSR(JLON)
 
-  write(*,*) "ZTAU = ",ZTAU(JLON)
-  write(*,*) "SENS = ",ZHF(JLON)
-  write(*,*) "LAT = ",ZEF(JLON)
+  WRITE(*,*) "ZTAU = ",ZTAU(JLON)
+  WRITE(*,*) "SENS = ",ZHF(JLON)
+  WRITE(*,*) "LAT = ",ZEF(JLON)
 
 !       4.2. Exchange coefficients PCD, PCH, PCE
@@ -634 +634 @@
   PCE(JLON) = (ZUSR(JLON)*ZQSR(JLON))/(ZDDU(JLON)*ZDDQ(JLON))
 
-  write(*,*) "ZUSR = ",ZUSR(JLON)
-  write(*,*) "ZTSR = ",ZTSR(JLON)
-  write(*,*) "ZQSR = ",ZQSR(JLON)
+  WRITE(*,*) "ZUSR = ",ZUSR(JLON)
+  WRITE(*,*) "ZTSR = ",ZTSR(JLON)
+  WRITE(*,*) "ZQSR = ",ZQSR(JLON)
 
 !       4.3. Stochastic perturbation of turbulent fluxes
@@ -753 +753 @@
 !ENDIF
 
-write(*,*) "JLON ",JLON
-write(*,*) "PTA ",klon,PTA
-write(*,*) "PCD ",SIZE(PCD),PCD
-write(*,*) "PCQ ",SIZE(PCE),PCE
-write(*,*) "PCH ",SIZE(PCH),PCH
+WRITE(*,*) "JLON ",JLON
+WRITE(*,*) "PTA ",klon,PTA
+WRITE(*,*) "PCD ",SIZE(PCD),PCD
+WRITE(*,*) "PCQ ",SIZE(PCE),PCE
+WRITE(*,*) "PCH ",SIZE(PCH),PCH
 
 coeffs = [PCD,&