Changeset 3620 for trunk/LMDZ.COMMON/libf/evolution/read_data_PCM_mod.F90

Timestamp:

Feb 12, 2025, 12:03:26 PM (3 days ago)

Author:

jbclement

Message:

PEM:
Few small adjustments in the code.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/read_data_PCM_mod.F90 (modified) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/read_data_PCM_mod.F90

@@ -52 +52 @@
 real, dimension(ngrid,nsoilmx,nslope,timelen), intent(out) :: watersoil_density_timeseries ! Water density timeseries in the soil layer [kg/m^3]
 ! Local variables
-integer                               :: i, j, l, islope                       ! Loop variables
+integer                               :: i, j, t, islope                       ! Loop variables
 real                                  :: A, B                                  ! Intermediate variables to compute the mean molar mass of the layer
 character(:), allocatable             :: num                                   ! For reading sloped variables
@@ -243 +243 @@
     !------
     call get_var3("tsurf"//num,var3_read_1)
-       write(*,*) "Data for tsurf"//num//" downloaded."
+    write(*,*) "Data for tsurf"//num//" downloaded."
 
     ! Compute average over the year for each point
@@ -263 +263 @@
         var3_read_3 = sum(var4_read,4)/timelen
 #ifndef CPP_1D
-        do l = 1,nsoilmx
-            call gr_dyn_fi(timelen,iim_input + 1,jjm_input + 1,ngrid,var4_read(:,:,l,:),tsoil_timeseries(:,l,islope,:))
+        do t = 1,timelen
+            call gr_dyn_fi(nsoilmx,iim_input + 1,jjm_input + 1,ngrid,var4_read(:,:,:,t),tsoil_timeseries(:,:,islope,t))
         enddo
         call gr_dyn_fi(nsoilmx,iim_input + 1,jjm_input + 1,ngrid,var3_read_3,tsoil_avg(:,:,islope))
@@ -277 +277 @@
         call get_var4("waterdensity_soil"//num,var4_read)
 #ifndef CPP_1D
-        do l = 1,nsoilmx
-            call gr_dyn_fi(timelen,iim_input + 1,jjm_input + 1,ngrid,var4_read(:,:,l,:),watersoil_density_timeseries(:,l,islope,:))
+        do t = 1,timelen
+            call gr_dyn_fi(nsoilmx,iim_input + 1,jjm_input + 1,ngrid,var4_read(:,:,:,t),watersoil_density_timeseries(:,:,islope,t))
         enddo
 #else