--- MITgcm_contrib/jscott/igsm/src/input.F	2009/09/01 22:03:56	1.7
+++ MITgcm_contrib/jscott/igsm/src/input.F	2009/09/17 15:48:38	1.8
@@ -1,4 +1,4 @@
-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 $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jscott/igsm/src/input.F,v 1.8 2009/09/17 15:48:38 jscott Exp $
 C $Name:  $
 
 #include "ctrparam.h"
@@ -89,7 +89,7 @@
      &  ,bc_data
       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
+     & ,flin_nep,last_clm,SEN_dat
 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
@@ -105,7 +105,7 @@
       common /sulfdata/sulf1986,sulf2050,sulfamp,SO2_EM
       common /o3data/o3_data
       common /solardata/S0C_data
-      common/aexpc/AEXP,ISTRT1,ISTRTCHEM,LYEAREM
+      common/aexpc/AEXP,ISTRT1,ISTRTCHEM
       common/qfl/QFLUX(JM0,0:13),ZOAV(JM0),QFLUXT(JM0)
       common/fixcld/cldssm(JM0,LM0,0:13),cldmcm(JM0,LM0,0:13)
      &  ,CLDSST(JM0,LM0),CLDMCT(JM0,LM0)
@@ -120,7 +120,6 @@
      *  KM,KINC,COEK,INDAY,IMONTH,LDAY,LMONTH,LYEAR,AEXP,               1519.5  
      *  READGHG,wr25,LFR,ISTRT1,PCLOUD,QFCOR,TRANSR,WRCLD,NWRCLD,CONTRR,
      *  ISTWRC,CLDFEED,OBSFOR,ALFFOR,YEARGT,CO2IN,ISTRTCHEM,
-     &  LYEAREM,
      *  AERFOR,AERF4BC,
      *  S0RATE,CFS0X,
      *  CFAEROSOL,CFBC,
@@ -136,7 +135,7 @@
      &  dirdat1,dirdat2
      &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss
      &  ,chem_init,chem_init2,chem_init4pfc,chemrstfl
-     &	  ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av
+     &	  ,oco2file,co2rfile,caruptfile,flrco2av
      &  ,flin_nep,last_clm,SEN_dat,nemdatdir
      &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer
      &  ,STRARFOR,GSOEQ,CO2FOR,CO2F,FORSULF,FORBC,S0FOR,FORVOL
@@ -167,7 +166,6 @@
       common/COMCLD/READGHG,PCLOUD,WRCLD,NWRCLD,NWRCL,INYEAR,JNDAY
      &,CFAEROSOL,ALFA,CFBC,cfvolaer
       common/ BACKGRGHG/GHGBGR(5)
-      common/CO2EM/emiss_data
       COMMON/CO2TRND/ALFFOR,CO2TR,YEARGT,CO2IN,INYRAD
       common/ S0XR/S0RATE,CFS0X
       common/cldfdb/coefcl(3),CLDFEED,SEN_dat
@@ -195,6 +193,7 @@
       logical odifcarbon,ocarcont,ocarindata
 #endif
 
+      common /SO2EMIN/SO2EM
 #if (!defined PREDICTED_GASES)                                               
 #if (defined CPL_TEM || defined CPL_OCEANCO2 )
       common /ATCO2/atm_co2(jm0),oco2file,co2rfile
@@ -259,7 +258,6 @@
       ALFFOR=0.
       S0RATE=0.
       YEARGT=1958.
-      LYEAREM=2100
       CO2IN=280.
       INYRAD=1.
       coefcl(1)=0.
@@ -354,14 +352,12 @@
       chem_init4pfc = 'undefined'
       chemrstfl = 'undefined'
       nemdatdir = 'TEMDATA'
-      emiss_data='edaily.dat'
       last_nep = 'undefined'
       last_clm = 'undefined'
       flin_nep = 'undefined'
       init_4nem = 'undefined'
       fnememiss = 'undefined'
       pov_deepo = 'undefined'
-      SO2ERATIO= 'undefined'
       SKIPSE=1.
       TAUT=6.                                                           1537.   
       TAUT=24.
@@ -628,8 +624,8 @@
       else
          print *,' Without changes in  stochastic precip'
       endif
-!#if ( defined PREDICTED_AEROSOL )
-#if ( defined CPL_CHEM )
+#if ( defined PREDICTED_AEROSOL )
+!#if ( defined CPL_CHEM )
 Cold  AFBYCF=0.6725
 Cold  SO2EREF=123.57
 Cigsm1AFBYCF=0.6054
@@ -644,10 +640,11 @@
 
 
       SO2EREF=147.375
-      open(664,file=SO2ERATIO,
-     &        form='formatted',
-     &        status='old')
-      read(664,'(f10.6)')SO2EM
+!     open(664,file=SO2ERATIO,
+!    &        form='formatted',
+!    &        status='old')
+!     read(664,'(f10.6)')SO2EM
+       print *,' SO2EM from input ',SO2EM
       SO2ER=SO2EM/SO2EREF
 !     F90BYF80=0.948
 !     AERFOR90=AERFOR*F90BYF80
@@ -735,11 +732,11 @@
 #endif
 #if ( defined CPL_TEM )
           open (876,file=last_nep,form='unformatted',status='new')
-          open (877,file=last_clm,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
 #if ( defined CPL_NEM )
-          open (368,file=init_4nem,form='unformatted',status='old')
+!         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='replace')
@@ -757,9 +754,9 @@
      *    status='OLD',form='unformatted')
 #if ( defined CPL_TEM )
           open (876,file=last_nep,form='unformatted',status='old')
-          open (877,file=last_clm,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')
+!         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='replace')