/[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.2 by jscott, Tue Aug 22 20:25:52 2006 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.COM"  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,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,
123       &  LYEAREM,       *  AERFOR,AERF4BC,
      *  AERFOR,  
124       *  S0RATE,CFS0X,       *  S0RATE,CFS0X,
125       *  CFAEROSOL,CFBC,       *  CFAEROSOL,CFBC,
126       &  cfocdif,rkv,diffcar0,       &  cfocdif,rkv,diffcar0,ocarcont,ocarindata,
127  !       Kvc=diffcar0+cfocdif*Kvh  !       Kvc=diffcar0+cfocdif*Kvh
128       *  file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf,clmsen,cfdif0,       *  file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf,clmsen,cfdif0,
129       *  t3file,tsfile,zmfile,ochemfile,deepco2in,       *  t3file,tsfile,zmfile,ochemfile,deepco2in,
130       &  fl_init_alkt,fl_init_salt,       &  fl_init_alkt,fl_init_salt,fl_dic_eq,
131       *  ghg_monthly,ghg_monthly2,co2_data,o3_data,       *  ghg_monthly,ghg_monthly2,co2_data,o3_data,
132       &    bgrghg_data,       &    bgrghg_data,
133       *  sulf1986,sulf2050,sulfamp,SO2_EM,       *  sulf1986,sulf2050,sulfamp,SO2_EM,
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 148  Cjrs     &  ,dtatm,dtocn Line 153  Cjrs     &  ,dtatm,dtocn
153        character * 120 vegfile        character * 120 vegfile
154       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir
155        common/wrcom/wr25,TRANSR,CONTRR,OBSFOR        common/wrcom/wr25,TRANSR,CONTRR,OBSFOR
156  c jrs     common/TIMESTEPS/dtatm,dtocn  c jrs      common/TIMESTEPS/dtatm,dtocn
157        LOGICAL LFR,NLFR,wr25,TRANSR,WRCLD,CONTRR,CLDFEED,OBSFOR        LOGICAL LFR,NLFR,wr25,TRANSR,WRCLD,CONTRR,CLDFEED,OBSFOR
158       &,GHSF,VEGCH,TRVEG,GSOEQ,OCNGEOM,GHSFALB,STRARFOR,CO2FOR       &,GHSF,VEGCH,TRVEG,GSOEQ,OCNGEOM,GHSFALB,STRARFOR,CO2FOR
159       & ,FORSULF,FORBC,S0FOR,FORVOL       & ,FORSULF,FORBC,S0FOR,FORVOL
# 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 171  c jrs     common/TIMESTEPS/dtatm,dtocn Line 175  c jrs     common/TIMESTEPS/dtatm,dtocn
175        common/atmos_lo/fland_atm(jm0)  !jrs not sure this does anything        common/atmos_lo/fland_atm(jm0)  !jrs not sure this does anything
176    
177  #if ( defined CLM )  #if ( defined CLM )
178  #include "CLM.COM"  cjrs alreadyin DRIVER.h#include "CLM.COM"
179        dimension clmlice(jm0),baresoil(jm0),        dimension clmlice(jm0),baresoil(jm0),
180       &  w1maxclm(jm0),w2maxclm(jm0),vmaskclm(jm0)       &  w1maxclm(jm0),w2maxclm(jm0),vmaskclm(jm0)
181        character * 120 lineclm        character * 120 lineclm
182  #endif  #endif
183  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
184        common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0)        common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0)
185    !     common/Garydiff/depthml(jm0),edzon(jm0,lmo),dzg(lmo),dzog(lmo-1),
186    !    &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)
187        common/Garydiff/depthml(jm0),edzon(jm0),dzg(lmo),dzog(lmo-1),        common/Garydiff/depthml(jm0),edzon(jm0),dzg(lmo),dzog(lmo-1),
188       &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)       &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)
189        common /Garychem/Hg(jm0)        common /Garychem/Hg(jm0)
190        common /Garyvdif/iyearocm,vdfocm,acvdfc,cfocdif,diffcar0        common /Garyvdif/iyearocm,vdfocm,acvdfc,cfocdif,diffcar0
191          common /Garyvlog/odifcarbon,ocarcont
192        real  Rco2in(jm0,lmo),Hgin(jm0)        real  Rco2in(jm0,lmo),Hgin(jm0)
193          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 192  c jrs     common/TIMESTEPS/dtatm,dtocn Line 201  c jrs     common/TIMESTEPS/dtatm,dtocn
201  #endif  #endif
202    
203        integer PCLOUD        integer PCLOUD
204        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TLANDD(JM0),TSURFD(JM0),DTSURF(JM0)
205    !     common/TSLD/TLANDC(JM0,0:13),TLANDT(JM0),TLANDD(JM0),DTLAND(JM0)
206    #include "TSRF.COM"
207        DATA DISK/'DISK'/,ANDEND/' &END   '/                              1520.          DATA DISK/'DISK'/,ANDEND/' &END   '/                              1520.  
208        DATA IDAYS0/0,1,32,60,91,121,152,182,213,244,274,305,335/         1521.          DATA IDAYS0/0,1,32,60,91,121,152,182,213,244,274,305,335/         1521.  
209        DATA NAMD60/'AUSD','MWST','SAHL','EPAC'/                          1522.          DATA NAMD60/'AUSD','MWST','SAHL','EPAC'/                          1522.  
# Line 247  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 267  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 277  C Line 288  C
288        cfocdif=0.42        cfocdif=0.42
289        cfocdif=1.375        cfocdif=1.375
290        cfocdif=0.6        cfocdif=0.6
291          cfocdif=3.0
292        diffcar0=2.85        diffcar0=2.85
293          diffcar0=1.00
294  cc    ALFA = 8.0*1.e3  cc    ALFA = 8.0*1.e3
295          AERF4BC=-0.35
296        TRANSR=.false.        TRANSR=.false.
297        WRCLD=.FALSE.        WRCLD=.FALSE.
298        clfile = 'undefined'        clfile = 'undefined'
# Line 301  cc    ALFA = 8.0*1.e3 Line 315  cc    ALFA = 8.0*1.e3
315        ISTWRC=0        ISTWRC=0
316        QFCOR=.FALSE.        QFCOR=.FALSE.
317        READGHG=0.        READGHG=0.
318  C      dtatm=1  cjrs      dtatm=1
319  C      dtocn=1  cjrs      dtocn=1
320        ghg_monthly =  'undefined'        ghg_monthly =  'undefined'
321        ghg_monthly2 =  'undefined'        ghg_monthly2 =  'undefined'
322        co2_data =  'undefined'        co2_data =  'undefined'
# Line 314  C      dtocn=1 Line 328  C      dtocn=1
328        deepco2in =  'undefined'        deepco2in =  'undefined'
329        fl_init_alkt = 'undefined'        fl_init_alkt = 'undefined'
330        fl_init_salt = 'undefined'        fl_init_salt = 'undefined'
331          fl_dic_eq = 'undefined'
332          ocarcont=.true.
333          ocarindata=.false.
334        vegfile =  'undefined'        vegfile =  'undefined'
335        fclmlice =  'undefined'        fclmlice =  'undefined'
336        fprratio =  'undefined'        fprratio =  'undefined'
# Line 332  C      dtocn=1 Line 349  C      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 454  C**** COPY INPUTZ NAMELIST ONTO CORE TAP Line 470  C**** COPY INPUTZ NAMELIST ONTO CORE TAP
470        IF(RECORD(1).NE.ANDEND) GO TO 60                                  1603.          IF(RECORD(1).NE.ANDEND) GO TO 60                                  1603.  
471        rewind 514        rewind 514
472        read (UNIT=514,NML=INPUTZ)        read (UNIT=514,NML=INPUTZ)
473        REWIND 514                                                        1606.        REWIND 514                                                        1606.  
474  C JRS ignore name.dat values, start Jan 1. with couple.nml years    C JRS ignore name.dat values, start Jan 1. with couple.nml years  
475        inyear = startYear        inyear = startYear
476        lyear = endYear +1        lyear = endYear +1
# Line 462  C JRS ignore name.dat values, start Jan Line 478  C JRS ignore name.dat values, start Jan
478        IMONTH   = 1        IMONTH   = 1
479        LDAY     = 1        LDAY     = 1
480        LMONTH   = 1        LMONTH   = 1
   
