Changeset 2274 for trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

Timestamp:

Apr 2, 2020, 5:45:44 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Some code cleanup:

• move profile.F to dyn1d since it is only used by the 1D model
• remove scatter.F, which is not used (and moreover a synonym of the scatter routine from the "mod_phys_lmdz_para" module)
• remove obsolete "multipl.F" routine, replace the calls to it in vdifc_mod by modern Fortran syntax.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (3 diffs)

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

@@ -444 +444 @@
 c       et /zkv/ = Ku
  
-      CALL multipl((nlay-1)*ngrid,zkv(1,2),zb0(1,2),zb(1,2))
-      CALL multipl(ngrid,zcdv,zb0,zb)
+      zb(1:ngrid,2:nlay)=zkv(1:ngrid,2:nlay)*zb0(1:ngrid,2:nlay)
+      zb(1:ngrid,1)=zcdv(1:ngrid)*zb0(1:ngrid,1)
 
       DO ig=1,ngrid
@@ -534 +534 @@
 
 c Mass variation scheme:
-      CALL multipl((nlay-1)*ngrid,zkh(1,2),zb0(1,2),zb(1,2))
-      CALL multipl(ngrid,zcdh,zb0,zb)
+      zb(1:ngrid,2:nlay)=zkh(1:ngrid,2:nlay)*zb0(1:ngrid,2:nlay)
+      zb(1:ngrid,1)=zcdh(1:ngrid)*zb0(1:ngrid,1)
 
 c on initialise dm c
@@ -784 +784 @@
 c       --------------------------------
         do iq=1,nq  !for all tracers including stormdust
-          CALL multipl((nlay-1)*ngrid,zkh(1,2),zb0(1,2),zb(1,2))
+          zb(1:ngrid,2:nlay)=zkh(1:ngrid,2:nlay)*zb0(1:ngrid,2:nlay)
 
           if ((water).and.(iq.eq.igcm_h2o_vap)) then 
 c            This line is required to account for turbulent transport 
 c            from surface (e.g. ice) to mid-layer of atmosphere:
-             CALL multipl(ngrid,zcdv,zb0,zb(1,1))
-             CALL multipl(ngrid,dryness,zb(1,1),zb(1,1)) 
+             zb(1:ngrid,1)=zcdv(1:ngrid)*zb0(1:ngrid,1)
+             zb(1:ngrid,1)=dryness(1:ngrid)*zb(1:ngrid,1)
           else ! (re)-initialize zb(:,1)
              zb(1:ngrid,1)=0