Changeset 3787 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jun 3, 2025, 3:59:23 PM (4 weeks ago)

Author:

mmaurice

Message:

Generic PCM

Add flag read_time_from_start1D to read start day and time from
start1D.nc's controle tab when restarting a 1D simulation. Default
behavior is to read start day and time from rcm1d.def file, but it is
convenient for chained simulations to read them from start1D.nc.

MM

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -85 +85 @@
       INTEGER ilayer,ilevel,isoil,idt,iq
       LOGICAl firstcall,lastcall
-      LOGICAL restart
+      LOGICAL restart, read_time_from_start1D
       INTEGER nid_restart1D, nid_restartfi
       INTEGER controleid, uid, vid, tempid, tsurfid, did, tid
@@ -526 +526 @@
         endif
         psurf = controle1D(1)
-        PRINT *,"In restart1D.nc, day0=",controle1D(2),
-     &          " and time=",controle1D(3)
-        PRINT*,"(These values are NOT used, ",
-     &         "those given in rcm1d.def are used)"
       else ! restart
         psurf = -99999.
@@ -574 +570 @@
 
       !else
+      call getin("read_time_from_start1D",read_time_from_start1D)
+      if (restart .and. read_time_from_start1D) then
+            PRINT *,"In restart1D.nc, day0=",controle1D(2),
+     &              " and time=",controle1D(3)
+            day0 = controle1D(2)
+            time = controle1D(3)
+      else
         day0 = 0 ! default value for day0
         write(*,*) 'Initial date (in martian sols ; =0 at Ls=0)?'
@@ -585 +588 @@
         write(*,*)" time = ",time
         time=time/24.E+0 ! convert time (hours) to fraction of sol
+      end if
       !endif