Changeset 7041 for branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90

Timestamp:

2016-10-19T12:55:07+02:00 (8 years ago)

Author:

cetlod

Message:

ROBUST5_CNRS : implementation of part I of new TOP interface - 1st step -, see ticket #1782

File:

• branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90 (modified) (4 diffs)

branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90

@@ -6 +6 @@
    !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) 
    !!----------------------------------------------------------------------
-#if defined key_cfc
-   !!----------------------------------------------------------------------
-   !!   'key_cfc'                                               CFC tracers
    !!----------------------------------------------------------------------
    !! trc_ini_cfc      : CFC model initialisation
@@ -46 +43 @@
       INTEGER  ::  iskip = 6   ! number of 1st descriptor lines
       REAL(wp) ::  zyy, zyd
+      CHARACTER(len = 20)  ::  cltra
       !!----------------------------------------------------------------------
 
@@ -52 +50 @@
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
 
+      
+      ! assign an index in trc arrays for each CFC prognostic variables
+      jp_cfc = 1 
+      IF( ln_cfc11 .AND. ln_cfc12 )  jp_cfc = 2
+
+      ! assign an index in trc arrays for each prognostic variables
+      DO jn = 1, jptra
+        cltra = ctrcnm(jn) 
+        IF( cltra == 'CFC11'  .OR. cltra == 'cfc11' )   jpc11 = jn 
+        IF( cltra == 'CFC12'  .OR. cltra == 'cfc12' )   jpc12 = jn 
+      ENDDO
+   
+      IF( jp_cfc == 1 ) THEN
+        IF( ln_cfc11 )  jp_cfc0  = jpc11
+        IF( ln_cfc12 )  jp_cfc0  = jpc12
+      ELSE
+        jp_cfc0 = MIN( jpc11, jpc12 ) 
+      ENDIF
+      jp_cfc1 = jp_cfc0 + jp_cfc - 1
+
+      IF( lwp ) THEN
+        WRITE(numout,*) ''
+        WRITE(numout,*) ' First index of CFC tracer in the passive tracer array   jp_cfc0 = ', jp_cfc0
+        WRITE(numout,*) ' Last  index of CFC tracer in the passive tracer array   jp_cfc1 = ', jp_cfc1
+        WRITE(numout,*) 
+      ENDIF
 
       IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112atm'
@@ -146 +170 @@
    END SUBROUTINE trc_ini_cfc
    
-#else
-   !!----------------------------------------------------------------------
-   !!   Dummy module                                         No CFC tracers
-   !!----------------------------------------------------------------------
-CONTAINS
-   SUBROUTINE trc_ini_cfc             ! Empty routine
-   END SUBROUTINE trc_ini_cfc
-#endif
-
    !!======================================================================
 END MODULE trcini_cfc