Changeset 7262 for branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/orchideedriver.f90

Timestamp:

2021-07-27T19:02:48+02:00 (3 years ago)

Author:

agnes.ducharne

Message:

As in r7065 and r76158 for the new driver. Runs OK on jean-zay with the old driver, but not yet tested with the new driver. The corresponding configurations need to be added following r7081, r7082, and r7719.

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/orchideedriver.f90 (modified) (10 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/orchideedriver.f90

@@ -50 +50 @@
   USE ioipslctrl
   USE xios_orchidee
+  
   !
   !-
@@ -317 +318 @@
   IF ( is_root_prc) THEN
      CALL globgrd_getdomsz(gridfilename, iim_glo, jjm_glo, nbindex_g, model_guess, file_id, forfilename, zoom_lon, zoom_lat)
+     write(numout,*) "nbindex_g after calling globgrd_getdomsz in orchideedriver", nbindex_g
      nbseg = 4
   ENDIF
@@ -342 +344 @@
           &               kindex_g, contfrac_glo, calendar)
   ENDIF
+  
   !
   CALL bcast(lon_glo)
@@ -369 +372 @@
   ENDDO
   !
+  
   WRITE(*,*) "Rank", mpi_rank, " Before parallel region All land points : ",  nbindex_g
   WRITE(*,*) "Rank", mpi_rank, " from ", iim_glo, " point in Lon. and ", jjm_glo, "in Lat."
@@ -382 +386 @@
   !-
   CALL grid_set_glo(iim_glo, jjm_glo, nbindex_g)
+  
+
   CALL grid_allocate_glo(nbseg)
   ! Copy the list of indexes of land points into index_g used by ORCHIDEE and then broacast to all
@@ -477 +483 @@
   !-
   CALL forcing_integration_time(date0, dt, nbdt)
+  write(numout, *) "orchideedriver date0", date0, dt, nbdt
   !
   !-
@@ -558 +565 @@
   !
   CALL ioipslctrl_restini(itau, date0, dt, rest_id, rest_id_stom, itau_offset, date0_shifted)
-  WRITE(numout,*) "itau_offset : ", itau_offset, date0, date0_shifted
-  WRITE(numout,*) "itau_offset diff = ", date0_shifted, date0, date0_shifted-date0
   !
   ! To ensure that itau starts with 0 at date0 for the restart, we have to set an off-set to achieve this. 
@@ -574 +579 @@
        date0,    year_end, month_end, day_end, julian_diff,  &
        lon,      lat,      znt)
-
   CALL sechiba_xios_initialize 
 
@@ -599 +603 @@
      timestep_interval(1) = julian_start
      timestep_interval(2) = julian_end
-     julian = julian_end
+     julian = (julian_start + julian_end) /2.0  !julian_end
+     
      !
      ! Get the forcing data
@@ -605 +610 @@
      CALL forcing_getvalues(timestep_interval, dt, zlev_tq, zlev_uv, temp_air, qair, &
           &                 precip_rain, precip_snow, swdown, lwdown, sinang, u, v, pb)
+     
      !-
      !