Changeset 7516

Timestamp:

2022-03-11T12:01:03+01:00 (2 years ago)

Author:

josefine.ghattas

Message:

Correct error in initalization of grid cells where no information in the soils_param.nc source grid were found (such as islands in the oceans) for the case using zobler map with standard interpolation method.
This changes slightly clayfraction, sandfraction and siltfraction and could also change njsc but it was not the case in the test using crujra 2deg grid. Old restart files should not be used.

See ticket #812

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90 (modified) (1 diff)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90

@@ -3429 +3429 @@
                    ENDDO
 
-                   IF ( sgn .LT. min_sechiba) THEN
+                   IF ( (sgn .LT. min_sechiba) .OR. (atext(ib) .LT. min_sechiba) ) THEN
                       ! Set default values if grid cells were only covered by glaciers or ocean
-                      ! or if now information on the source grid was found.
+                      ! or if now information on the source grid was found (atext(ib)=-1).
                       nbexp = nbexp + 1
                       njsc(ib) = usda_default ! 6 = Loam