/[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.3 by jscott, Mon Apr 23 21:20:18 2007 UTC revision 1.8 by jscott, Thu Sep 17 15:48:38 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "ctrparam.h"  #include "ctrparam.h"
5    
# Line 44  C**** Line 46  C****
46  #include "AGRID.h"  #include "AGRID.h"
47  #endif  #endif
48    
49  #if ( defined CPL_TEM )  cjrs done in driver.h #if ( defined CPL_TEM )
50  #include "TEM.h"  cjrs#include "TEM.h"
51  #endif  cjrs#endif
52  !  !
53  ! === Chien Wang 062904  ! === Chien Wang 062904
54  !  !
# Line 77  c     CHARACTER*8 RECORD,ANDEND,NLREC*80 Line 79  c     CHARACTER*8 RECORD,ANDEND,NLREC*80
79        character *120 file1,file2,plotfl,nwrfl        character *120 file1,file2,plotfl,nwrfl
80        character * 120  t3file,tsfile,zmfile,qffile,clfile,wrcldf        character * 120  t3file,tsfile,zmfile,qffile,clfile,wrcldf
81       &    ,ochemfile,deepco2in       &    ,ochemfile,deepco2in
82       &    ,oco2file,co2rfile,caruptfile,flrco2av  cjrs DRIVER.h caruptfile  & ,oco2file,co2rfile,caruptfile,flrco2av
83         &    ,oco2file,co2rfile,flrco2av
84       &    ,ghg_monthly,ghg_monthly2,co2_data,o3_data       &    ,ghg_monthly,ghg_monthly2,co2_data,o3_data
85       &    ,bgrghg_data       &    ,bgrghg_data
86        character * 120 sulf1986,sulf2050,sulfamp,SO2_EM,        character * 120 sulf1986,sulf2050,sulfamp,SO2_EM,
87       &  S0C_data,       &  S0C_data,
88       &  dirdat1,dirdat2       &  dirdat1,dirdat2
89       &  ,bc_data       &  ,bc_data
90        character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo        character * 120 chemdata,chemout,init_4nem,pov_deepo
91       & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat  cjrs DRIVER.h last_nep character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo
92       &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,       & ,flin_nep,last_clm,SEN_dat
93       & chem_init,chem_init2,chemrstfl  cjrs DRIVER.h fnememiss &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,
94         &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,
95         & chem_init,chem_init2,chem_init4pfc,chemrstfl
96        common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf        common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf
97       *,t3file,tsfile,zmfile,ochemfile,deepco2in       *,t3file,tsfile,zmfile,ochemfile,deepco2in
98        character * 120 ghostfile        character * 120 ghostfile
# Line 100  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 115  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 130  c==== 012201 Line 134  c==== 012201
134       &  S0C_data,cfvolaer,       &  S0C_data,cfvolaer,
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,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 139  c==== 012201 Line 143  c==== 012201
143       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3
144       &  ,OCNGEOM,ocngmfile,ocndata4atm       &  ,OCNGEOM,ocngmfile,ocndata4atm
145       &  ,bc_data       &  ,bc_data
146         &  ,PRTREND
147  Cjrs     &  ,dtatm,dtocn  Cjrs     &  ,dtatm,dtocn
148  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
149       &  ,init_co2       &  ,init_co2
# Line 161  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 189  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 253  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 273  C Line 277  C
277        GSOEQ=.false.        GSOEQ=.false.
278        VEGCH=.false.        VEGCH=.false.
279        TRVEG=.false.        TRVEG=.false.
280          PRTREND=.false.
281        OCNGEOM=.false.        OCNGEOM=.false.
282        do L=1,LM+1        do L=1,LM+1
283          ghostfv(L)=0.          ghostfv(L)=0.
# Line 344  cjrs      dtocn=1 Line 349  cjrs      dtocn=1
349        chemdata = 'DATA'        chemdata = 'DATA'
350        chem_init = 'init-data_46x11_1991'        chem_init = 'init-data_46x11_1991'
351        chem_init2 = 'init-data2_46x11_1991'        chem_init2 = 'init-data2_46x11_1991'
352          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 556  c Line 560  c
560  c     open file for carbon uptake  c     open file for carbon uptake
561  #if ( defined CPL_TEM || defined CPL_OCEANCO2 )  #if ( defined CPL_TEM || defined CPL_OCEANCO2 )
562  c     open(333,file=caruptfile,status='new',form='formatted')  c     open(333,file=caruptfile,status='new',form='formatted')
563        open(333,file=caruptfile,form='formatted')        open(333,file=caruptfile,status='replace',form='formatted')
564          close(333)
565  #endif  #endif
566    
567  #if ( defined CPL_OCEANCO2 && defined ML_2D)  #if ( defined CPL_OCEANCO2 && defined ML_2D)
# Line 614  C       OLD Line 619  C       OLD
619        else        else
620          print *,'No diffusion into deep ocean'          print *,'No diffusion into deep ocean'
621        endif        endif
622  !#if ( defined PREDICTED_AEROSOL )        if (PRTREND) then
623  #if ( defined CPL_CHEM )           print *,' With changes in  stochastic precip'
624          else
625             print *,' Without changes in  stochastic precip'
626          endif
627    #if ( defined PREDICTED_AEROSOL )
628    !#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 630  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 720  c     LMCM=LTM+2 Line 731  c     LMCM=LTM+2
731       *    status='OLD',form='unformatted')       *    status='OLD',form='unformatted')
732  #endif  #endif
733  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
734            open (367,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='new')            open (277,file=fnememiss,form='unformatted',status='replace')
743              close(277)
744  #endif  #endif
745  #endif  #endif
746  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
# Line 741  C      For restart of the run Line 753  C      For restart of the run
753            open( unit=547,file=nwrfl,            open( unit=547,file=nwrfl,
754       *    status='OLD',form='unformatted')       *    status='OLD',form='unformatted')
755  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
756            open (367,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='new')            open (277,file=fnememiss,form='unformatted',status='replace')
763              close(277)
764  #endif  #endif
765  #endif  #endif
766  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
# Line 1021  c     print *,' NWRGHG=',NWRGHG Line 1034  c     print *,' NWRGHG=',NWRGHG
1034         print *,TSCNTR,' form EXP=',EXPTSF         print *,TSCNTR,' form EXP=',EXPTSF
1035  !      read(576)EXPTSF  !      read(576)EXPTSF
1036  !      print *,' TSURF form EXP=',EXPTSF  !      print *,' TSURF form EXP=',EXPTSF
1037         read(576)TSURFC,TLANDC  #if ( defined CLM )
1038            read(576)TSURFC,TLANDC
1039    #else
1040            read(576)TSURFC
1041    #endif
1042  !      read(576)TLANDC  !      read(576)TLANDC
1043        endif        endif
1044          CORSR=1.0          CORSR=1.0
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22