Changeset 3992 for trunk/LMDZ.PLUTO/libf

Timestamp:

Dec 23, 2025, 4:10:55 PM (3 days ago)

Author:

tbertrand

Message:

PLUTO PCM :
Fixing bug of tracer index when asking a tendancy as output
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• condense_n2.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (1 diff)

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

@@ -666 +666 @@
               write(*,*) 'zmflux*DT, l=1' ,  zmflux(1)*ptimestep
               write(*,*) 'zmflux*DT, l=2' ,  zmflux(2)*ptimestep
-          write(*,*) 'pq, l=1,2,3' ,  pq(1,1,1), pq(1,2,1),pq(1,3,1)
-          write(*,*) 'zq, l=1,2,3' ,  zq(1,1), zq(2,1),zq(3,1)
-              write(*,*) 'dq*Dt l=1' ,  pdq(1,1,1)*ptimestep
-              write(*,*) 'dqcond*Dt l=1' ,  pdqc(1,1,1)*ptimestep
+          write(*,*) 'pq, l=1,2,3' ,  pq(1,1,igcm_n2), pq(1,2,igcm_n2),pq(1,3,igcm_n2)
+          write(*,*) 'zq, l=1,2,3' ,  zq(1,igcm_n2), zq(2,igcm_n2),zq(3,igcm_n2)
+              write(*,*) 'dq*Dt l=1' ,  pdq(1,1,igcm_n2)*ptimestep
+              write(*,*) 'dqcond*Dt l=1' ,  pdqc(1,1,igcm_n2)*ptimestep
      end if
 

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

@@ -2485 +2485 @@
             call write_output("zduc","tendancy U cond N2","m s-1",zduc)
             call write_output("zdvc","tendancy V cond N2","m s-1",zdvc)
-            call write_output("zdqc_n2","tendancy tracer cond N2","kg kg-1 s-1",zdqc(:,:,1))
-            call write_output("zdqsc_n2","tendancy tracer surf cond N2","kg kg-1 s-1",zdqsc(:,1))
-            call write_output("zdqdif_n2","tendancy tracer diff","kg kg-1 s-1",zdqdif(:,:,1))
-            call write_output("zdqsdif_n2","tendancy tracer surf diff","kg kg-1 s-1",zdqsdif(:,1))
-            call write_output("zdqadj_n2","tendancy tracer adj","K s-1",zdqadj(:,:,1))
+            call write_output("zdqc_n2","tendancy tracer cond N2","kg kg-1 s-1",zdqc(:,:,igcm_n2))
+            call write_output("zdqsc_n2","tendancy tracer surf cond N2","kg kg-1 s-1",zdqsc(:,igcm_n2))
+            call write_output("zdqdif_n2","tendancy tracer diff","kg kg-1 s-1",zdqdif(:,:,igcm_n2))
+            call write_output("zdqsdif_n2","tendancy tracer surf diff","kg kg-1 s-1",zdqsdif(:,igcm_n2))
+            call write_output("zdqadj_n2","tendancy tracer adj","K s-1",zdqadj(:,:,igcm_n2))
          endif