--- MITgcm/pkg/fizhi/fizhi_readparms.F	2004/08/06 21:12:14	1.6
+++ MITgcm/pkg/fizhi/fizhi_readparms.F	2005/05/24 16:04:10	1.11
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/fizhi/fizhi_readparms.F,v 1.6 2004/08/06 21:12:14 edhill Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/fizhi/fizhi_readparms.F,v 1.11 2005/05/24 16:04:10 molod Exp $
 C $Name:  $
 
 #include "FIZHI_OPTIONS.h"
@@ -19,6 +19,8 @@
 #include "chronos.h"
 #include "fizhi_io_comms.h"
 #include "EEPARAMS.h"
+#include "SIZE.h"
+#include "PARAMS.h"
 
 C     !INPUT PARAMETERS:
       integer myThid
@@ -26,19 +28,24 @@
 
 C     !LOCAL VARIABLES:
       character*(MAX_LEN_MBUF) msgBuf
-      integer ku
+      integer ku, ku2
       integer nymdbegin, nhmsbegin
       integer nymdcurrent, nhmscurrent
+      real runlength
+      integer  nincr
+      integer mmdd,hhmmss,nsecf2
 
       namelist / fizhi_list / 
-     &     nymdbegin, nhmsbegin,
-     &     fizhi_mnc_write_pickup, fizhi_mnc_read_pickup
+     .     nymdbegin, nhmsbegin,
+     .     fizhi_mnc_write_pickup, fizhi_mnc_read_pickup,
+     .     runlength
 
 C     Set defaults
-      fizhi_mdsio_read_pickup = .FALSE.
-      fizhi_mdsio_write_pickup = .FALSE.
-      fizhi_mnc_write_pickup = .TRUE.
-      fizhi_mnc_read_pickup = .TRUE.
+      fizhi_mdsio_read_pickup = .TRUE.
+      fizhi_mdsio_write_pickup = .TRUE.
+      fizhi_mnc_write_pickup = .FALSE.
+      fizhi_mnc_read_pickup = .FALSE.
+      runlength = 0.0
 
 C     Read Fizhi Namelist
       WRITE(msgBuf,'(A)') ' FIZHI_READPARMS: opening data.fizhi'
@@ -50,11 +57,28 @@
       close (ku)
 
 C     Read Supplemental Ascii File with Current Time Info
-      open(10,file='datetime',form='formatted')
-      read(10,1000)nymdcurrent,nhmscurrent
-      close (10)
+      CALL MDSFINDUNIT( ku2, myThid )
+      open(ku2,file='datetime0',form='formatted')
+      read(ku2,1000)nymdcurrent,nhmscurrent
+      close (ku2)
  1000 format(i8,2x,i6)
 
+C Change the length of the model run, ie, change ntimesteps
+C  if runlength has been set in the fizhi namelist
+
+      if(runlength.gt.0.) then
+       mmdd = int(runlength)
+       hhmmss = int((runlength - int(runlength))*1.e6)
+       nincr = nsecf2(hhmmss,mmdd,nymdcurrent)
+       ntimesteps = int(nincr/deltat)
+       endTime = startTime + deltat*float(ntimeSteps)
+       WRITE(msgBuf,'(A,I8)') 
+     .         ' CHANGING NUMBER OF MODEL TIMESTEPS TO',ntimesteps
+       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,SQUEEZE_RIGHT,1)
+       WRITE(msgBuf,'(A,I8)') ' CHANGING END TIME TO',endtime
+       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,SQUEEZE_RIGHT,1)
+      endif
+
 C     Fill chronos Common Block with Namelist and Other File Info
       nymd0 = nymdbegin
       nhms0 = nhmsbegin