Changeset 5546 for LMDZ6

Timestamp:

Feb 17, 2025, 3:05:27 PM (4 months ago)

Author:

aborella

Message:

Read airplane data

Location:

LMDZ6/branches/contrails

Files:

• DefLists/context_input_lmdz.xml (modified) (1 diff)
• libf/phylmd/limit_read_mod.f90 (modified) (2 diffs)
• libf/phylmd/lmdz_aviation.f90 (modified) (1 diff)
• libf/phylmd/physiq_mod.F90 (modified) (1 diff)

LMDZ6/branches/contrails/DefLists/context_input_lmdz.xml

@@ -391 +391 @@
 
 
+<!------------------------------------------------------------------------------------->
+                     <!--  read aviation traffic density  --> 
+<!------------------------------------------------------------------------------------->
+
+<file_definition>
+
+  <file id="aviation_file" name="aviation" mode="read" output_freq="1mo" type="one_file" time_counter_name="toto" >
+
+    <field id="KMFLOWN_id" name="KMFLOWN" operation="instant"  grid_ref="aviation_grid" freq_offset="1ts"  />
+    <field id="levaviation_id" name="LEV" axis_ref="aviation_lev" operation="instant" freq_offset="1ts" />
+    <field id="timeaviation_id" name="TIME" axis_ref="aviation_time" operation="instant" freq_offset="1ts" />
+
+  </file>
+
+</file_definition>
+
+ <domain_definition>                            
+
+  <domain id="aviation_dom" type="rectilinear"> 
+    <generate_rectilinear_domain   />           
+  </domain>                                     
+
+ </domain_definition>                           
+
+<grid_definition>
+
+    <grid id="aviation_grid">
+      <domain domain_ref="aviation_dom" /> 
+      <axis axis_ref="aviation_lev" />
+      <axis axis_ref="aviation_time"  />
+    </grid>
+
+</grid_definition>
+
+<axis_definition>
+
+  <axis id="aviation_lev" unit="hPa" axis_type="Z"/>
+  <axis id="aviation_time" name="TIME" /> 
+
+</axis_definition>
+
+<field_definition>
+
+  <field id ="KMFLOWN_read" field_ref="KMFLOWN_id"  read_access="true" />
+  <field id ="levaviation_read" field_ref="levaviation_id"  read_access="true" />
+  <field id ="timeaviation_read" field_ref="timeaviation_id"  read_access="true" />
+
+  <field id ="KMFLOWN_interp" field_ref="KMFLOWN_read"  read_access="true" grid_ref="grid_from_aviation" /> 
+
+</field_definition>
+
+<domain_definition>
+
+  <domain id="dom_from_aviation" domain_ref="dom_glo">
+      <interpolate_domain order="1" />
+  </domain>
+
+</domain_definition>
+
+<grid_definition>
+
+    <grid id="grid_from_aviation">
+      <domain domain_ref="dom_from_aviation" />
+      <axis axis_ref="klev" />
+      <axis axis_ref="aviation_time"  /> 
+    </grid>
+
+</grid_definition>
+
 
 

LMDZ6/branches/contrails/libf/phylmd/limit_read_mod.f90

@@ -41 +41 @@
     INTEGER, INTENT(IN) :: first_day
     
-    
     IF ( type_ocean /= 'couple') THEN
       IF (grid_type==unstructured) THEN
@@ -50 +49 @@
   END SUBROUTINE init_limit_read
   
+  SUBROUTINE init_aviation_read(first_day)
+  USE mod_grid_phy_lmdz
+  USE surface_data
+  USE mod_phys_lmdz_para
+  USE lmdz_xios
+
+  IMPLICIT NONE
+    INTEGER, INTENT(IN) :: first_day
+    LOGICAL :: ok_plane_contrail, ok_plane_h2o
+    
+    IF ( ok_plane_contrail .OR. ok_plane_h2o ) THEN
+      IF (grid_type==unstructured) THEN
+        IF (is_omp_master) CALL xios_set_file_attr("aviation_file", enabled=.TRUE.)
+      ENDIF  
+    ENDIF
+
+  END SUBROUTINE init_aviation_read
+
   SUBROUTINE limit_read_frac(itime, dtime, jour, pctsrf_new, is_modified)
 !

LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90

@@ -621 +621 @@
 END FUNCTION contrail_cross_section_onera
 
