/[MITgcm]/MITgcm/pkg/obcs/obcs_prescribe_read.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_prescribe_read.F

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

revision 1.20 by heimbach, Fri Mar 28 03:49:54 2008 UTC revision 1.26 by heimbach, Sun Dec 19 20:29:29 2010 UTC
# Line 54  c     == end of interface == Line 54  c     == end of interface ==
54       U     OBNv,   OBNv0,   OBNv1,   OBNvfile,       U     OBNv,   OBNv0,   OBNv1,   OBNvfile,
55       U     OBNt,   OBNt0,   OBNt1,   OBNtfile,       U     OBNt,   OBNt0,   OBNt1,   OBNtfile,
56       U     OBNs,   OBNs0,   OBNs1,   OBNsfile,       U     OBNs,   OBNs0,   OBNs1,   OBNsfile,
57    #   ifdef NONLIN_FRSURF
58         U     OBNeta, OBNeta0, OBNeta1, OBNetafile,
59    #   endif
60  #   ifdef ALLOW_SEAICE  #   ifdef ALLOW_SEAICE
61       I     siobNstartdate, siobNperiod,       I     siobNstartdate, siobNperiod,
62       U     OBNa,   OBNa0,   OBNa1,   OBNafile,       U     OBNa,   OBNa0,   OBNa1,   OBNafile,
# Line 78  c     == end of interface == Line 81  c     == end of interface ==
81       U     OBSv,   OBSv0,   OBSv1,   OBSvfile,       U     OBSv,   OBSv0,   OBSv1,   OBSvfile,
82       U     OBSt,   OBSt0,   OBSt1,   OBStfile,       U     OBSt,   OBSt0,   OBSt1,   OBStfile,
83       U     OBSs,   OBSs0,   OBSs1,   OBSsfile,       U     OBSs,   OBSs0,   OBSs1,   OBSsfile,
84    #   ifdef NONLIN_FRSURF
85         U     OBSeta, OBSeta0, OBSeta1, OBSetafile,
86    #   endif
87  #   ifdef ALLOW_SEAICE  #   ifdef ALLOW_SEAICE
88       I     siobSstartdate, siobSperiod,       I     siobSstartdate, siobSperiod,
89       U     OBSa,   OBSa0,   OBSa1,   OBSafile,       U     OBSa,   OBSa0,   OBSa1,   OBSafile,
# Line 102  c     == end of interface == Line 108  c     == end of interface ==
108       U     OBEv,   OBEv0,   OBEv1,   OBEvfile,       U     OBEv,   OBEv0,   OBEv1,   OBEvfile,
109       U     OBEt,   OBEt0,   OBEt1,   OBEtfile,       U     OBEt,   OBEt0,   OBEt1,   OBEtfile,
110       U     OBEs,   OBEs0,   OBEs1,   OBEsfile,       U     OBEs,   OBEs0,   OBEs1,   OBEsfile,
111    #   ifdef NONLIN_FRSURF
112         U     OBEeta, OBEeta0, OBEeta1, OBEetafile,
113    #   endif
114  #   ifdef ALLOW_SEAICE  #   ifdef ALLOW_SEAICE
115       I     siobEstartdate, siobEperiod,       I     siobEstartdate, siobEperiod,
116       U     OBEa,   OBEa0,   OBEa1,   OBEafile,       U     OBEa,   OBEa0,   OBEa1,   OBEafile,
# Line 126  c     == end of interface == Line 135  c     == end of interface ==
135       U     OBWv,   OBWv0,   OBWv1,   OBWvfile,       U     OBWv,   OBWv0,   OBWv1,   OBWvfile,
136       U     OBWt,   OBWt0,   OBWt1,   OBWtfile,       U     OBWt,   OBWt0,   OBWt1,   OBWtfile,
137       U     OBWs,   OBWs0,   OBWs1,   OBWsfile,       U     OBWs,   OBWs0,   OBWs1,   OBWsfile,
138    #   ifdef NONLIN_FRSURF
139         U     OBWeta, OBWeta0, OBWeta1, OBWetafile,
140    #   endif
141  #   ifdef ALLOW_SEAICE  #   ifdef ALLOW_SEAICE
142       I     siobEstartdate, siobEperiod,       I     siobWstartdate, siobWperiod,
143       U     OBWa,   OBWa0,   OBWa1,   OBWafile,       U     OBWa,   OBWa0,   OBWa1,   OBWafile,
144       U     OBWh,   OBWh0,   OBWh1,   OBWhfile,       U     OBWh,   OBWh0,   OBWh1,   OBWhfile,
145       U     OBWsl,  OBWsl0,  OBWsl1,  OBWslfile,       U     OBWsl,  OBWsl0,  OBWsl1,  OBWslfile,
# Line 167  cgg    calendar information as other bou Line 179  cgg    calendar information as other bou
179        call ctrl_getobcse ( mycurrenttime, mycurrentiter, mythid )        call ctrl_getobcse ( mycurrenttime, mycurrentiter, mythid )
180  # endif  # endif
181    
 #ifndef ALLOW_EXF  