481  Cjrs      dtatmo=dtatm  Cjrs      dtatmo=dtatm
482  C      dtocno=dtocn  cjrs      dtocno=dtocn
483  cb    open statments  cb    open statments
484  c  c
485  c     File which depend on resolution  c     File which depend on resolution
# Line 516  c Line 531  c
531        endif        endif
532  #if ( !defined CPL_CHEM )  #if ( !defined CPL_CHEM )
533  #if ( defined PREDICTED_BC)  #if ( defined PREDICTED_BC)
534          if(READGHG.eq.0)then
535    !      data for BC only
536         print *,'Data for black carbon'         print *,'Data for black carbon'
537        open(769,file=bc_data,         open(769,file=bc_data,
538       &  status='old',form='unformatted')       &  status='old',form='unformatted')
539          endif
540  #endif  #endif
541  #endif  #endif
542    
# Line 542  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 550  c     open(333,file=caruptfile,status='n Line 569  c     open(333,file=caruptfile,status='n
569       &   form='unformatted',status='old')       &   form='unformatted',status='old')
570        open(669,file=fl_init_salt,        open(669,file=fl_init_salt,
571       &   form='unformatted',status='old')       &   form='unformatted',status='old')
572    !     open(670,file=fl_dic_eq,
573    !    &   form='unformatted',status='old')
574        open(602,file=flrco2av,status='new',form='unformatted')        open(602,file=flrco2av,status='new',form='unformatted')
575  #endif  #endif
576    
577  ce    open statments  ce    open statments
 !#if ( defined CPL_CHEM )  
 !       YEARGT=1977.  
 !#endif  
 !     if(YEARGT.eq.1860)then  
 !        GHGBGR(1)=GHGBGR1860(1)  
 !        GHGBGR(2)=GHGBGR1860(2)  
 !        GHGBGR(3)=GHGBGR1860(3)  
 !        GHGBGR(4)=GHGBGR1860(4)  
 !        GHGBGR(5)=GHGBGR1860(5)  
 !     else if (YEARGT.eq.1958)then  
 !        GHGBGR(1)=GHGBGR1958(1)  
 !        GHGBGR(2)=GHGBGR1958(2)  
 !        GHGBGR(3)=GHGBGR1958(3)  
 !        GHGBGR(4)=GHGBGR1958(4)  
 !        GHGBGR(5)=GHGBGR1958(5)  
 !     else if (YEARGT.eq.1977)then  
 !        GHGBGR(1)=GHGBGR1977(1)  
 !        GHGBGR(2)=GHGBGR1977(2)  
 !        GHGBGR(3)=GHGBGR1977(3)  
 !        GHGBGR(4)=GHGBGR1977(4)  
 !        GHGBGR(5)=GHGBGR1977(5)  
 !     else if (YEARGT.eq.1980)then  
 !        GHGBGR(1)=GHGBGR1980(1)  
 !        GHGBGR(2)=GHGBGR1980(2)  
 !        GHGBGR(3)=GHGBGR1980(3)  
 !        GHGBGR(4)=GHGBGR1980(4)  
 !        GHGBGR(5)=GHGBGR1980(5)  
 !     else if (YEARGT.eq.1990)then  
 !        GHGBGR(1)=GHGBGR1990(1)  
 !        GHGBGR(2)=GHGBGR1990(2)  
 !        GHGBGR(3)=GHGBGR1990(3)  
 !        GHGBGR(4)=GHGBGR1990(4)  
 !        GHGBGR(5)=GHGBGR1990(5)  
 !     else if (YEARGT.eq.1991)then  
 !        GHGBGR(1)=GHGBGR1991(1)  
 !        GHGBGR(2)=GHGBGR1991(2)  
 !        GHGBGR(3)=GHGBGR1991(3)  
 !        GHGBGR(4)=GHGBGR1991(4)  
 !        GHGBGR(5)=GHGBGR1991(5)  
 !     else if (YEARGT.eq.2000)then  
 !        GHGBGR(1)=GHGBGR2000(1)  
 !        GHGBGR(2)=GHGBGR2000(2)  
 !        GHGBGR(3)=GHGBGR2000(3)  
 !        GHGBGR(4)=GHGBGR2000(4)  
 !        GHGBGR(5)=GHGBGR2000(5)  
 !     else if (YEARGT.eq.1765)then  
 !        GHGBGR(1)=GHGBGR1765(1)  
 !         if (CO2.le.-20.0) GHGBGR(1)=277.6  
 !        GHGBGR(2)=GHGBGR1765(2)  
 !        GHGBGR(3)=GHGBGR1765(3)  
 !        GHGBGR(4)=GHGBGR1765(4)  
 !        GHGBGR(5)=GHGBGR1765(5)  
 !     else if (YEARGT.eq.-1765.)then  
 !        GHGBGR(1)=4.*GHGBGR1765(1)  
 !        GHGBGR(2)=GHGBGR1765(2)  
 !        GHGBGR(3)=GHGBGR1765(3)  
 !        GHGBGR(4)=GHGBGR1765(4)  
 !        GHGBGR(5)=GHGBGR1765(5)  
 !     else  
 !        print *,' Wrong YEARGT ', YEARGT  
 !        stop  
 !     endif  
578        call bgrghg(YEARGT)        call bgrghg(YEARGT)
579  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
580        if(YEARGT.eq.1765)then  !     if(YEARGT.eq.1765)then
581         GHGBGR(1)=277.6  !      GHGBGR(1)=277.6
582         open (unit=861,file=init_co2,         open (unit=861,file=init_co2,
583       &     status='OLD',form='formatted')       &     status='OLD',form='formatted')
584         read (861,*)xco2init         read (861,*)xco2init
585         CO2=xco2init/GHGBGR(1)         CO2=xco2init/GHGBGR(1)
586         print *,'IPCC EMI CO2=',CO2         print *,'IPCC EMI CO2=',CO2
587        else  !     else
588           print *,' Wrong YEARGT ', YEARGT  !        print *,' Wrong YEARGT ', YEARGT
589           stop  !        stop
590        endif  !     endif
591  #endif  #endif
592        print *,'Background GHGs for year ',YEARGT        print *,'Background GHGs for year ',YEARGT
593        print '(5E12.4)',GHGBGR        print '(5E12.4)',GHGBGR
# Line 660  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 (PRTREND) then
623             print *,' With changes in  stochastic precip'
624          else
625             print *,' Without changes in  stochastic precip'
626          endif
627  #if ( defined PREDICTED_AEROSOL )  #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 674  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/ Line 639  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/
639  !      CFAEROSOL=(-AERFOR/AFBYCF)**1.035/(SO2ER**1.0391)  !      CFAEROSOL=(-AERFOR/AFBYCF)**1.035/(SO2ER**1.0391)
640    
641    
642        SO2EREF=149.07        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
651        AFBYCF90=1.043  !     AFBYCF90=1.04
652    !     print *,'AFBYCF90=',AFBYCF90
653    !     12/21/2006 (from runs 265x.06)
654    !     Adjustment for BC forcing
655           print *,' AERF4BC=',AERF4BC
656          AERFOR=AERFOR+AERF4BC
657    !     Adjustment for BC forcing
658          AFBYCF=1.05
659        print *,'SO2ER=',SO2ER        print *,'SO2ER=',SO2ER
660        print *,'AFBYCF90=',AFBYCF90        print *,'AFBYCF=',AFBYCF
661        CFAEROSOL=(-AERFOR90/AFBYCF90)**1.119        CFAEROSOL=(-AERFOR/AFBYCF)**1.21
662        CFAEROSOL=CFAEROSOL/(SO2ER**1.0391)        CFAEROSOL=CFAEROSOL/(SO2ER**1.01)
663        print *,'AERFOR90=',AERFOR90,'CFAEROSOL=',CFAEROSOL        print *,'AERFOR=',AERFOR,'CFAEROSOL=',CFAEROSOL
664  #if ( defined OCEAN_3D )  #if ( defined OCEAN_3D )
665         CFAEROSOL=CFAEROSOL/1.35         CFAEROSOL=CFAEROSOL/1.35
666         print *,'CFAEROSOL_3D=',CFAEROSOL         print *,'CFAEROSOL_3D=',CFAEROSOL
# Line 758  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 779  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 824  c      print *,' Form input ' Line 799  c      print *,' Form input '
799          READ (KDISK0,ERR=840)AEXPX,TAUX,JC1,C1,RC1,KEYNR,U,V,T,P,Q,          READ (KDISK0,ERR=840)AEXPX,TAUX,JC1,C1,RC1,KEYNR,U,V,T,P,Q,
800       &  ODATA,       &  ODATA,
801       *  GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,   1739.         *  GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,   1739.  
802       *  CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0       *  CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0
803         else         else
804  !       READ (KDISK0,ERR=840)AEXPX,TAUX,JC1,C1,RC1,KEYNR,U,V,T,P,Q,  !       READ (KDISK0,ERR=840)AEXPX,TAUX,JC1,C1,RC1,KEYNR,U,V,T,P,Q,
805          READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,          READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,
806       &  ODATA,       &  ODATA,
807       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.
808       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0,       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0,
809       *   TG3M,RTGO,STG3,DTG3       *   TG3M,RTGO,STG3,DTG3
810          endif          endif
811         if(.not.CONTRR)then         if(.not.CONTRR)then
# Line 879  c       print *,' WITH DEEP MIXED LAYER' Line 854  c       print *,' WITH DEEP MIXED LAYER'
854  #endif  #endif
855    
856  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
857           if(ocarcont) then
858            if(ocarindata)then
859              print *,'Wrong setting of ocarcont and ocarindata'
860              print *,ocarcont,ocarindata
861              stop
862            endif
863          open(116,file=deepco2in,          open(116,file=deepco2in,
864       * status='old',form='unformatted')       *   status='old',form='unformatted')
865          print *,' AFTER OPEN INIT. data for ocean chem.'          print *,' AFTER OPEN INIT. data for ocean chem.'
866          print *,deepco2in          print *,deepco2in
867          read(116)iyearocm,vdfocm          read(116)iyearocm,vdfocm
# Line 900  c       print *,' WITH DEEP MIXED LAYER' Line 881  c       print *,' WITH DEEP MIXED LAYER'
881           Rco2(j,k)=Rco2in(j,k)           Rco2(j,k)=Rco2in(j,k)
882          end do          end do
883          end do          end do
884           else
885            if(ocarindata)then
886             print *,' Reading initial data for ocean carbon'
887             open(116,file=deepco2in,
888         *    status='old',form='unformatted')
889             read(116)
890             read(116)Hgin
891             read(116)Rco2in
892             print *,(Rco2in(j,1),j=1,jm)
893             do k=1,lmo
894             do j=1,jm0
895    !         if(k.eq.1)Hg(j)=Hgin(j)
896              if(k.eq.1)Hg(j)=1.0e-8
897              Rco2(j,k)=Rco2in(j,k)
898             end do
899             end do
900             close(116)
901            else
902             do k=1,lmo
903             do j=1,jm0
904              if(k.eq.1)Hg(j)=1.0e-8
905              Rco2(j,k)=0.0
906             end do
907             end do
908            endif
909           endif
910  #endif  #endif
911    
912        else        else
# Line 909  C ******* Line 916  C *******
916          if(TRANSR)then          if(TRANSR)then
917           READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,ODATA, 1738.           READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,ODATA, 1738.
918       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.
919       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0,       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0,
920       *   TG3M,RTGO,STG3,DTG3       *   TG3M,RTGO,STG3,DTG3
921          else          else
922           READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,ODATA, 1738.           READ (KDISK0,ERR=840)AEXPX,TAUX,JC,C,RC,KEYNR,U,V,T,P,Q,ODATA, 1738.
923       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.       *   GDATA,BLDATA,RQT,SRHR,TRHR,(AJ(K,1),K=1,KACC),TAUY,TSSFC,CKS,  1739.
924       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0
925          endif          endif
926  C  C
927           if(abs(AEXPX-AEXP).gt.0.05)then           if(abs(AEXPX-AEXP).gt.0.05)then
# Line 943  c       if(iyearocm.ne.JYEAR) then Line 950  c       if(iyearocm.ne.JYEAR) then
950          endif          endif
951          read(369)Hgin          read(369)Hgin
952          read(369)Rco2in          read(369)Rco2in
953            print *,(Rco2in(j,1),j=1,jm)
954          do k=1,lmo          do k=1,lmo
955          do j=1,jm0          do j=1,jm0
956           if(k.eq.1)Hg(j)=Hgin(j)           if(k.eq.1)Hg(j)=Hgin(j)
# Line 952  c       if(iyearocm.ne.JYEAR) then Line 960  c       if(iyearocm.ne.JYEAR) then
960  #endif  #endif
961    
962        endif        ! endif for ISTRT1        endif        ! endif for ISTRT1
963    #if ( defined CPL_OCEANCO2 && defined ML_2D )
964           print *,'Rco2 from input'
965    !        print *,(Rco2(j,1),j=1,jm)
966             print *,Rco2
967    #endif
968    
969        if(KOCEAN.eq.1) then        if(KOCEAN.eq.1) then
970         print *,' T1 ocean'         print *,' T1 ocean'
# Line 980  c     endif Line 993  c     endif
993        TAUP=TAUX                                                         1746.          TAUP=TAUX                                                         1746.  
994  C**** UPDATE C ARRAY FROM INPUTZ                                        1747.    C**** UPDATE C ARRAY FROM INPUTZ                                        1747.  
995    500 READ (514,INPUTZ)                                                 1748.      500 READ (514,INPUTZ)                                                 1748.  
996    !      nrad=NDYN
997  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
998         CO2=xco2init/GHGBGR(1)         CO2=xco2init/GHGBGR(1)
999  #endif  #endif
# Line 1020  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  #if ( defined CLM )
1038            read(576)TSURFC,TLANDC
1039    #else
1040            read(576)TSURFC
1041    #endif
1042    !      read(576)TLANDC
1043        endif        endif
1044          CORSR=1.0          CORSR=1.0
1045        if(KOCEAN.eq.0)then        if(KOCEAN.eq.0)then
# Line 1063  c     PRINT *,' NLFR=',NLFR Line 1082  c     PRINT *,' NLFR=',NLFR
1082        if(READGHG.eq.2)then        if(READGHG.eq.2)then
1083          open( unit=569,file=ghg_monthly,          open( unit=569,file=ghg_monthly,
1084       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1085          open( unit=669,file=ghg_monthly2,          open( unit=679,file=ghg_monthly2,
1086       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1087          print *,' GHGs from ',ghg_monthly          print *,' GHGs from ',ghg_monthly
1088          if(ISTRT1.eq.1)then          if(ISTRT1.eq.1)then
1089            do 369 i=1,NWRGHG            do 369 i=1,NWRGHG
1090             do 369 ii=1,13             do 369 ii=1,13
1091             read(569)             read(569)
1092             if(ii.le.3)read(669)             if(ii.le.3)read(679)
1093    369     continue    369     continue
1094        endif        endif
1095        endif        endif
# Line 1079  c       open( unit=569,file=dirdat1(1:id Line 1098  c       open( unit=569,file=dirdat1(1:id
1098  c    *  status='OLD',form='unformatted')  c    *  status='OLD',form='unformatted')
1099          open( unit=569,file=ghg_monthly,          open( unit=569,file=ghg_monthly,
1100       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1101          open( unit=669,file=ghg_monthly2,          open( unit=679,file=ghg_monthly2,
1102       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1103          print *,' GHGs from ',ghg_monthly          print *,' GHGs from ',ghg_monthly
1104        endif        endif
# Line 1287  c          IF(FDATA(I,J,2).LT.0.01)FDATA Line 1306  c          IF(FDATA(I,J,2).LT.0.01)FDATA
1306          fland_temp(j)=FDATA(1,J,2)          fland_temp(j)=FDATA(1,J,2)
1307    284 continue    284 continue
1308  #if( defined OCEAN_3D)  #if( defined OCEAN_3D)
1309  Cjrs       if(jmocean.ne.jm0-2)then  Cjrs        if(jmocean.ne.jm0-2)then
1310  C          print *,"Wrong jm or jmocean"  C          print *,"Wrong jm or jmocean"
1311  C          stop  C          stop
1312  C       endif  C       endif

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8

  ViewVC Help
Powered by ViewVC 1.1.22