--- MITgcm_contrib/jscott/igsm/src/input.F 2007/10/04 20:40:27 1.4 +++ MITgcm_contrib/jscott/igsm/src/input.F 2009/09/17 15:48:38 1.8 @@ -1,3 +1,5 @@ +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" @@ -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 - & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat - & ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss, - & chem_init,chem_init2,chemrstfl + 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,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 common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf *,t3file,tsfile,zmfile,ochemfile,deepco2in character * 120 ghostfile @@ -100,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) @@ -115,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, @@ -130,8 +134,8 @@ & S0C_data,cfvolaer, & dirdat1,dirdat2 & ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss - & ,chem_init,chem_init2,chemrstfl - & ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av + & ,chem_init,chem_init2,chem_init4pfc,chemrstfl + & ,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 @@ -139,6 +143,7 @@ & ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3 & ,OCNGEOM,ocngmfile,ocndata4atm & ,bc_data + & ,PRTREND Cjrs & ,dtatm,dtocn #if ( defined IPCC_EMI ) & ,init_co2 @@ -161,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 @@ -189,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 @@ -253,7 +258,6 @@ ALFFOR=0. S0RATE=0. YEARGT=1958. - LYEAREM=2100 CO2IN=280. INYRAD=1. coefcl(1)=0. @@ -273,6 +277,7 @@ GSOEQ=.false. VEGCH=.false. TRVEG=.false. + PRTREND=.false. OCNGEOM=.false. do L=1,LM+1 ghostfv(L)=0. @@ -344,16 +349,15 @@ 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' 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. @@ -556,7 +560,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,8 +619,13 @@ else print *,'No diffusion into deep ocean' endif -!#if ( defined PREDICTED_AEROSOL ) -#if ( defined CPL_CHEM ) + 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 Cold SO2EREF=123.57 Cigsm1AFBYCF=0.6054 @@ -630,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 @@ -720,15 +731,16 @@ * status='OLD',form='unformatted') #endif #if ( defined CPL_TEM ) - open (367,file=last_nep,form='unformatted',status='new') - open (877,file=last_clm,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 #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='new') + open (277,file=fnememiss,form='unformatted',status='replace') + close(277) #endif #endif #if ( defined CPL_OCEANCO2 && defined ML_2D ) @@ -741,13 +753,14 @@ open( unit=547,file=nwrfl, * status='OLD',form='unformatted') #if ( defined CPL_TEM ) - open (367,file=last_nep,form='unformatted',status='old') - open (877,file=last_clm,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') +! 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 +1034,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