Changeset 3577 for trunk/LMDZ.PLUTO/libf

Timestamp:

Jan 15, 2025, 11:04:37 AM (5 days ago)

Author:

tbertrand

Message:

LMDZ.PLUTO:
update of routine for testing mass conservation
TB

File:

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

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

@@ -4 +4 @@
       use comcstfi_mod, only: g
       use geometry_mod, only: cell_area
-      USE comgeomfi_h, only: totarea
+      USE comgeomfi_h, only: totarea,totarea_planet
       implicit none
 
@@ -49 +49 @@
       pstot=0.0
 
-      do ig = 2, ngrid-1
+      do ig = 1, ngrid
 
          ! mass atm kg/m2
@@ -62 +62 @@
 
       enddo
-      atmmasstot= atmmasstot + ps(1)/g*cell_area(1)*iim
-      atmmasstot= atmmasstot + ps(ngrid)/g*cell_area(ngrid)*iim
-
-      surfmasstot= surfmasstot + surfmass(1)*cell_area(1)*iim
-      surfmasstot= surfmasstot + surfmass(ngrid)*cell_area(ngrid)*iim
-
-      pstot= pstot + ps(1)*cell_area(1)*iim
-      pstot= pstot + ps(ngrid)*cell_area(ngrid)*iim
                
-      atmmasstot=  atmmasstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
-      pstot=  pstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
-      surfmasstot = surfmasstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
+      atmmasstot=  atmmasstot/totarea_planet
+      pstot=  pstot/totarea_planet
+      surfmasstot = surfmasstot/totarea_planet
 
       totmass=surfmasstot+atmmasstot