Changeset 5637 for LMDZ6/branches/contrails/libf

Timestamp:

Apr 30, 2025, 9:59:06 AM (8 weeks ago)

Author:

aborella

Message:

Fix for activating aviation emissions data file

Location:

LMDZ6/branches/contrails/libf/phylmd

Files:

• lmdz_aviation.f90 (modified) (5 diffs)
• physiq_mod.F90 (modified) (2 diffs)

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

@@ -248 +248 @@
 
 
+SUBROUTINE init_read_aviation_emissions
+
+  USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured
+  USE mod_phys_lmdz_para, ONLY: is_omp_master
+  USE lmdz_xios, ONLY: xios_set_file_attr, xios_set_field_attr
+
+  IF (grid_type==unstructured) THEN
+    IF (is_omp_master) THEN
+      ! Activate aviation file
+      CALL xios_set_file_attr("aviation_file", enabled=.TRUE.)
+      ! Activate aviation fields
+      CALL xios_set_field_attr("KMFLOWN_read", enabled=.TRUE.)
+      CALL xios_set_field_attr("KMFLOWN_interp", enabled=.TRUE.)
+    ENDIF
+  ENDIF
+
+END SUBROUTINE init_read_aviation_emissions
+
+
 SUBROUTINE read_aviation_emissions(klon, klev)
     ! This subroutine allows to read the traffic density data read in the file aviation.nc
@@ -270 +289 @@
     REAL, ALLOCATABLE :: flight_dist_mpi(:,:,:), flight_h2o_mpi(:,:,:)
     INTEGER :: ierr
+    LOGICAL, SAVE :: first = .TRUE.
+    !$OMP THREADPRIVATE(first)
 
     ! FOR REGULAR LON LAT
@@ -288 +309 @@
 IF (grid_type==unstructured) THEN
 
-    IF (is_omp_master) THEN
-      ! Activate aviation file
-      CALL xios_set_file_attr("aviation", enabled=.TRUE.)
-      ! Activate aviation fields
-      CALL xios_set_field_attr("KMFLOWN_read", enabled=.TRUE.)
-      CALL xios_set_field_attr("KMFLOWN_interp", enabled=.TRUE.)
-
-      ! Get number of vertical levels and level values
-      CALL xios_get_axis_attr( "aviation_lev", n_glo=nleva )
-    ENDIF
-    CALL bcast_omp(nleva)
-
-    ! Allocation of arrays
-    ALLOCATE(flight_dist_read(klon, nleva,1), STAT=ierr)
-    IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate flight_dist',1)
-    ALLOCATE(flight_h2o_read(klon, nleva,1), STAT=ierr)
-    IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate flight_h2o',1)
-    ALLOCATE(aviation_lev(nleva), STAT=ierr)
-    IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate aviation_lev',1)
+    IF (first) THEN
+      IF (is_omp_master) THEN
+        ! Get number of vertical levels and level values
+        CALL xios_get_axis_attr("aviation_lev", n_glo=nleva)
+      ENDIF
+      CALL bcast_omp(nleva)
+
+      ! Allocation of arrays
+      ALLOCATE(flight_dist_read(klon, nleva, 1), STAT=ierr)
+      IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate flight_dist',1)
+      ALLOCATE(flight_h2o_read(klon, nleva, 1), STAT=ierr)
+      IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate flight_h2o',1)
+      ALLOCATE(aviation_lev(nleva), STAT=ierr)
+      IF (ierr /= 0) CALL abort_physic('lmdz_aviation', 'problem to allocate aviation_lev',1)
+
+      IF (is_omp_master) THEN
+        ! Get number of vertical levels and level values
+        CALL xios_get_axis_attr("aviation_lev", value=aviation_lev(:))
+      ENDIF
+      CALL bcast_omp(aviation_lev)
+      first = .FALSE.
+    ENDIF ! first
 
     ! Read the data from the file
@@ -318 +342 @@
         !CALL xios_recv_field("KGH2O_interp", flight_h2o_mpi(:,:,1))
         flight_h2o_mpi(:,:,:) = 0.
-        ! Get number of vertical levels and level values
-        CALL xios_get_axis_attr( "aviation_lev", value=aviation_lev(:))
     ENDIF
 
@@ -325 +347 @@
     CALL scatter_omp(flight_dist_mpi, flight_dist_read)
     CALL scatter_omp(flight_h2o_mpi, flight_h2o_read)
-    CALL bcast_omp(aviation_lev)
 
 ELSE

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

@@ -75 +75 @@
     USE write_field_phy
     use wxios_mod, ONLY: g_ctx, wxios_set_context
-    USE lmdz_aviation, ONLY : read_aviation_emissions, aviation_water_emissions, vertical_interpolation_aviation
+    USE lmdz_aviation, ONLY : init_read_aviation_emissions, read_aviation_emissions, &
+        aviation_water_emissions, vertical_interpolation_aviation
     USE lmdz_lscp, ONLY : lscp
     USE lmdz_call_cloud_optics_prop, ONLY : call_cloud_optics_prop
@@ -1393 +1394 @@
        IF (read_climoz>=1 .AND. create_etat0_limit .AND. grid_type==unstructured) &
           CALL regr_horiz_time_climoz(read_climoz,ok_daily_climoz)
+       CALL init_read_aviation_emissions
 
        print*, '================================================='