Changeset 2538

Timestamp:

Jun 22, 2021, 2:40:58 PM (4 years ago)

Author:

cmathe

Message:

MARS: fix mtotco2, icetotco2, vaptotco2 outputs

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3400 +3400 @@
 == 09/06/2021 == EM
 Bug (wrong usage of surface pressure) fix in vdifc (bug was introduced in r2515)
+
+== 22/06/2021 == CM
+co2clouds: fix mtotco2, icetotco2, vaptotco2 outputs

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -396 +396 @@
       REAL mdusttot(ngrid)      ! Total mass of dust tracer (kg/m2)
       REAL icetot(ngrid)        ! Total mass of water ice (kg/m2)
-      REAL mtotco2       ! Total mass of co2, including ice at the surface (kg/m2)
-      REAL vaptotco2     ! Total mass of co2 vapor (kg/m2)
-      REAL icetotco2     ! Total mass of co2 ice (kg/m2)
+      REAL mtotco2(ngrid)      ! Total mass of co2, including ice at the surface (kg/m2)
+      REAL vaptotco2(ngrid)     ! Total mass of co2 vapor (kg/m2)
+      REAL icetotco2(ngrid)     ! Total mass of co2 ice (kg/m2)
       REAL Nccntot(ngrid)       ! Total number of ccn (nbr/m2)
       REAL NccnCO2tot(ngrid)    ! Total number of ccnCO2 (nbr/m2)
@@ -2509 +2509 @@
 
            endif ! of if (water)
+
+
+          if (co2clouds) then
+            mtotco2(1:ngrid) = 0.
+            icetotco2(1:ngrid) = 0.
+            vaptotco2(1:ngrid) = 0.
+            do ig=1,ngrid
+              do l=1,nlayer
+                vaptotco2(ig) = vaptotco2(ig) + 
+     &                          zq(ig,l,igcm_co2) * 
+     &                          (zplev(ig,l) - zplev(ig,l+1)) / g
+                icetotco2(ig) = icetot(ig) + 
+     &                          zq(ig,l,igcm_co2_ice) * 
+     &                          (zplev(ig,l) - zplev(ig,l+1)) / g
+              end do
+              mtotco2(ig) = icetotco2(ig) + vaptotco2(ig)
+            end do
+          end if
         endif                   ! of if (tracer)
 #ifndef MESOSCALE
@@ -2633 +2651 @@
 
              endif ! of if (water)
-             
+
+             if (co2clouds) then
+               call wstats(ngrid,"mtotco2",
+     &                    "total mass atm of co2","kg/m2",
+     &                    2,mtotco2)
+               call wstats(ngrid,"icetotco2",
+     &                    "total mass atm of co2 ice","kg/m2",
+     &                    2,icetotco2)
+               call wstats(ngrid,"vaptotco2",
+     &                    "total mass atm of co2 vapor","kg/m2",
+     &                    2,icetotco2)
+             end if
              
            if (dustbin.ne.0) then
@@ -2929 +2958 @@
      &                     'kg.kg-1', 3, zq(:,:,igcm_co2_ice))
 
+         call WRITEDIAGFI(ngrid,"mtotco2","total mass atm of co2",
+     &                    "kg.m-2",2, mtotco2)
+         call WRITEDIAGFI(ngrid,"icetotco2","total mass atm of co2 ice",
+     &                    "kg.m-2", 2, icetotco2)
+         call WRITEDIAGFI(ngrid,"vaptotco2","total mass atm of co2
+     &                    vapor","kg.m-2", 2, vaptotco2)
          call WRITEDIAGFI(ngrid,"emis","Surface emissivity","w.m-1",2,
      &                  emis)
@@ -3820 +3855 @@
      &                     '%',0,(co2totA-co2totB)/co2totA)
       endif ! of if (igcm_co2_ice.ne.0)
+
 ! XIOS outputs
 #ifdef CPP_XIOS