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

Timestamp:

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

Author:

afalco

Message:

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

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (3 diffs)

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'