Changeset 1966

Timestamp:

Jul 6, 2018, 12:12:47 PM (6 years ago)

Author:

jvatant

Message:

Improve computation of saturation profiles that was done in planetary average
-> Now done for every grid point ( enables sensitivity to winter low temps )
-> firstcall "inicondens.F90" replaced by every-step-call-to "calc_ysat.F90"
--JVO

Location:

trunk/LMDZ.TITAN

Files:

• README (modified) (1 diff)
• libf/phytitan/calc_ysat.F90 (added)
• libf/phytitan/inicondens.F90 (deleted)
• libf/phytitan/physiq_mod.F90 (modified) (6 diffs)

trunk/LMDZ.TITAN/README

@@ -1488 +1488 @@
 Replace Numerical Recipes "solve.c" in C chemistry by "solve_lapack.c" call to Lapack dgesv
 -> ~60 times quicker ( quasi-instantaneous now )
+
+== 06/07/18 == JVO
+Improve computation of saturation profiles that was done in planetary average
+-> Now done for every grid point ( enables sensitivity to winter low temps )
+-> firstcall "inicondens.F90" replaced by every-step-call-to "calc_ysat.F90"

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

@@ -35 +35 @@
       use time_phylmdz_mod, only: daysec
       use logic_mod, only: moyzon_ch
-      use moyzon_mod, only: tmoy, playmoy, zphibar, zphisbar, zplevbar, &
-                            zplaybar, zzlevbar, zzlaybar, ztfibar, zqfibar
+      use moyzon_mod, only:  zphibar, zphisbar, zplevbar, zplaybar, &
+                             zzlevbar, zzlaybar, ztfibar, zqfibar
       use callkeys_mod
       use vertical_layers_mod, only: presnivs, pseudoalt
@@ -383 +383 @@
 
       ! Saturation profiles
-      real, dimension(:,:), allocatable, save     :: qysat ! (mol/mol)
-!$OMP THREADPRIVATE(qysat)
-      real temp_eq(nlayer), press_eq(nlayer) ! Planetary averages for the init. of saturation profiles (K,mbar)
+      real, dimension(ngrid,nlayer,nkim)          :: ysat ! (mol/mol)
 
       ! Surface methane
-      real, dimension(:), allocatable, save :: tankCH4    ! Depth of surface methane tank (m)
+      real, dimension(:), allocatable, save       :: tankCH4    ! Depth of surface methane tank (m)
 !$OMP THREADPRIVATE(tankCH4)
 
@@ -510 +508 @@
          endif
 
-!        Initialize names, timestep and saturation profiles for chemistry
-!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+!        Initialize names and timestep for chemistry
+!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
          if ( callchim ) then
@@ -522 +520 @@
 
             allocate(dycchi(ngrid,nlayer,nkim)) ! only for chemical tracers
-            allocate(qysat(nlayer,nkim))
             
             ! Chemistry timestep
             ctimestep = ptimestep*REAL(ichim)
 
-            ! qysat is taken at the equator ( small variations of t,p )
-            if (ngrid.ne.1) then ! TODO : a patcher (lecture d'un profil?) si jamais on a plus acces a moyzon_mod !
-              temp_eq(:)  = tmoy(:)
-              press_eq(:) = playmoy(:)/100. ! in mbar
-            else
-              temp_eq(:)  = pt(1,:)
-              press_eq(:) = pplay(1,:)/100.
-            endif
-            
-            call inicondens(press_eq,temp_eq,qysat)
-          
             zdqchi(:,:,:)  = 0.0
 
@@ -1146 +1132 @@
             ! i. Condensation of the 3D tracers after the transport
             ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+            
+            call calc_ysat(ngrid,nlayer,pplay,pt,ysat) ! Compute saturation profiles for every grid point
+
+            dyccond(:,:,:) = 0.0 ! Default value -> no condensation
+
             do iq=1,nkim
-               do l=1,nlayer
-                  do ig=1,ngrid
-                     if ( ychim(ig,l,iq).gt.qysat(l,iq) ) then
-                        dyccond(ig,l,iq+nmicro) = ( -ychim(ig,l,iq)+qysat(l,iq) ) / ptimestep
-                     else
-                        dyccond(ig,l,iq+nmicro) = 0.0
-                     endif
-                  enddo
-               enddo
+               where ( ychim(:,:,iq).gt.ysat(:,:,iq) )   &
+                     dyccond(:,:,iq+nmicro) = ( -ychim(:,:,iq)+ysat(:,:,iq) ) / ptimestep
             enddo
 
@@ -1171 +1155 @@
             if ( moyzon_ch .and. mod(icount-1,ichim).eq. 0 ) then
 
+              call calc_ysat(ngrid,nlayer,zzlaybar,ztfibar,ysat) ! Compute saturation profiles for every grid point for the zon-ave fields
+
+              dyccondbar(:,:,:) = 0.0 ! Default value -> no condensation
+              
               do iq = 1,nkim
-                do l=1,nlayer
-                  do ig=1,ngrid
-                    if ( ychimbar(ig,l,iq).gt.qysat(l,iq) ) then
-                       dyccondbar(ig,l,iq+nmicro) = ( -ychimbar(ig,l,iq)+qysat(l,iq) ) / ptimestep
-                    else
-                       dyccondbar(ig,l,iq+nmicro) = 0.0
-                    endif
-                  enddo
-                enddo
+                 where ( ychimbar(:,:,iq).gt.ysat(:,:,iq) )  &
+                       dyccondbar(:,:,iq+nmicro) = ( -ychimbar(:,:,iq)+ysat(:,:,iq) ) / ptimestep
               enddo