Changeset 278 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Sep 5, 2011, 5:52:58 PM (13 years ago)

Author:

aslmd

Message:

LMDZ.MARS and MESOSCALE: minor modification when printing diagnostics for (slightly) negative tracers. this is now only done if it is the case for the local quantity computed by thermals routines.

File:

• trunk/LMDZ.MARS/libf/phymars/thermcell_dqupdown.F90 (modified) (1 diff)

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

@@ -210 +210 @@
 
             if ((qa(ig,k).lt.0.) .and. (charvar .ne. 'momentum')) then
-                 print*,'qa<0!!!',charvar,ig,k,fm0(ig,k),qa(ig,k-1),entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
-                 print*,'---------> Cancelling qa'
-                 qa(ig,k)=q(ig,k)
+                 if ((q(ig,k) .gt. 0.) .and. ( q(ig,k) .gt. 1.e-15 )) then
+                   !! only print if it is the thermal scheme which makes qa<0
+                   print*,'qa<0 created by thermals!!!',charvar,ig,k,fm0(ig,k),qa(ig,k-1),entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
+                   print*,'---------> Cancelling qa'
+                 endif
+                 qa(ig,k)=q(ig,k)  !! no action of thermals in this case.
             endif
             else