Changeset 4026 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Jan 20, 2026, 9:50:44 AM (3 days ago)

Author:

tbertrand

Message:

PLUTO PCM:
Adding Paleo options for Triton from Pluto.old
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (6 diffs)

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

@@ -38 +38 @@
       use tabfi_mod, only: tab_cntrl_mod
       use wstats_mod, only: callstats, wstats, mkstats
-      use phyredem, only: physdem0, physdem1
+      use phyredem, only: physdem0, physdem1, physdem0pal, physdem1pal
       use planetwide_mod, only: planetwide_minval,planetwide_maxval,planetwide_sumval
       use mod_phys_lmdz_para, only : is_master
@@ -2174 +2174 @@
 
          if (paleo) then
-            ! time range for tendencies of ice flux qsurfyear
+            ! time range to be applied for tendencies of ice flux qsurfyear
             zdt_tot=year_day   ! Last year of simulation
 
+            ! update new reservoir of ice on the surface
             masslost(:)=0.
             massacc(:)=0.
-
             DO ig=1,ngrid
-               ! update new reservoir of ice on the surface
                DO iq=1,nq
                 ! kg/m2 to be sublimed or condensed during paleoyears
@@ -2187 +2186 @@
                            paleoyears*365.25/(zdt_tot*daysec/86400.)
 
-               ! special case if we sublime the entire reservoir
-               !! AF: TODO : fix following lines (real_area), using line below:
-            ! call planetwide_sumval((-qsurfyear(:,iq)-qsurf(:,iq))*cell_area(:),masslost)
-
-               !  IF (-qsurfyear(ig,iq).gt.qsurf(ig,iq)) THEN
-               !    masslost(iq)=masslost(iq)+real_area(ig)*   &
-               !          (-qsurfyear(ig,iq)-qsurf(ig,iq))
-               !    qsurfyear(ig,iq)=-qsurf(ig,iq)
-               !  ENDIF
-
-               !  IF (qsurfyear(ig,iq).gt.0.) THEN
-               !    massacc(iq)=massacc(iq)+real_area(ig)*qsurfyear(ig,iq)
-               !  ENDIF
-
+                ! special case if we sublime the entire reservoir
+                IF (-qsurfyear(ig,iq).gt.qsurf(ig,iq)) THEN
+                   masslost(iq)=masslost(iq)+cell_area(ig)* &
+                         (-qsurfyear(ig,iq)-qsurf(ig,iq))
+                   qsurfyear(ig,iq)=-qsurf(ig,iq)
+                ENDIF
+                ! massacc needed to redistribute masslost 
+                IF (qsurfyear(ig,iq).gt.0.) THEN
+                   massacc(iq)=massacc(iq)+cell_area(ig)*qsurfyear(ig,iq)
+                ENDIF
 
                ENDDO
@@ -2207 +2202 @@
             DO ig=1,ngrid
                DO iq=1,nq
+                 ! New reservoir paleo
                  qsurfpal(ig,iq)=qsurf(ig,iq)+qsurfyear(ig,iq)
+                 ! Redistribution of some of the mass lost when entire reservoir has sublimed
+                 ! Accumulated mass is less than calculated because sublimed mass is less
                  IF (qsurfyear(ig,iq).gt.0.) THEN
                   qsurfpal(ig,iq)=qsurfpal(ig,iq)- &
@@ -2214 +2212 @@
                ENDDO
             ENDDO
+
             ! Finally ensure conservation of qsurf
             DO iq=1,nq
@@ -2255 +2254 @@
             ! create restartfi
             if (ngrid.ne.1) then
-               print*, "physdem1pal not yet implemented"
-               stop
-               !TODO: import this routine from pluto.old
-               ! call physdem1pal("restartfi.nc",long,lati,nsoilmx,nq, &
+               call physdem0pal("restartfi.nc",longitude,latitude,nsoilmx,ngrid,nlayer,nq, &
+                         ptimestep,pdaypal,time_phys,cell_area,          &
+                         albedo_bareground,zmea,zstd,zsig,zgam,zthe,     &
+                         oblipal,eccpal,tpalnew,adjustnew,phisfipal,peri_daypal)  
+  
+                 !call physdem1pal("restartfi.nc",long,lati,nsoilmx,nq, &
                !      ptimestep,pdaypal, &
                !      ztime_restart,tsurf,tsoil,emis,q2,qsurfpal, &