Changeset 2372 for LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd

Timestamp:

Oct 13, 2015, 2:41:47 PM (9 years ago)

Author:

acozic

Message:

Change call to inca initialisation to fit with new sections dynamique/physic

File:

• LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/iniphysiq_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/iniphysiq_mod.F90

@@ -13 +13 @@
   USE dimphy, ONLY: init_dimphy
   USE mod_grid_phy_lmdz, ONLY: klon_glo,  & ! number of atmospheric columns (on full grid)
-                               regular_lonlat  ! regular longitude-latitude grid type
+                               regular_lonlat, &  ! regular longitude-latitude grid type
+                               nbp_lon, nbp_lat, nbp_lev
   USE mod_phys_lmdz_para, ONLY: klon_omp, & ! number of columns (on local omp grid)
                                 klon_omp_begin, & ! start index of local omp subgrid
@@ -32 +33 @@
   USE CHEM_REP, ONLY : Init_chem_rep_phys
 #endif
-  USE control_mod, ONLY: dayref,anneeref,day_step,nday,offline
+  USE control_mod, ONLY: dayref,anneeref,day_step,nday,offline, iphysiq
   USE inifis_mod, ONLY: inifis
   USE time_phylmdz_mod, ONLY: init_time
@@ -44 +45 @@
                                  south_west, south_east
   USE mod_interface_dyn_phys, ONLY :  init_interface_dyn_phys
+#ifdef INCA
+  USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic
+  USE parallel_lmdz, ONLY : mpi_size
+  USE mod_const_mpi, ONLY : COMM_LMDZ
+  USE bands, ONLY : distrib_phys
+  USE mod_phys_lmdz_omp_data, ONLY: klon_omp 
+#endif
   IMPLICIT NONE
 
@@ -282 +290 @@
 #endif
   ENDIF
+!$OMP END PARALLEL
+
+  IF (type_trac == 'inca') THEN
+#ifdef INCA
+     call init_const_lmdz( &
+          anneeref,dayref, &
+          iphysiq,day_step,nday,  &
+          nbsrf, is_oce,is_sic, &
+          is_ter,is_lic, calend)
+     call init_inca_para( &
+          nbp_lon,nbp_lat,nbp_lev,klon_glo,mpi_size, &
+          distrib_phys,COMM_LMDZ)
+#endif
+  END IF
+!$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/)
 
   ! Additional initializations for aquaplanets 
@@ -287 +310 @@
     CALL iniaqua(klon_omp,iflag_phys)
   END IF
+
+  IF (type_trac == 'inca') THEN
+#ifdef INCA
+     CALL init_inca_dim(klon_omp,nbp_lev,nbp_lon,nbp_lat - 1, &
+          rlonu,rlatu,rlonv,rlatv)
+#endif
+  END IF
+
 !$OMP END PARALLEL