Changeset 7991 for branches/ORCHIDEE_2_2

Timestamp:

2023-04-21T16:37:49+02:00 (19 months ago)

Author:

josefine.ghattas

Message:

Correction on lake_reservoir in routing simple, done by Y Meurdesoif

File:

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

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/routing_simple.f90

@@ -2669 +2669 @@
        !
        lake_reservoir(ig) = lake_reservoir(ig) + lakeinflow(ig)
-       !uptake in Kg/dt
-       refill = MAX(zero, maxevap_lake * (un - humrel(ig)) * dt_routing * total_area)
-       return_lakes(ig) = MIN(refill, lake_reservoir(ig))
-       lake_reservoir(ig) = lake_reservoir(ig) - return_lakes(ig)
-       !Return in Kg/m^2/dt
+
        IF ( doswamps ) THEN
+          ! uptake in Kg/dt
+          refill = MAX(zero, maxevap_lake * (un - humrel(ig)) * dt_routing * total_area)
+          return_lakes(ig) = MIN(refill, lake_reservoir(ig))
+          lake_reservoir(ig) = lake_reservoir(ig) - return_lakes(ig)
+          ! Return in Kg/m^2/dt
           return_lakes(ig) = return_lakes(ig)/total_area
        ELSE