Changeset 1787 for trunk

Timestamp:

Sep 25, 2017, 4:02:05 PM (7 years ago)

Author:

jvatant

Message:

Fixed a bug in calchim and physiq_mod on conversion of tendencies to 3D and molar ratio / mass ratio

-> we assume same 2D relative tendency on all lat in calchim and adjust longitudinal variations

--JVO

Location:

trunk/LMDZ.TITAN

Files:

• README (modified) (1 diff)
• libf/phytitan/calchim.F90 (modified) (2 diffs)
• libf/phytitan/physiq_mod.F90 (modified) (5 diffs)

trunk/LMDZ.TITAN/README

@@ -1322 +1322 @@
  - Added a tmpk variable
  - Save ~ 50% time on the RT, ~30% on the whole code on the tested simulations
+
+== 25/09/2017 == JVO
+Fixed a bug in calchim and physiq_mod on conversion of tendencies to 3D and molar ratio / mass ratio
+ ->  we assume same 2D relative tendency on all lat in calchim and adjust longitudinal variations

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

@@ -277 +277 @@
 !                BOUCLE SUR LES LATITUDES
 !
+!                * Permet de faire le calcul une seule fois par lat
+!
       DO j=1,nlon 
       
@@ -492 +494 @@
 !              FIN: BOUCLE SUR LES LATITUDES 
 
-      else      ! same latitude, we don't do calculations again
-        dqyc(j,:,:) = dqyc(jm1,:,:)
+      else      ! same latitude, we don't do calculations again, only adjust longitudinal variations
+        dqyc(j,:,:) = dqyc(jm1,:,:)/qy_c(jm1,:,:)*qy_c(j,:,:)
       endif
 

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

@@ -380 +380 @@
       
       real ychim(ngrid,nlayer,nq-nmicro)   
+
+      real rat_mmol(nq) ! Molar fraction ratio
 
       ! 2D vmr tendencies ( chemistry and condensation )
@@ -470 +472 @@
          endif
 
+         rat_mmol(:) = mmol(:) / mugaz
 
 !        Read 'startfi.nc' file.
@@ -1005 +1008 @@
          if (callchim) then
 
+            ! Utilisation de la moyenne zonale dans calchim
             zplev(:,:) = zplevbar(:,:)
             zplay(:,:) = zplaybar(:,:)
@@ -1013 +1017 @@
             if (nq.gt.nmicro) then
                do iq = nmicro+1,nq
-                  ychim(:,:,iq-nmicro) = pq(:,:,iq) 
+                  ychim(:,:,iq-nmicro) = pq(:,:,iq) * rat_mmol(iq) ! convert to molar fraction
                enddo
             endif
@@ -1040 +1044 @@
             
             if (nq.gt.nmicro) then
-               ! We convert tendencies back to 3D and mass mixing ratio
+               ! We convert tendencies back to mass mixing ratio
                do iq=nmicro+1,nq
-                  zdqchi(:,:,iq) = dycchi(:,:,iq-nmicro)*pq(:,:,iq) / ychim(:,:,iq-nmicro)
-                  zdqmph(:,:,iq) = dycmph(:,:,iq-nmicro)*pq(:,:,iq) / ychim(:,:,iq-nmicro)
+                  zdqchi(:,:,iq) = dycchi(:,:,iq-nmicro) / rat_mmol(iq)
+                  zdqmph(:,:,iq) = dycmph(:,:,iq-nmicro) / rat_mmol(iq)
                enddo