Changeset 1992 for LMDZ5/trunk/libf/phylmd/ustarhb.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (11 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/ustarhb.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/ustarhb.F) (1 diff)

LMDZ5/trunk/libf/phylmd/ustarhb.F90

@@ -1 @@
-!
+
 ! $Header$
-!
-      SUBROUTINE ustarhb(knon,u,v,cd_m, ustar)
-      use dimphy
-      IMPLICIT none
-c======================================================================
-c Laurent Li (LMD/CNRS), le 30 septembre 1998
-c Couche limite non-locale. Adaptation du code du CCM3.
-c Code non teste, donc a ne pas utiliser.
-c======================================================================
-c Nonlocal scheme that determines eddy diffusivities based on a
-c diagnosed boundary layer height and a turbulent velocity scale.
-c Also countergradient effects for heat and moisture are included.
-c
-c For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
-c Local versus nonlocal boundary-layer diffusion in a global climate
-c model. J. of Climate, vol. 6, 1825-1842.
-c======================================================================
-cym#include "dimensions.h"
-cym#include "dimphy.h"
-#include "YOMCST.h"
-c
-c Arguments:
-c
-      INTEGER knon ! nombre de points a calculer
-      REAL u(klon,klev) ! vitesse U (m/s)
-      REAL v(klon,klev) ! vitesse V (m/s)
-      REAL cd_m(klon) ! coefficient de friction au sol pour vitesse
-      REAL ustar(klon)
-c
-      INTEGER i, k
-      REAL zxt, zxq, zxu, zxv, zxmod, taux, tauy
-      REAL zx_alf1, zx_alf2 ! parametres pour extrapolation
-      LOGICAL unssrf(klon)  ! unstb pbl w/lvls within srf pbl lyr
-      LOGICAL unsout(klon)  ! unstb pbl w/lvls in outer pbl lyr
-      LOGICAL check(klon)   ! True=>chk if Richardson no.>critcal
-c
-#include "YOETHF.h"
-#include "FCTTRE.h"
-      DO i = 1, knon
-        zx_alf1 = 1.0
-        zx_alf2 = 1.0 - zx_alf1
-        zxu = u(i,1)*zx_alf1+u(i,2)*zx_alf2
-        zxv = v(i,1)*zx_alf1+v(i,2)*zx_alf2
-        zxmod = 1.0+SQRT(zxu**2+zxv**2)
-        taux = zxu *zxmod*cd_m(i)
-        tauy = zxv *zxmod*cd_m(i)
-        ustar(i) = SQRT(taux**2+tauy**2)
-c       print*,'Ust ',zxu,zxmod,taux,ustar(i)
-      ENDDO
-c
-      return
-      end
+
+SUBROUTINE ustarhb(knon, u, v, cd_m, ustar)
+  USE dimphy
+  IMPLICIT NONE
+  ! ======================================================================
+  ! Laurent Li (LMD/CNRS), le 30 septembre 1998
+  ! Couche limite non-locale. Adaptation du code du CCM3.
+  ! Code non teste, donc a ne pas utiliser.
+  ! ======================================================================
+  ! Nonlocal scheme that determines eddy diffusivities based on a
+  ! diagnosed boundary layer height and a turbulent velocity scale.
+  ! Also countergradient effects for heat and moisture are included.
+
+  ! For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
+  ! Local versus nonlocal boundary-layer diffusion in a global climate
+  ! model. J. of Climate, vol. 6, 1825-1842.
+  ! ======================================================================
+  ! ym#include "dimensions.h"
+  ! ym#include "dimphy.h"
+  include "YOMCST.h"
+
+  ! Arguments:
+
+  INTEGER knon ! nombre de points a calculer
+  REAL u(klon, klev) ! vitesse U (m/s)
+  REAL v(klon, klev) ! vitesse V (m/s)
+  REAL cd_m(klon) ! coefficient de friction au sol pour vitesse
+  REAL ustar(klon)
+
+  INTEGER i, k
+  REAL zxt, zxq, zxu, zxv, zxmod, taux, tauy
+  REAL zx_alf1, zx_alf2 ! parametres pour extrapolation
+  LOGICAL unssrf(klon) ! unstb pbl w/lvls within srf pbl lyr
+  LOGICAL unsout(klon) ! unstb pbl w/lvls in outer pbl lyr
+  LOGICAL check(klon) ! True=>chk if Richardson no.>critcal
+
+  include "YOETHF.h"
+  include "FCTTRE.h"
+  DO i = 1, knon
+    zx_alf1 = 1.0
+    zx_alf2 = 1.0 - zx_alf1
+    zxu = u(i, 1)*zx_alf1 + u(i, 2)*zx_alf2
+    zxv = v(i, 1)*zx_alf1 + v(i, 2)*zx_alf2
+    zxmod = 1.0 + sqrt(zxu**2+zxv**2)
+    taux = zxu*zxmod*cd_m(i)
+    tauy = zxv*zxmod*cd_m(i)
+    ustar(i) = sqrt(taux**2+tauy**2)
+    ! print*,'Ust ',zxu,zxmod,taux,ustar(i)
+  END DO
+
+  RETURN
+END SUBROUTINE ustarhb