Changeset 2197 for LMDZ5/trunk/libf/phylmd

Timestamp:

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

Author:

Ehouarn Millour

Message:

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• clift.F90 (modified) (2 diffs)
• convect3.F90 (modified) (2 diffs)
• cv3_inicp.F90 (modified) (2 diffs)
• hines_gwd.F90 (modified) (13 diffs)
• ini_wake.F90 (modified) (2 diffs)
• tilft43.F90 (modified) (1 diff)
• tlift.F90 (modified) (2 diffs)
• wake.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/clift.F90

@@ -3 +3 @@
 
 SUBROUTINE clift(p, t, rr, rs, plcl, dplcldt, dplcldq)
+IMPLICIT NONE
   ! ***************************************************************
   ! *                                                             *
@@ -41 +42 @@
 
   include "YOMCST.h"
+  real :: p,t,rr,rs,plcl,dplcldt,dplcldq,cpd,cpv,cl,cpvmcl,eps,alv0,a,b
+  real :: rh,chi,alv
 
   cpd = rcpd

LMDZ5/trunk/libf/phylmd/convect3.F90

@@ -17 +17 @@
   USE dimphy
   USE infotrac, ONLY: nbtr
-
+  IMPLICIT NONE
   include "dimensions.h"
   INTEGER na
@@ -73 +73 @@
 
 
-
+  REAL :: cpv,cl,cpvmcl,eps,alv0,rdcp,pbcrit,ptcrit,sigd,spfac
+  REAL :: tau,beta,alpha,dtcrit,dtovsh,ahm,rm,um,vm,dphinv
+  REAL :: a2,x,tvx,tvy,plcl,pden,dpbase,tvpbase,tvbase,tdif
+  REAL :: ath1,ath,delti,deltap,dcape,dlnp,sigold,dtmin,fac,w
+  REAL :: amu,rti,cpd,bf2,anum,denom,dei,altem,cwat,stemp,qp
+  REAL :: scrit,alt,smax,asij,wgh,sjmax,sjmin,smid,delp,delm
+  REAL :: asum,bsum,csum,wflux,tinv,wdtrain,awat,afac,afac1,afac2
+  REAL :: bfac,pr1,pr2,sigt,b6,c6,revap,tevap,delth,amfac,amp2
+  REAL :: xf,tf,af,bf,fac2,ur,sru,d,ampmax,dpinv,am,amde,cpinv
+  REAL :: amp1,ad,rat,ax,bx,cx,dx,ex,dsum
+  INTEGER :: nk,i,j,nopt,jn,k,im,jm,n
 
   REAL dnwd0(nd) !  precipitation driven unsaturated downdraft flux

LMDZ5/trunk/libf/phylmd/cv3_inicp.F90

@@ -9 +9 @@
   ! modified by :                                               *
   ! **************************************************************
-
+  IMPLICIT NONE
   include "YOMCST2.h"
 
@@ -19 +19 @@
 
   REAL qcoef1, qcoef2, qff, qfff, qmix, rmix, qmix1, rmix1, qmix2, rmix2, f
+  REAL :: sumcoef,sigma,aire,pdf,mu,df,ff
 
   qcoef1(f) = tanh(f/gammas)

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.

LMDZ5/trunk/libf/phylmd/ini_wake.F90

@@ -4 +4 @@
 SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, &
     alp_bl_prescr, ale_bl_prescr)
+  IMPLICIT NONE
   ! **************************************************************
   ! *
@@ -39 +40 @@
   include 'iniprint.h'
   ! declarations
+  REAL wape, fip, wape_prescr, fip_prescr
+  INTEGER it_wape_prescr
   REAL ale_bl_prescr
   REAL alp_bl_prescr
   REAL it
+  REAL w,f,alebl,alpbl
 
   ! FH A mettre si besoin dans physiq.def

LMDZ5/trunk/libf/phylmd/tilft43.F90

@@ -3 +3 @@
 
 SUBROUTINE tlift43(p, t, q, qs, gz, icb, nk, tvp, tpk, clw, nd, nl, kk)
+  IMPLICIT NONE
   REAL gz(nd), tpk(nd), clw(nd), p(nd)
   REAL t(nd), q(nd), qs(nd), tvp(nd), lv0
-
+  INTEGER icb, nk, nd, nl, kk
+  REAL cpd, cpv,  cl, g, rowl, gravity, cpvmcl, eps, epsi
+  REAL ah0, cpp, cpinv, tg, qg, alv, s, ahg, tc, denom, es
+  INTEGER i, nst, nsb, j 
   ! ***   ASSIGN VALUES OF THERMODYNAMIC CONSTANTS     ***
 

LMDZ5/trunk/libf/phylmd/tlift.F90

@@ -4 +4 @@
 SUBROUTINE tlift(p, t, rr, rs, gz, plcl, icb, nk, tvp, tpk, clw, nd, nl, &
     dtvpdt1, dtvpdq1)
-
+  IMPLICIT NONE
   ! Argument NK ajoute (jyg) = Niveau de depart de la
   ! convection
-
-  PARAMETER (na=60)
-  REAL gz(nd), tpk(nd), clw(nd)
+  INTEGER icb, nk, nd, nl
+  INTEGER,PARAMETER :: na=60
+  REAL gz(nd), tpk(nd), clw(nd), plcl
   REAL t(nd), rr(nd), rs(nd), tvp(nd), p(nd)
   REAL dtvpdt1(nd), dtvpdq1(nd) ! Derivatives of parcel virtual
@@ -17 +17 @@
   REAL dtpdt1(na), dtpdq1(na) ! Derivatives of parcel temperature
   ! wrt T1 and Q1
-
+  REAL gravity, cpd, cpv, cl, ci, cpvmcl, clmci, eps, alv0, alf0
+  REAL cpp, cpinv, ah0, alf, tg, s, ahg, tc, denom, alv, es, esi
+  REAL qsat_new, snew
+  INTEGER icbl, i, imin, j, icb1
 
   LOGICAL ice_conv

LMDZ5/trunk/libf/phylmd/wake.F90

@@ -1756 +1756 @@
   ! a une humidite positive dans la zone (x) et dans la zone (w).
   ! ------------------------------------------------------
-
+  IMPLICIT NONE
 
   ! Input
@@ -1772 +1772 @@
   REAL epsilon
   ! DATA epsilon/1.e-15/
+  INTEGER i,k
 
   DO k = 1, nl