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

Timestamp:

Apr 13, 2012, 3:15:09 PM (13 years ago)

Author:

acolaitis

Message:

THERMALS: corrected improper tracer conservation that would appear in sharp tracer gradient cases. Tracer conservation is now identical to convective adjustment

File:

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

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

@@ -271 +271 @@
 
       do l=2,nlayermx
-         rhobarz(:,l)=0.5*(rho(:,l)+rho(:,l-1))
+!         rhobarz(:,l)=0.5*(rho(:,l)+rho(:,l-1))
+          rhobarz(:,l)=pplev(:,l)/(r*0.5*(pt(:,l)+pt(:,l-1)))
       enddo
 
@@ -1533 +1534 @@
       modname='tracer'
       call thermcell_dqup(ngridmx,nlayermx,ptimestep     &
-     &     ,fm,entr,detr,  &
+     &     ,fm,entr,fraca,  &
      &    masse,pq(:,:,ico2),pdqadj(:,:,ico2),modname,zdzfull)
 !      endif