Changeset 1904 for trunk/LMDZ.TITAN/libf/dynphy_lonlat

Timestamp:

Feb 25, 2018, 7:44:55 PM (7 years ago)

Author:

jvatant

Message:

Fixed 2 critical bugs.
+ One about rat_mmol (* instead of /)
+ One about zonal means (or not) arrays sent in chemistry and mufi
In any case zonal means seems to be responsible for crashes of chem and mufi !
Stop using them until further informations (for mufi, no problem, full 3D is quick)
--JVO

File:

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

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

@@ -77 +77 @@
       ENDDO
       ! mass -> molar mixing ratio to be comparable to ykim_up later
-      avg_qtop(ichem,:)=avg_qtop(ichem,:)*rat_mmol(chimi_indx(ichem)) 
+      avg_qtop(ichem,:)=avg_qtop(ichem,:)/rat_mmol(chimi_indx(ichem)) 
     ENDDO
+
+  avg_qtop(:,:) = avg_qtop(:,:) / real(iim)
 
   ENDIF
@@ -150 +152 @@
         ! 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)/rat_mmol(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)) &
+        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))
 
@@ -161 +163 @@
             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)) / rat_mmol(chimi_indx(ichem))
+              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))
@@ -168 +170 @@
             ! Periodicity on longitude at 180 and -180
 
-            ykimlon = 0.5*(ykim_up(ichem,ng0+1,1)+ykim_up(ichem,ng0+iim,1)) / rat_mmol(chimi_indx(ichem))
+            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))