/[lmdze]/trunk/dyn3d/calfis.f

Diff of /trunk/dyn3d/calfis.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/dyn3d/calfis.f90
revision 70 by guez,
Mon Jun 24 15:39:52 2013 UTC
+trunk/dyn3d/calfis.f
revision 90 by guez,
Wed Mar 12 21:16:36 2014 UTC
 Line 4 
 module calfis_m
  contains
-   SUBROUTINE calfis(rdayvrai, time, ucov, vcov, teta, q, masse, ps, pk, phis, &
+   SUBROUTINE calfis(rdayvrai, time, ucov, vcov, teta, q, ps, pk, phis, phi, &
-        phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, lafin)
+        w, dufi, dvfi, dtetafi, dqfi, dpfi, lafin)
-     ! From dyn3d/calfis.F, version 1.3 2005/05/25 13:10:09
+     ! From dyn3d/calfis.F, version 1.3, 2005/05/25 13:10:09
      ! Authors: P. Le Van, F. Hourdin
-     ! 1. Réarrangement des tableaux et transformation des variables
+     ! 1. R\'earrangement des tableaux et transformation des variables
      ! dynamiques en variables physiques
      ! 2. Calcul des termes physiques
 Line 18 
 contains
      ! Remarques:
-     ! - Les vents sont donnés dans la physique par leurs composantes
+     ! - Les vents sont donn\'es dans la physique par leurs composantes
      ! naturelles.
      ! - La variable thermodynamique de la physique est une variable
      ! intensive : T.
      ! Pour la dynamique on prend T * (preff / p(l))**kappa
-     ! - Les deux seules variables dépendant de la géométrie
+     ! - Les deux seules variables d\'ependant de la g\'eom\'etrie
-     ! nécessaires pour la physique sont la latitude pour le
+     ! n\'ecessaires pour la physique sont la latitude (pour le
-     ! rayonnement et l'aire de la maille quand on veut intégrer une
+     ! rayonnement) et l'aire de la maille (quand on veut int\'egrer une
-     ! grandeur horizontalement.
+     ! grandeur horizontalement).
      use comconst, only: kappa, cpp, dtphys, g
      use comgeom, only: apoln, cu_2d, cv_2d, unsaire_2d, apols, rlonu, rlonv
 Line 41 
 contains
      use physiq_m, only: physiq
      use pressure_var, only: p3d, pls
-     ! Arguments :
+     REAL, intent(in):: rdayvrai
+     REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour
-     ! Input :
+     REAL, intent(in):: ucov(iim + 1, jjm + 1, llm)
      ! ucov covariant zonal velocity
-     ! vcov covariant meridional velocity
-     ! teta potential temperature
-     ! ps surface pressure
-     ! masse masse d'air dans chaque maille
-     ! pts surface temperature (K)
-     ! callrad clef d'appel au rayonnement
-     ! Output :
-     ! dufi tendency for the natural zonal velocity (ms-1)
-     ! dvfi tendency for the natural meridional velocity
-     ! dtetafi tendency for the potential temperature
-     ! pdtsfi tendency for the surface temperature
-     ! pdtrad radiative tendencies \ input and output
+     REAL, intent(in):: vcov(iim + 1, jjm, llm)
-     ! pfluxrad radiative fluxes / input and output
+     ! vcov covariant meridional velocity
-     REAL, intent(in):: rdayvrai
+     REAL, intent(in):: teta(iim + 1, jjm + 1, llm) ! teta potential temperature
-     REAL, intent(in):: time ! heure de la journée en fraction de jour
-     REAL, intent(in):: ucov(iim + 1, jjm + 1, llm)
-     REAL vcov(iim + 1, jjm, llm)
-     REAL, intent(in):: teta(iim + 1, jjm + 1, llm)
      REAL, intent(in):: q(iim + 1, jjm + 1, llm, nqmx)
-     ! (mass fractions of advected fields)
+     ! mass fractions of advected fields
+     REAL, intent(in):: ps(iim + 1, jjm + 1) ! ps surface pressure
-     REAL masse(iim + 1, jjm + 1, llm)
-     REAL ps(iim + 1, jjm + 1)
      REAL, intent(in):: pk(iim + 1, jjm + 1, llm)
+     ! Exner = cp * (p / preff)**kappa
      REAL, intent(in):: phis(iim + 1, jjm + 1)
      REAL, intent(in):: phi(iim + 1, jjm + 1, llm)
-     REAL dudyn(iim + 1, jjm + 1, llm)
-     REAL dv(iim + 1, jjm, llm)
-     REAL dq(iim + 1, jjm + 1, llm, nqmx)
      REAL, intent(in):: w(iim + 1, jjm + 1, llm)