-
-SUBROUTINE read_aviation_emissions( &
-        klon, klev, latitude_deg, longitude_deg, pplay, &
-        flight_dist, flight_h2o &
-        )
-
-IMPLICIT NONE
-!
-! Input
-!
-INTEGER, INTENT(IN)                        :: klon, klev       ! number of horizontal grid points and vertical levels
-REAL,    INTENT(IN),  DIMENSION(klon)      :: latitude_deg     ! latitude of the grid points [deg]
-REAL,    INTENT(IN),  DIMENSION(klon)      :: longitude_deg    ! longitude of the grid points [deg]
-REAL,    INTENT(IN),  DIMENSION(klon,klev) :: pplay            ! layer pressure [Pa]
-!
-! Output
-!
-REAL,    INTENT(OUT), DIMENSION(klon,klev) :: flight_dist      ! aviation distance flown within the mesh [m/s/mesh]
-REAL,    INTENT(OUT), DIMENSION(klon,klev) :: flight_h2o       ! aviation H2O emitted within the mesh [kgH2O/s/mesh]
-!
-! Local
-!
-INTEGER :: i, k
-
-!--Initialisation
-flight_dist(:,:) = 0.
-flight_h2o(:,:) = 0.
-
-DO i=1, klon
- IF ( ( latitude_deg(i) .GE. 42. ) .AND. ( latitude_deg(i) .LE. 48. ) ) THEN 
-   flight_dist(i,15) = 50000.  !--5000 m of flight/second in grid cell x 10 scaling
-   flight_h2o(i,15) = 100.  !--10 kgH2O/second in grid cell x 10 scaling
- ENDIF
-ENDDO
+SUBROUTINE read_aviation_emissions(klon, klev, flight_dist, flight_h2o)
+    ! This subroutine allows to read the traffic density data read in the file aviation.nc
+    ! This file is defined in ./COMP/lmdz.card
+    ! Field names in context_input_lmdz.xml should be the same as in the file.
+    USE netcdf
+    USE mod_phys_lmdz_para
+    USE mod_grid_phy_lmdz, ONLY: nbp_lon_=>nbp_lon, nbp_lat_=>nbp_lat, nbp_lev_=>nbp_lev, &
+                                 klon_glo, grid2Dto1D_glo, grid_type, unstructured
+    USE lmdz_xios
+    USE print_control_mod, ONLY: lunout
+    !USE dimphy
+    !USE phys_state_var_mod
+    IMPLICIT NONE
+
+    REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_dist ! Aviation distance flown within the mesh [m/s/mesh]
+    REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_h2o  ! Aviation H2O emitted within the mesh [kgH2O/s/mesh]
+
+    INTEGER, INTENT(IN)                        :: klon, klev       ! number of horizontal grid points and vertical levels
+
+    !----------------------------------------------------
+    ! Local variable
+    !----------------------------------------------------
+    INTEGER :: i, error, varid, iret
+    CHARACTER(len=*), PARAMETER :: filename='aviation.nc'
+    REAL, ALLOCATABLE :: flight_dist_mpi(:,:,:)
+
+    !--Initialisation
+    flight_dist(:,:) = 0.
+    flight_h2o(:,:) = 0.
+
+    ! is_omp_master is necessary to make XIOS works
+    IF (is_omp_master) THEN
+
+        ! Allocation of KMFLOWN for the local MPI process
+        ALLOCATE(flight_dist_mpi(klon_mpi,klev,1), STAT=error)
+        IF (error /= 0) CALL abort_physic('airplane_mod.F90', 'pb in allocation pointer flight_dist_mpi',1)
+
+        ! Read the data from the file
+        CALL xios_recv_field("KMFLOWN_interp", flight_dist_mpi(:,:,1))
+   ENDIF !- is_omp_master
+
+   ! Propagate to other OMP threads: flight_dist_mpi(klon_mpi,klev) to flight_dist(klon,klev)
+   ! (klon_mpi,klon) = (200,50) avec 80 MPI, 4 OMP, nbp40
+   CALL scatter_omp(flight_dist_mpi(:,:,1), flight_dist)
 
 END SUBROUTINE read_aviation_emissions
+
+!SUBROUTINE read_aviation_emissions( &
+!        klon, klev, latitude_deg, longitude_deg, pplay, &
+!        flight_dist, flight_h2o &
+!        )
+!
+!IMPLICIT NONE
+!!
+!! Input
+!!
+!INTEGER, INTENT(IN)                        :: klon, klev       ! number of horizontal grid points and vertical levels
+!REAL,    INTENT(IN),  DIMENSION(klon)      :: latitude_deg     ! latitude of the grid points [deg]
+!REAL,    INTENT(IN),  DIMENSION(klon)      :: longitude_deg    ! longitude of the grid points [deg]
+!REAL,    INTENT(IN),  DIMENSION(klon,klev) :: pplay            ! layer pressure [Pa]
+!!
+!! Output
+!!
+!REAL,    INTENT(OUT), DIMENSION(klon,klev) :: flight_dist      ! aviation distance flown within the mesh [m/s/mesh]
+!REAL,    INTENT(OUT), DIMENSION(klon,klev) :: flight_h2o       ! aviation H2O emitted within the mesh [kgH2O/s/mesh]
+!!
+!! Local
+!!
+!INTEGER :: i, k
+!
+!!--Initialisation
+!flight_dist(:,:) = 0.
+!flight_h2o(:,:) = 0.
+!
+!DO i=1, klon
+! IF ( ( latitude_deg(i) .GE. 42. ) .AND. ( latitude_deg(i) .LE. 48. ) ) THEN 
+!   flight_dist(i,15) = 50000.  !--5000 m of flight/second in grid cell x 10 scaling
+!   flight_h2o(i,15) = 100.  !--10 kgH2O/second in grid cell x 10 scaling
+! ENDIF
+!ENDDO
+!
+!END SUBROUTINE read_aviation_emissions
 
 END MODULE lmdz_aviation

LMDZ6/branches/contrails/libf/phylmd/physiq_mod.F90

@@ -3867 +3867 @@
     !--Read the aviation emissions
     IF ( ok_plane_h2o .OR. ok_plane_contrail ) THEN
-      CALL read_aviation_emissions(klon, klev, latitude_deg, longitude_deg, pplay, &
-                                   flight_dist, flight_h2o)
+      CALL read_aviation_emissions(klon, klev, flight_dist, flight_h2o)
     ENDIF