Changeset 3361

Timestamp:

Jun 10, 2024, 11:53:45 AM (7 months ago)

Author:

tbertrand

Message:

LMDZ.PLUTO
Fixing a bug in physiq.F when calling callcorrk_pluto (the old version of callcorrk for pluto)
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (1 diff)
• callcorrk_pluto_mod.F90 (modified) (4 diffs)
• callkeys_mod.F90 (modified) (1 diff)
• gfluxi.F (modified) (1 diff)
• physiq_mod.F90 (modified) (4 diffs)
• sfluxi_pluto_mod.F (modified) (3 diffs)
• sfluxv_pluto_mod.F (modified) (3 diffs)
• surfprop.F90 (modified) (1 diff)

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

@@ -856 +856 @@
 
       plevrad(1) = 0.
-!      plevrad(2) = 0.   !! JL18 enabling this line puts the radiative top at p=0 which was the idea before, but does not seem to perform best after all.
+      !!plevrad(2) = 0.   !! JL18 enabling this line puts the radiative top at p=0 which was the idea before, but does not seem to perform best after all.
 
       tlevrad(1) = tlevrad(2)

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

@@ -265 +265 @@
          Cmk= 0.01 * 1.0 / (g * mugaz * 1.672621e-27) ! q_main=1.0 assumed
 
-         epsi_ch4=mmol(igcm_ch4_gas)/mmol(igcm_n2)
+         if (methane) then
+           epsi_ch4=mmol(igcm_ch4_gas)/mmol(igcm_n2)
+         endif
 
          ! If fixed profile of CH4 gas
@@ -280 +282 @@
 
 !=======================================================================
-
 !     L_NSPECTV is the number of Visual(or Solar) spectral intervals
 !     how much light we get
@@ -653 +654 @@
          ELSE                          ! during the night, fluxes = 0
             nfluxtopv=0.0
+            fmnetv_nu(:,:)=0.0
+            nfluxtopv_nu(:)=0.0
             DO l=1,L_NLAYRAD
                fmnetv(l)=0.0
@@ -712 +715 @@
                end do
          END DO
-
          ! values at top of atmosphere
          dpp = g/(cpp*scalep*(plevrad(3)-plevrad(1)))

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

@@ -83 +83 @@
       logical,save :: convergeps,conservn2,condensn2,no_n2frost
 !$OMP THREADPRIVATE(convergeps,conservn2,condensn2,no_n2frost)
-      logical,save :: globmean1d,kmix_proffix
-!$OMP THREADPRIVATE(globmean1d,kmix_proffix)
+      logical,save :: kmix_proffix
+!$OMP THREADPRIVATE(kmix_proffix)
       logical,save :: ch4fix,vmrch4_proffix
 !$OMP THREADPRIVATE(ch4fix,vmrch4_proffix)

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

@@ -86 +86 @@
 !-----------------------------------------------------------------------
          if (W0(L).eq.1.D0) then
-            W0(L) = 0.99999D0
+               W0(L) = 0.99999D0
          endif
          

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

@@ -981 +981 @@
                                firstcall,lastcall)
                   albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
+                  fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+       &
+                               fluxsurf_sw(1:ngrid)*(1.-albedo(1:ngrid,1))
                else
                 call callcorrk(ngrid,nlayer,pq,nq,qsurf,  &
@@ -991 +993 @@
                               tau_col,cloudfrac,totcloudfrac,                     &
                               .false.,firstcall,lastcall)
+                  ! Radiative flux from the sky absorbed by the surface (W.m-2).
+                  GSR=0.0
+                  fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+        &
+                              fluxsurfabs_sw(1:ngrid)
                endif ! oldplutocorrk
                 !GG (feb2021): Option to "artificially" decrease the raditive time scale in
@@ -1010 +1016 @@
                 ! AF24: removed CLFvarying Option
 
-               ! Radiative flux from the sky absorbed by the surface (W.m-2).
-               GSR=0.0
-               fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+fluxsurfabs_sw(1:ngrid)
 
                             !if(noradsurf)then ! no lower surface; SW flux just disappears
@@ -2166 +2169 @@
 
          call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
-         !call writediagfi(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
+         call writediagfi(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
          call writediagfi(ngrid,"ISR","incoming stellar rad.","W m-2",2,fluxtop_dn)
          call writediagfi(ngrid,"ASR","absorbed stellar rad.","W m-2",2,fluxabs_sw)

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

@@ -14 +14 @@
       use radcommon_h, only: planckir, tlimit,sigma, gweight
       use comcstfi_mod, only: pi
-      use gfluxi_mod, only: gfluxi
+      use gfluxi_old_mod, only: gfluxi_old
 
 
@@ -136 +136 @@
 C         WITHIN EACH INTERVAL AT THE MIDPOINT WAVENUMBER
 
-          CALL GFLUXI(NLEVRAD,TLEV,NW,DWNI(NW),DTAUI(1,NW,NG),
+          CALL GFLUXI_OLD(NLEVRAD,TLEV,NW,DWNI(NW),DTAUI(1,NW,NG),
      *                TAUCUMI(1,NW,NG),
      *                WBARI(1,NW,NG),COSBI(1,NW,NG),UBARI,RSFI,BTOP,
@@ -193 +193 @@
 
 
-       CALL GFLUXI(NLEVRAD,TLEV,NW,DWNI(NW),DTAUI(1,NW,NG),
+       CALL GFLUXI_OLD(NLEVRAD,TLEV,NW,DWNI(NW),DTAUI(1,NW,NG),
      *                TAUCUMI(1,NW,NG),
      *                WBARI(1,NW,NG),COSBI(1,NW,NG),UBARI,RSFI,BTOP,

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

@@ -12 +12 @@
       use radinc_h, only: L_NLAYRAD, L_NLEVRAD
       use radcommon_h, only: tlimit, gweight
-      use gfluxv_mod, only: gfluxv
+      use gfluxv_old_mod, only: gfluxv_old
 
       implicit none
@@ -100 +100 @@
 
 
-          CALL GFLUXV(DTAUV(1,NW,NG),TAUV(1,NW,NG),TAUCUMV(1,NW,NG),
+          CALL GFLUXV_OLD(DTAUV(1,NW,NG),TAUV(1,NW,NG),TAUCUMV(1,NW,NG),
      *                WBARV(1,NW,NG),COSBV(1,NW,NG),UBAR0,F0PI,RSFV,
      *                BTOP,BSURF,FMUPV,FMDV,DIFFV,FLUXUP,FLUXDN)
@@ -155 +155 @@
 C       RETURN FLUXES FOR A GIVEN NT
 
-        CALL GFLUXV(DTAUV(1,NW,NG),TAUV(1,NW,NG),TAUCUMV(1,NW,NG),
+        CALL GFLUXV_OLD(DTAUV(1,NW,NG),TAUV(1,NW,NG),TAUCUMV(1,NW,NG),
      *              WBARV(1,NW,NG),COSBV(1,NW,NG),UBAR0,F0PI,RSFV,
      *              BTOP,BSURF,FMUPV,FMDV,DIFFV,FLUXUP,FLUXDN)

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

@@ -460 +460 @@
            enddo
         ENDDO
-        print*, emin
-        print*, tidat(1,:)
-        print*, mlayer(:)
 
       ELSE