/[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.17 by heimbach, Thu Jul 28 13:47:49 2005 UTC revision 1.18 by heimbach, Thu Jul 28 19:51:22 2005 UTC
# Line 125  cph( Line 125  cph(
125  cph    index  7-10 reserved for atmos. state,  cph    index  7-10 reserved for atmos. state,
126  cph    index 11-14 reserved for open boundaries,  cph    index 11-14 reserved for open boundaries,
127  cph    index 15-16 reserved for mixing coeff.  cph    index 15-16 reserved for mixing coeff.
128  cph    index 17    reserved for precip. (atmos. state)  cph    index 17    reserved for passive tracer TR1
129  cph    index 18,19 reserved for sst, sss  cph    index 18,19 reserved for sst, sss
130  cph    index 20             for hFacC  cph    index 20             for hFacC
131  cph    index 21-22          for efluxy, efluxp  cph    index 21-22          for efluxy, efluxp
# Line 134  cph    index 24 Line 134  cph    index 24
134  cph    index 25-26          for edtaux, edtauy  cph    index 25-26          for edtaux, edtauy
135  cph    index 27-29          for uvel0, vvel0, etan0  cph    index 27-29          for uvel0, vvel0, etan0
136  cph    index 30-31          for relax. SST, SSS  cph    index 30-31          for relax. SST, SSS
137  cph    index 32    reserved for passive tracer TR1  cph    index 32    reserved for precip (atmos. state)
138    cph    index 33    reserved for swflux (atmos. state)
139  cph)  cph)
140    
141  c----------------------------------------------------------------------  c----------------------------------------------------------------------
# Line 503  c-- Line 503  c--
503    
504  c----------------------------------------------------------------------  c----------------------------------------------------------------------
505  c--  c--
506  #ifdef ALLOW_PRECIP_CONTROL  #ifdef ALLOW_TR10_CONTROL
 c--   Atmos. precipitation  
   
 # ifdef ALLOW_CAL  
         call cal_FullDate( xx_precipstartdate1, xx_precipstartdate2,  
      &                     xx_precipstartdate , mythid )  
         call cal_TimePassed( xx_precipstartdate, modelstartdate,  
      &                       difftime, mythid )  
         call cal_ToSeconds ( difftime, diffsecs, mythid )  
         startrec = int((modelstart + startTime - diffsecs)/  
      &                 xx_precipperiod) + 1  
         endrec   = int((modelend + startTime - diffsecs + modelstep/2)/  
      &                 xx_precipperiod) + 2  
 # else  
         startrec = 1  
         endrec   = 1  
 # endif  
         diffrec  = endrec - startrec + 1  
507          call ctrl_init_ctrlvar (          call ctrl_init_ctrlvar (
508       &       xx_precip_file, 17, 117, diffrec, startrec, endrec,       &       xx_tr1_file, 17, 117, 1, 1, 1,
509       &       snx, sny, 1, 'c', 'xy', mythid )       &       snx, sny, nr, 'c', '3d', mythid )
510    #endif /* ALLOW_TR10_CONTROL */
 #endif /* ALLOW_PRECIP_CONTROL */  
511    
512  c----------------------------------------------------------------------  c----------------------------------------------------------------------
513  c--  c--
# Line 639  c-- Line 621  c--
621    
622  c----------------------------------------------------------------------  c----------------------------------------------------------------------
623  c--  c--
624  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_PRECIP_CONTROL
625    c--   Atmos. precipitation
626    
627    # ifdef ALLOW_CAL
628            call cal_FullDate( xx_precipstartdate1, xx_precipstartdate2,
629         &                     xx_precipstartdate , mythid )
630            call cal_TimePassed( xx_precipstartdate, modelstartdate,
631         &                       difftime, mythid )
632            call cal_ToSeconds ( difftime, diffsecs, mythid )
633            startrec = int((modelstart + startTime - diffsecs)/
634         &                 xx_precipperiod) + 1
635            endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
636         &                 xx_precipperiod) + 2
637    # else
638            startrec = 1
639            endrec   = 1
640    # endif
641            diffrec  = endrec - startrec + 1
642          call ctrl_init_ctrlvar (          call ctrl_init_ctrlvar (
643       &       xx_tr1_file, 32, 132, 1, 1, 1,       &       xx_precip_file, 32, 132, diffrec, startrec, endrec,
644       &       snx, sny, nr, 'c', '3d', mythid )       &       snx, sny, 1, 'c', 'xy', mythid )
645  #endif /* ALLOW_TR10_CONTROL */  
646    #endif /* ALLOW_PRECIP_CONTROL */
647    
648    c----------------------------------------------------------------------
649    c--
650    #ifdef ALLOW_SWFLUX_CONTROL
651    c--   Atmos. swflux
652    
653    # ifdef ALLOW_CAL
654            call cal_FullDate( xx_swfluxstartdate1, xx_swfluxstartdate2,
655         &                     xx_swfluxstartdate , mythid )
656            call cal_TimePassed( xx_swfluxstartdate, modelstartdate,
657         &                       difftime, mythid )
658            call cal_ToSeconds ( difftime, diffsecs, mythid )
659            startrec = int((modelstart + startTime - diffsecs)/
660         &                 xx_swfluxperiod) + 1
661            endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
662         &                 xx_swfluxperiod) + 2
663    # else
664            startrec = 1
665            endrec   = 1
666    # endif
667            diffrec  = endrec - startrec + 1
668            call ctrl_init_ctrlvar (
669         &       xx_swflux_file, 33, 133, diffrec, startrec, endrec,
670         &       snx, sny, 1, 'c', 'xy', mythid )
671    
672    #endif /* ALLOW_SWFLUX_CONTROL */
673    
674  c----------------------------------------------------------------------  c----------------------------------------------------------------------
675  c----------------------------------------------------------------------  c----------------------------------------------------------------------
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22