182        IF ( .NOT. useEXF ) THEN        IF ( .NOT. useEXF ) THEN
183    cph#ifndef ALLOW_AUTODIFF_TAMC
184         CALL OBCS_EXTERNAL_FIELDS_LOAD(         CALL OBCS_EXTERNAL_FIELDS_LOAD(
185       &     myCurrentTime, myCurrentIter, myThid )       &     myCurrentTime, myCurrentIter, myThid )
186    cph#else
187    cph       STOP 'PH HAS DISABLED THIS RUNTIME OPTION FOR ALLOW_EXF'
188    cph#endif
189        ENDIF        ENDIF
 #else  
       STOP 'PH HAS DISABLED THIS RUNTIME OPTION FOR ALLOW_EXF'  
 #endif  
190    
191  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE */  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE */
192    
# Line 192  C======================================= Line 204  C=======================================
204       U     OBv,   OBv0,   OBv1,   OBvfile,       U     OBv,   OBv0,   OBv1,   OBvfile,
205       U     OBt,   OBt0,   OBt1,   OBtfile,       U     OBt,   OBt0,   OBt1,   OBtfile,
206       U     OBs,   OBs0,   OBs1,   OBsfile,       U     OBs,   OBs0,   OBs1,   OBsfile,
207    #ifdef NONLIN_FRSURF
208         U     OBeta, OBeta0, OBeta1, OBetafile,
209    #endif
210  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
211       I     siobstartdate, siobperiod,       I     siobstartdate, siobperiod,
212       U     OBa,   OBa0,   OBa1,   OBafile,       U     OBa,   OBa0,   OBa1,   OBafile,
# Line 247  c     == routine arguments == Line 262  c     == routine arguments ==
262        _RL OBt1    (1-Olx:sNx+Olx,Nr,nSx,nSy)        _RL OBt1    (1-Olx:sNx+Olx,Nr,nSx,nSy)
263        _RL OBs1    (1-Olx:sNx+Olx,Nr,nSx,nSy)        _RL OBs1    (1-Olx:sNx+Olx,Nr,nSx,nSy)
264        CHARACTER*(MAX_LEN_FNAM) OBuFile,OBvFile,OBtFile,OBsFile        CHARACTER*(MAX_LEN_FNAM) OBuFile,OBvFile,OBtFile,OBsFile
265    #ifdef NONLIN_FRSURF
266          _RL OBeta   (1-Olx:sNx+Olx,nSx,nSy)
267          _RL OBeta0  (1-Olx:sNx+Olx,nSx,nSy)
268          _RL OBeta1  (1-Olx:sNx+Olx,nSx,nSy)
269          CHARACTER*(MAX_LEN_FNAM) OBetaFile
270    #endif
271  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
272        _RL     siobstartdate        _RL     siobstartdate
273        _RL     siobperiod        _RL     siobperiod
# Line 292  c     == local variables == Line 313  c     == local variables ==
313  # ifdef ALLOW_PTRACERS  # ifdef ALLOW_PTRACERS
314        integer iTracer        integer iTracer
315  # endif /* ALLOW_PTRACERS */  # endif /* ALLOW_PTRACERS */
   
316  c     == end of interface ==  c     == end of interface ==
317        if ( obcsperiod .eq. -12 ) then  
318          if ( obcsperiod .eq. -12. _d 0 ) then
319  c     obcsperiod=-12 means input file contains 12 monthly means  c     obcsperiod=-12 means input file contains 12 monthly means
320  c     record numbers are assumed 1 to 12 corresponding to  c     record numbers are assumed 1 to 12 corresponding to
321  c     Jan. through Dec.  c     Jan. through Dec.
# Line 304  c     Jan. through Dec. Line 325  c     Jan. through Dec.
325       I                        mycurrenttime, mycurrentiter, mythid       I                        mycurrenttime, mycurrentiter, mythid
326       &           )       &           )
327    
328        elseif ( obcsperiod .lt. 0 ) then        elseif ( obcsperiod .lt. 0. _d 0 ) then
329         print *, 'obcsperiod is out of range'         print *, 'obcsperiod is out of range'
330         STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_XZ'         STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_XZ'
331        else        else
# Line 316  c     get record numbers and interpolati Line 337  c     get record numbers and interpolati
337       O                       count0, count1, year0, year1,       O                       count0, count1, year0, year1,
338       I                       mycurrenttime, mycurrentiter, mythid       I                       mycurrenttime, mycurrentiter, mythid
339       &                      )       &                      )
        call exf_GetFFieldRec(  
      I                       siobstartdate, siobperiod,  
      I                       useYearlyFields,  
      O                       fac, first, changed,  
      O                       count0, count1, year0, year1,  
      I                       mycurrenttime, mycurrentiter, mythid  
      &                      )  
