--- MITgcm_contrib/jscott/igsm/src/input.F	2007/04/23 21:20:18	1.3
+++ MITgcm_contrib/jscott/igsm/src/input.F	2009/09/01 22:03:56	1.7
@@ -1,3 +1,5 @@
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jscott/igsm/src/input.F,v 1.7 2009/09/01 22:03:56 jscott Exp $
+C $Name:  $
 
 #include "ctrparam.h"
 
@@ -44,9 +46,9 @@
 #include "AGRID.h"
 #endif
 
-#if ( defined CPL_TEM )
-#include "TEM.h"
-#endif
+cjrs done in driver.h #if ( defined CPL_TEM )
+cjrs#include "TEM.h"
+cjrs#endif
 !
 ! === Chien Wang 062904
 !
@@ -77,17 +79,20 @@
       character *120 file1,file2,plotfl,nwrfl
       character * 120  t3file,tsfile,zmfile,qffile,clfile,wrcldf
      &	  ,ochemfile,deepco2in
-     &	  ,oco2file,co2rfile,caruptfile,flrco2av
+cjrs DRIVER.h caruptfile  & ,oco2file,co2rfile,caruptfile,flrco2av
+     &	  ,oco2file,co2rfile,flrco2av
      &    ,ghg_monthly,ghg_monthly2,co2_data,o3_data
      &    ,bgrghg_data
       character * 120 sulf1986,sulf2050,sulfamp,SO2_EM,
      &  S0C_data,
      &  dirdat1,dirdat2
      &  ,bc_data
-      character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo
+      character * 120 chemdata,chemout,init_4nem,pov_deepo
+cjrs DRIVER.h last_nep character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo
      & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat
-     &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,
-     & chem_init,chem_init2,chemrstfl
+cjrs DRIVER.h fnememiss &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,
+     &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,
+     & chem_init,chem_init2,chem_init4pfc,chemrstfl
       common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf
      *,t3file,tsfile,zmfile,ochemfile,deepco2in
       character * 120 ghostfile
@@ -130,7 +135,7 @@
      &  S0C_data,cfvolaer,
      &  dirdat1,dirdat2
      &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss
-     &  ,chem_init,chem_init2,chemrstfl
+     &  ,chem_init,chem_init2,chem_init4pfc,chemrstfl
      &	  ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av
      &  ,flin_nep,last_clm,SEN_dat,nemdatdir
      &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer
@@ -139,6 +144,7 @@
      &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3
      &  ,OCNGEOM,ocngmfile,ocndata4atm
      &  ,bc_data
+     &  ,PRTREND
 Cjrs     &  ,dtatm,dtocn
 #if ( defined IPCC_EMI )
      &  ,init_co2
@@ -273,6 +279,7 @@
       GSOEQ=.false.
       VEGCH=.false.
       TRVEG=.false.
+      PRTREND=.false.
       OCNGEOM=.false.
       do L=1,LM+1
         ghostfv(L)=0.
@@ -344,6 +351,7 @@
       chemdata = 'DATA'
       chem_init = 'init-data_46x11_1991'
       chem_init2 = 'init-data2_46x11_1991'
+      chem_init4pfc = 'undefined'
       chemrstfl = 'undefined'
       nemdatdir = 'TEMDATA'
       emiss_data='edaily.dat'
@@ -556,7 +564,8 @@
 c     open file for carbon uptake
 #if ( defined CPL_TEM || defined CPL_OCEANCO2 )
 c     open(333,file=caruptfile,status='new',form='formatted')
-      open(333,file=caruptfile,form='formatted')
+      open(333,file=caruptfile,status='replace',form='formatted')
+      close(333)
 #endif
 
 #if ( defined CPL_OCEANCO2 && defined ML_2D)
@@ -614,6 +623,11 @@
       else
         print *,'No diffusion into deep ocean'
       endif
+      if (PRTREND) then
+         print *,' With changes in  stochastic precip'
+      else
+         print *,' Without changes in  stochastic precip'
+      endif
 !#if ( defined PREDICTED_AEROSOL )
 #if ( defined CPL_CHEM )
 Cold  AFBYCF=0.6725
@@ -720,7 +734,7 @@
      *    status='OLD',form='unformatted')
 #endif
 #if ( defined CPL_TEM )
-          open (367,file=last_nep,form='unformatted',status='new')
+          open (876,file=last_nep,form='unformatted',status='new')
           open (877,file=last_clm,form='unformatted',status='new')
 c         file last_clm contains data for posible restart of NEM
 c         this file is writen at the end of the run
@@ -728,7 +742,8 @@
           open (368,file=init_4nem,form='unformatted',status='old')
 c         file init_4nem contains data for the restart of NEM
 c         from the results of a previous run
-          open (277,file=fnememiss,form='unformatted',status='new')
+          open (277,file=fnememiss,form='unformatted',status='replace')
+          close(277)
 #endif
 #endif
 #if ( defined CPL_OCEANCO2 && defined ML_2D )
@@ -741,13 +756,14 @@
           open( unit=547,file=nwrfl,
      *    status='OLD',form='unformatted')
 #if ( defined CPL_TEM )
-          open (367,file=last_nep,form='unformatted',status='old')
+          open (876,file=last_nep,form='unformatted',status='old')
           open (877,file=last_clm,form='unformatted',status='old')
 #if ( defined CPL_NEM )
           open (368,file=init_4nem,form='unformatted',status='old')
 c         file init_4nem contains data for the restart of NEM
 c         from the results of a previous run
-          open (277,file=fnememiss,form='unformatted',status='new')
+          open (277,file=fnememiss,form='unformatted',status='replace')
+          close(277)
 #endif
 #endif
 #if ( defined CPL_OCEANCO2 && defined ML_2D )
@@ -1021,7 +1037,11 @@
        print *,TSCNTR,' form EXP=',EXPTSF
 !      read(576)EXPTSF
 !      print *,' TSURF form EXP=',EXPTSF
-       read(576)TSURFC,TLANDC
+#if ( defined CLM )
+        read(576)TSURFC,TLANDC
+#else
+        read(576)TSURFC
+#endif
 !      read(576)TLANDC
       endif
         CORSR=1.0