Changeset 768 for branches/dev_001_GM/NEMO/TOP_SRC/CFC/trccfc.F90

Timestamp:

2007-12-16T15:46:18+01:00 (16 years ago)

Author:

gm

Message:

dev_001_GM - create 1 trclsm_ module by trc model (CFC, LOBSTER, PISCES..) + some bug corrections

File:

• branches/dev_001_GM/NEMO/TOP_SRC/CFC/trccfc.F90 (modified) (5 diffs)

branches/dev_001_GM/NEMO/TOP_SRC/CFC/trccfc.F90

@@ -25 +25 @@
    PUBLIC   trc_cfc       ! called in ???    
 
-   INTEGER , PARAMETER ::   jpyear = 100   ! temporal parameter 
-   INTEGER , PARAMETER ::   jphem  =   2   ! parameter for the 2 hemispheres
+   INTEGER , PUBLIC, PARAMETER ::   jpyear = 100   ! temporal parameter 
+   INTEGER , PUBLIC, PARAMETER ::   jphem  =   2   ! parameter for the 2 hemispheres
    INTEGER , PUBLIC    ::   ndate_beg      ! initial calendar date (aammjj) for CFC
    INTEGER , PUBLIC    ::   nyear_res      ! restoring time constant (year)
    INTEGER , PUBLIC    ::   nyear_beg      ! initial year (aa) 
    
-   REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, jpf_cfc:jpl_cfc) ::   p_cfc   ! partial hemispheric pressure for CFC          
+   REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, jp_cfc0:jp_cfc1) ::   p_cfc   ! partial hemispheric pressure for CFC          
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj)                       ::   xphem    ! spatial interpolation factor for patm
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj      ,jpf_cfc:jpl_cfc) ::   qtr      ! input function
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj      ,jpf_cfc:jpl_cfc) ::   qint     ! flux function
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj      ,jp_cfc0:jp_cfc1) ::   qtr      ! input function
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj      ,jp_cfc0:jp_cfc1) ::   qint     ! flux function
 
    REAL(wp), DIMENSION(jptra) ::   soa1, soa2, soa3, soa4   ! coefficient for solubility of CFC [mol/l/atm]
@@ -50 +50 @@
    !!----------------------------------------------------------------------
    !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
-   !! $Id:$ 
+   !! $Id$ 
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -107 +107 @@
 
       !                                                         !------------!
-      DO jn = jpf_cfc, jpl_cfc                                  !  CFC loop  !
+      DO jn = jp_cfc0, jp_cfc1                                  !  CFC loop  !
          !                                                      !------------!
          ! time interpolation at time kt
@@ -180 +180 @@
       !!---------------------------------------------------------------------
 
-      DO jn = 1, jptra
+      DO jn = jp_cfc0, jp_cfc1
          IF ( jn == jp11 ) THEN
             ! coefficient for solubility of CFC11 in  mol/l/atm
@@ -214 +214 @@
             sca4(jn) =    -0.067430
          ENDIF
-      ENDDO
-
-      DO jn = 1, jptra
+
          WRITE(numout,*) 'coefficient for solubility of tracer',ctrcnm(jn)
          WRITE(numout,*) soa1(jn), soa2(jn),soa3(jn), soa4(jn), &