/[MITgcm]/MITgcm_contrib/jscott/igsm/src/input.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/igsm/src/input.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.7 by jscott, Tue Sep 1 22:03:56 2009 UTC revision 1.8 by jscott, Thu Sep 17 15:48:38 2009 UTC
# Line 89  cjrs DRIVER.h caruptfile  & ,oco2file,co Line 89  cjrs DRIVER.h caruptfile  & ,oco2file,co
89       &  ,bc_data       &  ,bc_data
90        character * 120 chemdata,chemout,init_4nem,pov_deepo        character * 120 chemdata,chemout,init_4nem,pov_deepo
91  cjrs DRIVER.h last_nep character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo  cjrs DRIVER.h last_nep character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo
92       & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat       & ,flin_nep,last_clm,SEN_dat
93  cjrs DRIVER.h fnememiss &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,  cjrs DRIVER.h fnememiss &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,
94       &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,       &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,
95       & chem_init,chem_init2,chem_init4pfc,chemrstfl       & chem_init,chem_init2,chem_init4pfc,chemrstfl
# Line 105  c==== 012201 Line 105  c==== 012201
105        common /sulfdata/sulf1986,sulf2050,sulfamp,SO2_EM        common /sulfdata/sulf1986,sulf2050,sulfamp,SO2_EM
106        common /o3data/o3_data        common /o3data/o3_data
107        common /solardata/S0C_data        common /solardata/S0C_data
108        common/aexpc/AEXP,ISTRT1,ISTRTCHEM,LYEAREM        common/aexpc/AEXP,ISTRT1,ISTRTCHEM
109        common/qfl/QFLUX(JM0,0:13),ZOAV(JM0),QFLUXT(JM0)        common/qfl/QFLUX(JM0,0:13),ZOAV(JM0),QFLUXT(JM0)
110        common/fixcld/cldssm(JM0,LM0,0:13),cldmcm(JM0,LM0,0:13)        common/fixcld/cldssm(JM0,LM0,0:13),cldmcm(JM0,LM0,0:13)
111       &  ,CLDSST(JM0,LM0),CLDMCT(JM0,LM0)       &  ,CLDSST(JM0,LM0),CLDMCT(JM0,LM0)
# Line 120  c==== 012201 Line 120  c==== 012201
120       *  KM,KINC,COEK,INDAY,IMONTH,LDAY,LMONTH,LYEAR,AEXP,               1519.5         *  KM,KINC,COEK,INDAY,IMONTH,LDAY,LMONTH,LYEAR,AEXP,               1519.5  
121       *  READGHG,wr25,LFR,ISTRT1,PCLOUD,QFCOR,TRANSR,WRCLD,NWRCLD,CONTRR,       *  READGHG,wr25,LFR,ISTRT1,PCLOUD,QFCOR,TRANSR,WRCLD,NWRCLD,CONTRR,
122       *  ISTWRC,CLDFEED,OBSFOR,ALFFOR,YEARGT,CO2IN,ISTRTCHEM,       *  ISTWRC,CLDFEED,OBSFOR,ALFFOR,YEARGT,CO2IN,ISTRTCHEM,
      &  LYEAREM,  
123       *  AERFOR,AERF4BC,       *  AERFOR,AERF4BC,
124       *  S0RATE,CFS0X,       *  S0RATE,CFS0X,
125       *  CFAEROSOL,CFBC,       *  CFAEROSOL,CFBC,
# Line 136  c==== 012201 Line 135  c==== 012201
135       &  dirdat1,dirdat2       &  dirdat1,dirdat2
136       &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss       &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss
137       &  ,chem_init,chem_init2,chem_init4pfc,chemrstfl       &  ,chem_init,chem_init2,chem_init4pfc,chemrstfl
138       &    ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av       &    ,oco2file,co2rfile,caruptfile,flrco2av
139       &  ,flin_nep,last_clm,SEN_dat,nemdatdir       &  ,flin_nep,last_clm,SEN_dat,nemdatdir
140       &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer       &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer
141       &  ,STRARFOR,GSOEQ,CO2FOR,CO2F,FORSULF,FORBC,S0FOR,FORVOL       &  ,STRARFOR,GSOEQ,CO2FOR,CO2F,FORSULF,FORBC,S0FOR,FORVOL
# Line 167  c jrs      common/TIMESTEPS/dtatm,dtocn Line 166  c jrs      common/TIMESTEPS/dtatm,dtocn
166        common/COMCLD/READGHG,PCLOUD,WRCLD,NWRCLD,NWRCL,INYEAR,JNDAY        common/COMCLD/READGHG,PCLOUD,WRCLD,NWRCLD,NWRCL,INYEAR,JNDAY
167       &,CFAEROSOL,ALFA,CFBC,cfvolaer       &,CFAEROSOL,ALFA,CFBC,cfvolaer
168        common/ BACKGRGHG/GHGBGR(5)        common/ BACKGRGHG/GHGBGR(5)
       common/CO2EM/emiss_data  