340        endif        endif
   
341        call exf_set_obcs_xz(  OBu, OBu0, OBu1, OBufile, 'u'        call exf_set_obcs_xz(  OBu, OBu0, OBu1, OBufile, 'u'
342       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
343       I                     , obcsperiod, count0, count1, year0, year1       I                     , obcsperiod, count0, count1, year0, year1
# Line 341  c     get record numbers and interpolati Line 354  c     get record numbers and interpolati
354       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
355       I                     , obcsperiod, count0, count1, year0, year1       I                     , obcsperiod, count0, count1, year0, year1
356       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
357    # ifdef NONLIN_FRSURF
358           call exf_set_obcs_x ( OBeta, OBeta0, OBeta1, OBetaFile, 's'
359         I                     , fac, first, changed, useYearlyFields
360         I                     , obcsperiod, count0, count1, year0, year1
361         I                     , mycurrenttime, mycurrentiter, mythid )
362    # endif /* NONLIN_FRSURF */
363    # ifdef ALLOW_PTRACERS
364          if ( usePTRACERS ) then
365           do iTracer = 1, PTRACERS_numInUse
366            call exf_set_obcs_xz(  OBptr (1-Olx,1,1,1,iTracer)
367         I                       , OBptr0(1-Olx,1,1,1,iTracer)
368         I                       , OBptr1(1-Olx,1,1,1,iTracer)
369         I                       , OBptrFile(iTracer), 's'
370         I                       , fac, first, changed, useYearlyFields
371         I                       , obcsperiod, count0, count1, year0, year1
372         I                       , mycurrenttime, mycurrentiter, mythid )
373           enddo
374          endif
375    # endif /* ALLOW_PTRACERS */
376  # ifdef ALLOW_SEAICE  # ifdef ALLOW_SEAICE
377        IF (useSEAICE) THEN        IF (useSEAICE) THEN
378           if ( siobperiod .eq. -12. _d 0 ) then
379    c     siobperiod=-12 means input file contains 12 monthly means
380    c     record numbers are assumed 1 to 12 corresponding to
381    c     Jan. through Dec.
382            call cal_GetMonthsRec(
383         O                        fac, first, changed,
384         O                        count0, count1,
385         I                        mycurrenttime, mycurrentiter, mythid
386         &           )
387    
388           elseif ( siobperiod .lt. 0. _d 0 ) then
389            print *, 'siobperiod is out of range'
390            STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_XZ'
391           else
392    c     get record numbers and interpolation factor
393            call exf_GetFFieldRec(
394         I                       siobstartdate, siobperiod,
395         I                       useYearlyFields,
396         O                       fac, first, changed,
397         O                       count0, count1, year0, year1,
398         I                       mycurrenttime, mycurrentiter, mythid
399         &                      )
400           endif
401         call exf_set_obcs_x (  OBa, OBa0, OBa1, OBafile, 's'         call exf_set_obcs_x (  OBa, OBa0, OBa1, OBafile, 's'
402       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
403       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
# Line 359  c     get record numbers and interpolati Line 414  c     get record numbers and interpolati
414       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
415       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
416       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
417         call exf_set_obcs_x ( OBuice,OBuice0,OBuice1,OBuicefile,'s'         call exf_set_obcs_x ( OBuice,OBuice0,OBuice1,OBuicefile,'u'
418       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
419       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
420       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
421         call exf_set_obcs_x ( OBvice,OBvice0,OBvice1,OBvicefile,'s'         call exf_set_obcs_x ( OBvice,OBvice0,OBvice1,OBvicefile,'v'
422       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
423       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
424       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
425        ENDIF        ENDIF
426  # endif /* ALLOW_SEAICE */  # endif /* ALLOW_SEAICE */
 # ifdef ALLOW_PTRACERS  
       if ( usePTRACERS ) then  
        do iTracer = 1, PTRACERS_numInUse  
         call exf_set_obcs_xz(  OBptr (1-Olx,1,1,1,iTracer)  
      I                       , OBptr0(1-Olx,1,1,1,iTracer)  
      I                       , OBptr1(1-Olx,1,1,1,iTracer)  
      I                       , OBptrFile(iTracer), 's'  
      I                       , fac, first, changed, useYearlyFields  
      I                       , obcsperiod, count0, count1, year0, year1  
      I                       , mycurrenttime, mycurrentiter, mythid )  
        enddo  
       endif  
 # endif /* ALLOW_PTRACERS */  
427    
428  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE and ALLOW_EXF */  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE and ALLOW_EXF */
429        RETURN        RETURN
# Line 396  C======================================= Line 438  C=======================================
438       U     OBv,   OBv0,   OBv1,   OBvfile,       U     OBv,   OBv0,   OBv1,   OBvfile,
439       U     OBt,   OBt0,   OBt1,   OBtfile,       U     OBt,   OBt0,   OBt1,   OBtfile,
440       U     OBs,   OBs0,   OBs1,   OBsfile,       U     OBs,   OBs0,   OBs1,   OBsfile,
441    #ifdef NONLIN_FRSURF
442         U     OBeta, OBeta0, OBeta1, OBetafile,
443    #endif
444  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
445       I     siobstartdate, siobperiod,       I     siobstartdate, siobperiod,
446       U     OBa,   OBa0,   OBa1,   OBafile,       U     OBa,   OBa0,   OBa1,   OBafile,
# Line 451  c     == routine arguments == Line 496  c     == routine arguments ==
496        _RL OBt1    (1-Oly:sNy+Oly,Nr,nSx,nSy)        _RL OBt1    (1-Oly:sNy+Oly,Nr,nSx,nSy)
497        _RL OBs1    (1-Oly:sNy+Oly,Nr,nSx,nSy)        _RL OBs1    (1-Oly:sNy+Oly,Nr,nSx,nSy)
498        CHARACTER*(MAX_LEN_FNAM) OBuFile,OBvFile,OBtFile,OBsFile        CHARACTER*(MAX_LEN_FNAM) OBuFile,OBvFile,OBtFile,OBsFile
499    #ifdef NONLIN_FRSURF
500          _RL OBeta   (1-Olx:sNx+Olx,nSx,nSy)
501          _RL OBeta0  (1-Olx:sNx+Olx,nSx,nSy)
502          _RL OBeta1  (1-Olx:sNx+Olx,nSx,nSy)
503          CHARACTER*(MAX_LEN_FNAM) OBetaFile
504    #endif
505  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
506        _RL     siobstartdate        _RL     siobstartdate
507        _RL     siobperiod        _RL     siobperiod
# Line 498  c     == local variables == Line 549  c     == local variables ==
549  # endif /* ALLOW_PTRACERS */  # endif /* ALLOW_PTRACERS */
550    
551  c     == end of interface ==  c     == end of interface ==
552        if ( obcsperiod .eq. -12 ) then        if ( obcsperiod .eq. -12. _d 0 ) then
553  c     obcsperiod=-12 means input file contains 12 monthly means  c     obcsperiod=-12 means input file contains 12 monthly means
554  c     record numbers are assumed 1 to 12 corresponding to  c     record numbers are assumed 1 to 12 corresponding to
555  c     Jan. through Dec.  c     Jan. through Dec.
# Line 508  c     Jan. through Dec. Line 559  c     Jan. through Dec.
559       I                        mycurrenttime, mycurrentiter, mythid       I                        mycurrenttime, mycurrentiter, mythid
560       &           )       &           )
561    
562        elseif ( obcsperiod .lt. 0 ) then        elseif ( obcsperiod .lt. 0. _d 0 ) then
563         print *, 'obcsperiod is out of range'         print *, 'obcsperiod is out of range'
564         STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_YZ'         STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_YZ'
565        else        else
# Line 521  c     get record numbers and interpolati Line 572  c     get record numbers and interpolati
572       I                       mycurrenttime, mycurrentiter, mythid       I                       mycurrenttime, mycurrentiter, mythid
573       &                      )       &                      )
574        endif        endif
   
575        call exf_set_obcs_yz(  OBu, OBu0, OBu1, OBufile, 'u'        call exf_set_obcs_yz(  OBu, OBu0, OBu1, OBufile, 'u'
576       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
577       I                     , obcsperiod, count0, count1, year0, year1       I                     , obcsperiod, count0, count1, year0, year1
# Line 538  c     get record numbers and interpolati Line 588  c     get record numbers and interpolati
588       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
589       I                     , obcsperiod, count0, count1, year0, year1       I                     , obcsperiod, count0, count1, year0, year1
590       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
591    # ifdef NONLIN_FRSURF
592           call exf_set_obcs_y ( OBeta, OBeta0, OBeta1, OBetaFile, 's'
593         I                     , fac, first, changed, useYearlyFields
594         I                     , obcsperiod, count0, count1, year0, year1
595         I                     , mycurrenttime, mycurrentiter, mythid )
596    # endif /* NONLIN_FRSURF */
597    # ifdef ALLOW_PTRACERS
598          if ( usePTRACERS ) then
599           do iTracer = 1, PTRACERS_numInUse
600            call exf_set_obcs_yz(  OBptr (1-Olx,1,1,1,iTracer)
601         I                       , OBptr0(1-Olx,1,1,1,iTracer)
602         I                       , OBptr1(1-Olx,1,1,1,iTracer)
603         I                       , OBptrFile(iTracer), 's'
604         I                       , fac, first, changed, useYearlyFields
605         I                       , obcsperiod, count0, count1, year0, year1
606         I                       , mycurrenttime, mycurrentiter, mythid )
607           enddo
608          endif
609    # endif /* ALLOW_PTRACERS */
610  # ifdef ALLOW_SEAICE  # ifdef ALLOW_SEAICE
611        IF (useSEAICE) THEN        IF (useSEAICE) THEN
612           if ( siobperiod .eq. -12. _d 0 ) then
613    c     siobperiod=-12 means input file contains 12 monthly means
614    c     record numbers are assumed 1 to 12 corresponding to
615    c     Jan. through Dec.
616            call cal_GetMonthsRec(
617         O                        fac, first, changed,
618         O                        count0, count1,
619         I                        mycurrenttime, mycurrentiter, mythid
620         &           )
621    
622           elseif ( siobperiod .lt. 0. _d 0 ) then
623            print *, 'siobperiod is out of range'
624            STOP 'ABNORMAL END: S/R OBCS_PRESCIBE_EXF_XZ'
625           else
626    c     get record numbers and interpolation factor
627            call exf_GetFFieldRec(
628         I                       siobstartdate, siobperiod,
629         I                       useYearlyFields,
630         O                       fac, first, changed,
631         O                       count0, count1, year0, year1,
632         I                       mycurrenttime, mycurrentiter, mythid
633         &                      )
634           endif
635         call exf_set_obcs_y (  OBa, OBa0, OBa1, OBafile, 's'         call exf_set_obcs_y (  OBa, OBa0, OBa1, OBafile, 's'
636       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
637       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
# Line 556  c     get record numbers and interpolati Line 648  c     get record numbers and interpolati
648       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
649       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
650       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
651         call exf_set_obcs_y ( OBuice,OBuice0,OBuice1,OBuicefile,'s'         call exf_set_obcs_y ( OBuice,OBuice0,OBuice1,OBuicefile,'u'
652       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
653       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
654       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
655         call exf_set_obcs_y ( OBvice,OBvice0,OBvice1,OBvicefile,'s'         call exf_set_obcs_y ( OBvice,OBvice0,OBvice1,OBvicefile,'v'
656       I                     , fac, first, changed, useYearlyFields       I                     , fac, first, changed, useYearlyFields
657       I                     , siobperiod, count0, count1, year0, year1       I                     , siobperiod, count0, count1, year0, year1
658       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
659        ENDIF        ENDIF
660  # endif /* ALLOW_SEAICE */  # endif /* ALLOW_SEAICE */
 # ifdef ALLOW_PTRACERS  
       if ( usePTRACERS ) then  
        do iTracer = 1, PTRACERS_numInUse  
         call exf_set_obcs_yz(  OBptr (1-Olx,1,1,1,iTracer)  
      I                       , OBptr0(1-Olx,1,1,1,iTracer)  
      I                       , OBptr1(1-Olx,1,1,1,iTracer)  
      I                       , OBptrFile(iTracer), 's'  
      I                       , fac, first, changed, useYearlyFields  
      I                       , obcsperiod, count0, count1, year0, year1  
      I                       , mycurrenttime, mycurrentiter, mythid )  
        enddo  
       endif  
 # endif /* ALLOW_PTRACERS */  
661    
662  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE and ALLOW_EXF */  #endif /* ALLOW_OBCS and ALLOW_OBCS_PRESCRIBE and ALLOW_EXF */
663        RETURN        RETURN
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22