Changeset 7509 for branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/thermosoil.f90

Timestamp:

2022-03-04T18:31:10+01:00 (2 years ago)

Author:

josefine.ghattas

Message:

• Bug correction when interpolating reftemp.nc file with XIOS in context_input_orchidee.xml : the latitudes needs to be inverted.
• Added output of reftemp just after it has been read.

See ticket #812

See ticket #812

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/thermosoil.f90 (modified) (2 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/thermosoil.f90

@@ -1881 +1881 @@
     IF (xios_interpolation) THEN
 
-       CALL xios_orchidee_recv_field('reftemp_interp',reftemp(:,1))
+       CALL xios_orchidee_recv_field('reftemp_interp',reftemp_file)
 
        DO ib=1, kjpindex
-           reftemp(ib,:) = reftemp(ib,1) + ZeroCelsius 
+           reftemp(ib,:) = reftemp_file(ib) + ZeroCelsius 
        END DO
        areftemp = 1.0
@@ -1933 +1933 @@
     ! Write diagnostics
     CALL xios_orchidee_send_field("areftemp",areftemp)
+    CALL xios_orchidee_send_field("interp_diag_reftemp",reftemp_file)
     
   END SUBROUTINE thermosoil_read_reftempfile