Changeset 1899 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/vert_regrid_kim.F90

Timestamp:

Jan 26, 2018, 3:02:55 PM (7 years ago)

Author:

jvatant

Message:

Fixing an incoherence - the upper chemistry fields
are in molar mixing ratio not mass mixing ratio
--JVO

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/vert_regrid_kim.F90 (modified) (8 diffs)

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/vert_regrid_kim.F90

@@ -16 +16 @@
   ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-  USE comchem_h, ONLY: cnames, nlaykim_up, preskim, ykim_up
+  USE comchem_h, ONLY: cnames, rat_mmol, nlaykim_up, preskim, ykim_up
   USE comchem_newstart_h
   USE tracer_h
@@ -34 +34 @@
   REAL,DIMENSION(iim+1,jjm+1,llm,nq), INTENT(INOUT) :: q   ! Advected fields (kg/kg) on 3D dyn. grid
 
-  REAL, DIMENSION(:,:), ALLOCATABLE   :: avg_qtop  ! Zonally averaged q on top layer
+  REAL, DIMENSION(:,:), ALLOCATABLE   :: avg_qtop  ! Zonally averaged q (mol/mol) on top layer
 
   REAL :: coef, ykimlon
@@ -52 +52 @@
   ngridmx = size(ykim_up,DIM=2)
 
-  ! -----------------------------------------------------------------------------------------
-  ! 1. Compute zonal mean of last layer for every chem and lat and convert to physics grid 
+  ! ------------------------------------------------------------
+  ! 1. Compute zonal mean of last layer for every chem and lat
+  ! and convert it to molar mixing fraction and to physics grid 
   ! Preliminary, only in case ceiling has been lowered
-  ! -----------------------------------------------------------------------------------------
+  ! ------------------------------------------------------------
 
   lowered = .FALSE.
@@ -75 +76 @@
         ENDDO
       ENDDO
+      ! mass -> molar mixing ratio to be comparable to ykim_up later
+      avg_qtop(ichem,:)=avg_qtop(ichem,:)*rat_mmol(chimi_indx(ichem)) 
     ENDDO
 
@@ -125 +128 @@
   ! 3. Correct the 3D advected chem. tracer fields if model ceiling is highered
   ! In this case we interpolate between tracers below and upper_chemistry values
+  ! Doing this we convert via rat_mmol ykim_up from molar to mass mixing ratio
   ! ----------------------------------------------------------------------------
  
@@ -146 +150 @@
         ! so we deal with mono-gridpoints for North and South Poles   
 
-        q(:,1,ilay,chimi_indx(ichem)) = (1.0-coef)*ykim_up(ichem,1,1) + coef*q(:,1,isup,chimi_indx(ichem))
-        q(:,jjm+1,ilay,chimi_indx(ichem)) = (1.0-coef)*ykim_up(ichem,ngridmx,1) + coef*q(:,jjm+1,isup,chimi_indx(ichem))
+        q(:,1,ilay,chimi_indx(ichem)) = (1.0-coef)*ykim_up(ichem,1,1)/rat_mmol(chimi_indx(ichem)) &
+                                      + coef*q(:,1,isup,chimi_indx(ichem))
+        q(:,jjm+1,ilay,chimi_indx(ichem)) = (1.0-coef)*ykim_up(ichem,ngridmx,1)/rat_mmol(chimi_indx(ichem)) &
+                                          + coef*q(:,jjm+1,isup,chimi_indx(ichem))
 
           DO ilat=2,jjm
@@ -155 +161 @@
             DO ilon=2,iim
               ! ykim_up and q are shifted one to the other on longitudinal grid
-              ykimlon = 0.5*(ykim_up(ichem,ng0+ilon-1,1)+ykim_up(ichem,ng0+ilon,1)) 
+              ykimlon = 0.5*(ykim_up(ichem,ng0+ilon-1,1)+ykim_up(ichem,ng0+ilon,1)) / rat_mmol(chimi_indx(ichem))
 
               q(ilon,ilat,ilay,chimi_indx(ichem)) = (1.0-coef)*ykimlon + coef*q(ilon,ilat,isup,chimi_indx(ichem))
@@ -162 +168 @@
             ! Periodicity on longitude at 180 and -180
 
-            ykimlon = 0.5*(ykim_up(ichem,ng0+1,1)+ykim_up(ichem,ng0+iim,1)) 
+            ykimlon = 0.5*(ykim_up(ichem,ng0+1,1)+ykim_up(ichem,ng0+iim,1)) / rat_mmol(chimi_indx(ichem))
 
             q(1,ilat,ilay,chimi_indx(ichem)) = (1.0-coef)*ykimlon + coef*q(1,ilat,isup,chimi_indx(ichem))