New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 35 for trunk/NEMO/OPA_SRC/OBC – NEMO

Changeset 35 for trunk/NEMO/OPA_SRC/OBC


Ignore:
Timestamp:
2004-02-18T12:06:07+01:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE002 : Concerns Open Boundaries

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/OBC/obcdta.F90

    r26 r35  
    88   !!   'key_obc'         :                                Open Boundary Conditions 
    99   !!------------------------------------------------------------------------------ 
    10    !!   obc_dta_uvt       : read u, v, t, s data along each open boundary 
     10   !!   obc_dta           : read u, v, t, s data along each open boundary 
    1111   !!   obc_dta_psi       : read psi data along each open boundary (rigid lid only) 
    1212   !!------------------------------------------------------------------------------ 
     
    2020   USE in_out_manager  ! I/O logical units 
    2121   USE lib_mpp         ! distribued memory computing 
     22   USE dynspg_rl       !  
    2223 
    2324 
     
    3031 
    3132   !! * Accessibility 
    32    PUBLIC obc_dta_uvt                     ! routines called by step.F90 
    33    PUBLIC obc_dta_psi                     ! routines called by step.F90 
     33   PUBLIC obc_dta                         ! routines called by step.F90 
    3434    
    3535   !! * Substitutions 
     
    4141CONTAINS 
    4242 
    43    SUBROUTINE obc_dta_uvt ( kt ) 
     43   SUBROUTINE obc_dta( kt ) 
    4444      !!--------------------------------------------------------------------------- 
    45       !!                      ***  SUBROUTINE obc_dta_uvt  *** 
     45      !!                      ***  SUBROUTINE obc_dta  *** 
    4646      !!                     
    4747      !! ** Purpose :   Find the climatological  boundary arrays for the specified date,  
     
    7575      !!--------------------------------------------------------------------- 
    7676 
     77 
     78      IF( lk_dynspg_rl )   CALL obc_dta_psi( kt )     ! update bsf data at open boundaries 
     79 
     80 
    7781      ! 0. Initialization of date 
    7882      !    imois is the index (1 to 12) of the first month to be used in the  
     
    9195      IF( zxy > 1.01 .OR. zxy < 0. ) THEN 
    9296         IF(lwp) WRITE(numout,*)'           ' 
    93          IF(lwp) WRITE(numout,*)'obc_dta_uvt: Pbm with the the weight of the after field zxy ' 
     97         IF(lwp) WRITE(numout,*)'obc_dta: Pbm with the the weight of the after field zxy ' 
    9498         IF(lwp) WRITE(numout,*)'~~~~~~~~~~~' 
    9599         nstop = nstop + 1 
     
    101105      IF( kt == nit000 ) THEN 
    102106         IF(lwp) WRITE(numout,*)'           ' 
    103          IF(lwp) WRITE(numout,*)'obcdta: initial step in obc_dta_uvt' 
     107         IF(lwp) WRITE(numout,*)'obcdta: initial step in obc_dta' 
    104108         IF(lwp) WRITE(numout,*)'~~~~~~  months ',imois,' and', imoisp1,' read' 
    105109 
     
    171175                    RECL   = 4096 ) 
    172176               IF( ios > 0 ) THEN 
    173                   IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obceast.dta file ' 
    174                   IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     177                  IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obceast.dta file ' 
     178                  IF(lwp) WRITE(numout,*) '~~~~~~~' 
    175179                  nstop = nstop + 1 
    176180               END IF 
     
    186190                    RECL   = irecl ) 
    187191               IF( ios > 0 ) THEN 
    188                   IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obceast.dta file ' 
    189                   IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     192                  IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obceast.dta file ' 
     193                  IF(lwp) WRITE(numout,*) '~~~~~~~' 
    190194                  nstop = nstop + 1 
    191195               END IF 
     
    268272                 RECL   = 4096 ) 
    269273               IF( ios > 0 ) THEN 
    270                   IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcwest.dta file ' 
    271                   IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     274                  IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcwest.dta file ' 
     275                  IF(lwp) WRITE(numout,*) '~~~~~~~' 
    272276                  nstop = nstop + 1 
    273277               END IF 
     
    283287                 RECL   = irecl ) 
    284288               IF( ios > 0 ) THEN 
    285                   IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcwest.dta file ' 
    286                   IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     289                  IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcwest.dta file ' 
     290                  IF(lwp) WRITE(numout,*) '~~~~~~~' 
    287291                  nstop = nstop + 1 
    288292               END IF 
     
    338342                 RECL   = 4096 ) 
    339343            IF( ios > 0 ) THEN 
    340                IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcnorth.dta file ' 
    341                IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     344               IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcnorth.dta file ' 
     345               IF(lwp) WRITE(numout,*) '~~~~~~~' 
    342346               nstop = nstop + 1 
    343347            END IF 
     
    353357                 RECL   = irecl ) 
    354358            IF( ios > 0 ) THEN 
    355                IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcnorth.dta file ' 
    356                IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     359               IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcnorth.dta file ' 
     360               IF(lwp) WRITE(numout,*) '~~~~~~~' 
    357361               nstop = nstop + 1 
    358362            END IF 
     
    410414                 RECL   = 4096 ) 
    411415            IF( ios > 0 ) THEN 
    412                IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcsouth.dta file ' 
    413                IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     416               IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcsouth.dta file ' 
     417               IF(lwp) WRITE(numout,*) '~~~~~~~' 
    414418               nstop = nstop + 1 
    415419            END IF 
     
    425429                 RECL   = irecl ) 
    426430            IF( ios > 0 ) THEN 
    427                IF(lwp) WRITE(numout,*) 'obc_dta_uvt: Pbm to OPEN the obcsouth.dta file ' 
    428                IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     431               IF(lwp) WRITE(numout,*) 'obc_dta: Pbm to OPEN the obcsouth.dta file ' 
     432               IF(lwp) WRITE(numout,*) '~~~~~~~' 
    429433               nstop = nstop + 1 
    430434            END IF 
     
    518522      END IF 
    519523 
    520    END SUBROUTINE obc_dta_uvt 
     524   END SUBROUTINE obc_dta 
    521525 
    522526# if defined key_dynspg_fsc 
     
    700704 
    701705   END SUBROUTINE obc_dta_psi 
     706 
    702707# endif 
    703708 
    704709#else 
     710   !!------------------------------------------------------------------------------ 
     711   !!   default option:           Dummy module          NO Open Boundary Conditions 
     712   !!------------------------------------------------------------------------------ 
     713CONTAINS 
     714   SUBROUTINE obc_dta( kt )             ! Dummy routine 
     715      INTEGER, INTENT (in) :: kt 
     716      WRITE(*,*) 'obc_dta: You should not have seen this print! error?', kt 
     717   END SUBROUTINE obc_dta 
     718#endif 
     719 
    705720   !!============================================================================== 
    706    !!                            ***  MODULE obcdta  *** 
    707    !!  Empty module : no open boundary conditions 
    708    !!============================================================================== 
    709 CONTAINS 
    710    SUBROUTINE obc_dta_uvt( kt )             ! Empty routine 
    711       INTEGER, INTENT (in) :: kt 
    712       WRITE(*,*) 'obc_dta_uvt: You should not have seen this print! error?', kt 
    713    END SUBROUTINE obc_dta_uvt 
    714  
    715    SUBROUTINE obc_dta_psi( kt )             ! Empty routine 
    716       INTEGER, INTENT (in) :: kt 
    717       WRITE(*,*) 'obc_dta_psi: You should not have seen this print! error?', kt 
    718    END SUBROUTINE obc_dta_psi 
    719  
    720 #endif 
    721  
    722721END MODULE obcdta 
Note: See TracChangeset for help on using the changeset viewer.