Changeset 6075 for LMDZ6

Timestamp:

Feb 12, 2026, 10:54:21 PM (2 weeks ago)

Author:

evignon

Message:

petits changements pour pouvoir prescrire la concentration en nombre de gouttelettes et cristaux

Location:

LMDZ6/trunk/libf/phylmd

Files:

• conf_phys_m.f90 (modified) (1 diff)
• lmdz_lscp_ini.f90 (modified) (3 diffs)
• lmdz_lscp_phase.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -2490 +2490 @@
 
     ! flag_aerosol need to be different to zero if ok_cdnc is activated
-    IF (ok_cdnc .AND. flag_aerosol .EQ. 0) THEN
-       CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if ok_cdnc is activated ', 1)
-    ENDIF
+    !IF (ok_cdnc .AND. flag_aerosol .EQ. 0) THEN
+    !   CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if ok_cdnc is activated ', 1)
+    !ENDIF
 
     ! ok_cdnc must be set to y if ok_aie is activated

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90

@@ -244 +244 @@
    REAL, SAVE, PROTECTED :: capa_crystal = 1.                  ! Crystal capacitance (shape factor) for lscp_icefrac_turb [-]
    !$OMP THREADPRIVATE(capa_crystal)
+
+   LOGICAL, SAVE, PROTECTED :: fix_icnc = .FALSE.             ! fix ice crystal number concentration for cloud phase determination  
+   !$OMP THREADPRIVATE(fix_icnc)
+
+   REAL, SAVE, PROTECTED :: icnc0 = 10.                       ! prescribed ice crystal number concentration for cloud phase determination  [L-1]
+   !$OMP THREADPRIVATE(icnc0)
 
    REAL, SAVE, PROTECTED :: naero5 = 0.5                       ! Number concentration of aerosol larger than 0.5 microns [scm-3]
@@ -494 +500 @@
       CALL getin_p('dist_liq', dist_liq)
       CALL getin_p('tresh_cl', tresh_cl)
+      CALL getin_p('fix_icnc', fix_icnc)
+      CALL getin_p('icnc0', icnc0)
       CALL getin_p('capa_crystal', capa_crystal)
       CALL getin_p('naero5', naero5)
@@ -611 +619 @@
       WRITE (lunout, *) 'lscp_ini, capa_crystal', capa_crystal
       WRITE (lunout, *) 'lscp_ini, naero5', naero5
+      WRITE (lunout, *) 'lscp_ini, fix_icnc', fix_icnc
+      WRITE (lunout, *) 'lscp_ini, icnc0', icnc0
       WRITE (lunout, *) 'lscp_ini, gamma_snwretro', gamma_snwretro
       WRITE (lunout, *) 'lscp_ini, gamma_mixth', gamma_mixth

LMDZ6/trunk/libf/phylmd/lmdz_lscp_phase.f90

@@ -163 +163 @@
       USE lmdz_lscp_ini, ONLY: seuil_neb, temp_nowater
       USE lmdz_lscp_ini, ONLY: naero5, gamma_snwretro, gamma_taud, capa_crystal, rho_ice
-      USE lmdz_lscp_ini, ONLY: eps, snow_fallspeed
+      USE lmdz_lscp_ini, ONLY: eps, snow_fallspeed, fix_icnc, icnc0
       USE lmdz_lscp_tools, ONLY: calc_qsat_ecmwf
 
@@ -317 +317 @@
                   !--We assume an exponential ice PSD whose parameters
                   !--are computed following Morrison&Gettelman 2008
-                  !--Ice number density is assumed equals to INP density
+                  !--Ice number density is assumed equals to a constant value
+                  !--if fix_icnc = True
+                  !--or to the INP density
                   !--which is for naero5>0 a function of temperature (DeMott 2010)
                   !--bi and B0 are microphysical function characterizing
@@ -329 +331 @@
                   !--the Southern Ocean and the Antarctic
 
-                  IF (naero5 .LE. 0) THEN
-                     IF (temp(i) .GT. tempvig1) THEN
-                        nb_crystals = 1.e3*10**(-0.14*(temp(i) - tempvig1) - 2.88)
-                     ELSE IF (temp(i) .GT. tempvig2) THEN
-                        nb_crystals = 1.e3*10**(-0.31*(temp(i) - tempvig1) - 2.88)
+                  IF (fix_icnc) THEN 
+                  
+                     nb_crystals = icnc0 * 1000.! from L-1 to m-3
+                     
+                  ELSE
+
+                     IF (naero5 .LE. 0) THEN
+                        IF (temp(i) .GT. tempvig1) THEN
+                           nb_crystals = 1.e3*10**(-0.14*(temp(i) - tempvig1) - 2.88)
+                        ELSE IF (temp(i) .GT. tempvig2) THEN
+                           nb_crystals = 1.e3*10**(-0.31*(temp(i) - tempvig1) - 2.88)
+                        ELSE
+                           nb_crystals = 1.e3*10**(0.)
+                        END IF
                      ELSE
-                        nb_crystals = 1.e3*10**(0.)
-                     END IF
-                  ELSE
-                     nb_crystals = 1.e3*5.94e-5*(RTT - temp(i))**3.33*naero5**(0.0264*(RTT - temp(i)) + 0.0033)
-                  END IF
+                        nb_crystals = 1.e3*5.94e-5*(RTT - temp(i))**3.33*naero5**(0.0264*(RTT - temp(i)) + 0.0033)
+                     END IF
+                  
+                  ENDIF
+
                   lambda_PSD = ((RPI*rho_ice*nb_crystals)/(rho_air*MAX(qiceini_incl, eps)))**(1./3.)
                   N0_PSD = nb_crystals*lambda_PSD