Changeset 4127 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Apr 20, 2022, 11:01:09 AM (2 years ago)

Author:

acozic

Message:

modify lmdz code in order to transfer to inca model all infomations necessary for the coupling with dynamico

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -34 +34 @@
     USE FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
     USE fonte_neige_mod, ONLY  : fonte_neige_get_vars
-    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
+    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg, longitude,latitude, &
+         boundslon,boundslat, dx, dy, ind_cell_glo
     USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, &
          histwrite, ju2ymds, ymds2ju, getin
@@ -55 +56 @@
     USE phystokenc_mod, ONLY: offline, phystokenc
     USE phys_cal_mod, only: year_len, mth_len, days_elapsed, jh_1jan, &
-         year_cur, mth_cur,jD_cur, jH_cur, jD_ref, day_cur, hour
+         year_cur, mth_cur,jD_cur, jH_cur, jD_ref, day_cur, hour, calend
 !!  USE phys_local_var_mod, ONLY : a long list of variables
 !!              ==> see below, after "CPP Keys" section
@@ -70 +71 @@
     USE surface_data,     ONLY : type_ocean, ok_veget, landice_opt
     USE time_phylmdz_mod, only: annee_ref, current_time, day_ini, day_ref, &
-          day_step_phy, itau_phy, pdtphys, raz_date, start_time, update_time
+          day_step_phy, itau_phy, pdtphys, raz_date, start_time, update_time, ndays
     USE tracinca_mod, ONLY: config_inca
     USE tropopause_m,     ONLY: dyn_tropopause
@@ -2028 +2029 @@
           WRITE(lunout,*) 'initial time chemini', days_elapsed, calday
 
-          CALL chemini(  &
-               rg, &
-               ra, &
-               cell_area, &
+          call init_const_lmdz( &
+          ndays, nbsrf, is_oce,is_sic, is_ter,is_lic, calend, &
+          config_inca)
+
+          CALL init_inca_geometry( & 
+               longitude, latitude, &
+               boundslon, boundslat, &
+               dx, dy, cell_area, ind_cell_glo) 
+
+
+          CALL chemini(  pplay, &
+               nbp_lon, nbp_lat, &
                latitude_deg, &
                longitude_deg, &
@@ -5083 +5092 @@
             pphis, &
             zx_rh, &
-            aps, bps, ap, bp)
+            aps, bps, ap, bp, lafin)
 
        CALL VTe(VTinca)
@@ -5345 +5354 @@
 #ifdef CPP_XIOS
        IF (is_omp_master) CALL xios_context_finalize
+
+#ifdef INCA
+       if (type_trac == 'inca' ) then 
+          IF (is_omp_master .and. grid_type==unstructured) THEN 
+             CALL finalize_inca
+          ENDIF
+       endif
+#endif
+
 #endif
        WRITE(lunout,*) ' physiq fin, nombre de steps ou cvpas = 1 : ', Ncvpaseq1