Changeset 3844

Timestamp:

Jul 13, 2025, 10:11:41 PM (4 days ago)

Author:

tbertrand

Message:

Pluto PCM
CO cycle now possible with the option no_n2frost (CO condensation and sublimation is limited to within the N2 ice deposits only).
NLTE factor 0.33 on the cooling rates now removed in callcorrk.
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (1 diff)
• cocloud.F90 (modified) (1 diff)
• cosat.F90 (modified) (3 diffs)
• cosurf.F (modified) (1 diff)
• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (1 diff)
• vdifc_pluto_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

@@ -1100 +1100 @@
             dtlw_nu(l,nw)=dtlw_nu(l,nw)*eps_nlte_lw(ig,l)
           else
-            !dtlw_nu(l,nw)=1.*dtlw_nu(l,nw) ! no CO correction (Strobbel 1996)
-            dtlw_nu(l,nw)=0.33*dtlw_nu(l,nw) ! CO correction (Strobbel 1996)
-!               dtlw_co(ig,l)=dtlw_co(ig,l)+ dtlw_nu(l,nw) ! diagnostic
+            dtlw_nu(l,nw)=1.*dtlw_nu(l,nw) ! no CO correction (Strobbel 1996)
+            !dtlw_nu(l,nw)=0.33*dtlw_nu(l,nw) ! CO correction (Strobbel 1996)
+            !dtlw_co(ig,l)=dtlw_co(ig,l)+ dtlw_nu(l,nw) ! diagnostic
           end if
           dtlw(ig,l)=dtlw(ig,l)+ dtlw_nu(l,nw) !average now on each wavelength

trunk/LMDZ.PLUTO/libf/phypluto/cocloud.F90

@@ -159 +159 @@
 !         enddo
 
-         call cosat(ngrid*nlay,zt,pplay,zqsat,vecnull,vecnull)
+         call cosat(ngrid*nlay,zt,pplay,zqsat,vecnull)
 !        TEMPORAIRE :
 !        test sans condensation atmospherique

trunk/LMDZ.PLUTO/libf/phypluto/cosat.F90

@@ -1 @@
-subroutine cosat(nsize,t,p,qsat,qsurf_n2,qsurf_ch4)
+subroutine cosat(nsize,t,p,qsat,qsurf_n2)
 
   IMPLICIT NONE
@@ -15 +15 @@
       real t(nsize) , p(nsize)
 !   OUTPUT
-      real qsat(nsize),qsurf_n2(nsize),qsurf_ch4(nsize)
+      real qsat(nsize),qsurf_n2(nsize)
       INTEGER i
 
@@ -21 +21 @@
       do i=1,nsize
 
-!          from Fray and schmitt fit by formulation L = 2.74e5 J/kg
-           qsat(i)=0.1537*exp((28*274./8.314)*(1/68.1-1/t(i)))*100000*28.0/(28.0*p(i))
-
-!          Raoult law if mixte in N2 ice and CH4 ice : assume 0.3% CO in N2
-!          (Merlin) and 0.3% in CH4
-!          if ((qsurf_n2(i).gt.0.001).or.(qsurf_ch4(i).gt.0.001)) then
+!       from Fray and schmitt fit by formulation L = 2.74e5 J/kg
+        qsat(i)=0.1537*exp((28*274./8.314)*(1/68.1-1/t(i)))*100000*28.0/(28.0*p(i))
 
 !       Raoult law if mixte in N2 ice : assume 0.3% CO in N2

trunk/LMDZ.PLUTO/libf/phypluto/cosurf.F

@@ -74 +74 @@
       ENDDO
 
-      call cosat(ngrid,tsurf,zpsrf,qsat,pqsurf(:,igcm_n2),
-     &                                       pqsurf(:,igcm_ch4_ice))
+      call cosat(ngrid,tsurf,zpsrf,qsat,pqsurf(:,igcm_n2))
 
       DO ig=1,ngrid

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -1364 +1364 @@
          call abort_physic(rname, 'for now, haze/rad_proffix only works with optichaze=T', 1)
      endif
-      if (carbox.and.condcosurf.and.no_n2frost) then
-        call abort_physic(rname, "CO surface condensation and no_n2frost are both active which may not be relevant", 1)
-      end if
+     ! if (carbox.and.condcosurf.and.no_n2frost) then
+     !   call abort_physic(rname, "CO surface condensation and no_n2frost are both active which may not be relevant", 1)
+     ! end if
 
      if ((cpp_mugaz_mode >= 1).and.(is_master).and.(ngrid>1)) then

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -2110 +2110 @@
       endif ! end of 'lastcall'
 
+      if(mod(icount,diagfi_output_rate).eq.0) then
+         write_restartfi = .true.
+      endif
       if(startfi_output_rate.ne.0)then
          if (mod(icount,startfi_output_rate).eq.0) then

trunk/LMDZ.PLUTO/libf/phypluto/vdifc_pluto_mod.F90

@@ -15 +15 @@
       use comgeomfi_h
       use callkeys_mod, only: carbox, methane, condcosurf, condensn2, condmetsurf,&
-                              kmix_proffix, vertdiff, tracer, kmixmin
+                              kmix_proffix, vertdiff, tracer, kmixmin, no_n2frost
       use datafile_mod, only: datadir
       use surfdat_h, only: phisfi
@@ -100 +100 @@
       REAL zcst1
       REAL zu2
-
+      REAL qsat_co_factor(ngrid)   ! factor to prevent co frost formation if no n2 frost
       EXTERNAL SSUM,SCOPY
       REAL SSUM
@@ -159 +159 @@
             close(114)
         ENDIF
+
+        ! If fixed distribution of N2, then no CO frost either
+        qsat_co_factor(:)=1.
+        IF (no_n2frost) then
+           DO ig=1,ngrid
+              if (pqsurf(ig,igcm_n2).eq.0.) then
+                   qsat_co_factor(ig) = 1.e6
+              endif
+           ENDDO
+        ENDIF
+
       ENDIF
 
@@ -645 +656 @@
           ELSE IF (carbox.and.(iq.eq.igcm_co_gas)) then
 
-           ! calcul de la valeur de q a la surface :
+            !! Calculating saturation mixing ratio at surface 
             call cosat(ngrid,ptsrf,pplev(1,1),qsat_co, &
-                   pqsurf(:,igcm_n2),pqsurf(:,igcm_ch4_ice))
+                   pqsurf(:,igcm_n2))
 
             !! Prevent CO condensation at the surface
             if (.not.condcosurf) then
                qsat_co=qsat_co*1.e6
+            endif
+            if (no_n2frost) then
+               qsat_co=qsat_co*qsat_co_factor     
             endif