Changeset 4025 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto

Timestamp:

Jan 20, 2026, 9:49:28 AM (38 hours ago)

Author:

tbertrand

Message:

PLUTO PCM:
Fix in newstart to allow topography to remain unchanged when using startarchive.nc
TB

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F

@@ -544 +544 @@
         write(*,*) 'Please choose the relevant file',
      &  ' (e.g. "surface_pluto.nc")'
-        write(*,*) ' or "none" to not use any (and set planetary'
-        write(*,*) '  albedo and surface thermal inertia)'
+        write(*,*) ' or "none" to not use any'
         read(*,fmt='(a50)') surfacefile
 
@@ -552 +551 @@
           CALL datareadnc(relief,surfacefile,phis,alb,surfith,
      &          zmeaS,zstdS,zsigS,zgamS,ztheS)
+
+          CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phis,phisfi)
+
         else
-        ! specific case when not using a "surface.nc" file
-          phis(:,:)=0
+          ! specific case when not using a "surface.nc" file
           zmeaS(:,:)=0
           zstdS(:,:)=0
@@ -562 +563 @@
           n2fracfi(:)=1.
 
-          write(*,*) "Enter value of albedo of the bare ground:"
-          write(*,*) "Might not count"
-          read(*,*) alb(1,1)
-          alb(:,:)=alb(1,1)
-
-          write(*,*) "Enter value of thermal inertia of soil:"
-          write(*,*) "Might not count"
-          read(*,*) surfith(1,1)
-          surfith(:,:)=surfith(1,1)
+          write(*,*)'OK : skipping topography change'
 
         endif ! of if (surfacefile.ne."none")
 
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phis,phisfi)
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,surfith,surfithfi)
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi) ! erased at next step
+        !CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,surfith,surfithfi)
+        !CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi) ! erased at next step
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zmeaS,zmea)
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zstdS,zstd)
@@ -598 +590 @@
 
         write(*,*) "OK, read start_archive file"
-        ! copy soil thermal inertia
+
+        ! copy topography 
+        phis(:,:)=phisold_newgrid(:,:)
+        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phis,phisfi)
+
+        ! copy soil thermal inertia 
         ithfi(:,:)=inertiedat(:,:)
-        ! copy topo
-        ! CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phisold_newgrid,phisfi)
-        ! Initialize n2frac
-        n2fracfi(1:ngridmx)=1.
 
         ierr= NF_CLOSE(nid)