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

Timestamp:

Jun 11, 2020, 9:35:53 PM (5 years ago)

Author:

jvatant

Message:

Titan GCM:Fix missing MESOSCALE cpp key in calchim

File:

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

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

@@ -67 +67 @@
   USE comcstfi_mod, ONLY: g, rad, pi, r, kbol
   USE geometry_mod, ONLY: latitude
+#ifndef MESOSCALE  
   USE logic_mod, ONLY: moyzon_ch
   USE moyzon_mod, ONLY: tmoy, playmoy
+#endif
 
   IMPLICIT NONE
@@ -223 +225 @@
      PRINT*,'Init chemistry : pressure, density, temperature ... :'
      PRINT*,'level, press_c (mbar), nb (cm-3), temp_c (K)'
-     
+
+#ifndef MESOSCALE
+!! JVO 20 : I put this CPP key because mesoscale cannot access tmoy and playmoy,
+!  you should fix this when you will propeerly use chemistry !
+
      IF (ngrid.NE.1) THEN
        DO l=1,klev
@@ -232 +238 @@
        ENDDO
      ELSE
+#endif
        DO l=1,klev
           temp_c(l)  = ctemp(1,l)                             ! K 
@@ -238 +245 @@
           PRINT*, l, press_c(l), nb(l), temp_c(l)
        ENDDO
+#ifndef MESOSCALE
      ENDIF
+#endif
 
      ! b. Extension in upper atmosphere with Vervack profile
@@ -349 +358 @@
      ! -------------------------------------------------------------------------------------------------------
 
+#ifndef MESOSCALE
      IF ( ( moyzon_ch .AND. ( ig.EQ.1 .OR. (ABS(latitude(ig)-latitude(igm1)).GT.0.1*pi/180.0)) ) .OR. (.NOT. moyzon_ch) ) THEN
-
+#endif
         ! 1. Compute altitude for the grid point with hydrostat. equilib.
         ! ---------------------------------------------------------------
@@ -365 +375 @@
            nb(l)      = 1.e-4*press_c(l) / (kbol*temp_c(l))  ! cm-3
         ENDDO
-        rinter(klev+1)=(czlev(ig,klev+1)+rad)/1000.
+        rinter(klev+1)=( (czlay(ig,klev) + ( czlay(ig,klev) - czlev(ig,klev))) +rad )/1000. ! You shall not use czlev(zlev+1) whose value is 1.e7 !
 
         ! b. Extension in upper atmosphere with Vervack profile
@@ -524 +534 @@
         ! NB: The full vertical composition grid will be created only for the outputs
 
-
+#ifndef MESOSCALE
      ELSE ! In 2D chemistry, if following grid point at same latitude, same zonal mean so don't do calculations again !
         dqyc(ig,:,:)    = dqyc(igm1,:,:) ! will be put back in 3D with longitudinal variations assuming same relative tendencies within a lat band
         ykim_up(:,ig,:) = ykim_up(:,igm1,:) ! no horizontal mixing in upper layers -> no longitudinal variations
      ENDIF
+#endif
 
   ENDDO