/[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.51 by gforget, Tue Sep 4 14:58:15 2012 UTC revision 1.52 by heimbach, Sat Jan 26 14:45:56 2013 UTC
# Line 6  C $Name$ Line 6  C $Name$
6  # include "EXF_OPTIONS.h"  # include "EXF_OPTIONS.h"
7  #endif  #endif
8    
 C--  File ctrl_init.F:  
 C--   Contents  
 C--   o CTRL_INIT  
 C--   o CTRL_INIT_REC  
   
9  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
10    
11        subroutine ctrl_init( mythid )        subroutine ctrl_init( mythid )
# Line 802  c-- Line 797  c--
797  #endif /* ALLOW_GENARR3D_CONTROL */  #endif /* ALLOW_GENARR3D_CONTROL */
798    
799  c----------------------------------------------------------------------  c----------------------------------------------------------------------
800    c--
801    #ifdef ALLOW_GENTIM2D_CONTROL
802           do iarr = 1, maxCtrlTim2D
803            call ctrl_init_rec ( xx_gentim2d_file(iarr)(1:MAX_LEN_FNAM),
804         I       xx_gentim2d_startdate1(iarr),
805         I       xx_gentim2d_startdate2(iarr),
806         I       xx_gentim2d_period(iarr),
807         I       1,
808         O       xx_gentim2d_startdate(1,iarr),
809         O       diffrec, startrec, endrec,
810         I       mythid )
811    C
812            call ctrl_init_ctrlvar (
813         &       xx_gentim2d_file(iarr)(1:MAX_LEN_FNAM),
814         &       300+iarr, 400+iarr,
815         &       diffrec, startrec, endrec,
816         &       snx, sny, 1, 'c', 'xy', mythid )
817           enddo
818    #endif /* ALLOW_GENTIM2D_CONTROL */
819    
820    c----------------------------------------------------------------------
821  c----------------------------------------------------------------------  c----------------------------------------------------------------------
822    
823        call ctrl_init_wet( mythid )        call ctrl_init_wet( mythid )
# Line 851  c--------------------------------------- Line 867  c---------------------------------------
867    
868        return        return
869        end        end
   
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  
   
       subroutine ctrl_init_rec(  
      I     fldname,  
      I     fldstartdate1, fldstartdate2, fldperiod, nfac,  
      O     fldstartdate, diffrec, startrec, endrec,  
      I     mythid )  
   
 c     ==================================================================  
 c     SUBROUTINE ctrl_init_rec  
 c     ==================================================================  
 c  
 c     helper routine to compute the first and last record of a  
 c     time dependent control variable  
 c  
 c     Martin.Losch@awi.de, 2011-Mar-15  
 c  
 c     ==================================================================  
 c     SUBROUTINE ctrl_init_rec  
 c     ==================================================================  
   
       implicit none  
   
 c     == global variables ==  
 #include "SIZE.h"  
 #include "EEPARAMS.h"  
 #include "PARAMS.h"  
 #ifdef ALLOW_CAL  
 # include "cal.h"  
 #endif  
   
 c     == input variables ==  
 c     fldstartdate1/2 : start time (date/time) of fld  
 c     fldperod        : sampling interval of fld  
 c     nfac            : factor for the case that fld is an obcs variable  
 c                       in this case nfac = 4, otherwise nfac = 1  
 c     mythid          : thread ID of this instance  
       character*(*) fldname  
       integer fldstartdate1  
       integer fldstartdate2  
       _RL     fldperiod  
       integer nfac  
       integer mythid  
   
 c     == output variables ==  
 c     fldstartdate : full date from fldstartdate1 and 2  
 c     startrec     : first record of ctrl variable  
 c     startrec     : last record of ctrl variable  
 c     diffrec      : difference between first and last record of ctrl variable  
       integer fldstartdate(4)  
       integer startrec  
       integer endrec  
       integer diffrec  
   
 c     == local variables ==  
       integer i  
 #ifdef ALLOW_CAL  
       integer difftime(4)  
       _RL     diffsecs  
 #endif /* ALLOW_CAL */  
       character*(max_len_mbuf) msgbuf  
       integer il  
   
 c     == functions ==  
       integer  ilnblnk  
       external ilnblnk  
   
       if ( debugLevel .GE. debLevB ) then  
        il=ilnblnk(fldname)  
        WRITE( msgBuf,'(A,A)')  
      &     'CTRL_INIT_REC: Getting record indices for ',fldname(1:il)  
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &      SQUEEZE_RIGHT , myThid )  
       endif  
   
 c     initialise some output  
       do i = 1,4  
        fldstartdate(i) = 0  
       end do  
       startrec = 0  
       endrec   = 0  
       diffrec  = 0  
       if ( fldperiod .EQ. -12. ) then  
        startrec = 1  
        endrec   = 12*nfac  
       elseif ( fldperiod .EQ. 0. ) then  
        startrec = 1  
        endrec   = 1*nfac  
       else  
 # ifdef ALLOW_CAL  
        call cal_FullDate( fldstartdate1, fldstartdate2,  
      &                        fldstartdate , mythid )  
        call cal_TimePassed( fldstartdate, modelstartdate,  
      &                           difftime, mythid )  
        call cal_ToSeconds ( difftime, diffsecs, mythid )  
        startrec = int((modelstart + startTime - diffsecs)/  
      &                fldperiod) + 1  
        endrec   = int((modelend + startTime - diffsecs + modelstep/2)/  
      &                fldperiod) + 2  
        if ( nfac .ne. 1 ) then  
 c     This is the case of obcs.  
         startrec = (startrec - 1)*nfac + 1  
         endrec   = endrec*nfac  
        endif  
 # else /* ndef ALLOW_CAL */  
        startrec = 1  
        endrec   = (int((endTime - startTime)/fldperiod) + 1)*nfac  
 #endif /* ALLOW_CAL */  
       endif  
       diffrec  = endrec - startrec + 1  
   
       if ( debugLevel .GE. debLevB ) then  
        WRITE( msgBuf,'(A,A,A)')  
      &      'CTRL_INIT_REC: Record indices for ',fldname(1:il),':'  
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &      SQUEEZE_RIGHT , myThid )  
        WRITE( msgBuf,'(A,I10,A,I10)')  
      &      'CTRL_INIT_REC: startrec = ',startrec,', endrec = ',endrec  
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &      SQUEEZE_RIGHT , myThid )  
       endif  
   
       return  
       end  
Legend:
Removed from v.1.51  
changed lines
  Added in v.1.52
  ViewVC Help
Powered by ViewVC 1.1.22