Changeset 3853 for trunk/LMDZ.PLUTO/libf

Timestamp:

Jul 17, 2025, 4:36:29 AM (10 days ago)

Author:

tbertrand

Message:

Pluto PCM:
Implementing the CO cycle in the mode no_n2frost in the VTM and GCM
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• cosurf.F (modified) (2 diffs)
• vdifc_pluto_mod.F90 (modified) (3 diffs)

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

@@ -2 +2 @@
      &  tsurf,pplev,pdpsurf,pq,pdq,pqsurf,pdqsurf,pdqco,pdqsco)
        
-      use callkeys_mod, only: dayfrac
+      use callkeys_mod, only: dayfrac,no_n2frost
       use comcstfi_mod, only: g, r
       use comgeomfi_h 
@@ -52 +52 @@
       DATA vonk/0.4/
       
+      LOGICAL firstcall
+      SAVE firstcall
+      real,dimension(:),save,allocatable :: qsat_co_factor ! factor to prevent co frost formation if no n2 frost
+!$OMP THREADPRIVATE(qsat_co_factor)
+
+      IF (firstcall) THEN
+        firstcall=.false.
+        ! If fixed distribution of N2, then no CO frost either
+        ALLOCATE(qsat_co_factor(ngrid))
+        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
+
       ! Calculation of turbulent flux : F=rho*cdrag*uv*(qsat-zq)
    

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

@@ -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
       LOGICAL firstcall
       SAVE firstcall
+      real,dimension(:),save,allocatable :: qsat_co_factor ! factor to prevent co frost formation if no n2 frost
+!$OMP THREADPRIVATE(qsat_co_factor)
 
       !!read fixed profile for kmix
@@ -161 +162 @@
 
         ! If fixed distribution of N2, then no CO frost either
+        ALLOCATE(qsat_co_factor(ngrid))
         qsat_co_factor(:)=1.
         IF (no_n2frost) then
@@ -662 +664 @@
             !! Prevent CO condensation at the surface
             if (.not.condcosurf) then
-               qsat_co=qsat_co*1.e6
+               qsat_co(:)=qsat_co(:)*1.e6
             endif
             if (no_n2frost) then
-               qsat_co=qsat_co*qsat_co_factor     
+               qsat_co(:)=qsat_co(:)*qsat_co_factor(:)     
             endif