169        COMMON/CO2TRND/ALFFOR,CO2TR,YEARGT,CO2IN,INYRAD        COMMON/CO2TRND/ALFFOR,CO2TR,YEARGT,CO2IN,INYRAD
170        common/ S0XR/S0RATE,CFS0X        common/ S0XR/S0RATE,CFS0X
171        common/cldfdb/coefcl(3),CLDFEED,SEN_dat        common/cldfdb/coefcl(3),CLDFEED,SEN_dat
# Line 195  cjrs alreadyin DRIVER.h#include "CLM.COM Line 193  cjrs alreadyin DRIVER.h#include "CLM.COM
193        logical odifcarbon,ocarcont,ocarindata        logical odifcarbon,ocarcont,ocarindata
194  #endif  #endif
195    
196          common /SO2EMIN/SO2EM
197  #if (!defined PREDICTED_GASES)                                                #if (!defined PREDICTED_GASES)                                              
198  #if (defined CPL_TEM || defined CPL_OCEANCO2 )  #if (defined CPL_TEM || defined CPL_OCEANCO2 )
199        common /ATCO2/atm_co2(jm0),oco2file,co2rfile        common /ATCO2/atm_co2(jm0),oco2file,co2rfile
# Line 259  C Line 258  C
258        ALFFOR=0.        ALFFOR=0.
259        S0RATE=0.        S0RATE=0.
260        YEARGT=1958.        YEARGT=1958.
       LYEAREM=2100  
261        CO2IN=280.        CO2IN=280.
262        INYRAD=1.        INYRAD=1.
263        coefcl(1)=0.        coefcl(1)=0.
# Line 354  cjrs      dtocn=1 Line 352  cjrs      dtocn=1
352        chem_init4pfc = 'undefined'        chem_init4pfc = 'undefined'
353        chemrstfl = 'undefined'        chemrstfl = 'undefined'
354        nemdatdir = 'TEMDATA'        nemdatdir = 'TEMDATA'
       emiss_data='edaily.dat'  
355        last_nep = 'undefined'        last_nep = 'undefined'
356        last_clm = 'undefined'        last_clm = 'undefined'
357        flin_nep = 'undefined'        flin_nep = 'undefined'
358        init_4nem = 'undefined'        init_4nem = 'undefined'
359        fnememiss = 'undefined'        fnememiss = 'undefined'
360        pov_deepo = 'undefined'        pov_deepo = 'undefined'
       SO2ERATIO= 'undefined'  
361        SKIPSE=1.        SKIPSE=1.
362        TAUT=6.                                                           1537.          TAUT=6.                                                           1537.  
363        TAUT=24.        TAUT=24.
# Line 628  C       OLD Line 624  C       OLD
624        else        else
625           print *,' Without changes in  stochastic precip'           print *,' Without changes in  stochastic precip'
626        endif        endif
627  !#if ( defined PREDICTED_AEROSOL )  #if ( defined PREDICTED_AEROSOL )
628  #if ( defined CPL_CHEM )  !#if ( defined CPL_CHEM )
629  Cold  AFBYCF=0.6725  Cold  AFBYCF=0.6725
630  Cold  SO2EREF=123.57  Cold  SO2EREF=123.57
631  Cigsm1AFBYCF=0.6054  Cigsm1AFBYCF=0.6054
# Line 644  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/ Line 640  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/
640    
641    
642        SO2EREF=147.375        SO2EREF=147.375
643        open(664,file=SO2ERATIO,  !     open(664,file=SO2ERATIO,
644       &        form='formatted',  !    &        form='formatted',
645       &        status='old')  !    &        status='old')
646        read(664,'(f10.6)')SO2EM  !     read(664,'(f10.6)')SO2EM
647           print *,' SO2EM from input ',SO2EM
648        SO2ER=SO2EM/SO2EREF        SO2ER=SO2EM/SO2EREF
649  !     F90BYF80=0.948  !     F90BYF80=0.948
650  !     AERFOR90=AERFOR*F90BYF80  !     AERFOR90=AERFOR*F90BYF80
# Line 735  c     LMCM=LTM+2 Line 732  c     LMCM=LTM+2
732  #endif  #endif
733  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
734            open (876,file=last_nep,form='unformatted',status='new')            open (876,file=last_nep,form='unformatted',status='new')
735            open (877,file=last_clm,form='unformatted',status='new')  !         open (877,file=last_clm,form='unformatted',status='new')
736  c         file last_clm contains data for posible restart of NEM  c         file last_clm contains data for posible restart of NEM
737  c         this file is writen at the end of the run  c         this file is writen at the end of the run
738  #if ( defined CPL_NEM )  #if ( defined CPL_NEM )
739            open (368,file=init_4nem,form='unformatted',status='old')  !         open (368,file=init_4nem,form='unformatted',status='old')
740  c         file init_4nem contains data for the restart of NEM  c         file init_4nem contains data for the restart of NEM
741  c         from the results of a previous run  c         from the results of a previous run
742            open (277,file=fnememiss,form='unformatted',status='replace')            open (277,file=fnememiss,form='unformatted',status='replace')
# Line 757  C      For restart of the run Line 754  C      For restart of the run
754       *    status='OLD',form='unformatted')       *    status='OLD',form='unformatted')
755  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
756            open (876,file=last_nep,form='unformatted',status='old')            open (876,file=last_nep,form='unformatted',status='old')
757            open (877,file=last_clm,form='unformatted',status='old')  !         open (877,file=last_clm,form='unformatted',status='old')
758  #if ( defined CPL_NEM )  #if ( defined CPL_NEM )
759            open (368,file=init_4nem,form='unformatted',status='old')  !         open (368,file=init_4nem,form='unformatted',status='old')
760  c         file init_4nem contains data for the restart of NEM  c         file init_4nem contains data for the restart of NEM
761  c         from the results of a previous run  c         from the results of a previous run
762            open (277,file=fnememiss,form='unformatted',status='replace')            open (277,file=fnememiss,form='unformatted',status='replace')
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22