/[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.4 by jscott, Thu Oct 4 20:40:27 2007 UTC revision 1.7 by jscott, Tue Sep 1 22:03:56 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "ctrparam.h"  #include "ctrparam.h"
5    
# 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    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,emiss_data,SO2ERATIO,SEN_dat
93       &  ,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       & chem_init,chem_init2,chemrstfl       &  ,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 130  c==== 012201 Line 135  c==== 012201
135       &  S0C_data,cfvolaer,       &  S0C_data,cfvolaer,
136       &  dirdat1,dirdat2       &  dirdat1,dirdat2
137       &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss       &  ,chemdata,chemout,last_nep,init_4nem,pov_deepo,fnememiss
138       &  ,chem_init,chem_init2,chemrstfl       &  ,chem_init,chem_init2,chem_init4pfc,chemrstfl
139       &    ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av       &    ,oco2file,co2rfile,caruptfile,emiss_data,SO2ERATIO,flrco2av
140       &  ,flin_nep,last_clm,SEN_dat,nemdatdir       &  ,flin_nep,last_clm,SEN_dat,nemdatdir
141       &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer       &  ,GHSFALB,GHSF,ALBCF,FVOLADD,ghostfile,fl_volaer
# Line 139  c==== 012201 Line 144  c==== 012201
144       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3       &  ,fclmlice,fbaresoil,fwmax,fprratio,o3datadir,CLIMO3
145       &  ,OCNGEOM,ocngmfile,ocndata4atm       &  ,OCNGEOM,ocngmfile,ocndata4atm
146       &  ,bc_data       &  ,bc_data
147         &  ,PRTREND
148  Cjrs     &  ,dtatm,dtocn  Cjrs     &  ,dtatm,dtocn
149  #if ( defined IPCC_EMI )  #if ( defined IPCC_EMI )
150       &  ,init_co2       &  ,init_co2
# Line 273  C Line 279  C
279        GSOEQ=.false.        GSOEQ=.false.
280        VEGCH=.false.        VEGCH=.false.
281        TRVEG=.false.        TRVEG=.false.
282          PRTREND=.false.
283        OCNGEOM=.false.        OCNGEOM=.false.
284        do L=1,LM+1        do L=1,LM+1
285          ghostfv(L)=0.          ghostfv(L)=0.
# Line 344  cjrs      dtocn=1 Line 351  cjrs      dtocn=1
351        chemdata = 'DATA'        chemdata = 'DATA'
352        chem_init = 'init-data_46x11_1991'        chem_init = 'init-data_46x11_1991'
353        chem_init2 = 'init-data2_46x11_1991'        chem_init2 = 'init-data2_46x11_1991'
354          chem_init4pfc = 'undefined'
355        chemrstfl = 'undefined'        chemrstfl = 'undefined'
356        nemdatdir = 'TEMDATA'        nemdatdir = 'TEMDATA'
357        emiss_data='edaily.dat'        emiss_data='edaily.dat'
# Line 556  c Line 564  c
564  c     open file for carbon uptake  c     open file for carbon uptake
565  #if ( defined CPL_TEM || defined CPL_OCEANCO2 )  #if ( defined CPL_TEM || defined CPL_OCEANCO2 )
566  c     open(333,file=caruptfile,status='new',form='formatted')  c     open(333,file=caruptfile,status='new',form='formatted')
567        open(333,file=caruptfile,form='formatted')        open(333,file=caruptfile,status='replace',form='formatted')
568          close(333)
569  #endif  #endif
570    
571  #if ( defined CPL_OCEANCO2 && defined ML_2D)  #if ( defined CPL_OCEANCO2 && defined ML_2D)
# Line 614  C       OLD Line 623  C       OLD
623        else        else
624          print *,'No diffusion into deep ocean'          print *,'No diffusion into deep ocean'
625        endif        endif
626          if (PRTREND) then
627             print *,' With changes in  stochastic precip'
628          else
629             print *,' Without changes in  stochastic precip'
630          endif
631  !#if ( defined PREDICTED_AEROSOL )  !#if ( defined PREDICTED_AEROSOL )
632  #if ( defined CPL_CHEM )  #if ( defined CPL_CHEM )
633  Cold  AFBYCF=0.6725  Cold  AFBYCF=0.6725
# Line 720  c     LMCM=LTM+2 Line 734  c     LMCM=LTM+2
734       *    status='OLD',form='unformatted')       *    status='OLD',form='unformatted')
735  #endif  #endif
736  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
737            open (367,file=last_nep,form='unformatted',status='new')            open (876,file=last_nep,form='unformatted',status='new')
738            open (877,file=last_clm,form='unformatted',status='new')            open (877,file=last_clm,form='unformatted',status='new')
739  c         file last_clm contains data for posible restart of NEM  c         file last_clm contains data for posible restart of NEM
740  c         this file is writen at the end of the run  c         this file is writen at the end of the run
# Line 728  c         this file is writen at the end Line 742  c         this file is writen at the end
742            open (368,file=init_4nem,form='unformatted',status='old')            open (368,file=init_4nem,form='unformatted',status='old')
743  c         file init_4nem contains data for the restart of NEM  c         file init_4nem contains data for the restart of NEM
744  c         from the results of a previous run  c         from the results of a previous run
745            open (277,file=fnememiss,form='unformatted',status='new')            open (277,file=fnememiss,form='unformatted',status='replace')
746              close(277)
747  #endif  #endif
748  #endif  #endif
749  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
# Line 741  C      For restart of the run Line 756  C      For restart of the run
756            open( unit=547,file=nwrfl,            open( unit=547,file=nwrfl,
757       *    status='OLD',form='unformatted')       *    status='OLD',form='unformatted')
758  #if ( defined CPL_TEM )  #if ( defined CPL_TEM )
759            open (367,file=last_nep,form='unformatted',status='old')            open (876,file=last_nep,form='unformatted',status='old')
760            open (877,file=last_clm,form='unformatted',status='old')            open (877,file=last_clm,form='unformatted',status='old')
761  #if ( defined CPL_NEM )  #if ( defined CPL_NEM )
762            open (368,file=init_4nem,form='unformatted',status='old')            open (368,file=init_4nem,form='unformatted',status='old')
763  c         file init_4nem contains data for the restart of NEM  c         file init_4nem contains data for the restart of NEM
764  c         from the results of a previous run  c         from the results of a previous run
765            open (277,file=fnememiss,form='unformatted',status='new')            open (277,file=fnememiss,form='unformatted',status='replace')
766              close(277)
767  #endif  #endif
768  #endif  #endif
769  #if ( defined CPL_OCEANCO2 && defined ML_2D )  #if ( defined CPL_OCEANCO2 && defined ML_2D )
# Line 1021  c     print *,' NWRGHG=',NWRGHG Line 1037  c     print *,' NWRGHG=',NWRGHG
1037         print *,TSCNTR,' form EXP=',EXPTSF         print *,TSCNTR,' form EXP=',EXPTSF
1038  !      read(576)EXPTSF  !      read(576)EXPTSF
1039  !      print *,' TSURF form EXP=',EXPTSF  !      print *,' TSURF form EXP=',EXPTSF
1040         read(576)TSURFC,TLANDC  #if ( defined CLM )
1041            read(576)TSURFC,TLANDC
1042    #else
1043            read(576)TSURFC
1044    #endif
1045  !      read(576)TLANDC  !      read(576)TLANDC
1046        endif        endif
1047          CORSR=1.0          CORSR=1.0
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22