Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/coare_cp_mod.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

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

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

@@ -6 +6 @@
 contains
 
-  real function psit_30(zet)
+  REAL function psit_30(zet)
     IMPLICIT NONE
-    real, intent(in) :: zet
+    REAL, INTENT(IN) :: zet
 
     REAL :: x, psik, psic, f, c
@@ -26 +26 @@
   END FUNCTION psit_30
 
-  real function psiuo(zet)
+  REAL function psiuo(zet)
     IMPLICIT NONE
-    real, intent(in) :: zet
+    REAL, INTENT(IN) :: zet
 
     REAL :: x, psik, psic, f, c
 
-    if (zet<0) THEN
+    IF (zet<0) THEN
        x=(1.-15.*zet)**.25
        psik=2.*log((1.+x)/2.)+log((1.+x*x)/2.)-2.*atan(x)+2.*atan(1.)
@@ -42 +42 @@
      c=min(50.,.35*zet)
      psiuo=-((1.+1.0*zet)**1.0+.667*(zet-14.28)/exp(c)+8.525)
-  endif
+  ENDIF
 
 END FUNCTION psiuo
@@ -67 +67 @@
   IMPLICIT NONE
 
-  real, intent(in) :: du,dt,dq,t,q
-  real, intent(in) :: zu,zt,zq,p,zi
-  integer, intent(in) :: nits
-  ! real, dimension (nits), intent(out) :: zo,tau,hsb,hlb
-  ! real, dimension(nits,3), intent(out) :: var
-  real, dimension(3), intent(out) :: coeffs
-  real, intent(out) :: rugosm
-  real, intent(out) :: rugosh
-
-  real, parameter ::  beta=1.2, von=.4, fdg = 1. ,&
+  REAL, INTENT(IN) :: du,dt,dq,t,q
+  REAL, INTENT(IN) :: zu,zt,zq,p,zi
+  INTEGER, INTENT(IN) :: nits
+  ! REAL, dimension (nits), INTENT(OUT) :: zo,tau,hsb,hlb
+  ! REAL, DIMENSION(nits,3), INTENT(OUT) :: var
+  REAL, DIMENSION(3), INTENT(OUT) :: coeffs
+  REAL, INTENT(OUT) :: rugosm
+  REAL, INTENT(OUT) :: rugosh
+
+  REAL, parameter ::  beta=1.2, von=.4, fdg = 1. ,&
        tdk = 273.16, pi = 3.141593, grav = 9.82,&
        rgas = 287.1
 
-  integer, dimension(3) :: shape_input
-
-  real bf, cc, rhoa, visa,&
+  INTEGER, DIMENSION(3) :: shape_input
+
+  REAL bf, cc, rhoa, visa,&
        u10, ut, uts, ut0, ug,&
        cd10, ch10, ct, ct10,&
@@ -98 +98 @@
 
 
-  real old_usr, old_tsr, old_qsr,tmp
-
-  real, external :: grv
-  integer i,j,k
+  REAL old_usr, old_tsr, old_qsr,tmp
+
+  REAL, external :: grv
+  INTEGER i,j,k
 !---------------- Rajout pour prendre en compte différent Z0 --------------------------------!
 !  INTEGER :: NGRVWAVES        ! Pour le choix du z0
@@ -153 +153 @@
   Ribu=grav*zu/t*(dt+.61*t*dq)/ut**2
 
-  if (Ribu < 0) THEN
+  IF (Ribu < 0) THEN
      zetu=CC*Ribu/(1+Ribu/Ribcu)
   else
      zetu=CC*Ribu*(1+27/9*Ribu/CC)
-  endif
+  ENDIF
 
   L10=zu/zetu
@@ -163 +163 @@
   ! if (zetu .GT. 50) THEN
   !    nits=1
-  ! endif
+  ! END IF
 
   usr=ut*von/(log(zu/zom10)-psiuo(zetu))
@@ -171 +171 @@
 
   ! charnock constant - lin par morceau - constant
-  if ( ut <= 10. ) THEN
+  IF ( ut <= 10. ) THEN
      charn=0.011
   else
-     if (ut > 18) THEN
+     IF (ut > 18) THEN
         charn=0.018
      else
         charn=0.011+(ut-10)/(18-10)*(0.018-0.011)
      endif
-  endif
+  ENDIF
 
 !  ZHWAVE = 0.018*ut*ut*(1.+0.015*ut)
@@ -220 +220 @@
      Bf=-grav/t*usr*(tsr+.61*t*qsr)
 
-     if (Bf > 0) THEN
+     IF (Bf > 0) THEN
         ug=Beta*(Bf*zi)**.333
      else