-     REAL dufi(iim + 1, jjm + 1, llm)
-     REAL dvfi(iim + 1, jjm, llm)
+     REAL, intent(out):: dufi(iim + 1, jjm + 1, llm)
+     ! tendency for the covariant zonal velocity (m2 s-2)
+     REAL, intent(out):: dvfi(iim + 1, jjm, llm)
+     ! tendency for the natural meridional velocity
      REAL, intent(out):: dtetafi(iim + 1, jjm + 1, llm)
-     REAL dqfi(iim + 1, jjm + 1, llm, nqmx)
+     ! tendency for the potential temperature
-     REAL dpfi(iim + 1, jjm + 1)
+     REAL, intent(out):: dqfi(iim + 1, jjm + 1, llm, nqmx)
+     REAL, intent(out):: dpfi(iim + 1, jjm + 1) ! tendance sur la pression
      LOGICAL, intent(in):: lafin
-     ! Local variables :
+     ! Local:
-     INTEGER i, j, l, ig0, ig, iq, iiq
+     INTEGER i, j, l, ig0, iq, iiq
      REAL zpsrf(klon)
      REAL paprs(klon, llm+1), play(klon, llm)
+     ! paprs defini aux (llm +1) interfaces des couches
+     ! play defini aux (llm) milieux des couches
      REAL pphi(klon, llm), pphis(klon)
      REAL u(klon, llm), v(klon, llm)
-Line 99 
 contains
+Line 94 
 contains
      real qx(klon, llm, nqmx) ! mass fractions of advected fields
      REAL omega(klon, llm)
-     REAL d_u(klon, llm), d_v(klon, llm)
+     REAL d_u(klon, llm), d_v(klon, llm) ! tendances physiques du vent (m s-2)
      REAL d_t(klon, llm), d_qx(klon, llm, nqmx)
      REAL d_ps(klon)
      REAL z1(iim)
      REAL pksurcp(iim + 1, jjm + 1)
-     ! I. Musat: diagnostic PVteta, Amip2
-     INTEGER, PARAMETER:: ntetaSTD=3
-     REAL:: rtetaSTD(ntetaSTD) = (/350., 380., 405./)
-     REAL PVteta(klon, ntetaSTD)
      !-----------------------------------------------------------------------
      !!print *, "Call sequence information: calfis"
-     ! 1. Initialisations :
-     ! latitude, longitude et aires des mailles pour la physique:
      ! 40. transformation des variables dynamiques en variables physiques:
-     ! 41. pressions au sol (en Pascals)
-     zpsrf(1) = ps(1, 1)
-     ig0 = 2
-     DO j = 2, jjm
-        CALL SCOPY(iim, ps(1, j), 1, zpsrf(ig0), 1)
-        ig0 = ig0+iim
-     ENDDO
-     zpsrf(klon) = ps(1, jjm + 1)
      ! 42. pression intercouches :
-     ! paprs defini aux (llm +1) interfaces des couches
-     ! play defini aux (llm) milieux des couches
-     ! Exner = cp * (p(l) / preff) ** kappa
      forall (l = 1: llm+1) paprs(:, l) = pack(p3d(:, :, l), dyn_phy)
      ! 43. temperature naturelle (en K) et pressions milieux couches
-Line 187 
 contains
+Line 158 
 contains
      DO l=1, llm
         DO j=2, jjm
            ig0 = 1+(j-2)*iim
-           u(ig0+1, l)= 0.5 * &
+           u(ig0+1, l)= 0.5 &
-                (ucov(iim, j, l)/cu_2d(iim, j) + ucov(1, j, l)/cu_2d(1, j))
+                * (ucov(iim, j, l) / cu_2d(iim, j) + ucov(1, j, l) / cu_2d(1, j))
            DO i=2, iim
-              u(ig0+i, l)= 0.5 * &
+              u(ig0+i, l)= 0.5 * (ucov(i-1, j, l)/cu_2d(i-1, j) &
-                   (ucov(i-1, j, l)/cu_2d(i-1, j) &
                    + ucov(i, j, l)/cu_2d(i, j))
            end DO
         end DO
-Line 204 
 contains
+Line 174 
 contains
           + vcov(:iim, j, l) / cv_2d(:iim, j))
      zvfi(iim + 1, 2:jjm, :) = zvfi(1, 2:jjm, :)
-     ! 47. champs de vents au pôle nord
+     ! 47. champs de vents au p\^ole nord
      ! U = 1 / pi * integrale [ v * cos(long) * d long ]
      ! V = 1 / pi * integrale [ v * sin(long) * d long ]
-Line 218 
 contains
+Line 188 
 contains
         zvfi(:, 1, l) = SUM(SIN(rlonv(:iim)) * z1) / pi
      ENDDO
-     ! 48. champs de vents au pôle sud:
+     ! 48. champs de vents au p\^ole sud:
      ! U = 1 / pi * integrale [ v * cos(long) * d long ]
      ! V = 1 / pi * integrale [ v * sin(long) * d long ]
-Line 235 
 contains
+Line 205 
 contains
      forall(l= 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy)
