Changeset 4596

Timestamp:

Jun 30, 2023, 5:03:07 PM (11 months ago)

Author:

Laurent Fairhead

Message:

Resolve an unitialised variable issue
ACaubel

Location:

LMDZ6/trunk/libf/phylmd

Files:

• oasis.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/oasis.F90

@@ -141 +141 @@
     CHARACTER (len = 80)               :: abort_message 
     LOGICAL, SAVE                      :: cpl_current_omp
+    INTEGER, DIMENSION(klon_mpi)       :: ind_cell_glo_mpi
 
 !*    1. Initializations
@@ -175 +176 @@
     cpl_current = cpl_current_omp
     WRITE(lunout,*) 'Couple ocean currents, cpl_current = ',cpl_current 
+
+!************************************************************************************
+! Gather global index to be used for oasis decomposition
+!************************************************************************************
+    CALL gather_omp(ind_cell_glo,ind_cell_glo_mpi)
 
 !************************************************************************************
@@ -291 +297 @@
 
       DO jf=1, klon_mpi_para_nb(mpi_rank)
-        ig_paral(2+jf) = ind_cell_glo(jf)
+        ig_paral(2+jf) = ind_cell_glo_mpi(jf)
       ENDDO
 

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1315 +1315 @@
        irneb= strIdx(tracers(:)%name, addPhase('H2O', 'r'))
        ibs  = strIdx(tracers(:)%name, addPhase('H2O', 'b'))
-       CALL init_etat0_limit_unstruct
-       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
+!       CALL init_etat0_limit_unstruct
+!       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
        !CR:nvelles variables convection/poches froides
 
@@ -1339 +1339 @@
             read_climoz, &
             alp_offset)
+       CALL init_etat0_limit_unstruct
+       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
        CALL phys_state_var_init(read_climoz)
        CALL phys_output_var_init