Changeset 3901 for trunk/LMDZ.MARS/libf/phymars/vlz_fi.F

Timestamp:

Aug 20, 2025, 4:25:12 PM (3 days ago)

Author:

emillour

Message:

Mars PCM:
Some code tidying:

• turn aeroptproperties, albedocaps, cvmgp and convadj into modules
• remove useless check in callradite
• clean module vlz_fi (remove obsolete #ifdef CRAY cpp alternatives)
• remove function cvmgp since it was only called under #ifdef CRAY in vlz_fi

EM

File:

• trunk/LMDZ.MARS/libf/phymars/vlz_fi.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/vlz_fi.F

@@ -10 +10 @@
 c
 c    ********************************************************************
-c     Shema  d'advection " pseudo amont " dans la verticale
-c    pour appel dans la physique (sedimentation)
+c    "pseudo upstream" Advection scheme along the vertical
+c    to be used in the physics (sedimentation)
 c    ********************************************************************
-c    q rapport de melange (kg/kg)...
-c    masse : masse de la couche Dp/g
-c    w : masse d'atm ``transferee'' a chaque pas de temps (kg.m-2)
-c    pente_max = 2 conseillee
-c
-c
-c   --------------------------------------------------------------------
+
       IMPLICIT NONE
 c
@@ -29 +23 @@
       integer,intent(in) :: ngrid ! number of atmospheric columns
       integer,intent(in) :: nlay ! number of atmospheric layers
-      real masse(ngrid,nlay),pente_max
-      REAL q(ngrid,nlay)
-      REAL w(ngrid,nlay)
-      REAL wq(ngrid,nlay+1)
+      real,intent(in) :: masse(ngrid,nlay) ! mass of atmospheric layer delta(P)/g
+      real,intent(in) :: pente_max ! maximum slope for the scheme (2 is recommended)
+      real,intent(inout) :: q(ngrid,nlay) ! tracer mixing ratio (kg/kg)
+      real,intent(inout) :: w(ngrid,nlay) ! mass of atmosphere "transfered" over the time step (kg.m-2)
+      real,intent(out) :: wq(ngrid,nlay+1) ! trancer increment due to advection (kg)
 c
 c      Local 
@@ -45 +40 @@
       integer m
 
-      REAL      SSUM,CVMGP,CVMGT
-      integer ismax,ismin
 
-
-c    On oriente tout dans le sens de la pression c'est a dire dans le
-c    sens de W
+c    Orientation follows pressure, i.e. follows W
 
       do l=2,nlay
@@ -61 +52 @@
       do l=2,nlay-1
          do ij=1,ngrid
-#ifdef CRAY
-            dzq(ij,l)=0.5*
-     ,      cvmgp(dzqw(ij,l)+dzqw(ij,l+1),0.,dzqw(ij,l)*dzqw(ij,l+1))
-#else
             if(dzqw(ij,l)*dzqw(ij,l+1).gt.0.) then
                 dzq(ij,l)=0.5*(dzqw(ij,l)+dzqw(ij,l+1))
@@ -70 +57 @@
                 dzq(ij,l)=0.
             endif
-#endif
+
             dzqmax=pente_max*min(adzqw(ij,l),adzqw(ij,l+1))
             dzq(ij,l)=sign(min(abs(dzq(ij,l)),dzqmax),dzq(ij,l))
@@ -81 +68 @@
       enddo
 c ---------------------------------------------------------------
-c   .... calcul des termes d'advection verticale  .......
+c   .... compute vertical advection terms   .......
 c ---------------------------------------------------------------
 
-c calcul de  - d( q   * w )/ d(sigma)    qu'on ajoute a  dq pour calculer dq
+c compute  - d( q   * w )/ d(sigma) , later added to dq to compute dq
 c
 c      No flux at the model top: