/[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.3 by jscott, Mon Apr 23 21:20:18 2007 UTC
# Line 45  C**** Line 45  C****
45  #endif  #endif
46    
47  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
48  #include "TEM.COM"  #include "TEM.h"
49  #endif  #endif
50  !  !
51  ! === Chien Wang 062904  ! === Chien Wang 062904
# Line 86  c     CHARACTER*8 RECORD,ANDEND,NLREC*80 Line 86  c     CHARACTER*8 RECORD,ANDEND,NLREC*80
86       &  ,bc_data       &  ,bc_data
87        character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo        character * 120 chemdata,chemout,last_nep,init_4nem,pov_deepo
88       & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat       & ,flin_nep,last_clm,emiss_data,SO2ERATIO,SEN_dat
89       &  ,fl_init_alkt,fl_init_salt,fnememiss,       &  ,fl_init_alkt,fl_init_salt,fl_dic_eq,fnememiss,
90       & chem_init,chem_init2,chemrstfl       & chem_init,chem_init2,chemrstfl
91        common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf        common/files/file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf
92       *,t3file,tsfile,zmfile,ochemfile,deepco2in       *,t3file,tsfile,zmfile,ochemfile,deepco2in
# Line 116  c==== 012201 Line 116  c==== 012201
116       *  READGHG,wr25,LFR,ISTRT1,PCLOUD,QFCOR,TRANSR,WRCLD,NWRCLD,CONTRR,       *  READGHG,wr25,LFR,ISTRT1,PCLOUD,QFCOR,TRANSR,WRCLD,NWRCLD,CONTRR,
117       *  ISTWRC,CLDFEED,OBSFOR,ALFFOR,YEARGT,CO2IN,ISTRTCHEM,       *  ISTWRC,CLDFEED,OBSFOR,ALFFOR,YEARGT,CO2IN,ISTRTCHEM,
118       &  LYEAREM,       &  LYEAREM,
119       *  AERFOR,       *  AERFOR,AERF4BC,
120       *  S0RATE,CFS0X,       *  S0RATE,CFS0X,
121       *  CFAEROSOL,CFBC,       *  CFAEROSOL,CFBC,
122       &  cfocdif,rkv,diffcar0,       &  cfocdif,rkv,diffcar0,ocarcont,ocarindata,
123  !       Kvc=diffcar0+cfocdif*Kvh  !       Kvc=diffcar0+cfocdif*Kvh
124       *  file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf,clmsen,cfdif0,       *  file1,file2,plotfl,nwrfl,qffile,clfile,wrcldf,clmsen,cfdif0,
125       *  t3file,tsfile,zmfile,ochemfile,deepco2in,       *  t3file,tsfile,zmfile,ochemfile,deepco2in,
126       &  fl_init_alkt,fl_init_salt,       &  fl_init_alkt,fl_init_salt,fl_dic_eq,
127       *  ghg_monthly,ghg_monthly2,co2_data,o3_data,       *  ghg_monthly,ghg_monthly2,co2_data,o3_data,
128       &    bgrghg_data,       &    bgrghg_data,
129       *  sulf1986,sulf2050,sulfamp,SO2_EM,       *  sulf1986,sulf2050,sulfamp,SO2_EM,
# Line 148  Cjrs     &  ,dtatm,dtocn Line 148  Cjrs     &  ,dtatm,dtocn
148        character * 120 vegfile        character * 120 vegfile
149       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir
150        common/wrcom/wr25,TRANSR,CONTRR,OBSFOR        common/wrcom/wr25,TRANSR,CONTRR,OBSFOR
151  c jrs     common/TIMESTEPS/dtatm,dtocn  c jrs      common/TIMESTEPS/dtatm,dtocn
152        LOGICAL LFR,NLFR,wr25,TRANSR,WRCLD,CONTRR,CLDFEED,OBSFOR        LOGICAL LFR,NLFR,wr25,TRANSR,WRCLD,CONTRR,CLDFEED,OBSFOR
153       &,GHSF,VEGCH,TRVEG,GSOEQ,OCNGEOM,GHSFALB,STRARFOR,CO2FOR       &,GHSF,VEGCH,TRVEG,GSOEQ,OCNGEOM,GHSFALB,STRARFOR,CO2FOR
154       & ,FORSULF,FORBC,S0FOR,FORVOL       & ,FORSULF,FORBC,S0FOR,FORVOL
# Line 171  c jrs     common/TIMESTEPS/dtatm,dtocn Line 171  c jrs     common/TIMESTEPS/dtatm,dtocn
171        common/atmos_lo/fland_atm(jm0)  !jrs not sure this does anything        common/atmos_lo/fland_atm(jm0)  !jrs not sure this does anything
172    
173  #if ( defined CLM )  #if ( defined CLM )
174  #include "CLM.COM"  cjrs alreadyin DRIVER.h#include "CLM.COM"
175        dimension clmlice(jm0),baresoil(jm0),        dimension clmlice(jm0),baresoil(jm0),
176       &  w1maxclm(jm0),w2maxclm(jm0),vmaskclm(jm0)       &  w1maxclm(jm0),w2maxclm(jm0),vmaskclm(jm0)
177        character * 120 lineclm        character * 120 lineclm
178  #endif  #endif
179  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
180        common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0)        common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0)
181    !     common/Garydiff/depthml(jm0),edzon(jm0,lmo),dzg(lmo),dzog(lmo-1),
182    !    &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)
183        common/Garydiff/depthml(jm0),edzon(jm0),dzg(lmo),dzog(lmo-1),        common/Garydiff/depthml(jm0),edzon(jm0),dzg(lmo),dzog(lmo-1),
184       &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)       &Rco2(jm0,lmo),edohd(lmo),zg(lmo),focean(jm0)
185        common /Garychem/Hg(jm0)        common /Garychem/Hg(jm0)
186        common /Garyvdif/iyearocm,vdfocm,acvdfc,cfocdif,diffcar0        common /Garyvdif/iyearocm,vdfocm,acvdfc,cfocdif,diffcar0
187          common /Garyvlog/odifcarbon,ocarcont
188        real  Rco2in(jm0,lmo),Hgin(jm0)        real  Rco2in(jm0,lmo),Hgin(jm0)
189          logical odifcarbon,ocarcont,ocarindata
190  #endif  #endif
191    
192  #if (!defined PREDICTED_GASES)                                                #if (!defined PREDICTED_GASES)                                              
# Line 192  c jrs     common/TIMESTEPS/dtatm,dtocn Line 196  c jrs     common/TIMESTEPS/dtatm,dtocn
196  #endif  #endif
197    
198        integer PCLOUD        integer PCLOUD
199        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TLANDD(JM0),TSURFD(JM0),DTSURF(JM0)
200    !     common/TSLD/TLANDC(JM0,0:13),TLANDT(JM0),TLANDD(JM0),DTLAND(JM0)
201    #include "TSRF.COM"
202        DATA DISK/'DISK'/,ANDEND/' &END   '/                              1520.          DATA DISK/'DISK'/,ANDEND/' &END   '/                              1520.  
203        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.  
204        DATA NAMD60/'AUSD','MWST','SAHL','EPAC'/                          1522.          DATA NAMD60/'AUSD','MWST','SAHL','EPAC'/                          1522.  
# Line 277  C Line 283  C
283        cfocdif=0.42        cfocdif=0.42
284        cfocdif=1.375        cfocdif=1.375
285        cfocdif=0.6        cfocdif=0.6
286          cfocdif=3.0
287        diffcar0=2.85        diffcar0=2.85
288          diffcar0=1.00
289  cc    ALFA = 8.0*1.e3  cc    ALFA = 8.0*1.e3
290          AERF4BC=-0.35
291        TRANSR=.false.        TRANSR=.false.
292        WRCLD=.FALSE.        WRCLD=.FALSE.
293        clfile = 'undefined'        clfile = 'undefined'
# Line 301  cc    ALFA = 8.0*1.e3 Line 310  cc    ALFA = 8.0*1.e3
310        ISTWRC=0        ISTWRC=0
311        QFCOR=.FALSE.        QFCOR=.FALSE.
312        READGHG=0.        READGHG=0.
313  C      dtatm=1  cjrs      dtatm=1
314  C      dtocn=1  cjrs      dtocn=1
315        ghg_monthly =  'undefined'        ghg_monthly =  'undefined'
316        ghg_monthly2 =  'undefined'        ghg_monthly2 =  'undefined'
317        co2_data =  'undefined'        co2_data =  'undefined'
# Line 314  C      dtocn=1 Line 323  C      dtocn=1
323        deepco2in =  'undefined'        deepco2in =  'undefined'
324        fl_init_alkt = 'undefined'        fl_init_alkt = 'undefined'
325        fl_init_salt = 'undefined'        fl_init_salt = 'undefined'
326          fl_dic_eq = 'undefined'
327          ocarcont=.true.
328          ocarindata=.false.
329        vegfile =  'undefined'        vegfile =  'undefined'
330        fclmlice =  'undefined'        fclmlice =  'undefined'
331        fprratio =  'undefined'        fprratio =  'undefined'
# Line 454  C**** COPY INPUTZ NAMELIST ONTO CORE TAP Line 466  C**** COPY INPUTZ NAMELIST ONTO CORE TAP
466        IF(RECORD(1).NE.ANDEND) GO TO 60                                  1603.          IF(RECORD(1).NE.ANDEND) GO TO 60                                  1603.  
467        rewind 514        rewind 514
468        read (UNIT=514,NML=INPUTZ)        read (UNIT=514,NML=INPUTZ)
469        REWIND 514                                                        1606.        REWIND 514                                                        1606.  
470  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  
471        inyear = startYear        inyear = startYear
472        lyear = endYear +1        lyear = endYear +1
# Line 462  C JRS ignore name.dat values, start Jan Line 474  C JRS ignore name.dat values, start Jan
474        IMONTH   = 1        IMONTH   = 1
475        LDAY     = 1        LDAY     = 1
476        LMONTH   = 1        LMONTH   = 1
   
477  Cjrs      dtatmo=dtatm  Cjrs      dtatmo=dtatm
478  C      dtocno=dtocn  cjrs      dtocno=dtocn
479  cb    open statments  cb    open statments
480  c  c
481  c     File which depend on resolution  c     File which depend on resolution
# Line 516  c Line 527  c
527        endif        endif
528  #if ( !defined CPL_CHEM )  #if ( !defined CPL_CHEM )
529  #if ( defined PREDICTED_BC)  #if ( defined PREDICTED_BC)
530          if(READGHG.eq.0)then
531    !      data for BC only
532         print *,'Data for black carbon'         print *,'Data for black carbon'
533        open(769,file=bc_data,         open(769,file=bc_data,
534       &  status='old',form='unformatted')       &  status='old',form='unformatted')
535          endif
536  #endif  #endif
537  #endif  #endif
538    
# Line 550  c     open(333,file=caruptfile,status='n Line 564  c     open(333,file=caruptfile,status='n
564       &   form='unformatted',status='old')       &   form='unformatted',status='old')
565        open(669,file=fl_init_salt,        open(669,file=fl_init_salt,
566       &   form='unformatted',status='old')       &   form='unformatted',status='old')
567    !     open(670,file=fl_dic_eq,
568    !    &   form='unformatted',status='old')
569        open(602,file=flrco2av,status='new',form='unformatted')        open(602,file=flrco2av,status='new',form='unformatted')
570  #endif  #endif
571    
572  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  
573        call bgrghg(YEARGT)        call bgrghg(YEARGT)
574  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
575        if(YEARGT.eq.1765)then  !     if(YEARGT.eq.1765)then
576         GHGBGR(1)=277.6  !      GHGBGR(1)=277.6
577         open (unit=861,file=init_co2,         open (unit=861,file=init_co2,
578       &     status='OLD',form='formatted')       &     status='OLD',form='formatted')
579         read (861,*)xco2init         read (861,*)xco2init
580         CO2=xco2init/GHGBGR(1)         CO2=xco2init/GHGBGR(1)
581         print *,'IPCC EMI CO2=',CO2         print *,'IPCC EMI CO2=',CO2
582        else  !     else
583           print *,' Wrong YEARGT ', YEARGT  !        print *,' Wrong YEARGT ', YEARGT
584           stop  !        stop
585        endif  !     endif
586  #endif  #endif
587        print *,'Background GHGs for year ',YEARGT        print *,'Background GHGs for year ',YEARGT
588        print '(5E12.4)',GHGBGR        print '(5E12.4)',GHGBGR
# Line 660  C       OLD Line 614  C       OLD
614        else        else
615          print *,'No diffusion into deep ocean'          print *,'No diffusion into deep ocean'
616        endif        endif
617  #if ( defined PREDICTED_AEROSOL )  !#if ( defined PREDICTED_AEROSOL )
618    #if ( defined CPL_CHEM )
619  Cold  AFBYCF=0.6725  Cold  AFBYCF=0.6725
620  Cold  SO2EREF=123.57  Cold  SO2EREF=123.57
621  Cigsm1AFBYCF=0.6054  Cigsm1AFBYCF=0.6054
# Line 674  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/ Line 629  Cold  CFAEROSOL=(-AERFOR/AFBYCF)**1.035/
629  !      CFAEROSOL=(-AERFOR/AFBYCF)**1.035/(SO2ER**1.0391)  !      CFAEROSOL=(-AERFOR/AFBYCF)**1.035/(SO2ER**1.0391)
630    
631    
632        SO2EREF=149.07        SO2EREF=147.375
633        open(664,file=SO2ERATIO,        open(664,file=SO2ERATIO,
634       &        form='formatted',       &        form='formatted',
635       &        status='old')       &        status='old')
636        read(664,'(f10.6)')SO2EM        read(664,'(f10.6)')SO2EM
637        SO2ER=SO2EM/SO2EREF        SO2ER=SO2EM/SO2EREF
638        F90BYF80=0.948  !     F90BYF80=0.948
639        AERFOR90=AERFOR*F90BYF80  !     AERFOR90=AERFOR*F90BYF80
640        AFBYCF90=1.043  !     AFBYCF90=1.04
641    !     print *,'AFBYCF90=',AFBYCF90
642    !     12/21/2006 (from runs 265x.06)
643    !     Adjustment for BC forcing
644           print *,' AERF4BC=',AERF4BC
645          AERFOR=AERFOR+AERF4BC
646    !     Adjustment for BC forcing
647          AFBYCF=1.05
648        print *,'SO2ER=',SO2ER        print *,'SO2ER=',SO2ER
649        print *,'AFBYCF90=',AFBYCF90        print *,'AFBYCF=',AFBYCF
650        CFAEROSOL=(-AERFOR90/AFBYCF90)**1.119        CFAEROSOL=(-AERFOR/AFBYCF)**1.21
651        CFAEROSOL=CFAEROSOL/(SO2ER**1.0391)        CFAEROSOL=CFAEROSOL/(SO2ER**1.01)
652        print *,'AERFOR90=',AERFOR90,'CFAEROSOL=',CFAEROSOL        print *,'AERFOR=',AERFOR,'CFAEROSOL=',CFAEROSOL
653  #if ( defined OCEAN_3D )  #if ( defined OCEAN_3D )
654         CFAEROSOL=CFAEROSOL/1.35         CFAEROSOL=CFAEROSOL/1.35
655         print *,'CFAEROSOL_3D=',CFAEROSOL         print *,'CFAEROSOL_3D=',CFAEROSOL
# Line 824  c      print *,' Form input ' Line 786  c      print *,' Form input '
786          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,
787       &  ODATA,       &  ODATA,
788       *  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.  
789       *  CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0       *  CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0
790         else         else
791  !       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,
792          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,
793       &  ODATA,       &  ODATA,
794       *   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.
795       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0,       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0,
796       *   TG3M,RTGO,STG3,DTG3       *   TG3M,RTGO,STG3,DTG3
797          endif          endif
798         if(.not.CONTRR)then         if(.not.CONTRR)then
# Line 879  c       print *,' WITH DEEP MIXED LAYER' Line 841  c       print *,' WITH DEEP MIXED LAYER'
841  #endif  #endif
842    
843  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
844           if(ocarcont) then
845            if(ocarindata)then
846              print *,'Wrong setting of ocarcont and ocarindata'
847              print *,ocarcont,ocarindata
848              stop
849            endif
850          open(116,file=deepco2in,          open(116,file=deepco2in,
851       * status='old',form='unformatted')       *   status='old',form='unformatted')
852          print *,' AFTER OPEN INIT. data for ocean chem.'          print *,' AFTER OPEN INIT. data for ocean chem.'
853          print *,deepco2in          print *,deepco2in
854          read(116)iyearocm,vdfocm          read(116)iyearocm,vdfocm
# Line 900  c       print *,' WITH DEEP MIXED LAYER' Line 868  c       print *,' WITH DEEP MIXED LAYER'
868           Rco2(j,k)=Rco2in(j,k)           Rco2(j,k)=Rco2in(j,k)
869          end do          end do
870          end do          end do
871           else
872            if(ocarindata)then
873             print *,' Reading initial data for ocean carbon'
874             open(116,file=deepco2in,
875         *    status='old',form='unformatted')
876             read(116)
877             read(116)Hgin
878             read(116)Rco2in
879             print *,(Rco2in(j,1),j=1,jm)
880             do k=1,lmo
881             do j=1,jm0
882    !         if(k.eq.1)Hg(j)=Hgin(j)
883              if(k.eq.1)Hg(j)=1.0e-8
884              Rco2(j,k)=Rco2in(j,k)
885             end do
886             end do
887             close(116)
888            else
889             do k=1,lmo
890             do j=1,jm0
891              if(k.eq.1)Hg(j)=1.0e-8
892              Rco2(j,k)=0.0
893             end do
894             end do
895            endif
896           endif
897  #endif  #endif
898    
899        else        else
# Line 909  C ******* Line 903  C *******
903          if(TRANSR)then          if(TRANSR)then
904           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.
905       *   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.
906       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0,       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0,
907       *   TG3M,RTGO,STG3,DTG3       *   TG3M,RTGO,STG3,DTG3
908          else          else
909           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.
910       *   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.
911       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TSURFT,TSURFD,DWAV0       *   CKN,WMGE,TPRIM2,MRCHT,TRSURF,SRSURF,TLANDD,TSURFD,DWAV0
912          endif          endif
913  C  C
914           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 937  c       if(iyearocm.ne.JYEAR) then
937          endif          endif
938          read(369)Hgin          read(369)Hgin
939          read(369)Rco2in          read(369)Rco2in
940            print *,(Rco2in(j,1),j=1,jm)
941          do k=1,lmo          do k=1,lmo
942          do j=1,jm0          do j=1,jm0
943           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 947  c       if(iyearocm.ne.JYEAR) then
947  #endif  #endif
948    
949        endif        ! endif for ISTRT1        endif        ! endif for ISTRT1
950    #if ( defined CPL_OCEANCO2 && defined ML_2D )
951           print *,'Rco2 from input'
952    !        print *,(Rco2(j,1),j=1,jm)
953             print *,Rco2
954    #endif
955    
956        if(KOCEAN.eq.1) then        if(KOCEAN.eq.1) then
957         print *,' T1 ocean'         print *,' T1 ocean'
# Line 980  c     endif Line 980  c     endif
980        TAUP=TAUX                                                         1746.          TAUP=TAUX                                                         1746.  
981  C**** UPDATE C ARRAY FROM INPUTZ                                        1747.    C**** UPDATE C ARRAY FROM INPUTZ                                        1747.  
982    500 READ (514,INPUTZ)                                                 1748.      500 READ (514,INPUTZ)                                                 1748.  
983    !      nrad=NDYN
984  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
985         CO2=xco2init/GHGBGR(1)         CO2=xco2init/GHGBGR(1)
986  #endif  #endif
# Line 1020  c     print *,' NWRGHG=',NWRGHG Line 1021  c     print *,' NWRGHG=',NWRGHG
1021         print *,TSCNTR,' form EXP=',EXPTSF         print *,TSCNTR,' form EXP=',EXPTSF
1022  !      read(576)EXPTSF  !      read(576)EXPTSF
1023  !      print *,' TSURF form EXP=',EXPTSF  !      print *,' TSURF form EXP=',EXPTSF
1024         read(576)TSURFC         read(576)TSURFC,TLANDC
1025    !      read(576)TLANDC
1026        endif        endif
1027          CORSR=1.0          CORSR=1.0
1028        if(KOCEAN.eq.0)then        if(KOCEAN.eq.0)then
# Line 1063  c     PRINT *,' NLFR=',NLFR Line 1065  c     PRINT *,' NLFR=',NLFR
1065        if(READGHG.eq.2)then        if(READGHG.eq.2)then
1066          open( unit=569,file=ghg_monthly,          open( unit=569,file=ghg_monthly,
1067       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1068          open( unit=669,file=ghg_monthly2,          open( unit=679,file=ghg_monthly2,
1069       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1070          print *,' GHGs from ',ghg_monthly          print *,' GHGs from ',ghg_monthly
1071          if(ISTRT1.eq.1)then          if(ISTRT1.eq.1)then
1072            do 369 i=1,NWRGHG            do 369 i=1,NWRGHG
1073             do 369 ii=1,13             do 369 ii=1,13
1074             read(569)             read(569)
1075             if(ii.le.3)read(669)             if(ii.le.3)read(679)
1076    369     continue    369     continue
1077        endif        endif
1078        endif        endif
# Line 1079  c       open( unit=569,file=dirdat1(1:id Line 1081  c       open( unit=569,file=dirdat1(1:id
1081  c    *  status='OLD',form='unformatted')  c    *  status='OLD',form='unformatted')
1082          open( unit=569,file=ghg_monthly,          open( unit=569,file=ghg_monthly,
1083       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1084          open( unit=669,file=ghg_monthly2,          open( unit=679,file=ghg_monthly2,
1085       *  status='OLD',form='unformatted')       *  status='OLD',form='unformatted')
1086          print *,' GHGs from ',ghg_monthly          print *,' GHGs from ',ghg_monthly
1087        endif        endif
# Line 1287  c          IF(FDATA(I,J,2).LT.0.01)FDATA Line 1289  c          IF(FDATA(I,J,2).LT.0.01)FDATA
1289          fland_temp(j)=FDATA(1,J,2)          fland_temp(j)=FDATA(1,J,2)
1290    284 continue    284 continue
1291  #if( defined OCEAN_3D)  #if( defined OCEAN_3D)
1292  Cjrs       if(jmocean.ne.jm0-2)then  Cjrs        if(jmocean.ne.jm0-2)then
1293  C          print *,"Wrong jm or jmocean"  C          print *,"Wrong jm or jmocean"
1294  C          stop  C          stop
1295  C       endif  C       endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22