/[MITgcm]/MITgcm/pkg/ctrl/ctrl_init.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.7 by heimbach, Fri Mar 7 02:45:48 2003 UTC revision 1.8 by heimbach, Mon Jun 23 22:29:05 2003 UTC
# Line 98  c     == external == Line 98  c     == external ==
98    
99  c     == end of interface ==  c     == end of interface ==
100    
 c--   Read the namelist input.  
       namelist /ctrl_nml/  
      &  xx_theta_file,  
      &  xx_salt_file,  
      &  xx_hflux_file,  
      &  xx_hfluxstartdate1, xx_hfluxstartdate2, xx_hfluxperiod,  
      &  xx_sflux_file,  
      &  xx_sfluxstartdate1, xx_sfluxstartdate2, xx_sfluxperiod,  
      &  xx_tauu_file,  
      &  xx_tauustartdate1,  xx_tauustartdate2,  xx_tauuperiod,  
      &  xx_tauv_file,  
      &  xx_tauvstartdate1,  xx_tauvstartdate2,  xx_tauvperiod,  
      &  xx_atemp_file,  
      &  xx_atempstartdate1, xx_atempstartdate2, xx_atempperiod,  
      &  xx_aqh_file,  
      &  xx_aqhstartdate1, xx_aqhstartdate2, xx_aqhperiod,  
      &  xx_uwind_file,  
      &  xx_uwindstartdate1, xx_uwindstartdate2, xx_uwindperiod,  
      &  xx_vwind_file,  
      &  xx_vwindstartdate1, xx_vwindstartdate2, xx_vwindperiod,  
      &  xx_obcsn_file,  
      &  xx_obcsnstartdate1,  xx_obcsnstartdate2,  xx_obcsnperiod,  
      &  xx_obcss_file,  
      &  xx_obcssstartdate1,  xx_obcssstartdate2,  xx_obcssperiod,  
      &  xx_obcsw_file,  
      &  xx_obcswstartdate1,  xx_obcswstartdate2,  xx_obcswperiod,  
      &  xx_obcse_file,  
      &  xx_obcsestartdate1,  xx_obcsestartdate2,  xx_obcseperiod,  
      &  xx_diffkr_file,  
      &  xx_kapgm_file,  
      &  xx_tr1_file,  
      &  xx_sst_file,  
      &  xx_sss_file,  
      &  xx_hfacc_file,  
      &  xx_efluxy_file,  
      &  xx_efluxp_file,  
      &  xx_bottomdrag_file  
   
       namelist /ctrl_packnames/  
      &  yadmark, yctrlid,  
      &  ctrlname, costname, scalname, maskname, metaname  
   
