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

Timestamp:

Mar 11, 2024, 5:26:53 PM (9 months ago)

Author:

afalco

Message:

Pluto PCM:
Methane/CO taken into account in callcorrk
AF

File:

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

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

@@ -55 +55 @@
                               n2cond,nearn2cond,noseason_day,conservn2, &
                               kbo,triton,paleo,paleoyears, &
-                              carbox, methane, UseVdifcPlutold, &
+                              carbox, methane, oldplutovdifc, oldplutocorrk, &
                               aerohaze,haze_proffix,source_haze, &
                               season, sedimentation,generic_condensation, &
@@ -67 +67 @@
       use phys_state_var_mod
       use callcorrk_mod, only: callcorrk
+    !   use callcorrk_pluto_mod, only: callcorrk_pluto
       use vdifc_mod, only: vdifc
       use vdifc_pluto_mod, only: vdifc_pluto
@@ -920 +921 @@
                ! standard callcorrk
                ! clearsky=.false.
-               call callcorrk(ngrid,nlayer,pq,nq,qsurf,                           &
+            !    if (oldplutocorrk) then
+            !       call callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,          &
+            !                    albedo,emis,mu0,pplev,pplay,pt,                   &
+            !                    tsurf,fract,dist_star,aerosol, &
+            !                    zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,   &
+            !                    fluxabs_sw,fluxtop_dn,reffrad,tau_col,ptime,pday, &
+            !                    firstcall,lastcall,zzlay)
+            !    else
+                call callcorrk(ngrid,nlayer,pq,nq,qsurf,                          &
                               albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,   &
                               tsurf,fract,dist_star,aerosol,muvar,                &
@@ -929 +938 @@
                               tau_col,cloudfrac,totcloudfrac,                     &
                               .false.,firstcall,lastcall)
-
+            !    endif ! oldplutocorrk
                 !GG (feb2021): Option to "artificially" decrease the raditive time scale in
                 !the deep atmosphere  press > 0.1 bar. Suggested by MT.
@@ -1042 +1051 @@
          zflubid(1:ngrid)=fluxrad(1:ngrid)+fluxgrd(1:ngrid)
 
-         if (UseVdifcPlutold) then
+         if (oldplutovdifc) then
             zflubid(1:ngrid)=fluxrad(1:ngrid)+fluxgrd(1:ngrid)
             zdum1(:,:) = 0
@@ -1098 +1107 @@
                           taux,tauy)
 
-         else ! if (UseVdifcPlutold) .and. (UseTurbDiff)
+         else ! if (oldplutovdifc) .and. (UseTurbDiff)
 
             zdh(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)/zpopsk(1:ngrid,1:nlayer)