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

Timestamp:

Oct 29, 2025, 9:06:42 PM (7 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

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• inifis_mod.F90 (modified) (1 diff)
• orbitetriton.F (modified) (1 diff)
• soil.F90 (modified) (2 diffs)

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

@@ -462 +462 @@
      if (is_master) write(*,*)trim(rname)//": alpha_soil = ",alpha_soil
 
+     if (is_master) write(*,*)trim(rname)//&
+       "Geothermal flux (W) at the bottom layer"
+     fluxgeo=0. ! default value
+     call getin_p("fluxgeo",fluxgeo)
+     if (is_master) write(*,*)trim(rname)//" fluxgeo = ",fluxgeo
+
+     if (is_master) write(*,*)trim(rname)//&
+       "Assymetry flux (W) at the bottom layer"
+     assymflux=.false. ! default value
+     call getin_p("assymflux",assymflux)
+     if (is_master) write(*,*)trim(rname)//" assymflux = ",assymflux
+
+     if (is_master) write(*,*)trim(rname)//&
+       "Geothermal flux (W) for assymetry"
+     fluxgeo2=fluxgeo ! default value
+     call getin_p("fluxgeo2",fluxgeo2)
+     if (is_master) write(*,*)trim(rname)//" fluxgeo2 = ",fluxgeo2
+
+     if (is_master) write(*,*)trim(rname)//&
+       "Warm patch of flux"
+     patchflux=0 ! default value
+     call getin_p("patchflux",patchflux)
+     if (is_master) write(*,*)trim(rname)//" patchflux = ",patchflux
+
 ! Chemistry in the thermosphere
      if (is_master) write(*,*) trim(rname)//": Use deposition ?"

trunk/LMDZ.PLUTO/libf/phypluto/orbitetriton.F

@@ -62 +62 @@
       zint = zint + BB*sin(zalpha0-zbeta0+(zbeta1-zalpha1)*zan)
       zint = zint + CC*sin(zalpha0+zbeta0-(zalpha1+zbeta1)*zan)
-      pdecli = -asin (zint)
+      pdecli = asin (zint)
 
       RETURN

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)