101        jtlo = mybylo(mythid)        jtlo = mybylo(mythid)
102        jthi = mybyhi(mythid)        jthi = mybyhi(mythid)
103        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
# Line 152  c--   Read the namelist input. Line 110  c--   Read the namelist input.
110        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
111    
112  c--     Set default values.  c--     Set default values.
         xx_theta_file      = ' '  
         xx_salt_file       = ' '  
         xx_hfluxstartdate1 =   0  
         xx_hfluxstartdate2 =   0  
         xx_hfluxperiod     =   0. _d 0  
         xx_hflux_file      = ' '  
         xx_sfluxstartdate1 =   0  
         xx_sfluxstartdate2 =   0  
         xx_sfluxperiod     =   0. _d 0  
         xx_sflux_file      = ' '  
         xx_tauustartdate1  =   0  
         xx_tauustartdate2  =   0  
         xx_tauuperiod      =   0. _d 0  
         xx_tauu_file       = ' '  
         xx_tauvstartdate1  =   0  
         xx_tauvstartdate2  =   0  
         xx_tauvperiod      =   0. _d 0  
         xx_tauv_file       = ' '  
         xx_atempstartdate1  =   0  
         xx_atempstartdate2  =   0  
         xx_atempperiod      =   0. _d 0  
         xx_atemp_file       = ' '  
         xx_aqhstartdate1    =   0  
         xx_aqhstartdate2    =   0  
         xx_aqhperiod        =   0. _d 0  
         xx_aqh_file         = ' '  
         xx_uwindstartdate1  =   0  
         xx_uwindstartdate2  =   0  
         xx_uwindperiod      =   0. _d 0  
         xx_uwind_file       = ' '  
         xx_vwindstartdate1  =   0  
         xx_vwindstartdate2  =   0  
         xx_vwindperiod      =   0. _d 0  
         xx_vwind_file       = ' '  
         xx_obcsnstartdate1  =   0  
         xx_obcsnstartdate2  =   0  
         xx_obcsnperiod      =   0. _d 0  
         xx_obcsn_file       = ' '  
         xx_obcssstartdate1  =   0  
         xx_obcssstartdate2  =   0  
         xx_obcssperiod      =   0. _d 0  
         xx_obcss_file       = ' '  
         xx_obcswstartdate1  =   0  
         xx_obcswstartdate2  =   0  
         xx_obcswperiod      =   0. _d 0  
         xx_obcsw_file       = ' '  
         xx_obcsestartdate1  =   0  
         xx_obcsestartdate2  =   0  
         xx_obcseperiod      =   0. _d 0  
         xx_obcse_file       = ' '  
         xx_diffkr_file      = ' '  
         xx_kapgm_file       = ' '  
         xx_tr1_file         = ' '  
         xx_sst_file         = ' '  
         xx_sss_file         = ' '  
         xx_hfacc_file       = ' '  
         xx_efluxy_file      = ' '  
         xx_efluxp_file      = ' '  
         xx_bottomdrag_file  = ' '  
   
         yadmark             = 'ad'  
         yctrlid             = 'MIT_CE_000'  
         ctrlname            = ' '  
         costname            = ' '  
         scalname            = ' '  
         maskname            = ' '  
         metaname            = ' '  
   
 c--     Check versions.  
   
         open(unit=scrunit1,status='scratch')  
   
 c--     Next, read the ecco data file.  
         open(unit = modeldataunit,file = 'data.ctrl',  
      &       status = 'old', iostat = errio)  
         if ( errio .lt. 0 ) then  
           stop ' stopped in ctrl_init'  
         endif  
   
         do while ( .true. )  
           read(modeldataunit, fmt='(a)', end=1001) record  
           il = max(ilnblnk(record),1)  
           if ( record(1:1) .ne. commentcharacter )  
      &        write(unit=scrunit1, fmt='(a)') record(:il)  
         enddo  
  1001   continue  
         close( modeldataunit )  
   
         rewind( scrunit1 )  
         read(unit = scrunit1, nml = ctrl_nml)  
         read(unit = scrunit1, nml = ctrl_packnames)  
         close( scrunit1 )  
   
 #ifdef ALLOW_CALENDAR  
   
 c--     Get the complete dates of the control variables.  
 #if     (defined  (ALLOW_HFLUX_CONTROL))  
 c--     The heat flux contribution.  
         call cal_FullDate( xx_hfluxstartdate1, xx_hfluxstartdate2,  
      &                     xx_hfluxstartdate , mythid )  
 #elif   (defined  (ALLOW_ATEMP_CONTROL))  
 c--     Atmos. temperature contribution.  
         call cal_FullDate( xx_atempstartdate1, xx_atempstartdate2,  
      &                     xx_atempstartdate , mythid )  
 #endif  
   
 #if     (defined  (ALLOW_SFLUX_CONTROL))  
 c--     The salt flux contribution.  
         call cal_FullDate( xx_sfluxstartdate1, xx_sfluxstartdate2,  
      &                     xx_sfluxstartdate , mythid )  
 #elif   (defined  (ALLOW_AQH_CONTROL))  
 c--     Atmospheric humidity contribution.  
         call cal_FullDate( xx_aqhstartdate1, xx_aqhstartdate2,  
      &                     xx_aqhstartdate , mythid )  
 #endif  
   
 #if     (defined  (ALLOW_USTRESS_CONTROL))  
 c--     The zonal wind stress contribution.  
         call cal_FullDate( xx_tauustartdate1,  xx_tauustartdate2,  
      &                     xx_tauustartdate,   mythid )  
 #elif   (defined  (ALLOW_UWIND_CONTROL))  
 c--     Zonal wind speed contribution.  
         call cal_FullDate( xx_uwindstartdate1, xx_uwindstartdate2,  
      &                     xx_uwindstartdate , mythid )  
 #endif  
   
 #if     (defined  (ALLOW_VSTRESS_CONTROL))  
 c--     The merid. wind stress contribution.  
         call cal_FullDate( xx_tauvstartdate1,  xx_tauvstartdate2,  
      &                     xx_tauvstartdate,   mythid )  
 #elif   (defined  (ALLOW_VWIND_CONTROL))  
 c--     Merid. wind speed contribution.  
         call cal_FullDate( xx_vwindstartdate1, xx_vwindstartdate2,  
      &                     xx_vwindstartdate , mythid )  
 #endif  
   
 #ifdef ALLOW_OBCS_CONTROL  
         call cal_FullDate( xx_obcsnstartdate1,  xx_obcsnstartdate2,  
      &                     xx_obcsnstartdate,   mythid )  
         call cal_FullDate( xx_obcssstartdate1,  xx_obcssstartdate2,  
      &                     xx_obcssstartdate,   mythid )  
         call cal_FullDate( xx_obcswstartdate1,  xx_obcswstartdate2,  
      &                     xx_obcswstartdate,   mythid )  
         call cal_FullDate( xx_obcsestartdate1,  xx_obcsestartdate2,  
      &                     xx_obcsestartdate,   mythid )  
 #endif  
   
 #endif /* ALLOW_CALENDAR */  
   
 c--     Set default values.  
113          do ivarindex = 1,maxcvars          do ivarindex = 1,maxcvars
114            ncvarindex(ivarindex) = -1            ncvarindex(ivarindex) = -1
115            ncvarrecs(ivarindex)  =  0            ncvarrecs(ivarindex)  =  0
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22