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

Timestamp:

Aug 20, 2025, 4:25:12 PM (4 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/swr_fouquart.F (modified) (4 diffs)

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

@@ -15 +15 @@
       use callkeys_mod, only: rayleigh
       use swrayleigh_mod, only: swrayleigh
+      use cvmgt_mod, only: cvmgt
 
       IMPLICIT NONE
@@ -112 +113 @@
      S     ,  ZTRA1(NDLON,NFLEV+1), ZTRA2(NDLON,NFLEV+1)
 
-c     Function
-c     --------
-      real CVMGT
-
 C    --------------------------------
 C     OPTICAL PARAMETERS FOR AEROSOLS 
@@ -138 +135 @@
             ZTAUAZ(JL,JK) = 0.
          END DO
-         DO 106 JAE=1,naerkind
-            DO 105 JL = 1 , KDLON
+         DO JAE=1,naerkind
+            DO JL = 1 , KDLON
 c              Mean Extinction optical depth in the spectral band
 c              ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -154 +151 @@
      S           QVISsQREF3d(JL,JK,KNU,JAE)*
      &           omegaVIS3d(JL,JK,KNU,JAE)*gVIS3d(JL,JK,KNU,JAE)
- 105        CONTINUE
- 106     CONTINUE
+            ENDDO
+         ENDDO
       END DO
 C     
       DO JK = 1 , nlaylte
          DO JL = 1 , KDLON
+           ! NB: function CVMGT(x1,x2,l) returns x1 if l==.true.
+           !     or x2 otherwise. Maybe worth inlining someday? 
             ZCGAZ(JL,JK) = CVMGT( 0., ZCGAZ(JL,JK) / ZPIZAZ(JL,JK),
      S            (ZPIZAZ(JL,JK).EQ.0) )