Changeset 1911 for trunk/LMDZ.MARS/libf/phymars/improvedCO2clouds.F

Timestamp:

Mar 9, 2018, 6:19:27 PM (7 years ago)

Author:

mvals

Message:

Mars GCM
"cosmetics" changes in co2cloud and improvedCO2clouds to better understand the role of the variables:

• in co2cloud : - change of variables subpdq and subpdt in sum_subpdq and sum_subpdt to make clear we ae making the sum of the tendencies in the microphysics loop.
  • change of variables pdqsed in subpdqsed as these are tendencies inside the microphysics loop
  • flag "sedimentation" has been added for the sedimentation block
  • variables names in the sedimentation block have been changed (tempo_traceurs becomes zqsed, sav_trac becomes zqsed0)
  • variable sum_subpdqs_sedco2 was added in the microphysics loop to make clear we make the sum of the surface sedimentation flux, output variable is still pdqs_sedco2
  • variable sum_subpdqs_sedco2 has been initialized to zero
  • zteff has been changed to pteff to be logical with the other variables names
• in improvedCO2clouds : change of the names of input and output variables according to their names in co2cloud in order to not confuse them (in improvedCO2clouds "ptimestep" corresponds actually to

"microtimestep").

MV

File:

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

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

@@ -1 @@
-      subroutine improvedCO2clouds(ngrid,nlay,ptimestep,
-     &             pplay,pplev,pt,pdt,
-     &             pq,pdq,pdqcloudco2,pdtcloudco2,
+      subroutine improvedCO2clouds(ngrid,nlay,microtimestep,
+     &             pplay,pplev,pteff,sum_subpdt,
+     &             pqeff,sum_subpdq,subpdqcloudco2,subpdtcloudco2,
      &             nq,tauscaling,
      &             memdMMccn,memdMMh2o,memdNNccn)
@@ -61 +61 @@
       INTEGER,INTENT(in) :: ngrid,nlay
       integer,intent(in) :: nq         ! number of tracers
-      REAL,INTENT(in) :: ptimestep     ! physics time step (s)
+      REAL,INTENT(in) :: microtimestep     ! physics time step (s)
       REAL,INTENT(in) :: pplay(ngrid,nlay)     ! mid-layer pressure (Pa)
       REAL,INTENT(in) :: pplev(ngrid,nlay+1)   ! inter-layer pressure (Pa)
-      REAL,INTENT(in) :: pt(ngrid,nlay) ! temperature at the middle of the
+      REAL,INTENT(in) :: pteff(ngrid,nlay) ! temperature at the middle of the
                                  !   layers (K)
-      REAL,INTENT(in) :: pdt(ngrid,nlay) ! tendency on temperature from
+      REAL,INTENT(in) :: sum_subpdt(ngrid,nlay) ! tendency on temperature from
                                  !  previous physical parametrizations
-      REAL,INTENT(in) :: pq(ngrid,nlay,nq) ! tracers (kg/kg)
-      REAL,INTENT(in) :: pdq(ngrid,nlay,nq) ! tendencies on tracers 
+      REAL,INTENT(in) :: pqeff(ngrid,nlay,nq) ! tracers (kg/kg)
+      REAL,INTENT(in) :: sum_subpdq(ngrid,nlay,nq) ! tendencies on tracers 
                                  !  before condensation (kg/kg.s-1)
       REAL,INTENT(in) :: tauscaling(ngrid) ! Convertion factor for qdust and Ndust
 c     Outputs:
-      REAL,INTENT(out) :: pdqcloudco2(ngrid,nlay,nq) ! tendency on tracers
+      REAL,INTENT(out) :: subpdqcloudco2(ngrid,nlay,nq) ! tendency on tracers
                                    ! due to CO2 condensation (kg/kg.s-1)
       ! condensation si igcm_co2_ice 
-      REAL,INTENT(out) :: pdtcloudco2(ngrid,nlay)  ! tendency on temperature due
+      REAL,INTENT(out) :: subpdtcloudco2(ngrid,nlay)  ! tendency on temperature due
                                    ! to latent heat
 
@@ -302 +302 @@
 
 c     Initialize the tendencies
-      pdqcloudco2(1:ngrid,1:nlay,1:nq)=0.
-      pdtcloudco2(1:ngrid,1:nlay)=0.
-      
-c pt temperature layer; pdt dT.s-1
-c pq traceur kg/kg; pdq tendance idem .s-1
+      subpdqcloudco2(1:ngrid,1:nlay,1:nq)=0.
+      subpdtcloudco2(1:ngrid,1:nlay)=0.
+      
+c pteff temperature layer; sum_subpdt dT.s-1
+c pqeff traceur kg/kg; sum_subpdq tendance idem .s-1
       zt(1:ngrid,1:nlay) = 
-     &      pt(1:ngrid,1:nlay) + 
-     &      pdt(1:ngrid,1:nlay) * ptimestep 
+     &      pteff(1:ngrid,1:nlay) + 
+     &      sum_subpdt(1:ngrid,1:nlay) * microtimestep 
       zq(1:ngrid,1:nlay,1:nq) = 
-     &      pq(1:ngrid,1:nlay,1:nq) + 
-     &      pdq(1:ngrid,1:nlay,1:nq) * ptimestep
+     &      pqeff(1:ngrid,1:nlay,1:nq) + 
+     &      sum_subpdq(1:ngrid,1:nlay,1:nq) * microtimestep
       WHERE( zq(1:ngrid,1:nlay,1:nq) < 1.e-30 )
      &     zq(1:ngrid,1:nlay,1:nq) = 1.e-30
@@ -443 +443 @@
               dMh2o = 0.
               do i = 1, nbinco2_cld
-                 Proba=1.0-dexp(-1.*ptimestep*rate(i))
-                 Probah2o=coeffh2o*(1.0-dexp(-1.*ptimestep*rateh2o(i))) !if co2useh2o=.false., this is =0
+                 Proba=1.0-dexp(-1.*microtimestep*rate(i))
+                 Probah2o=coeffh2o*
+     &            (1.0-dexp(-1.*microtimestep*rateh2o(i))) !if co2useh2o=.false., this is =0
                  dNh2o    = dNh2o + n_aer_h2oice(i) * Probah2o
                  dMh2o    = dMh2o + m_aer_h2oice(i) * Probah2o
@@ -517 +518 @@
                  Ic_rice=0.
                  flag_pourri=1
-                 pdtcloudco2(ig,l)=-pdt(ig,l)
+                 subpdtcloudco2(ig,l)=-sum_subpdt(ig,l)
                  dMice=0
                  
               else
-                 dMice=zq(ig,l,igcm_ccnco2_number)*Ic_rice*ptimestep
+                 dMice=zq(ig,l,igcm_ccnco2_number)*Ic_rice*microtimestep
      &                *tauscaling(ig) ! Kg par kg d'air, >0 si croissance !
                  !kg.s-1 par particule * nb particule par kg air*s
@@ -530 +531 @@
 !facteurmax maximum quantity of CO2 that can sublime/condense according to available thermal energy
 ! latent heat release       >0 if growth i.e. if dMice >0
-              pdtcloudco2(ig,l)=dMice*lw/cpp/ptimestep
+              subpdtcloudco2(ig,l)=dMice*lw/cpp/microtimestep
 ! kgco2/kgair* J/kgco2 * 1/(J.kgair-1.K-1)/s= K par seconde
               !Now update tracers 
@@ -588 +589 @@
 
           ! Get cloud tendencies
-        pdqcloudco2(1:ngrid,1:nlay,igcm_co2) =
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_co2) =
      &       (zq(1:ngrid,1:nlay,igcm_co2) - 
-     &       zq0(1:ngrid,1:nlay,igcm_co2))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_co2_ice) =
+     &       zq0(1:ngrid,1:nlay,igcm_co2))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_co2_ice) =
      &       (zq(1:ngrid,1:nlay,igcm_co2_ice) -
-     &       zq0(1:ngrid,1:nlay,igcm_co2_ice))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_h2o_ice) =
+     &       zq0(1:ngrid,1:nlay,igcm_co2_ice))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_h2o_ice) =
      &       (zq(1:ngrid,1:nlay,igcm_h2o_ice) - 
-     &       zq0(1:ngrid,1:nlay,igcm_h2o_ice))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_ccn_mass) =
+     &       zq0(1:ngrid,1:nlay,igcm_h2o_ice))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_ccn_mass) =
      &       (zq(1:ngrid,1:nlay,igcm_ccn_mass) -
-     &       zq0(1:ngrid,1:nlay,igcm_ccn_mass))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_ccn_number) =
+     &       zq0(1:ngrid,1:nlay,igcm_ccn_mass))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_ccn_number) =
      &       (zq(1:ngrid,1:nlay,igcm_ccn_number) -
-     &       zq0(1:ngrid,1:nlay,igcm_ccn_number))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_ccnco2_mass) =
+     &       zq0(1:ngrid,1:nlay,igcm_ccn_number))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_ccnco2_mass) =
      &       (zq(1:ngrid,1:nlay,igcm_ccnco2_mass) -
-     &       zq0(1:ngrid,1:nlay,igcm_ccnco2_mass))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_ccnco2_number) =
+     &       zq0(1:ngrid,1:nlay,igcm_ccnco2_mass))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_ccnco2_number) =
      &       (zq(1:ngrid,1:nlay,igcm_ccnco2_number) -
-     &       zq0(1:ngrid,1:nlay,igcm_ccnco2_number))/ptimestep
-        pdqcloudco2(1:ngrid,1:nlay,igcm_dust_mass) =
+     &       zq0(1:ngrid,1:nlay,igcm_ccnco2_number))/microtimestep
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_dust_mass) =
      &       (zq(1:ngrid,1:nlay,igcm_dust_mass) -
-     &       zq0(1:ngrid,1:nlay,igcm_dust_mass))/ptimestep 
-        pdqcloudco2(1:ngrid,1:nlay,igcm_dust_number) =
+     &       zq0(1:ngrid,1:nlay,igcm_dust_mass))/microtimestep 
+        subpdqcloudco2(1:ngrid,1:nlay,igcm_dust_number) =
      &       (zq(1:ngrid,1:nlay,igcm_dust_number) -
-     &       zq0(1:ngrid,1:nlay,igcm_dust_number))/ptimestep
+     &       zq0(1:ngrid,1:nlay,igcm_dust_number))/microtimestep
 
         end