-     !IM calcul PV a teta=350, 380, 405K
-     CALL PVtheta(klon, llm, ucov, vcov, teta, t, play, paprs, &
-          ntetaSTD, rtetaSTD, PVteta)
      ! Appel de la physique :
      CALL physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, u, &
-          v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn, PVteta)
+          v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps)
      ! transformation des tendances physiques en tendances dynamiques:
-     ! tendance sur la pression :
      dpfi = gr_fi_dyn(d_ps)
      ! 62. enthalpie potentielle
-Line 254 
 contains
+Line 218 
 contains
         dtetafi(:, :, l) = cpp * gr_fi_dyn(d_t(:, l)) / pk(:, :, l)
      end do
-     ! 62. humidite specifique
-     DO iq=1, nqmx
-        DO l=1, llm
-           DO i=1, iim + 1
-              dqfi(i, 1, l, iq) = d_qx(1, l, iq)
-              dqfi(i, jjm + 1, l, iq) = d_qx(klon, l, iq)
-           ENDDO
-           DO j=2, jjm
-              ig0=1+(j-2)*iim
-              DO i=1, iim
-                 dqfi(i, j, l, iq) = d_qx(ig0+i, l, iq)
-              ENDDO
-              dqfi(iim + 1, j, l, iq) = dqfi(1, j, l, iq)
-           ENDDO
-        ENDDO
-     ENDDO
      ! 63. traceurs
      ! initialisation des tendances
-Line 297 
 contains
+Line 243 
 contains
      ! 65. champ u:
      DO l=1, llm
         DO i=1, iim + 1
            dufi(i, 1, l) = 0.
            dufi(i, jjm + 1, l) = 0.
-Line 306 
 contains
+Line 251 
 contains
         DO j=2, jjm
            ig0=1+(j-2)*iim
            DO i=1, iim-1
-              dufi(i, j, l)= &
+              dufi(i, j, l)= 0.5*(d_u(ig0+i, l)+d_u(ig0+i+1, l))*cu_2d(i, j)
-.5*(d_u(ig0+i, l)+d_u(ig0+i+1, l))*cu_2d(i, j)
            ENDDO
-           dufi(iim, j, l)= &
+           dufi(iim, j, l)= 0.5*(d_u(ig0+1, l)+d_u(ig0+iim, l))*cu_2d(iim, j)
-.5*(d_u(ig0+1, l)+d_u(ig0+iim, l))*cu_2d(iim, j)
            dufi(iim + 1, j, l)=dufi(1, j, l)
         ENDDO
      ENDDO
      ! 67. champ v:
      DO l=1, llm
         DO j=2, jjm-1
            ig0=1+(j-2)*iim
            DO i=1, iim
-              dvfi(i, j, l)= &
+              dvfi(i, j, l)= 0.5*(d_v(ig0+i, l)+d_v(ig0+i+iim, l))*cv_2d(i, j)
-.5*(d_v(ig0+i, l)+d_v(ig0+i+iim, l))*cv_2d(i, j)
            ENDDO
            dvfi(iim + 1, j, l) = dvfi(1, j, l)
         ENDDO
      ENDDO
-     ! 68. champ v pres des poles:
+     ! 68. champ v pr\`es des p\^oles:
      ! v = U * cos(long) + V * SIN(long)
      DO l=1, llm
         DO i=1, iim
-           dvfi(i, 1, l)= &
+           dvfi(i, 1, l)= d_u(1, l)*COS(rlonv(i))+d_v(1, l)*SIN(rlonv(i))
-                d_u(1, l)*COS(rlonv(i))+d_v(1, l)*SIN(rlonv(i))
+           dvfi(i, jjm, l)=d_u(klon, l)*COS(rlonv(i)) +d_v(klon, l)*SIN(rlonv(i))
-           dvfi(i, jjm, l)=d_u(klon, l)*COS(rlonv(i)) &
+           dvfi(i, 1, l)= 0.5*(dvfi(i, 1, l)+d_v(i+1, l))*cv_2d(i, 1)
-                +d_v(klon, l)*SIN(rlonv(i))
+           dvfi(i, jjm, l)= 0.5 &
-           dvfi(i, 1, l)= &
+                * (dvfi(i, jjm, l) + d_v(klon - iim - 1 + i, l)) * cv_2d(i, jjm)
-.5*(dvfi(i, 1, l)+d_v(i+1, l))*cv_2d(i, 1)
-           dvfi(i, jjm, l)= &
-.5*(dvfi(i, jjm, l)+d_v(klon-iim-1+i, l))*cv_2d(i, jjm)
         ENDDO
         dvfi(iim + 1, 1, l) = dvfi(1, 1, l)

 Legend:



Removed from v.70
 


changed lines


 
Added in v.90
 Legend:



Removed from v.70
 


changed lines


 
Added in v.90
-Removed from v.70
+Added in v.90

	ViewVC Help
Powered by ViewVC 1.1.21