Changeset 2601 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jan 4, 2022, 1:20:22 PM (3 years ago)

Author:

abierjon

Message:

Mars GCM:
Some further cleaning of co2condens following r2599 and r2600
that fixed the bug appearing when the scavenging by CO2 ice
was activated.

AB

File:

• trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (7 diffs)

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

@@ -132 +132 @@
       real :: emisref(ngrid)
       
-      REAL zdq_scav(ngrid,nlayer,nq) ! tendancy due to scavenging by co2
+      REAL zdq_scav(ngrid,nlayer,nq) ! tendency due to scavenging by co2
       REAL zq(ngrid,nlayer,nq)
 
@@ -508 +508 @@
              zv(l)   =pv(ig,l)   +pdv( ig,l)  *ptimestep
             do iq=1,nq
-             zqc(l,iq)=zq(ig,l,iq)+zdq_scav(ig,l,iq)*ptimestep ! zdq_scav=0 if watercloud=false
+             zqc(l,iq)=zq(ig,l,iq)+zdq_scav(ig,l,iq)*ptimestep ! zdq_scav=0 if co2clouds=true
             enddo
            enddo
@@ -678 +678 @@
      &             + condens_layer(ig,l)*zqc(l,iq) )
 
-                    pdqc(ig,l,iq)=pdqc(ig,l,iq)+zdq_scav(ig,l,iq) ! zdq_scav=0 if watercloud=false
+                    pdqc(ig,l,iq)=pdqc(ig,l,iq)+zdq_scav(ig,l,iq) ! zdq_scav=0 if co2clouds=true
                   END DO
                 end if
@@ -886 +886 @@
       REAL,INTENT(IN) :: pq(ngrid,nlayer,nq)
       REAL,INTENT(IN) :: rdust(ngrid,nlayer) ! dust effective radius
-      REAL,INTENT(IN) :: pcondicea(ngrid,nlayer) ! condensation rate in layer  l  (kg/m2/s)
+      REAL,INTENT(IN) :: pcondicea(ngrid,nlayer) ! condensation rate in layer l (kg/m2/s)
       REAL,INTENT(IN) :: pfallice(ngrid,nlayer+1) ! amount of ice falling from layer l (kg/m2/s)
       
-      REAL,INTENT(OUT) :: pdq_scav(ngrid,nlayer,nq) ! tendancy due to scavenging by co2
+      REAL,INTENT(OUT) :: pdq_scav(ngrid,nlayer,nq) ! tendency due to scavenging by co2
       REAL,INTENT(OUT) :: pdqsc(ngrid,nq) ! tendency on surface tracers
       
@@ -906 +906 @@
       scav_ratio_wice = scav_ratio_dust
       pdq_scav(:,:,:)=0.
-      
+      pdqsc(:,:)=0.
+             
       DO ig=1,ngrid
         scav_dust_mass(nlayer+1)=0.
@@ -983 +984 @@
        ENDDO
        ! Calculation of the surface tendencies
-       pdqsc(ig,igcm_dust_mass)=0.
-       pdqsc(ig,igcm_dust_number)=0.
-       
        if (freedust) then
          pdqsc(ig,igcm_dust_mass)=pdqsc(ig,igcm_dust_mass)
@@ -1001 +999 @@
          pdqsc(ig,igcm_h2o_ice)=scav_h2o_ice(1)
        endif
-      ENDDO
+       
+      ENDDO ! loop on ngrid
       
       END SUBROUTINE scavenging_by_co2