Changeset 3421 for trunk/LMDZ.PLUTO/libf/phypluto

Timestamp:

Aug 26, 2024, 6:11:23 PM (16 months ago)

Author:

afalco

Message:

Pluto PCM: Added vertical mixing in no gcm, ch4/co surf, no_n2frost.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• ch4surf.F (added)
• condense_n2.F90 (modified) (6 diffs)
• cosurf.F (added)
• physiq_mod.F90 (modified) (3 diffs)

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

@@ -116 +116 @@
   real zdtlatent (klon,klev)
   real zdt (klon,klev)
-  REAL albediceF(klon)
+!  REAL albediceF(klon)
 !   SAVE albediceF
   INTEGER nsubtimestep,itsub    !number of subtimestep when calling vl1d
@@ -130 +130 @@
   REAL globzplevnew
 
-  REAL vmrn2(klon)
-!   SAVE vmrn2
+  real,dimension(:),save,allocatable :: vmrn2
+!$OMP THREADPRIVATE(vmrn2)
   REAL stephan
   DATA stephan/5.67e-08/  ! Stephan Boltzman constant
@@ -179 +179 @@
      ENDIF
 
+     ALLOCATE(vmrn2(klon))
      vmrn2(:) = 1.
      !IF (ch4lag) then
@@ -187 +188 @@
      !   ENDDO
      !ENDIF
-     !IF (no_n2frost) then
-     !   DO ig=1,klon
-     !      if (picen2(ig).eq.0.) then
-     !         vmrn2(ig) = 1.e-15
-     !      endif
-     !   ENDDO
-     !ENDIF
+     IF (no_n2frost) then
+        DO ig=1,klon
+           if (picen2(ig).eq.0.) then
+              vmrn2(ig) = 1.e-15
+           endif
+        ENDDO
+     ENDIF
      firstcall=.false.
   ENDIF
@@ -344 +345 @@
    DO ig=1,klon
      ! forecast of frost temperature ztcondsol
-     ztcondsol(ig) = tcond_n2(zplev(ig),zqn2(ig,1))
-     !ztcondsol(ig) = tcond_n2(zplev(ig),vmrn2(ig))
+     !ztcondsol(ig) = tcond_n2(zplev(ig),zqn2(ig,1))
+     ztcondsol(ig) = tcond_n2(zplev(ig),vmrn2(ig))
 
 !     Loop over where we have condensation / sublimation
@@ -898 +899 @@
    return
    end  subroutine vl1d
-

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

@@ -33 +33 @@
       USE tracer_h, only: noms, mmol, radius, rho_q, qext, &
                           igcm_n2,igcm_ch4_gas,igcm_ch4_ice,igcm_haze,&
-                          igcm_co_gas,igcm_co_ice,igcm_prec_haze,lw_n2,&
+                          igcm_co_gas,igcm_co_ice,igcm_prec_haze,lw_n2,lw_ch4,lw_co,&
                           alpha_lift, alpha_devil, qextrhor, &
                           nesp, is_chim, is_condensable,constants_epsi_generic
@@ -55 +55 @@
                               n2cond,nearn2cond,noseason_day,conservn2, &
                               convergeps,kbo,triton,paleo,paleoyears,glaflow, &
-                              carbox, methane,&
+                              carbox, methane,condmetsurf,condcosurf,&
                               oldplutovdifc,oldplutocorrk,oldplutosedim, &
                               aerohaze,haze_proffix,source_haze, tsurfmax, &
@@ -1248 +1248 @@
 
          ! if(.not.newtonian)then
-            zdtsurf(1:ngrid) = zdtsurf(1:ngrid) + (fluxrad(1:ngrid) + fluxgrd(1:ngrid))/capcal(1:ngrid)
+         zdtsurf(1:ngrid) = zdtsurf(1:ngrid) + (fluxrad(1:ngrid) + fluxgrd(1:ngrid))/capcal(1:ngrid)
+
+!        ------------------------------------------------------------------
+!        Methane surface sublimation and condensation in fast model (nogcm)
+!        ------------------------------------------------------------------
+         if ((methane).and.(fast).and.condmetsurf) THEN
+
+            call ch4surf(ngrid,nlayer,nq,ptimestep, &
+               tsurf,zdtsurf,pplev,pdpsrf,pq,pdq,qsurf,dqsurf, &
+               zdqch4fast,zdqsch4fast)
+
+            dqsurf(1:ngrid,igcm_ch4_ice)= dqsurf(1:ngrid,igcm_ch4_ice) + &
+                                         zdqsch4fast(1:ngrid)
+            pdq(1:ngrid,1,igcm_ch4_gas)= pdq(1:ngrid,1,igcm_ch4_gas) + &
+                                         zdqch4fast(1:ngrid)
+            zdtsurf(1:ngrid)=zdtsurf(1:ngrid)+lw_ch4*zdqsch4fast(1:ngrid)/capcal(1:ngrid)
+            end if
+!        ------------------------------------------------------------------
+!        CO surface sublimation and condensation in fast model (nogcm)
+!        ------------------------------------------------------------------
+         if ((carbox).and.(fast).and.condcosurf) THEN
+
+            call cosurf(ngrid,nlayer,nq,ptimestep, &
+               tsurf,pplev,pdpsrf,pq,pdq,qsurf,dqsurf, &
+               zdqcofast,zdqscofast)
+
+            dqsurf(1:ngrid,igcm_co_ice)= dqsurf(1:ngrid,igcm_co_ice) + &
+                                         zdqscofast(1:ngrid)
+            pdq(1:ngrid,1,igcm_co_gas)= pdq(1:ngrid,1,igcm_co_gas) + &
+                                        zdqcofast(1:ngrid)
+            zdtsurf(1:ngrid)=zdtsurf(1:ngrid)+lw_co*zdqscofast(1:ngrid)/capcal(1:ngrid)
+         end if
+
 
       endif ! end of 'calldifv'