Ignore:
Timestamp:
Jul 17, 2025, 4:36:29 AM (12 days ago)
Author:
tbertrand
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.PLUTO/libf/phypluto/vdifc_pluto_mod.F90

    r3844 r3853  
    100100      REAL zcst1
    101101      REAL zu2
    102       REAL qsat_co_factor(ngrid)   ! factor to prevent co frost formation if no n2 frost
    103102      EXTERNAL SSUM,SCOPY
    104103      REAL SSUM
    105104      LOGICAL firstcall
    106105      SAVE firstcall
     106      real,dimension(:),save,allocatable :: qsat_co_factor ! factor to prevent co frost formation if no n2 frost
     107!$OMP THREADPRIVATE(qsat_co_factor)
    107108
    108109      !!read fixed profile for kmix
     
    161162
    162163        ! If fixed distribution of N2, then no CO frost either
     164        ALLOCATE(qsat_co_factor(ngrid))
    163165        qsat_co_factor(:)=1.
    164166        IF (no_n2frost) then
     
    662664            !! Prevent CO condensation at the surface
    663665            if (.not.condcosurf) then
    664                qsat_co=qsat_co*1.e6
     666               qsat_co(:)=qsat_co(:)*1.e6
    665667            endif
    666668            if (no_n2frost) then
    667                qsat_co=qsat_co*qsat_co_factor     
     669               qsat_co(:)=qsat_co(:)*qsat_co_factor(:)     
    668670            endif
    669671
Note: See TracChangeset for help on using the changeset viewer.