Changeset 2197 for LMDZ5/trunk/libf/phylmd/hines_gwd.F90

Timestamp:

Feb 9, 2015, 8:13:05 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Added 'implicit none' statements and proper variable definitions where they were missing.
EM

File:

• LMDZ5/trunk/libf/phylmd/hines_gwd.F90 (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/hines_gwd.F90

@@ -625 +625 @@
     mmin_alpha, kstar, slope, f1, f2, f3, naz, levbot, levtop, il1, il2, &
     nlons, nlevs, nazmth, sigsqmcw, sigmatm, lorms, sigalpmc, f2mod)
-
   ! Smooth cutoff wavenumbers and total rms velocity in the vertical
   ! direction NSMAX times, using FLUX_U as temporary work array.
@@ -715 +714 @@
     i_alpha, mmin_alpha, kstar, slope, f1, f2, f3, naz, levbot, levtop, il1, &
     il2, nlons, nlevs, nazmth, sigsqmcw, sigmatm, lorms, sigalpmc, f2mod)
-
+  IMPLICIT NONE
   ! This routine calculates the cutoff vertical wavenumber and velocity
   ! variances on a longitude by altitude grid for the Hines' Doppler
@@ -766 +765 @@
 
   INTEGER naz, levbot, levtop, il1, il2, nlons, nlevs, nazmth
-  REAL slope, kstar(nlons), f1, f2, f3
+  REAL slope, kstar(nlons), f1, f2, f3, f2mfac
   REAL m_alpha(nlons, nlevs, nazmth)
   REAL sigma_alpha(nlons, nlevs, nazmth)
@@ -938 +937 @@
 SUBROUTINE hines_wind(v_alpha, vel_u, vel_v, naz, il1, il2, lev1, lev2, &
     nlons, nlevs, nazmth)
-
+    IMPLICIT NONE
   ! This routine calculates the azimuthal horizontal background wind
   ! components
@@ -1034 +1033 @@
     m_alpha, ak_alpha, k_alpha, slope, naz, il1, il2, lev1, lev2, nlons, &
     nlevs, nazmth, lorms)
-
+    IMPLICIT NONE
   ! Calculate zonal and meridional components of the vertical flux
   ! of horizontal momentum and corresponding wave drag (force per unit mass)
@@ -1089 +1088 @@
   ! Internal variables.
 
-  INTEGER i, l, lev1p, lev2m
+  INTEGER i, l, lev1p, lev2m, lev2p
   REAL cos45, prod2, prod4, prod6, prod8, dendz, dendz2
   DATA cos45/0.7071068/
@@ -1234 +1233 @@
     bvfreq, density, densb, sigma_t, visc_mol, kstar, slope, f2, f3, f5, f6, &
     naz, il1, il2, lev1, lev2, nlons, nlevs, nazmth)
-
+  IMPLICIT NONE
   ! This routine calculates the gravity wave induced heating and
   ! diffusion coefficient on a longitude by altitude grid for
@@ -1355 +1354 @@
 SUBROUTINE hines_sigma(sigma_t, sigma_alpha, sigsqh_alpha, naz, lev, il1, &
     il2, nlons, nlevs, nazmth)
-
+  IMPLICIT NONE
   ! This routine calculates the total rms and azimuthal rms horizontal
   ! velocities at a given level on a longitude by altitude grid for
@@ -1450 +1449 @@
 SUBROUTINE hines_intgrl(i_alpha, v_alpha, m_alpha, bvfb, slope, naz, lev, &
     il1, il2, nlons, nlevs, nazmth, lorms)
-
+  IMPLICIT NONE
   ! This routine calculates the vertical wavenumber integral
   ! for a single vertical level at each azimuth on a longitude grid
@@ -1643 +1642 @@
     alt_cutoff, smco, nsmax, iheatcal, k_alpha, ierror, nmessg, nlons, &
     nazmth, coslat)
-
+  IMPLICIT NONE
   ! This routine specifies various parameters needed for the
   ! the Hines' Doppler spread gravity wave drag parameterization scheme.
@@ -1772 +1771 @@
     sigma_alpha, v_alpha, m_alpha, iu_print, iv_print, nmessg, ilprt1, &
     ilprt2, levprt1, levprt2, naz, nlons, nlevs, nazmth)
-
+  IMPLICIT NONE
   ! Print out altitude profiles of various quantities from
   ! Hines' Doppler spread gravity wave drag parameterization scheme.
@@ -1864 +1863 @@
 SUBROUTINE hines_exp(data, data_zmax, alt, alt_exp, iorder, il1, il2, lev1, &
     lev2, nlons, nlevs)
-
+  IMPLICIT NONE
   ! This routine exponentially damps a longitude by altitude array
   ! of data above a specified altitude.
@@ -1941 +1940 @@
 SUBROUTINE vert_smooth(data, work, coeff, nsmooth, il1, il2, lev1, lev2, &
     nlons, nlevs)
-
+  IMPLICIT NONE
   ! Smooth a longitude by altitude array in the vertical over a
   ! specified number of levels using a three point smoother.