Changeset 616 for trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

Timestamp:

Apr 5, 2012, 3:43:29 PM (13 years ago)

Author:

acolaitis

Message:

Forgot to mask some diagnostic outputs from thermals in previous commit

File:

• trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

@@ -754 +754 @@
       do ig=1,ngridmx
             if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
+      IF (lwrite) THEN
              print*,'On tombe sur le cas particulier de thermcell_plume'
+      ENDIF
                 zw2(ig,l+1)=0.
                 linter(ig)=l+1
@@ -1029 +1031 @@
                ncorecfm2=ncorecfm2+1
                else
+               IF (lwrite) THEN
           print*,'fm(l+1)<0 : ig, l+1,lmax :',ig,l+1,lmax(ig),fm(ig,l+1)
+               ENDIF
                ncorecfm1=ncorecfm1+1
                fm(ig,l+1)=fm(ig,l)
@@ -1229 +1233 @@
      &   fm(ig,k)*ztv(ig,k)+fm(ig,k)*ztva(ig,k)-fm(ig,k+1)*ztva(ig,k+1))
             if (ztv(ig,k) + ptimestep*zdthladj(ig,k) .le. 0.) then
+      IF (lwrite) THEN
               print*,'Teta<0 in thermcell_dTeta up: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj(ig,k)
+      ENDIF
               if(ztv(ig,k) .gt. 0.) then
               zdthladj(ig,k)=0.
@@ -1315 +1321 @@
      & fm_down(ig,k)*ztv(ig,k)+fm_down(ig,k)*ztvd(ig,k)-fm_down(ig,k+1)*ztvd(ig,k+1))
             if (ztv(ig,k) + ptimestep*zdthladj_down(ig,k) .le. 0.) then
+      IF (lwrite) THEN
               print*,'q<0 in thermcell_dTeta down: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj_down(ig,k)
+      ENDIF
               if(ztv(ig,k) .gt. 0.) then
               zdthladj(ig,k)=0.