Changeset 4058 for branches/2013/dev_r3996_CMCC6_topbc

Timestamp:

2013-10-11T11:51:22+02:00 (11 years ago)

Author:

vichi

Message:

Finalize the structure of generic boundary conditions for trc

Include the subroutine trcbc.F90 that implements a structure
similar to trcdta for reading surface, coastal (river) and
open boundary conditions for tracers.
A new namelist trc_bc can be included in namelist_top
where the files to be read during stepping are specified.
At the moment this structure is only used by the BFM but it
is harmless for PISCES as this model still reads those
data internally.
NOTE: OBC reading has not been implementetested as it will be made
available only for BDY.
--This line, and th se below, will be ignored--

M trcdta.F90
A trcbc.F90
M par_trc.F90

Location:

branches/2013/dev_r3996_CMCC6_topbc/NEMOGCM/NEMO/TOP_SRC

Files:

• par_trc.F90 (modified) (1 diff)
• trcbc.F90 (added)
• trcdta.F90 (modified) (4 diffs)

branches/2013/dev_r3996_CMCC6_topbc/NEMOGCM/NEMO/TOP_SRC/par_trc.F90

@@ -17 +17 @@
 
    IMPLICIT NONE
+
+   ! Passive tracers : Maximum number of tracers. Needed to define data structures
+   ! --------------- 
+   INTEGER, PUBLIC,  PARAMETER ::   jpmaxtrc = 100
 
    ! Passive tracers : Total size

branches/2013/dev_r3996_CMCC6_topbc/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

@@ -29 +29 @@
    PUBLIC   trc_dta_init    ! called in trcini.F90 
 
-   INTEGER  , PARAMETER, PUBLIC                        :: MAXTRC=100  ! maximum number of tracers 
    INTEGER  , SAVE, PUBLIC                             :: nb_trcdta   ! number of tracers to be initialised with data
    INTEGER  , SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:)  :: n_trc_index ! indice of tracer which is initialised with data
@@ -55 +54 @@
       !!----------------------------------------------------------------------
       !
-      INTEGER,INTENT(IN) :: ntrc
-      INTEGER            :: jl, jn                   ! dummy loop indicies
+      INTEGER,INTENT(IN) :: ntrc                             ! number of tracers
+      INTEGER            :: jl, jn                           ! dummy loop indices
       INTEGER            :: ierr0, ierr1, ierr2, ierr3       ! temporary integers
       CHARACTER(len=100) :: clndta, clntrc
@@ -62 +61 @@
       !
       CHARACTER(len=100) :: cn_dir
-      TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) :: slf_i     ! array of namelist informations on the fields to read
-      TYPE(FLD_N), DIMENSION(MAXTRC) :: sn_trcdta
-      REAL(wp)   , DIMENSION(MAXTRC) :: rn_trfac    ! multiplicative factor for tracer values
+      TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) :: slf_i ! array of namelist informations on the fields to read
+      TYPE(FLD_N), DIMENSION(jpmaxtrc) :: sn_trcdta
+      REAL(wp)   , DIMENSION(jpmaxtrc) :: rn_trfac    ! multiplicative factor for tracer values
       !!
       NAMELIST/namtrc_dta/ sn_trcdta, cn_dir, rn_trfac 
@@ -150 +149 @@
       ENDIF
       !
+      DEALLOCATE( slf_i )          ! deallocate local field structure
       IF( nn_timing == 1 )  CALL timing_stop('trc_dta_init')
       !