Changeset 2383 for trunk

Timestamp:

Jul 2, 2020, 3:42:24 PM (4 years ago)

Author:

cmathe

Message:

physiq_mod: line 658 misleading => no

Try another place

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (2 diffs)

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

@@ -656 +656 @@
          endif         
 
-         if(callnlte.and.nltemodel.eq.2) then
-            call nlte_setup
-         end if
-         if(callnirco2.and.nircorr.eq.1) then
-          call NIR_leedat
-         end if
+         if(callnlte.and.nltemodel.eq.2) call nlte_setup
+         if(callnirco2.and.nircorr.eq.1) call NIR_leedat
+
 
         IF (tracer.AND.water.AND.(ngrid.NE.1)) THEN
@@ -3725 +3722 @@
         ENDIF ! hdo
 
+      if (co2clouds) then
+        call WRITEDIAGFI(ngrid,'mtotco2',
+     &                   'total mass of CO2, including ice surf',
+     &                   'kg/m2',0,mtotco2)
+
+        call WRITEDIAGFI(ngrid,'vaptotco2',
+     &                   'total mass of CO2 vapor',
+     &                   'kg/m2',0,vaptotco2)
+
+        call WRITEDIAGFI(ngrid,'icetotco2',
+     &                   'total mass of CO2 ice',
+     &                   'kg/m2',0,icetotco2)
+      end if
+      co2totB = 0. ! added by C.M.
+      do ig=1,ngrid
+        do l=1,nlayer
+          co2totB = co2totB + (zplev(ig,l)-zplev(ig,l+1))/g*
+     &           (pq(ig,l,igcm_co2)+pq(ig,l,igcm_co2_ice)
+     &      +(pdq(ig,l,igcm_co2)+pdq(ig,l,igcm_co2_ice))*ptimestep)
+        enddo
+        co2totB = co2totB + co2ice(ig)
+      enddo
+
+      call WRITEDIAGFI(ngrid,'co2conservation',
+     &                   'Total CO2 mass conservation in physic',
+     &                   '%',0,(co2totA-co2totB)/co2totA)
+      call test_vmr(pq(:,:,igcm_co2), pdq(:,:,igcm_co2), ptimestep,
+     &               ngrid, nlayer, 'end physiq_mod')
+
         call WRITEDIAGFI(ngrid,"rice","ice radius","m",1,
      &                    rice)