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

Timestamp:

Nov 23, 2022, 11:11:26 AM (2 years ago)

Author:

romain.vande

Message:

Mars GCM:
The variable co2ice is deleted. All the co2 ice on surface is now in qsurf(:,igcm_co2).
CO2 tracer is now mandatory. diagfi output is unchanged.
RV

File:

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

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

@@ -5 +5 @@
       CONTAINS
       
-      SUBROUTINE vdifc(ngrid,nlay,nq,co2ice,ppopsk,
+      SUBROUTINE vdifc(ngrid,nlay,nq,ppopsk,
      $                ptimestep,pcapcal,lecrit,
      $                pplay,pplev,pzlay,pzlev,pz0,
@@ -17 +17 @@
       use tracer_mod, only: noms, igcm_dust_mass, igcm_dust_number,
      &                      igcm_dust_submicron, igcm_h2o_vap,
-     &                      igcm_h2o_ice, alpha_lift,
+     &                      igcm_h2o_ice, alpha_lift, igcm_co2,
      &                      igcm_hdo_vap, igcm_hdo_ice,
      &                      igcm_stormdust_mass, igcm_stormdust_number
@@ -74 +74 @@
 
 c    Argument added for condensation:
-      REAL,INTENT(IN) :: co2ice (ngrid), ppopsk(ngrid,nlay)
+      REAL,INTENT(IN) :: ppopsk(ngrid,nlay)
       logical,INTENT(IN) :: lecrit
       REAL,INTENT(IN) :: pcondicea_co2microp(ngrid,nlay)! tendency due to CO2 condensation (kg/kg.s-1)
@@ -721 +721 @@
            if (doubleq.AND.submicron) then
              do ig=1,ngrid
-c              if(co2ice(ig).lt.1) then
+c              if(qsurf(ig,igcm_co2).lt.1) then
                  pdqsdif(ig,igcm_dust_mass) =
      &             -alpha_lift(igcm_dust_mass)  
@@ -734 +734 @@
                                           !or 2 (injection in CL)
               do ig=1,ngrid
-               if(co2ice(ig).lt.1) then ! pas de soulevement si glace CO2
+               if(pqsurf(ig,igcm_co2).lt.1) then ! pas de soulevement si glace CO2
                  pdqsdif(ig,igcm_dust_mass) =
      &             -alpha_lift(igcm_dust_mass)  
@@ -743 +743 @@
              elseif(dustinjection.eq.1)then ! dust injection scheme = 1 injection from surface
               do ig=1,ngrid
-               if(co2ice(ig).lt.1) then ! pas de soulevement si glace CO2
+               if(pqsurf(ig,igcm_co2).lt.1) then ! pas de soulevement si glace CO2
                 IF((ti_injection_sol.LE.local_time(ig)).and.
      &                  (local_time(ig).LE.tf_injection_sol)) THEN
@@ -775 +775 @@
                 ENDIF
                
-                end if ! of if(co2ice(ig).lt.1)
+                end if ! of if(qsurf(ig,igcm_co2).lt.1)
               end do
              endif ! end if dustinjection
@@ -785 +785 @@
            else
 #endif
-            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zcdh,co2ice,
-     &                    pdqsdif)
+            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zcdh,
+     &                    pqsurf(:,igcm_co2),pdqsdif)
 #ifndef MESOSCALE
            endif !doubleq.AND.submicron