Changeset 2562 for trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

Timestamp:

Sep 21, 2021, 9:02:49 AM (3 years ago)

Author:

cmathe

Message:

GCM MARS: CO2 clouds microphysics improvements

File:

• trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F (modified) (5 diffs)

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

@@ -159 +159 @@
       INTEGER,SAVE :: iccnco2_number ! index of tracer containing CCN number
       INTEGER,SAVE :: iccnco2_mass ! index of tracer containing CCN number
+      INTEGER,SAVE :: iccnco2_h2o_number ! index of tracer containing CCN number
+      INTEGER,SAVE :: iccnco2_h2o_mass_ice ! index of tracer containing CCN number
+      INTEGER,SAVE :: iccnco2_h2o_mass_ccn ! index of tracer containing CCN number
       INTEGER,SAVE :: ico2_ice ! index of tracer containing CCN number
 
@@ -241 +244 @@
           iccnco2_mass=0
           iccnco2_number=0
+          iccnco2_h2o_mass_ice=0
+          iccnco2_h2o_mass_ccn=0
+          iccnco2_h2o_number=0
           ico2_ice=0
           do iq=1,nq
@@ -254 +260 @@
               iccnco2_number=iq
               write(*,*)"callsedim: iccnco2_number=",iccnco2_number
+            endif
+            if (noms(iq).eq."ccnco2_h2o_number") then
+              iccnco2_h2o_number=iq
+              write(*,*)"callsedim: iccnco2_h2o_number=",
+     &                   iccnco2_h2o_number
+            endif
+            if (noms(iq).eq."ccnco2_h2o_mass_ice") then
+              iccnco2_h2o_mass_ice=iq
+              write(*,*)"callsedim: iccnco2_h2o_mass_ice=",
+     &                 iccnco2_h2o_mass_ice
+            endif
+            if (noms(iq).eq."ccnco2_h2o_mass_ccn") then
+              iccnco2_h2o_mass_ccn=iq
+              write(*,*)"callsedim: iccnco2_h2o_mass_ccn=",
+     &                 iccnco2_h2o_mass_ccn
             endif
           enddo
@@ -263 +284 @@
             call abort_physic(modname,"missing co2 ccn tracers",1)
           endif
+          if (co2useh2o) then
+            if((iccnco2_h2o_mass_ccn.eq.0).or.
+     &         (iccnco2_h2o_mass_ice.eq.0).or.
+     &         (iccnco2_h2o_number.eq.0)) then
+            write(*,*) 'callsedim: error! could not identify'
+            write(*,*) ' tracers for ccn co2 mass and number mixing'
+            write(*,*) ' ratio and co2clouds are activated!'
+            call abort_physic(modname,"missing co2 ccn tracers",1)
+           end if
+          end if
        ENDIF                    !of if (co2clouds)
 
@@ -398 +429 @@
         if(radius(iq).gt.1.e-9 .and.(iq.ne.ico2_ice) .and.
      &        (iq .ne. iccnco2_mass) .and. (iq .ne. 
-     &        iccnco2_number) .and. ! no sedim for gaz or CO2 clouds  (done in microtimestep)
+     &        iccnco2_number) .and.
+     &        (iq.ne.iccnco2_h2o_number).and.
+     &        (iq.ne.iccnco2_h2o_mass_ice).and.
+     &        (iq.ne.iccnco2_h2o_mass_ccn).and.  ! no sedim for gaz or CO2 clouds  (done in microtimestep)
      &        iq .ne. igcm_hdo_ice) then !MVals: hdo is transported by h2o
 c -----------------------------------------------------------------