Changeset 2109 for trunk/LMDZ.TITAN/libf/phytitan

Timestamp:

Mar 4, 2019, 4:03:08 PM (6 years ago)

Author:

jvatant

Message:

Fix some problems for the microphysics :
+ Altitude of the last level at 1e7m from physics was certainly source of divergence
+ Sanity check for negative is moved from within mm_microphysic to the end of calmufi avoiding rounding pbs
--JVO

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• calmufi.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/calmufi.F90

@@ -31 +31 @@
   REAL(kind=8), DIMENSION(:,:), INTENT(IN) :: temp  !! Temperature at the center of each layer (K).
 
-  REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: pq    !! Tracers (\(kg.kg^{-1}}\)).
-  REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: zdqfi !! Tendency from former processes for tracers (\(kg.kg^{-1}}\)).
-  REAL(kind=8), DIMENSION(:,:,:), INTENT(OUT) :: zdq   !! Microphysical tendency for tracers (\(kg.kg^{-1}}\)).
+  REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: pq    !! Tracers (\(X.kg^{-1}}\)).
+  REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: zdqfi !! Tendency from former processes for tracers (\(X.kg^{-1}}\)).
+  REAL(kind=8), DIMENSION(:,:,:), INTENT(OUT) :: zdq   !! Microphysical tendency for tracers (\(X.kg^{-1}}\)).
   
-  REAL(kind=8), DIMENSION(:,:,:), ALLOCATABLE :: zq !! Local tracers updated from former processes (\(kg.kg^{-1}}\)).
+  REAL(kind=8), DIMENSION(:,:,:), ALLOCATABLE :: zq !! Local tracers updated from former processes (\(X.kg^{-1}}\)).
   
   REAL(kind=8), DIMENSION(:), ALLOCATABLE :: m0as !! 0th order moment of the spherical mode (\(m^{-2}\)).
@@ -56 +56 @@
   REAL(kind=8), DIMENSION(:,:), ALLOCATABLE :: dgazs !! Tendencies of each condensible gaz species !(\(mol.mol^{-1}\)). 
 
-  REAL(kind=8), DIMENSION(:,:), ALLOCATABLE ::  int2ext
+  REAL(kind=8), DIMENSION(:,:), ALLOCATABLE ::  int2ext !! (\(m^{-2}\)).
   TYPE(error) :: err
 
@@ -165 +165 @@
       enddo
     endif
+
+    ! Sanity check ( way safer to be done here rather than within YAMMS )
+    WHERE (zq+zdq < 0.0) ; zdq = -zq ; END WHERE
+
   END DO ! loop on ilon
 

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -1082 +1082 @@
  
          if (callmufi) then
+            zzlev(:,nlayer+1)=zzlay(:,nlayer)+(zzlay(:,nlayer)-zzlev(:,nlayer)) ! JVO 19 : We assume zzlev isn't reused later on
 #ifdef USE_QTEST
             dtpq(:,:,:) = 0.0 ! we want tpq to go only through mufi
@@ -1087 +1088 @@
             tpq(:,:,:) = tpq(:,:,:) + zdqmufi(:,:,:)*ptimestep ! only manipulation of tpq->*ptimestep here
 #else
-            call calmufi(ptimestep,pplev,zzlev,pplay,zzlay,gzlat,pt,pq,pdq,zdqmufi) ! JVO 19 : To be fixed, what altitude do we need ?
+            call calmufi(ptimestep,pplev,zzlev,pplay,zzlay,gzlat,pt,pq,pdq,zdqmufi)
             pdq(:,:,:) = pdq(:,:,:) + zdqmufi(:,:,:)
 #endif