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

Timestamp:

Nov 21, 2025, 2:37:14 PM (4 weeks ago)

Author:

tbertrand

Message:

PLUTO PCM:
Cleaning and adding conservation tests
TB

File:

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

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

@@ -285 +285 @@
       REAL flusurfold(ngrid,nq)  ! old flux cond/sub kg.m-2.s-1
       REAL zplev(ngrid,nlayer+1),zplay(ngrid,nlayer)
+      REAL vecnull(ngrid,nq)      ! null vector used to check conservation of tracer
 
       REAL,SAVE :: ptime0    ! store the first time
@@ -531 +532 @@
 !        Variables set to 0
 !        ~~~~~~~~~~~~~~~~~~
+         vecnull(:,:) = 0.0
          dtrad(:,:) = 0.0
          fluxrad(:) = 0.0
@@ -1182 +1184 @@
 
       if (calldifv) then
-
+         zdqdif(:,:,:)=0.
+         zdqsdif(:,:)=0.
          zflubid(1:ngrid)=fluxrad(1:ngrid)+fluxgrd(1:ngrid)
 
@@ -1361 +1364 @@
          zdvadj(1:ngrid,1:nlayer)=0.0
          zdhadj(1:ngrid,1:nlayer)=0.0
-
+         zdqadj(1:ngrid,1:nlayer,1:nq)=0.0
 
          call convadj(ngrid,nlayer,nq,ptimestep,            &
@@ -1385 +1388 @@
          endif
 
-         ! ! Test water conservation !AF24: removed
+         ! Test tracer conservation 
+         if (methane.and.conservch4) then
+           write(*,*) 'conservch4 calladj'
+           call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdqadj,vecnull,'CH4','calladj')
+         endif
 
       endif ! end of 'calladj'
@@ -1430 +1439 @@
            pdq(1:ngrid,1:nlayer,1:nq)=pdq(1:ngrid,1:nlayer,1:nq) + &
                                       zdq_mix(1:ngrid,1:nlayer,1:nq)
+
+           ! Test tracer conservation
+           if (methane.and.conservch4) then
+             write(*,*) 'conservch4 nonoro'
+             call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdq_mix,vecnull,'CH4','nonorog')
+           endif
+
           ENDIF
-  
   
         ENDIF ! of IF (calllott_nonoro)
@@ -1527 +1544 @@
          zdqch4cloud(:,:,:)=0.
          if ((methane).and.(metcloud).and.(.not.fast)) THEN
+            zdqch4cloud(:,:,:)=0.
             call ch4cloud(ngrid,nlayer,ptimestep,  &
                       pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
@@ -1553 +1571 @@
                ENDDO
             ENDDO
+
+            ! Test tracer conservation
+            if (methane.and.conservch4) then
+             write(*,*) 'conservch4 ch4clouds'
+             call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdqch4cloud,zdqsch4cloud,'CH4',' cloud ')
+            endif
+
          end if
 
@@ -1558 +1585 @@
          ! CO ice condensation in the atmosphere
          !--------------------------------------
-         zdqcocloud(:,:,:)=0.
          IF ((carbox).and.(monoxcloud).and.(.not.fast)) THEN
+            zdqcocloud(:,:,:)=0.
             call cocloud(ngrid,nlayer,ptimestep,   &
                       pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
@@ -1614 +1641 @@
                               zdqphot_prec,zdqphot_ch4,zdqconv_prec,declin)
                pdq(:,:,:) = pdq(:,:,:) + pdqmufi_prod(:,:,:) ! Should be updated
+
+               ! Test tracer conservation
+               if (methane.and.conservch4) then
+                 write(*,*) 'conservch4 hazeprod'
+                 call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,pdqmufi_prod,vecnull,'CH4','hazepro')
+               endif
+
             ENDIF ! end call_haze_prod_pCH4
 
@@ -1673 +1709 @@
             ENDIF ! end fast.and.fasthaze
 
+            ! Test tracer conservation
+            if (methane.and.conservch4) then
+                 write(*,*) 'conservch4 zdqhaze'
+                 call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdqhaze,vecnull,'CH4','zdqhaze')
+            endif
+
             ! Sedimentation.
             if (sedimentation) then
@@ -1691 +1735 @@
             end if ! end of 'sedimentation'
 
+            ! Test tracer conservation
+            if (methane.and.conservch4) then
+                 write(*,*) 'conservch4 sedim'
+                 call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdqsed,zdqssed,'CH4',' sedim ')
+            endif
+
          ENDIF ! end callmufi
 
@@ -1726 +1778 @@
             pdpsrf(1:ngrid)            = pdpsrf(1:ngrid) + zdpsrfmr(1:ngrid)
             zdtsurf(1:ngrid)           = zdtsurf(1:ngrid) + zdtsurfmr(1:ngrid)
+
+            ! Test tracer conservation
+            if (methane.and.conservch4) then
+                 write(*,*) 'conservch4 mass redistrib'
+                 call testconserv(ngrid,nlayer,nq,pq,pdq,qsurf,dqsurf, &
+                   igcm_ch4_gas,igcm_ch4_ice, &
+                   ptimestep,pplev,zdqmr,zdqsurfmr,'CH4','massred')
+            endif
 
          endif
@@ -2295 +2355 @@
       ! Total energy balance diagnostics
       if(callrad)then
+         call write_output("mu0","ycosinus incident angle"," ",mu0)
          call write_output("ALB","Surface albedo"," ",albedo_equivalent)
          call write_output("ASR","absorbed stellar rad.","W m-2",fluxabs_sw)