Changeset 3661 for trunk/LMDZ.TITAN/libf

Timestamp:

Feb 27, 2025, 2:53:11 PM (9 months ago)

Author:

emoisan

Message:

Titan CRM:
Add Titan interface in INTERFACES_V4
Adapt module_model_constants.F to Titan
Add new tracer_mode for Titan (CH4 scalar)
Add new communication of variables between LMDZ.TITAN and WRF
Allow microphysics for Mesoscale in physiq_mod.F90
EMo

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• comm_wrf.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (17 diffs)

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

@@ -18 +18 @@
   REAL,SAVE,ALLOCATABLE :: comm_FLUXSURF_LW(:)
   REAL,SAVE,ALLOCATABLE :: comm_FLXGRD(:)
+  REAL,SAVE,ALLOCATABLE :: comm_zqfi_omp(:,:,:)
+  REAL,SAVE,ALLOCATABLE :: comm_zdtlc(:,:)
 
 contains
 
-  subroutine allocate_comm_wrf(ngrid,nlayer)
+  subroutine allocate_comm_wrf(ngrid,nlayer,nq)
   implicit none
-  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: ngrid  ! number of atmospheric columns
   integer,intent(in) :: nlayer ! number of atmospheric layers
+  integer,intent(in) :: nq     ! number of tracers
   allocate(comm_HR_SW(ngrid,nlayer))
   allocate(comm_HR_LW(ngrid,nlayer))
@@ -38 +41 @@
   allocate(comm_FLUXSURF_LW(ngrid))
   allocate(comm_FLXGRD(ngrid))
+  allocate(comm_zqfi_omp(ngrid,nlayer,nq))
+  allocate(comm_zdtlc(ngrid,nlayer))
 
   end subroutine allocate_comm_wrf
@@ -56 +61 @@
   deallocate(comm_FLUXSURF_LW)
   deallocate(comm_FLXGRD)
+  deallocate(comm_zqfi_omp)
+  deallocate(comm_zdtlc)
 
   end subroutine deallocate_comm_wrf

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

@@ -46 +46 @@
 #else
 use comm_wrf, only : comm_HR_SW, comm_HR_LW, &
-                     comm_FLUXTOP_DN,comm_FLUXABS_SW,&
-                     comm_FLUXTOP_LW,comm_FLUXSURF_SW,&
-                     comm_FLUXSURF_LW,comm_FLXGRD
+                     comm_FLUXTOP_DN, comm_FLUXABS_SW,&
+                     comm_FLUXTOP_LW, comm_FLUXSURF_SW,&
+                     comm_FLUXSURF_LW, comm_FLXGRD,&
+                     comm_zqfi_omp, comm_zdtlc
 #endif
 #ifdef CPP_XIOS      
@@ -138 +139 @@
 !    pdv(ngrid,nlayer)         \  Temporal derivative of the corresponding
 !    pdt(ngrid,nlayer)         /  variables due to physical processes.
-!    pdq(ngrid,nlayer)        /
+!    pdq(ngrid,nlayer,nq)     /
 !    pdpsrf(ngrid)           /
 !
@@ -331 +332 @@
       character*2 :: str2
 
-#ifndef MESOSCALE
+!#ifndef MESOSCALE
 
 ! Local variables for Titan chemistry and microphysics
@@ -363 +364 @@
 #endif
 
-      logical file_ok
+#ifdef MESOSCALE
+      LOGICAL, SAVE ::  moyzon_ch ! used for zonal averages in Titan
+      REAL, ALLOCATABLE :: zplevbar(:,:)
+      REAL, ALLOCATABLE :: zplaybar(:,:)
+      REAL, ALLOCATABLE :: ztfibar(:,:)
+      REAL, ALLOCATABLE :: zqfibar(:,:,:)
+      REAL, ALLOCATABLE :: zphibar(:,:)
+      REAL, ALLOCATABLE :: zphisbar(:)
+      REAL, ALLOCATABLE :: zzlevbar(:,:)
+      REAL, ALLOCATABLE :: zzlaybar(:,:)
+!     REAL, DIMENSION(:,:)   :: zplevbar(:,:)
+!     REAL, DIMENSION(:,:)   :: zplaybar(:,:)
+!     REAL, DIMENSION(:,:)   :: ztfibar(:,:)
+!     REAL, DIMENSION(:,:,:) :: zqfibar(:,:,:)
+!     REAL, DIMENSION(:,:)   :: zphibar(:,:)
+!     REAL, DIMENSION(:)     :: zphisbar(:)
+!     REAL, DIMENSION(:,:)   :: zzlevbar(:,:)
+!     REAL, DIMENSION(:,:)   :: zzlaybar(:,:)
+#endif
 
 !-----------------------------------------------------------------------------
@@ -385 +404 @@
     END INTERFACE
 
-#endif
-      
+!#endif
+      
+      logical file_ok
+
 !==================================================================================================
 
@@ -447 +468 @@
            int_dtauv(:,:,:) = 0.D0
           
-#ifndef MESOSCALE 
+!#ifndef MESOSCALE 
            IF (callmufi .AND. (.NOT. uncoupl_optic_haze)) THEN
              haze_opt_file=trim(datadir)//'/optical_tables/HAZE_OPTIC_'//trim(adjustl(tmp1))//'x'//trim(adjustl(tmp2))//'.DAT'
@@ -459 +480 @@
               endif
            ENDIF
-#endif           
+!#endif           
 
          endif
+
+#ifdef MESOSCALE
+         moyzon_ch = .false. !no zonal mean for mesoscale
+#endif
 
 #ifndef MESOSCALE
 !        Initialize names and timestep for chemistry
 !        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-
          if (callchim) then
 
@@ -479 +503 @@
 
          endif
+#endif
 
 !        Initialize microphysics.
@@ -491 +516 @@
 
          ENDIF
-#endif
 
 #ifdef CPP_XIOS
@@ -717 +741 @@
       ! JVO 19 : We shall always have correct altitudes in chemistry no matter what's in physics
       ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-#ifndef MESOSCALE
+!#ifndef MESOSCALE
       if (moyzon_ch) then ! Zonal averages
          
@@ -748 +772 @@
 
       else !  if not moyzon
-#endif
+!#endif
       
         DO ig=1,ngrid
@@ -763 +787 @@
         ENDDO
 
-#ifndef MESOSCALE
+!#ifndef MESOSCALE
       endif  ! moyzon
-#endif
+!#endif
 
       ! -------------------------------------------------------------------------------------
@@ -1065 +1089 @@
       if (tracer) then
 
-#ifndef MESOSCALE
    ! -------------------
    !   V.1 Microphysics
@@ -1072 +1095 @@
          ! We must call microphysics before chemistry, for condensation !
          if (callmufi) then
-
             zzlev(:,nlayer+1)=zzlay(:,nlayer)+(zzlay(:,nlayer)-zzlev(:,nlayer)) ! JVO 19 : We assume zzlev isn't reused later on (could be done cleaner)
 
@@ -1148 +1170 @@
          endif ! callmufi
       
+#ifndef MESOSCALE
   ! -----------------
   !   V.2. Chemistry
@@ -1292 +1315 @@
          endif ! end of 'callchim'
 
-! END MESOSCALE
+!! END ifndef MESOSCALE
 #endif
 
@@ -2010 +2033 @@
       comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid)
       sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid)
+      comm_zqfi_omp(1:ngrid,1:nlayer,1:nq) = zq(1:ngrid,1:nlayer,1:nq)
+      comm_zdtlc(1:ngrid,1:nlayer) = zdtlc(1:ngrid,1:nlayer)
 #endif