Changeset 3940 for trunk/LMDZ.PLUTO/libf/phypluto/soil.F90

Timestamp:

Oct 29, 2025, 9:06:42 PM (8 weeks ago)

Author:

tbertrand

Message:

PLUTO PCM :
Preparing for Triton simulations, implementation of internal heat flux and new North-South convention for the subsolar point
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/soil.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/soil.F90

@@ -9 +9 @@
       use planete_mod, only: year_day
       use geometry_mod, only: longitude, latitude ! in radians
-      use callkeys_mod, only: changeti
+      use callkeys_mod, only: changeti,fluxgeo,fluxgeo2,patchflux,assymflux
 
       implicit none
@@ -141 +141 @@
       enddo
 
+! Flux at the bottom
+      do ig=1,ngrid
+
+       if (assymflux) then
+         if ( (latitude(ig)*180./pi.le.90.).and. &
+                (latitude(ig)*180./pi.ge.0.)  ) then
+            tsoil(ig,nsoil) = tsoil(ig,nsoil)  &
+            + timestep*fluxgeo2/(volcapa*(layer(nsoil)-layer(nsoil-1)))
+         else   
+            tsoil(ig,nsoil) = tsoil(ig,nsoil) &
+            + timestep*fluxgeo/(volcapa*(layer(nsoil)-layer(nsoil-1)))
+         endif
+
+       elseif(patchflux.eq.1) then
+         if ( (longitude(ig)*180./pi.le.180.).and.(longitude(ig)*180./pi.ge.174.) &
+         .and.(((latitude(ig)*180./pi.le.46.).and.(latitude(ig)*180./pi.ge.42.)) &
+         .or.  ((latitude(ig)*180./pi.le.36.).and.(latitude(ig)*180./pi.ge.32.)) &
+         .or.  ((latitude(ig)*180./pi.le.26.).and.(latitude(ig)*180./pi.ge.22.)) &
+         .or.  ((latitude(ig)*180./pi.le.16.).and.(latitude(ig)*180./pi.ge.12.)) &
+         .or.  ((latitude(ig)*180./pi.le.6.).and.(latitude(ig)*180./pi.ge.2.)) &
+                ) ) then
+            tsoil(ig,nsoil) = tsoil(ig,nsoil) &
+            + timestep*fluxgeo2/(volcapa*(layer(nsoil)-layer(nsoil-1)))
+         else   
+            tsoil(ig,nsoil) = tsoil(ig,nsoil) &
+             + timestep*fluxgeo/(volcapa*(layer(nsoil)-layer(nsoil-1)))
+         endif
+
+       else 
+         tsoil(ig,nsoil) = tsoil(ig,nsoil) &
+         + timestep*fluxgeo/(volcapa*(layer(nsoil)-layer(nsoil-1)))
+
+       endif
+      enddo
+
       endif! of if (firstcall)