Changeset 2659 for trunk/LMDZ.MARS/libf/phymars/initracer.F

Timestamp:

Apr 7, 2022, 4:16:37 PM (3 years ago)

Author:

cmathe

Message:

ticket 100 resolved

File:

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

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

@@ -47 +47 @@
 c  rho_dust          ! Mars dust density
 c  rho_ice           ! Water ice density
-c  nuice_ref         ! Effective variance nueff of the
-c                    !   water-ice size distributions
 c  doubleq           ! if method with mass (iq=1) and number(iq=2) mixing ratio
 c  varian            ! Characteristic variance of log-normal distribution
@@ -520 +518 @@
       rho_ice=920.    ! Water ice density (kg.m-3)
       rho_ice_co2=1650. 
-      !Mangan et al., Icarus 2017 :CO2 density = 1.72391-2.53×10−4T – 2.87×10−6T^2
-      nuice_ref=0.1   ! Effective variance nueff of the
-                      ! water-ice size distribution
-      !!!nuice_sed=0.45   ! Sedimentation effective variance
-                      ! of the water-ice size distribution
-
-      nuiceco2_ref=0.2    !C.M. Effective variance nueff of the
-                          ! co2-ice size distribution
+
       if (doubleq) then
 c       "doubleq" technique 
@@ -720 +711 @@
 
       
-! Initialisation for CO2 clouds
-      if (co2clouds ) then 
-        radius(igcm_ccnco2_mass) = radius(igcm_dust_mass)
-        alpha_lift(igcm_ccnco2_mass) = 1e-30
-        alpha_devil(igcm_ccnco2_mass) = 1e-30
-        rho_q(igcm_ccnco2_mass) = rho_dust
-        radius(igcm_ccnco2_number) = radius(igcm_ccnco2_mass)
-        alpha_lift(igcm_ccnco2_number) = alpha_lift(igcm_ccnco2_mass)
-        alpha_devil(igcm_ccnco2_number) = alpha_devil(igcm_ccnco2_mass)
-        rho_q(igcm_ccnco2_number) = rho_q(igcm_ccnco2_mass)
-     
-        radius(igcm_co2)=0.
-        alpha_lift(igcm_co2) =0.
-        alpha_devil(igcm_co2)=0.
-        radius(igcm_co2_ice)=1.e-8
-        rho_q(igcm_co2_ice)=rho_ice_co2
-        alpha_lift(igcm_co2_ice) =0.
-        alpha_devil(igcm_co2_ice)=0.
-
-      endif 
-      
+
 c     Output for records:
 c     ~~~~~~~~~~~~~~~~~~
@@ -840 +811 @@
              call abort_physic("initracer","co2 clouds issue",1)
           end if
-          ! TODO: ajouter les tests avec co2useh2o, et meteo_flux, ainsi que les traceurs ccnco2_number
+          if (igcm_ccnco2_number .eq. 0) then
+             write(*,*) "initracer: error !!"
+             write(*,*) "  cannot use co2 clouds option without ",
+     &            "a ccnco2_number tracer !"
+             call abort_physic("initracer","co2 clouds issue",1)
+          end if
+          if (igcm_ccnco2_mass .eq. 0) then
+             write(*,*) "initracer: error !!"
+             write(*,*) "  cannot use co2 clouds option without ",
+     &            "a ccnco2_mass tracer !"
+             call abort_physic("initracer","co2 clouds issue",1)
+          end if
+          if (co2useh2o) then
+            if (igcm_ccnco2_h2o_number .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_h2o_number tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+            if (igcm_ccnco2_h2o_mass_ice .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_h2o_mass_ice tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+            if (igcm_ccnco2_h2o_mass_ccn .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_h2o_mass_ccn tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+          end if
+          if (meteo_flux) then
+            if (igcm_ccnco2_meteor_number .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_meteor_number tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+            if (igcm_ccnco2_meteor_mass .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_h2o_mass_ice tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+            if (igcm_ccnco2_meteor_mass .eq. 0) then
+               write(*,*) "initracer: error !!"
+               write(*,*) "  cannot use co2 clouds option without ",
+     &              "a ccnco2_meteor_mass tracer !"
+               call abort_physic("initracer","co2 clouds issue",1)
+            end if
+          end if
        endif
  
@@ -999 +1021 @@
        endif ! of if (photochem .or. callthermos)
 
+! Initialisation for CO2 clouds
+      if (co2clouds) then
+        radius(igcm_ccnco2_mass) = radius(igcm_dust_mass)
+        alpha_lift(igcm_ccnco2_mass) = 1e-30
+        alpha_devil(igcm_ccnco2_mass) = 1e-30
+        rho_q(igcm_ccnco2_mass) = rho_dust
+        radius(igcm_ccnco2_number) = radius(igcm_ccnco2_mass)
+        alpha_lift(igcm_ccnco2_number) = alpha_lift(igcm_ccnco2_mass)
+        alpha_devil(igcm_ccnco2_number) = alpha_devil(igcm_ccnco2_mass)
+        rho_q(igcm_ccnco2_number) = rho_q(igcm_ccnco2_mass)
+
+        radius(igcm_co2)=0.
+        alpha_lift(igcm_co2) =0.
+        alpha_devil(igcm_co2)=0.
+        radius(igcm_co2_ice)=1.e-8
+        rho_q(igcm_co2_ice)=rho_ice_co2
+        alpha_lift(igcm_co2_ice) =0.
+        alpha_devil(igcm_co2_ice)=0.
+
+      endif
+
       end