Changeset 2385 for trunk/LMDZ.MARS/libf

Timestamp:

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

Author:

cmathe

Message:

physiq_mod: writediagfi 1 => 0

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (1 diff)

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

@@ -3721 +3721 @@
 
         ENDIF ! hdo
-      write(*,*)'ok'
-      if (co2clouds) then
-       write(*,*)'i am here'
-        call WRITEDIAGFI(ngrid,'mtotco2',
-     &                   'total mass of CO2, including ice surf',
-     &                   'kg/m2',1,mtotco2)
-
-        call WRITEDIAGFI(ngrid,'vaptotco2',
-     &                   'total mass of CO2 vapor',
-     &                   'kg/m2',1,vaptotco2)
-
-        call WRITEDIAGFI(ngrid,'icetotco2',
-     &                   'total mass of CO2 ice',
-     &                   'kg/m2',1,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',
-     &                   '%',1,(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,