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

Timestamp:

Dec 17, 2018, 12:34:13 PM (6 years ago)

Author:

jvatant

Message:

Slight adjustements of the eddy diffusion profile for chemistry.

File:

• trunk/LMDZ.TITAN/libf/phytitan/calchim.F90 (modified) (3 diffs)

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

@@ -264 +264 @@
      ! 6. Eddy mixing coefficients (constant with time and space)
      ! ----------------------------------------------------------
-
+     
      kedd(:) = 1.e3 ! Default value =/= zero
 
@@ -271 +271 @@
 
      ! First calculate kedd for upper chemistry layers
-     DO l=klev+1,nlaykim_tot
+     DO l=klev-4,nlaykim_tot
         logp=-log10(press_c(l))
-     ! 1E6 at 300 km ~ 10-1 mbar
-        IF     ( logp.ge.1.0 .and. logp.le.4.0 ) THEN 
-              kedd(l) = 10.**(6.0+1.3*(logp-1.0)/3.0)
-     ! 2E7 above 600 km ~ 10-4 mbar
+     ! 2E6 at 400 km ~ 10-2 mbar
+        IF     ( logp.ge.2.0 .and. logp.le.3.0 ) THEN 
+              kedd(l) = 2.e6 * 5.0**(logp-2.0)
+     ! 1E7 at 500 km ~ 10-3 mbar
+        ELSE IF     ( logp.ge.3.0 .and. logp.le.4.0 ) THEN 
+              kedd(l) = 1.e7 * 3.0**(logp-3.0)
+     ! 3E7 above 700 km ~ 10-4 mbar
         ELSEIF ( logp.gt.4.0                   ) THEN 
-             kedd(l) = 2.e7
+             kedd(l) = 3.e7
         ENDIF
      ENDDO
      
-     ! Then adjust 15 last layers profile fading to default value depending on kedd(ptop)
-     DO l=klev-15,klev
-        temp1   = ( log10(press_c(l)/press_c(klev-15)) ) / ( log10(press_c(klev+1)/press_c(klev-15)) )
-        kedd(l) = 10.**( 3.0 + log10(kedd(klev+1)/1.e3) * temp1 )
+     ! JVO 18 : We keep the nominal profile in the GCM 5 upper layers 
+     !          to have  a correct vertical mixing in the sponge layer
+     
+     ! Then adjust 10 layers profile fading to default value depending on kedd(ptop)
+     DO l=klev-15,klev-5
+        temp1   = ( log10(press_c(l)/press_c(klev-15)) ) / ( log10(press_c(klev-4)/press_c(klev-15)) )
+        kedd(l) = 10.**( 3.0 + log10(kedd(klev-4)/1.e3) * temp1 )
      ENDDO
-     
+    
      firstcall = .FALSE.
   ENDIF  ! firstcall
@@ -311 +317 @@
      ! -------------------------------------------------------------------------------------------------------
 
-     IF ( ( moyzon_ch .AND. ( ig.EQ.1 .OR. (ABS(latitude(ig)-latitude(igm1)).GT.0.1*180.0/pi)) ) .OR. (.NOT. moyzon_ch) ) THEN
+     IF ( ( moyzon_ch .AND. ( ig.EQ.1 .OR. (ABS(latitude(ig)-latitude(igm1)).GT.0.1*pi/180)) ) .OR. (.NOT. moyzon_ch) ) THEN
 
         ! 1. Compute altitude for the grid point with hydrostat. equilib.